重磁数据处理解释技术在华南地区岩体圈定与形态反演中的应用研究
|
摘要
随着国民经济的高速发展，对矿产资源的消耗量迅速增长，过去探明矿床的资源储量满足不了日益增长的需求，迫使地质工作者转变思路，另辟蹊径，开拓新的找矿空间。向地层深部进军，寻找更多、更好的隐伏盲矿体，成为地质工作不可忽视的新领域。
区域重磁资料解释是揭示深部岩体最为快速、有效的途径。区域重磁资料在划分大地构造及区域地质构造，研究岩浆岩带分布与圈定隐伏岩体等方面独特的作用。
深部和隐伏岩体对于寻找深部隐伏矿具有极为重要的价值，通常在岩体顶盖、边部以及岩体外带、内带等都是金属矿床的赋存空间。所以隐伏岩体的发现与空间范围的定量化研究，不仅可以作为直接或间接的找矿手段，可以为深部找矿、矿产资源调查评价的选区和工作部署提供重要的地球物理依据，对于我国开展下一个深度的金属矿产潜力调查具有重要意义。
本文开展并完成了以下研究工作：
第一、分别从理论模型处理和实际资料处理结果两个方面对常用岩体的定性解释方法技术进行了评价，提出了一套深部岩体的圈定、反演和综合解释的技术流程。
第二、基于国内现有的区域重力、航磁资料进行二次开发，基于传统的重磁位场转换方法，结合斜导数法等场源边界识别新技术，参考地质图、地球化学和其它资料，完成了华南地区深部隐伏岩体圈定。
第三、在重磁数据处理技术方面，提出了空间域和频率域数据扩边、网格化数据空区插值技术，实现了2.5维剖面反演模型数据和三维自动反演模型数据转换为三维人机交互程序模型文件，并将这些技术应用到华南地区深部隐伏岩体圈定与形态研究工作中，取得了令人满意的应用效果。
第四、基于研究区内岩石的密度和磁性资料，充分参考利用地质、钻探、化探资料和其它地球物理资料，对华南地区的部分岩体异常进行了三维定量反演，推断出了这些岩体的埋深、大小和空间形态。
With the rapid development of the national economy, the demand for mineralresources is growing quickly. The existing resources can not meet the increasingdemand which forces the geological researcher to change thoughts and develop newprospecting area. In order to look for more and better blind ore, exploring the deepstrata becomes a new field that can not be ignored by geological researcher.
Regional gravity and magnetic data interpretation is the most efficient andeffective way to reveal the deep rock masses. It has a unique effect on the division oftectonics and regional geological structure, the distribution of igneous rocks and thedelineation of buried rock, etc.
Delimitation of deep rock masses is extremely important in finding concealedore. Metallic mineral usually deposits in the top and edge of the rock or the outer andinner zone of the rock. So finding the buried rock and determining its scope can beused not only as a direct or indirect prospecting means, but also as an importantgeophysical basis for the deep ore prospecting, the selection of district and workdeployment of investigation and assessment of mineral resources. which hasimportant significance for further works in the potential investigation of deep metalmineral.
The research contents of this paper are as follows:
Firstly, the qualitative interpretation techniques has been evaluated from theprocessing results of theoretical model and practical data, A set of technologicalprocess has been provided to delineate, inverse and integrated interpretation of deeprock masses.
Secondly, based on the existing regional gravity, aeromagnetic data and thetraditional potential field conversion method, combined with field source boundaryidentification technology such as the oblique derivative method etc. and referencinggeological map, geochemical and other data, the deep concealed rock mass inSouthern China.has been delineated.
Thirdly, in the aspects of gravity and magnetic data processing technology, I putforward the space domain and frequency domain data expansion and grid data spaceinterpolation technique, and realize the technology of transforming the2.5dimensional profile inversion model data and3D inversion model data into3Dhuman-computer interaction program model file. Finally, these technologies areapplied to deep concealed rock delineation and morphological study in the southernChina and have achieved the good application effect.
Fourthly, based on the rock density and magnetic data in the study area, I makefull use of geology, drilling and other geophysical data and finish thethree-dimensional quantitative inversion of a part of local gravity anomaly inSouthern China and infer the depth, size and spatial morphology of these concealedrock masses.
区域重磁资料解释是揭示深部岩体最为快速、有效的途径。区域重磁资料在划分大地构造及区域地质构造，研究岩浆岩带分布与圈定隐伏岩体等方面独特的作用。
深部和隐伏岩体对于寻找深部隐伏矿具有极为重要的价值，通常在岩体顶盖、边部以及岩体外带、内带等都是金属矿床的赋存空间。所以隐伏岩体的发现与空间范围的定量化研究，不仅可以作为直接或间接的找矿手段，可以为深部找矿、矿产资源调查评价的选区和工作部署提供重要的地球物理依据，对于我国开展下一个深度的金属矿产潜力调查具有重要意义。
本文开展并完成了以下研究工作：
第一、分别从理论模型处理和实际资料处理结果两个方面对常用岩体的定性解释方法技术进行了评价，提出了一套深部岩体的圈定、反演和综合解释的技术流程。
第二、基于国内现有的区域重力、航磁资料进行二次开发，基于传统的重磁位场转换方法，结合斜导数法等场源边界识别新技术，参考地质图、地球化学和其它资料，完成了华南地区深部隐伏岩体圈定。
第三、在重磁数据处理技术方面，提出了空间域和频率域数据扩边、网格化数据空区插值技术，实现了2.5维剖面反演模型数据和三维自动反演模型数据转换为三维人机交互程序模型文件，并将这些技术应用到华南地区深部隐伏岩体圈定与形态研究工作中，取得了令人满意的应用效果。
第四、基于研究区内岩石的密度和磁性资料，充分参考利用地质、钻探、化探资料和其它地球物理资料，对华南地区的部分岩体异常进行了三维定量反演，推断出了这些岩体的埋深、大小和空间形态。
With the rapid development of the national economy, the demand for mineralresources is growing quickly. The existing resources can not meet the increasingdemand which forces the geological researcher to change thoughts and develop newprospecting area. In order to look for more and better blind ore, exploring the deepstrata becomes a new field that can not be ignored by geological researcher.
Regional gravity and magnetic data interpretation is the most efficient andeffective way to reveal the deep rock masses. It has a unique effect on the division oftectonics and regional geological structure, the distribution of igneous rocks and thedelineation of buried rock, etc.
Delimitation of deep rock masses is extremely important in finding concealedore. Metallic mineral usually deposits in the top and edge of the rock or the outer andinner zone of the rock. So finding the buried rock and determining its scope can beused not only as a direct or indirect prospecting means, but also as an importantgeophysical basis for the deep ore prospecting, the selection of district and workdeployment of investigation and assessment of mineral resources. which hasimportant significance for further works in the potential investigation of deep metalmineral.
The research contents of this paper are as follows:
Firstly, the qualitative interpretation techniques has been evaluated from theprocessing results of theoretical model and practical data, A set of technologicalprocess has been provided to delineate, inverse and integrated interpretation of deeprock masses.
Secondly, based on the existing regional gravity, aeromagnetic data and thetraditional potential field conversion method, combined with field source boundaryidentification technology such as the oblique derivative method etc. and referencinggeological map, geochemical and other data, the deep concealed rock mass inSouthern China.has been delineated.
Thirdly, in the aspects of gravity and magnetic data processing technology, I putforward the space domain and frequency domain data expansion and grid data spaceinterpolation technique, and realize the technology of transforming the2.5dimensional profile inversion model data and3D inversion model data into3Dhuman-computer interaction program model file. Finally, these technologies areapplied to deep concealed rock delineation and morphological study in the southernChina and have achieved the good application effect.
Fourthly, based on the rock density and magnetic data in the study area, I makefull use of geology, drilling and other geophysical data and finish thethree-dimensional quantitative inversion of a part of local gravity anomaly inSouthern China and infer the depth, size and spatial morphology of these concealedrock masses.
引文
[1] Ahmed Salem, Dananjay Ravat. A combined analytic signal and Euler method(AN-EUL) forautomatic interpretation of magnetic data. Geophysics,2003,68(6):1952～1961
[2] An YuLin, Huang JinMing, Chen YuDong. Forward computation of the magnetic field of a3Dbody with arbitrary boundary and continually varying magnetization. Geophysical Prospecting,2004,52(1):49-63
[3] Atchuta Rao D., Ram Babu H.V., Sanker Narayan P.V.. Relationship of magnetic anomaliesdue to subsurface features and the interpretation of sloping contacts. Geophysics,1980,45(1):32～36
[4] Valeria Cristina F. Barbosa and Joao B. C. Silva. Generalized compact gravity inversion.Geophysics,1994,59(1):57～68
[5] Bean R.I.. A rapid graphical solution for the aeromagnetic anomaly of the two-dimensionaltabular body. Geophysics,1966,31(5):963～970
[6] Beer.K.E, Burley.A.J., Tombs.J.M. The concealed granite roof in south-west Cornwall. MinerReconnaissance Programme Rep,1975, Inst.Geol. Sci.No.1
[7] Bhattacharyya B.K. Two-dimensional harmonic analysis as a tool for magnetic interpretation.Geophysics,1965,30(5):829～857
[8] Blakely R.J, Simpson R.W. Approximating edges of source bodies from magnetic or grayityanomalies. Geophysics,1986,51(7):1494～1498
[9] Bruckshaw J.M., Kunaratnam K. The interpretation of magnetic anomalies due to dikes.Geophysical Prospecting,1963,11(4):509～522
[10] C. T. Barnett. Theoretical modeling of the magnetic and gravitational fields of an arbitrarilyshaped three-dimensional body. Geophysics,1976,41(6):1353～1364
[11] Camacho A, Montesinos F.G, Vieira R. Gravity inversion by means of growing bodies.Geophysics,2000,65(1):95～101
[12] Chokri Jallouli, Kevin Mickus, et al.Gravity and aeromagnetic constraints on the extent ofCenozoic volcanic rocks within the Nefza-Tabarka region, northwestern Tunisia. Journal ofVolcanology and Geothermal Research,2003,122(1):51～68
[13] Chris Wijns,Carlos Perez, Peter KowaJczyk. Theta map: Edge detection in magnetic data.Geophysics,2005,70(4):39～43
[14] Cordell L., Grauch V.J.S. Mapping basement magnetization zones from aeromagnetic data inthe San Juan basin, New Mexico. The Utility of Regional Gravity and Magnetic Anomaly Maps.Society of Exploration Geophysicists,1985:181～197
[15] D.Atchuta Rao, H.V.Ram Babu, P.V.Sanker Narayan. Interpretation of magnetic anomaliesdue to dikes: The complex gradient method. Geophysics,1981,46(11):1572～1578
[16] Davis Carvalho de Oliveira, Roberto Dall’Agnol, et al. Gravimetric, radiometric, andmagnetic susceptibility study of the Paleoproterozoic Redenc a o and Bannach plutons, easternAmazonian Craton, Brazil: Implications for architecture and zoning of A-type granites. Journal ofSouth American Earth Sciences,2008,25(1):100～115
[17] Derek J, Airhead F, Chris Green M. New insights into magnetic derivatives for structuralmapping. The Leading Edge,2004,23(2):116～119
[18] E.Lagios, R.G.Hipkin. A Geophysical Approach to the Granite Batholith under the EasternSouthern Uplands, Scotland. Pageoph,1982.120(2):375～388
[19] Elena Galadi-Enriquez, Jesus Galindo-Zaldivar. Diapiric emplacement in the upper crust of agranitic body: the La Bazana granite (SW Spain). Tectonophysics,2003,361(1):83～96
[20] F.J.Condi, C.A.Zelt. Gravity inversion for rifted margin deep structure using extension andisostatic constraints. Geophysical Journal International,1999,138(2):435～446
[21] Fairhead J.D., Bennett K.J., Gordon D.R.H, et al. Euler: beyond the “black box”. SEGExpanded Abstracts,1994,64(13):422～424
[22] Fedi M., Rapolla A..3-D inversion of gravity and magnetic data with depth resolution.Geophysics,1999,64(2):452～460
[23] G. Florio, M. Fedi, R. Pasteka. On the application of Euler deconvolution to the analyticsignal. Geophysics,2006,71(6):87～93
[24] Green R., Stanley J.M.. Application of a Hilbert transform method to the interpretation ofsurface vehicle magnetic data. Geophysics Prospecting,1975,23(1):18～27
[25] Hansen R.O., Powlowski R.S.,Wang X. Joint use of analytic signal and amplitude ofhorizonal gradient maxima for three-dimensional gravity data interpretation:57th Ann InternatMtg, SEG Expanded Abstract,1987,100～102
[26] Hood P.. Gradient measurement in aeromagnetic surveying. Geophysics,1965,30(5):891～902
[27] Hood, P., McClure, P.J.. Gradient measurements in ground magnetic prospecting. Geophysics,1965,30(3):403～410
[28] Hornby P., Boschetti F., Horowitz F.G.. Analysis of potential field data in the wavelet domain.Geophysical Journal International,1999,137(1):175～196
[29] Hsu S.K., Coppens D., Shyu C.T.. Depth to magnetic source using the generalized analyticsignal. Geophysics,1998,63(6):1947～1957
[30] Hsu S.K., Sibuet J.C., Shyu C.T.. High-resolution detection of geologic boundaries frompotential anomalies: An enhanced analytic signal technique. Geophysics,1996,61(2):373～386
[31] Hualin Zeng, Tianfeng Wanz, Christian Teyssier, et al. The3-D geometry of the Linglonggranitic complex from2-D gravity forward modeling, Shandong Province, east China. Geophysics,2000,65(2):421～425
[32] Huang Jin Ming,Liu Kun Liang.Design and Implement on Volume Data Visualization System.Advanced Materials Research,2012,Volumes433-440:5680～5685
[33] Hugh G, Miller, Vijay Singh. Potential filed tilt-a new concept for location of potential fìeldsources. Journa1of Applied Geophysics,1994,32(2),213～217
[34] I.Romeo, R.Tejero, R.Capote, et al.3D gravity modelling of the Aguablanca Stock, tectoniccontrol and emplacement of a Variscan gabbronorite bearing a Ni–Cu–PGE ore, SW Iberia.Geological Magazine.2008,145(3):345~359
[35] J.D Cornwell, I.F Smith. A possible concealed granite beneath part of Anglesey, North Wales.Journal of the Geological Society,1993,150(1):83～87
[36] J.M.C Tombs, Results of a gravity survey of the south-west margin of Dartmoor, Devon.Mineral Reconnaissance Programme Rep,1980, Inst.Geol.Sci., No.34
[37] Joao B.C., Silva, Walter E.. Gravity inversion using convexity constraint. Geophysics,2000,65(1):102～112
[38] John A. Donato, Joan B. Megson. A buried granite batholith beneath the East Midland Shelfof the Southern North Sea Basin. Journal of the Geological Sociely, London,1990,147(1):133～140
[39] K.E. Rollin, J.M.C. Tombs. Seismic and gravity surveys over the concealed granite ridge atBosworgy, Cornwall. Mineral Reconnaissance Programme Rep.1982, Inst.Geol.Sci., No.49
[40] Kimbell. A gravity investigation of the Middleton Granite, near Inverurie, Aberdeenshire.British Geological Survey Technical Report WF/91/6(BGS Mineral Reconnaissance ProgrammeReport120),1991
[41] Koulomzine T., Lamontagne Y., Nadeau A.. New methods for the direct interpretation ofmagnetic anomalies caused by inclined dikes of infinite length. Geophysics,1970,35(5):812～830
[42] L. Hecht, J.L. Vigneresse, G. Morteani. Constraints on the origin of zonation of the granitecomplexes in the Fichtelgebirge(Germany and Czech Republic): evidence from a gravity andgeochemical study. Geol Rundsch,1994,86(Suppl): S93～S109
[43] Li Y, Oldenburg D.W..3-D inversion of gravity data. Geophysics,1998,63(1):109～119
[44] Li Y, Oldenburg D.W..3-D inversion of magnetic data. Geophysics,1996,61(2):1678～1688
[45] M. Alva-Valdivia, J. Urrutia-Fucugauchi, et al. Aeromagnetic anomalies and paleomagnetismin Jalisco and Michoacan, southern Mexico continental margin. Tectonophysics,1991,192(1):169～190
[46] M.Fedi, G.Florio. Detection of potential fields source boundaries by enhanced horizontalderivative method. Geophysical Prospecting,2001,49(1):40～58
[47] Marson I., Klingele E.E.. Advantages of using the vertical gradient of gravity for3-Dinterpretation. Geophysics,1993,58(11):1588～1595
[48] Martelet G., Sailhac P., Moreau F., et al. Characterization of geological boundaries using1-Dwavelet transform on gravity data: Theory and application to the Himalayas. Geophysics,2001,66(4):1116～1129
[49] Medeiros W.E., Silva J.B.C. Geophysical inversion using approxing equality constraints.Geophysics,1996,61(12):1678～1688
[50] Moreau F., Gibert D., Holschneider M., et al. Identification sources of potential fields withthe continuous wavelet transform Basic theory. Journal of Geophysical Research,1999,104(b3):5003～5013
[51] Nabighian M.N.. Toward a three-dimensional automatic interpretation of potential field datavia generalized Hilbert transforms: fundamental relations. Geophysics,1984,49(6):780～786
[52] Nabighian, M.N.. The analytic signal of two-dimensional magnetic bodies with polygonalcross-section: Its properties and use for automated anomaly interpretation. Geophysics,1972,37(3):507～517
[53] Olivier Boulanger, Michel Chouteau. Contstraints in3D gravity inversion. Geophysicalprospecting,2001,49(5):265～280
[54] P. John Kurianl, M. Radha Krishnal, C.G. Nambiar, et al. Gravity field and subsurfacegeometry of the Kalpatta granite, South India and the Tectonic Significance. Gondwana Research,2001,4(1):105～111
[55] Paolo Mauriello, Domenico Patella. Localization of maximum-depth gravity anomaly sourcesby a distribution of equivalent point masses. Geophysics,2001,66(5):1431～1437
[56] Paolo Mauriello, Domenico Patella. Principles of probability tomography for natural-sourceelectromagnetic induction fields. Geophysics,1999,64(5):1404～1417
[57] Patella D.. Introduction to ground surface self-potential tomography. GeophysicalProspecting,1997,45(4):653～681
[58] Pierre B. Keating. Weighted Euler deconvolution of gravity data. Geophysics,1998,63(5):1595～1603
[59] Porniaguine O., Zhdanov M. Focusing geophysical inversion images. Geophysics,1999,64(3):874～887
[60] Portniaguin O., Zhdanov M.3-D magnetic inversion with data compression and imagefocusing. Geophysics,2002,67(5):1532～1541
[61] Rao B.S.R., Prakasa Rao T.K.S., Gopala Rao D., et al. Derivatives and dike anomalyinterpretation. Pure and Applied Geophysics,1972,99(1):120～129
[62] Reid A.B., Allsop J.M., Granser H., et al. Magnetic interpretation in three dimensions usingEuler deconvolution.Geophysics,1990,55(1):80～91
[63] Roest W.R., Verhoef J., Pilkington M.. Magnetic interpretation using the3-D analytic signal.Geophysics,1992,57(1):116～125
[64] Sailhac D., Galdeano A., Gibert D., et al. Identification of source of potential field swith thecontinuous wavelet transform: Comple wavelets and application to aeromagnetic profile in FrenchGuiana. Journal of Geophysical Research,2000,105(b8):5003～5013
[65] Smith R.S., J.B.Thurston, T. Dai, I.N. MacLeod. iSPITM-The improve source parameterimaging method. Geophysical Prospecting,1998,46(2):141～15
[66] Stanley J.M., Green, R.. Gravity gradients and the interpretation of the truncated plate.Geophysics,1976,41(6):1370～1376
[67] Teresa Iuliano, Paolo Mauriello, Domenico Patella. Looking inside Mount Vesuvius bypotential fields integrated probability tomographies. Journal of Volcanology and GeothermalResearch,2002,113(2):363～378
[68] Thompson D.T.. EULDPH A new technique for making computer-assisted depth estimatesfrom magnetic data. Geophysics,1982,47(1):31～37
[69] Thurston J.B., Smith R.S.. Automatic conversion of magnetic data to depth, dip andsusceptibility contrast using the SPI method. Geophysics,1997,62(3):807～813
[70] V. Pohanka. Optimum Expression for computation of the gravity field of a homogeneouspolyhedral body. Geophysical Prospecting,1988,36(7):733～751
[71] V.J.S. Grauch, L. Cordell. Limitations of determining density or magnetic boundaries fromthe horizontal gradient of gravity or pseudogravity data. Geophysics,1987,52(1):118～121
[72] Val’eria C.E.. Gravity inversion of a discontinuous relief stabilized by weighted smoothnessconstraints on depth. Geophysics,1999,64(5):1429～1437
[73] Vallee M.A., Keating P., Smith R.S., et al. Estimating depth and model type using thecontinuous wavelet transform of magnetic data. Geophysics,2004,69(1):191～199
[74] Xiong Li. Understanding3D analytic signal amplitude. Geophysics,2006,71(2):13～16
[75]安玉林,陈玉东,黄金明.重磁勘探正反演理论方法研究的新进展.地学前缘，2003,10(1):141~149
[76]安玉林,黄金明.局部重磁场源全方位成像（续）.北京:地震出版社,2003
[77]安玉林,谭保华.局部重磁场源全方位成像（计算地球物理丛书第三集），北京:地质出版社，1997
[78]柏冠军,吴汉宁,赵希刚,等．重力资料识别鄂尔多斯盆地线性构造方法研究．地球物理学进展,2007,22(5):1386~1392
[79]班丽.相关约束重磁三维定量反演方法研究:[博士学位论文].山东:中国石油大学（华东）,2009
[80]陈国能,洛尼·格拉佩斯.花岗岩成因:原地重熔与地壳演化.湖北:中国地质大学出版社,2009
[81]陈玉东.复小波变换反演重力异常.物探与化探,2003,27(5):354～361
[82]陈玉东.利用位场连续复小波变换识辨磁场源(上).物探化探计算技术,2003,25(2):113～118
[83]陈玉东.利用位场连续复小波变换识辨磁场源(下).物探化探计算技术,2003,25(3):220～225
[84]陈玉东.应用连续复小波变换反演位场.地学前缘,2003,10(1):212～212
[85]戴明刚,曲寿利.位场梯度模法及其在碳酸盐岩地区断裂识别中的应用.地球物理学进展,2009,24(3):951~958
[86]邓晋福,滕吉文,彭聪,等.中国地球物理场特征及深部地质与成矿.北京:地质出版社,2008
[87]段本春,徐世浙．磁（重力）异常局部场与区域场分离处理中的扩边方法研究．物化探计算技术,1997,19(4):298~304
[88]范美宁.欧拉反褶积方法的研究及应用:[博士学位论文].吉林:吉林大学,2006
[89]郭华,吴燕冈,高铁.重力斜导数方法在时间域中的理论模型与研究.吉林大学学报,2006,36(9):11～14
[90]郭华,于长春,吴燕冈.改进的斜导数方法与应用.物探与化探,2009,33(2):212～216
[91]郭华,于长春,吴燕冈.非绝对值斜导数水平导数方法及在重力资料中的应用.吉林:中国地球物理,2008:195
[92]郭良辉,孟小红,石磊.磁异常ΔT三维相关成像.地球物理学报,2010,53(2):435~441
[93]郭良辉,孟小红,石磊等.重力和重力梯度数据三维相关成像.地球物理学报,2009,52(4):1098~1106
[94]何昌礼.复杂形体的高精度重力异常正演方法.物化探计算技术,1988,10(2):121~128
[95]胡浩.内蒙古道朗和都格多金属矿区重磁资料处理解释:[硕士学位论文].吉林:吉林大学,2012
[96]胡中栋,余钦范,楼海.三维解析信号法.物探化探计算技术,1995,17(3):36～42
[97]黄临平,管志宁.利用磁异常总梯度模确定磁源边界位置.华东地质学院学报,1998,21(2):143～150
[98]黄永丽.基于断层数据的三维重建理论与技术的研究:[硕士学位论文].河南:郑州大学,2004
[99]李春芳.空间域位场分离方法研究:[硕士学位论文].陕西:长安大学,2011
[100]林乃雄,汪方展,林光锋,等.福建省1:20万区域重力调查总结报告.北京:全国地质资料馆内部资料,1984
[101]林振民,陈少强.计算机上的橡皮模技术.物探化探计算技术,1996,18(1):6～16
[102]刘金兰,李庆春,赵斌.位场场源边界识别新技术及其在山西古构造带与断裂探测中的应用研究.工程地质学报,2007,15(04):569～574
[103]刘军.重力异常频率域分离方法及应用.物探与化探,1998,22(6):446~452
[104]刘银萍.利用重力异常研究虎林盆地的构造分区和基底形态:[硕士学位论文].吉林:吉林大学,2009
[105]刘英会,余学中,黎津.中国东南沿海不同类型花岗岩的磁性特征.物探与化探,2007,31(6):526~528
[106]楼海.重力方法在地壳结构研究中的应用:[博士学位论文].北京:中国地震局地球物理研究所,2001
[107]毛立峰,王绪本,高永才.大地电磁成像的效果评价.地球物理学报,2005,48(2):429～433
[108]潘玉.重力异常线性构造信号在识别断裂上的理论及应用:[硕士学位论文].陕西:长安大学,2005
[109]孙建龙.滑动平均“窗口”选择初探.物探与化探,1985,9(3):235～238
[110]孙文珂,金宜生,涂承林.勘查地球物理勘查地球化学文集第22集区域重力调查专辑.北京:地质出版社，1997
[111]孙文珂,黄崇轲,丁鹏飞.重点成矿区带的区域构造和成矿构造文集.北京:地质出版社，2001
[112]唐泽圣.三维数据场可视化.北京:清华大学出版社,1992
[113]田黔宁,吴文鹂,管志宁.任意形状重磁异常三度体人机联作反演.物探化探计算技术,2001,23(2):125～129
[114]王邦华,林盛表,邓一谦.均匀磁化多面体的磁场.地球物理学报,1980,23(4):415~426
[115]王纪恒.关于磁异常化极的若干体会.物探化探计算技术,1993,15(4):333～338
[116]王亮,张应文,刘盛光.区域重磁资料圈定贵州境内侵入岩体及局部地质构造.物探与化探,2009,33(3):245~249
[117]王懋基,程家印,刘洪学,等.中国东部区域地球物理对地质构造和矿产预测研究-华南地区研究报告.北京:全国地质资料馆内部资料,1984
[118]王乃东.有关1:25万航磁系列图的几个问题.物探与化探,2007,31(5):459~464
[119]王万银.位场总水平导数极值位置空间变化规律研究.地球物理学报,2010,53(9):2257~2270
[120]王万银,邱之云,刘金兰,等.位场数据处理中的最小曲率扩边和补空方法研究.地球物理学进展,2009,24(4):1327~1338
[121]王万银,邱之云,杨永,等.位场边缘识别方法研究进展.地球物理学进展,2010,25(1):196~210
[122]王想,李桐林. Tilt梯度及其水平导数提取重磁源边界位置.地球物理学进展,2004,19(3):625～630
[123]王绪本,毛立峰,高永才.电磁导数场概率成像方法研究.成都理工大学学报(自然科学版),2004,31(6):679～683
[124]韦昌山,蔡锦辉,李林,等.华南成矿区成矿规律和找矿方向综合研究成果报告.北京:全国地质资料馆内部资料,2003
[125]吴宣志,刘光海,薛光琦,等.富里叶变换和位场谱分析方法及其应用.北京:测绘出版社,1987
[126]吴中海,叶澄清,潘云鹤．一个改进的简单多边形凸包算法.计算机辅助设计与图形学学报,1997,9(1):9～13
[127]夏玲燕,吴汉宁,柏冠军,等.柴达木盆地航磁资料做弱信息增强技术研究及在线性构造识别中的应用.地球物型学进展,2008,23(4):1058～1062
[128]肖锋.重力数据处理方法的研究及其在钾盐矿勘探中的应用:[博士学位论文].吉林:吉林大学,2009
[129]谢勇,张虎生,张立为,等.江西省全省重磁综合编图与找矿靶区优选成果报告.北京:全国地质资料馆内部资料,2007
[130]许令周,关继腾,房文静.高次导数的概率成像原理.青岛大学学报,2003,16(4):32～36
[131]许令周,关继腾,房文静.减小概率成像分辨宽度的一种方法.物探化探计算技术,2005,27(2):138～140
[132]严加永,吕庆田,孟贵祥,等．铜陵矿集区中酸性岩体三维成像及对深部找矿方向的指示.矿床地质,2009,28(6):838～849
[133]严良俊,胡文宝,姚长利．重磁资料面积处理中的滤波增强技术与应用．勘探地球物理进展，2006,29(2):102~103
[134]杨高印．位场数据处理的一项新技术—小子域滤波法．石油地球物理勘探,1995,30(2):240~244
[135]姚长利，管志宁，吴其斌，等．欧拉反演方法分析及实用技术改进．物探与化探，2004，28(2)：150～155
[136]姚长利,郝天珧等.重磁反演约束条件及三维物性反演技术策略.物探与化探,2002,26(4):253～257
[137]姚长利,郝天珧等.重磁遗传算法三维反演中高速及有效存储方法技术.地球物理学报,2003,46(2):252～258
[138]姚长利,黎益仕,管志宁.重磁异常正反演可视化实时方法技术改进.现代地质,1998,12(1):115～122
[139]姚长利,郑元满等.重力遗传算法三维反演中动态数组优化方法.物探化探计算技术,2002.24(3):240～245
[140]余钦范,楼海.水平梯度法提取重磁源边界位置.物探化探计算技术,1994,16(4):363～367
[141]袁奎荣.隐伏花岗岩预测及深部找矿.北京:科学出版,1990
[142]翟丽娜.重磁数据处理方法在东北地区构造识别中的应用:[博士学位论文].吉林:吉林大学,2012
[143]张贵宾.重力异常的正则化等维反演研究.重磁数据处理解释应用研讨会论文集,2008:138～141
[144]张季生.用三维解析信号技术处理华南航磁异常.物探与化探,2000,24(3):190～196
[145]张季生.增强解析信号技术的应用条件.物探与化探,1999,23(4):296～300
[146]张明华,乔计花,黄金明,等.重磁电数据处理解释软件RGIS.北京:地质出版社,2011
[147]张维宸,刘建芬,谢连文.利用航磁数据推断隐伏(半隐伏)岩体.东华理工大学学报(自然科学版),2008,41(4):349～356
[148]赵希刚,吴汉宁,柏冠军,等.重磁异常解释断裂构造的处理方法及图示技术．地球物理学进展,2008,23(2):414~421
[149]郑广如,乔春贵,刘英会.高精度航磁资料圈定隐伏岩体的效果.物探与化探,2003,27(1):18~22
[150]郑元满.重磁三维物性实时可视化反演关键技术研究:[博士学位论文].北京:中国地质大学（北京）,2011
[151]钟清.区域重力资料在地质填图中的边界定位问题研究:[博士学位论文].北京:中国地质大学（北京）,2006
[152]钟清,孟小红,刘士毅．重力资料定位地质体边界问题的探讨.物探化探计算技术,2007,29(增刊):35~38
[153]周风桐,赵永贵,曹洛华,等.勘查地球物理勘查地球化学文集第24集隐伏矿地球物理研究.北京:地质出版社，2004
[154]朱国器,黎海龙,黎旭勤,等.广西桂西地区1:20万区域重力综合报告,北京:全国地质资料馆内部资料,2005
[2] An YuLin, Huang JinMing, Chen YuDong. Forward computation of the magnetic field of a3Dbody with arbitrary boundary and continually varying magnetization. Geophysical Prospecting,2004,52(1):49-63
[3] Atchuta Rao D., Ram Babu H.V., Sanker Narayan P.V.. Relationship of magnetic anomaliesdue to subsurface features and the interpretation of sloping contacts. Geophysics,1980,45(1):32～36
[4] Valeria Cristina F. Barbosa and Joao B. C. Silva. Generalized compact gravity inversion.Geophysics,1994,59(1):57～68
[5] Bean R.I.. A rapid graphical solution for the aeromagnetic anomaly of the two-dimensionaltabular body. Geophysics,1966,31(5):963～970
[6] Beer.K.E, Burley.A.J., Tombs.J.M. The concealed granite roof in south-west Cornwall. MinerReconnaissance Programme Rep,1975, Inst.Geol. Sci.No.1
[7] Bhattacharyya B.K. Two-dimensional harmonic analysis as a tool for magnetic interpretation.Geophysics,1965,30(5):829～857
[8] Blakely R.J, Simpson R.W. Approximating edges of source bodies from magnetic or grayityanomalies. Geophysics,1986,51(7):1494～1498
[9] Bruckshaw J.M., Kunaratnam K. The interpretation of magnetic anomalies due to dikes.Geophysical Prospecting,1963,11(4):509～522
[10] C. T. Barnett. Theoretical modeling of the magnetic and gravitational fields of an arbitrarilyshaped three-dimensional body. Geophysics,1976,41(6):1353～1364
[11] Camacho A, Montesinos F.G, Vieira R. Gravity inversion by means of growing bodies.Geophysics,2000,65(1):95～101
[12] Chokri Jallouli, Kevin Mickus, et al.Gravity and aeromagnetic constraints on the extent ofCenozoic volcanic rocks within the Nefza-Tabarka region, northwestern Tunisia. Journal ofVolcanology and Geothermal Research,2003,122(1):51～68
[13] Chris Wijns,Carlos Perez, Peter KowaJczyk. Theta map: Edge detection in magnetic data.Geophysics,2005,70(4):39～43
[14] Cordell L., Grauch V.J.S. Mapping basement magnetization zones from aeromagnetic data inthe San Juan basin, New Mexico. The Utility of Regional Gravity and Magnetic Anomaly Maps.Society of Exploration Geophysicists,1985:181～197
[15] D.Atchuta Rao, H.V.Ram Babu, P.V.Sanker Narayan. Interpretation of magnetic anomaliesdue to dikes: The complex gradient method. Geophysics,1981,46(11):1572～1578
[16] Davis Carvalho de Oliveira, Roberto Dall’Agnol, et al. Gravimetric, radiometric, andmagnetic susceptibility study of the Paleoproterozoic Redenc a o and Bannach plutons, easternAmazonian Craton, Brazil: Implications for architecture and zoning of A-type granites. Journal ofSouth American Earth Sciences,2008,25(1):100～115
[17] Derek J, Airhead F, Chris Green M. New insights into magnetic derivatives for structuralmapping. The Leading Edge,2004,23(2):116～119
[18] E.Lagios, R.G.Hipkin. A Geophysical Approach to the Granite Batholith under the EasternSouthern Uplands, Scotland. Pageoph,1982.120(2):375～388
[19] Elena Galadi-Enriquez, Jesus Galindo-Zaldivar. Diapiric emplacement in the upper crust of agranitic body: the La Bazana granite (SW Spain). Tectonophysics,2003,361(1):83～96
[20] F.J.Condi, C.A.Zelt. Gravity inversion for rifted margin deep structure using extension andisostatic constraints. Geophysical Journal International,1999,138(2):435～446
[21] Fairhead J.D., Bennett K.J., Gordon D.R.H, et al. Euler: beyond the “black box”. SEGExpanded Abstracts,1994,64(13):422～424
[22] Fedi M., Rapolla A..3-D inversion of gravity and magnetic data with depth resolution.Geophysics,1999,64(2):452～460
[23] G. Florio, M. Fedi, R. Pasteka. On the application of Euler deconvolution to the analyticsignal. Geophysics,2006,71(6):87～93
[24] Green R., Stanley J.M.. Application of a Hilbert transform method to the interpretation ofsurface vehicle magnetic data. Geophysics Prospecting,1975,23(1):18～27
[25] Hansen R.O., Powlowski R.S.,Wang X. Joint use of analytic signal and amplitude ofhorizonal gradient maxima for three-dimensional gravity data interpretation:57th Ann InternatMtg, SEG Expanded Abstract,1987,100～102
[26] Hood P.. Gradient measurement in aeromagnetic surveying. Geophysics,1965,30(5):891～902
[27] Hood, P., McClure, P.J.. Gradient measurements in ground magnetic prospecting. Geophysics,1965,30(3):403～410
[28] Hornby P., Boschetti F., Horowitz F.G.. Analysis of potential field data in the wavelet domain.Geophysical Journal International,1999,137(1):175～196
[29] Hsu S.K., Coppens D., Shyu C.T.. Depth to magnetic source using the generalized analyticsignal. Geophysics,1998,63(6):1947～1957
[30] Hsu S.K., Sibuet J.C., Shyu C.T.. High-resolution detection of geologic boundaries frompotential anomalies: An enhanced analytic signal technique. Geophysics,1996,61(2):373～386
[31] Hualin Zeng, Tianfeng Wanz, Christian Teyssier, et al. The3-D geometry of the Linglonggranitic complex from2-D gravity forward modeling, Shandong Province, east China. Geophysics,2000,65(2):421～425
[32] Huang Jin Ming,Liu Kun Liang.Design and Implement on Volume Data Visualization System.Advanced Materials Research,2012,Volumes433-440:5680～5685
[33] Hugh G, Miller, Vijay Singh. Potential filed tilt-a new concept for location of potential fìeldsources. Journa1of Applied Geophysics,1994,32(2),213～217
[34] I.Romeo, R.Tejero, R.Capote, et al.3D gravity modelling of the Aguablanca Stock, tectoniccontrol and emplacement of a Variscan gabbronorite bearing a Ni–Cu–PGE ore, SW Iberia.Geological Magazine.2008,145(3):345~359
[35] J.D Cornwell, I.F Smith. A possible concealed granite beneath part of Anglesey, North Wales.Journal of the Geological Society,1993,150(1):83～87
[36] J.M.C Tombs, Results of a gravity survey of the south-west margin of Dartmoor, Devon.Mineral Reconnaissance Programme Rep,1980, Inst.Geol.Sci., No.34
[37] Joao B.C., Silva, Walter E.. Gravity inversion using convexity constraint. Geophysics,2000,65(1):102～112
[38] John A. Donato, Joan B. Megson. A buried granite batholith beneath the East Midland Shelfof the Southern North Sea Basin. Journal of the Geological Sociely, London,1990,147(1):133～140
[39] K.E. Rollin, J.M.C. Tombs. Seismic and gravity surveys over the concealed granite ridge atBosworgy, Cornwall. Mineral Reconnaissance Programme Rep.1982, Inst.Geol.Sci., No.49
[40] Kimbell. A gravity investigation of the Middleton Granite, near Inverurie, Aberdeenshire.British Geological Survey Technical Report WF/91/6(BGS Mineral Reconnaissance ProgrammeReport120),1991
[41] Koulomzine T., Lamontagne Y., Nadeau A.. New methods for the direct interpretation ofmagnetic anomalies caused by inclined dikes of infinite length. Geophysics,1970,35(5):812～830
[42] L. Hecht, J.L. Vigneresse, G. Morteani. Constraints on the origin of zonation of the granitecomplexes in the Fichtelgebirge(Germany and Czech Republic): evidence from a gravity andgeochemical study. Geol Rundsch,1994,86(Suppl): S93～S109
[43] Li Y, Oldenburg D.W..3-D inversion of gravity data. Geophysics,1998,63(1):109～119
[44] Li Y, Oldenburg D.W..3-D inversion of magnetic data. Geophysics,1996,61(2):1678～1688
[45] M. Alva-Valdivia, J. Urrutia-Fucugauchi, et al. Aeromagnetic anomalies and paleomagnetismin Jalisco and Michoacan, southern Mexico continental margin. Tectonophysics,1991,192(1):169～190
[46] M.Fedi, G.Florio. Detection of potential fields source boundaries by enhanced horizontalderivative method. Geophysical Prospecting,2001,49(1):40～58
[47] Marson I., Klingele E.E.. Advantages of using the vertical gradient of gravity for3-Dinterpretation. Geophysics,1993,58(11):1588～1595
[48] Martelet G., Sailhac P., Moreau F., et al. Characterization of geological boundaries using1-Dwavelet transform on gravity data: Theory and application to the Himalayas. Geophysics,2001,66(4):1116～1129
[49] Medeiros W.E., Silva J.B.C. Geophysical inversion using approxing equality constraints.Geophysics,1996,61(12):1678～1688
[50] Moreau F., Gibert D., Holschneider M., et al. Identification sources of potential fields withthe continuous wavelet transform Basic theory. Journal of Geophysical Research,1999,104(b3):5003～5013
[51] Nabighian M.N.. Toward a three-dimensional automatic interpretation of potential field datavia generalized Hilbert transforms: fundamental relations. Geophysics,1984,49(6):780～786
[52] Nabighian, M.N.. The analytic signal of two-dimensional magnetic bodies with polygonalcross-section: Its properties and use for automated anomaly interpretation. Geophysics,1972,37(3):507～517
[53] Olivier Boulanger, Michel Chouteau. Contstraints in3D gravity inversion. Geophysicalprospecting,2001,49(5):265～280
[54] P. John Kurianl, M. Radha Krishnal, C.G. Nambiar, et al. Gravity field and subsurfacegeometry of the Kalpatta granite, South India and the Tectonic Significance. Gondwana Research,2001,4(1):105～111
[55] Paolo Mauriello, Domenico Patella. Localization of maximum-depth gravity anomaly sourcesby a distribution of equivalent point masses. Geophysics,2001,66(5):1431～1437
[56] Paolo Mauriello, Domenico Patella. Principles of probability tomography for natural-sourceelectromagnetic induction fields. Geophysics,1999,64(5):1404～1417
[57] Patella D.. Introduction to ground surface self-potential tomography. GeophysicalProspecting,1997,45(4):653～681
[58] Pierre B. Keating. Weighted Euler deconvolution of gravity data. Geophysics,1998,63(5):1595～1603
[59] Porniaguine O., Zhdanov M. Focusing geophysical inversion images. Geophysics,1999,64(3):874～887
[60] Portniaguin O., Zhdanov M.3-D magnetic inversion with data compression and imagefocusing. Geophysics,2002,67(5):1532～1541
[61] Rao B.S.R., Prakasa Rao T.K.S., Gopala Rao D., et al. Derivatives and dike anomalyinterpretation. Pure and Applied Geophysics,1972,99(1):120～129
[62] Reid A.B., Allsop J.M., Granser H., et al. Magnetic interpretation in three dimensions usingEuler deconvolution.Geophysics,1990,55(1):80～91
[63] Roest W.R., Verhoef J., Pilkington M.. Magnetic interpretation using the3-D analytic signal.Geophysics,1992,57(1):116～125
[64] Sailhac D., Galdeano A., Gibert D., et al. Identification of source of potential field swith thecontinuous wavelet transform: Comple wavelets and application to aeromagnetic profile in FrenchGuiana. Journal of Geophysical Research,2000,105(b8):5003～5013
[65] Smith R.S., J.B.Thurston, T. Dai, I.N. MacLeod. iSPITM-The improve source parameterimaging method. Geophysical Prospecting,1998,46(2):141～15
[66] Stanley J.M., Green, R.. Gravity gradients and the interpretation of the truncated plate.Geophysics,1976,41(6):1370～1376
[67] Teresa Iuliano, Paolo Mauriello, Domenico Patella. Looking inside Mount Vesuvius bypotential fields integrated probability tomographies. Journal of Volcanology and GeothermalResearch,2002,113(2):363～378
[68] Thompson D.T.. EULDPH A new technique for making computer-assisted depth estimatesfrom magnetic data. Geophysics,1982,47(1):31～37
[69] Thurston J.B., Smith R.S.. Automatic conversion of magnetic data to depth, dip andsusceptibility contrast using the SPI method. Geophysics,1997,62(3):807～813
[70] V. Pohanka. Optimum Expression for computation of the gravity field of a homogeneouspolyhedral body. Geophysical Prospecting,1988,36(7):733～751
[71] V.J.S. Grauch, L. Cordell. Limitations of determining density or magnetic boundaries fromthe horizontal gradient of gravity or pseudogravity data. Geophysics,1987,52(1):118～121
[72] Val’eria C.E.. Gravity inversion of a discontinuous relief stabilized by weighted smoothnessconstraints on depth. Geophysics,1999,64(5):1429～1437
[73] Vallee M.A., Keating P., Smith R.S., et al. Estimating depth and model type using thecontinuous wavelet transform of magnetic data. Geophysics,2004,69(1):191～199
[74] Xiong Li. Understanding3D analytic signal amplitude. Geophysics,2006,71(2):13～16
[75]安玉林,陈玉东,黄金明.重磁勘探正反演理论方法研究的新进展.地学前缘，2003,10(1):141~149
[76]安玉林,黄金明.局部重磁场源全方位成像（续）.北京:地震出版社,2003
[77]安玉林,谭保华.局部重磁场源全方位成像（计算地球物理丛书第三集），北京:地质出版社，1997
[78]柏冠军,吴汉宁,赵希刚,等．重力资料识别鄂尔多斯盆地线性构造方法研究．地球物理学进展,2007,22(5):1386~1392
[79]班丽.相关约束重磁三维定量反演方法研究:[博士学位论文].山东:中国石油大学（华东）,2009
[80]陈国能,洛尼·格拉佩斯.花岗岩成因:原地重熔与地壳演化.湖北:中国地质大学出版社,2009
[81]陈玉东.复小波变换反演重力异常.物探与化探,2003,27(5):354～361
[82]陈玉东.利用位场连续复小波变换识辨磁场源(上).物探化探计算技术,2003,25(2):113～118
[83]陈玉东.利用位场连续复小波变换识辨磁场源(下).物探化探计算技术,2003,25(3):220～225
[84]陈玉东.应用连续复小波变换反演位场.地学前缘,2003,10(1):212～212
[85]戴明刚,曲寿利.位场梯度模法及其在碳酸盐岩地区断裂识别中的应用.地球物理学进展,2009,24(3):951~958
[86]邓晋福,滕吉文,彭聪,等.中国地球物理场特征及深部地质与成矿.北京:地质出版社,2008
[87]段本春,徐世浙．磁（重力）异常局部场与区域场分离处理中的扩边方法研究．物化探计算技术,1997,19(4):298~304
[88]范美宁.欧拉反褶积方法的研究及应用:[博士学位论文].吉林:吉林大学,2006
[89]郭华,吴燕冈,高铁.重力斜导数方法在时间域中的理论模型与研究.吉林大学学报,2006,36(9):11～14
[90]郭华,于长春,吴燕冈.改进的斜导数方法与应用.物探与化探,2009,33(2):212～216
[91]郭华,于长春,吴燕冈.非绝对值斜导数水平导数方法及在重力资料中的应用.吉林:中国地球物理,2008:195
[92]郭良辉,孟小红,石磊.磁异常ΔT三维相关成像.地球物理学报,2010,53(2):435~441
[93]郭良辉,孟小红,石磊等.重力和重力梯度数据三维相关成像.地球物理学报,2009,52(4):1098~1106
[94]何昌礼.复杂形体的高精度重力异常正演方法.物化探计算技术,1988,10(2):121~128
[95]胡浩.内蒙古道朗和都格多金属矿区重磁资料处理解释:[硕士学位论文].吉林:吉林大学,2012
[96]胡中栋,余钦范,楼海.三维解析信号法.物探化探计算技术,1995,17(3):36～42
[97]黄临平,管志宁.利用磁异常总梯度模确定磁源边界位置.华东地质学院学报,1998,21(2):143～150
[98]黄永丽.基于断层数据的三维重建理论与技术的研究:[硕士学位论文].河南:郑州大学,2004
[99]李春芳.空间域位场分离方法研究:[硕士学位论文].陕西:长安大学,2011
[100]林乃雄,汪方展,林光锋,等.福建省1:20万区域重力调查总结报告.北京:全国地质资料馆内部资料,1984
[101]林振民,陈少强.计算机上的橡皮模技术.物探化探计算技术,1996,18(1):6～16
[102]刘金兰,李庆春,赵斌.位场场源边界识别新技术及其在山西古构造带与断裂探测中的应用研究.工程地质学报,2007,15(04):569～574
[103]刘军.重力异常频率域分离方法及应用.物探与化探,1998,22(6):446~452
[104]刘银萍.利用重力异常研究虎林盆地的构造分区和基底形态:[硕士学位论文].吉林:吉林大学,2009
[105]刘英会,余学中,黎津.中国东南沿海不同类型花岗岩的磁性特征.物探与化探,2007,31(6):526~528
[106]楼海.重力方法在地壳结构研究中的应用:[博士学位论文].北京:中国地震局地球物理研究所,2001
[107]毛立峰,王绪本,高永才.大地电磁成像的效果评价.地球物理学报,2005,48(2):429～433
[108]潘玉.重力异常线性构造信号在识别断裂上的理论及应用:[硕士学位论文].陕西:长安大学,2005
[109]孙建龙.滑动平均“窗口”选择初探.物探与化探,1985,9(3):235～238
[110]孙文珂,金宜生,涂承林.勘查地球物理勘查地球化学文集第22集区域重力调查专辑.北京:地质出版社，1997
[111]孙文珂,黄崇轲,丁鹏飞.重点成矿区带的区域构造和成矿构造文集.北京:地质出版社，2001
[112]唐泽圣.三维数据场可视化.北京:清华大学出版社,1992
[113]田黔宁,吴文鹂,管志宁.任意形状重磁异常三度体人机联作反演.物探化探计算技术,2001,23(2):125～129
[114]王邦华,林盛表,邓一谦.均匀磁化多面体的磁场.地球物理学报,1980,23(4):415~426
[115]王纪恒.关于磁异常化极的若干体会.物探化探计算技术,1993,15(4):333～338
[116]王亮,张应文,刘盛光.区域重磁资料圈定贵州境内侵入岩体及局部地质构造.物探与化探,2009,33(3):245~249
[117]王懋基,程家印,刘洪学,等.中国东部区域地球物理对地质构造和矿产预测研究-华南地区研究报告.北京:全国地质资料馆内部资料,1984
[118]王乃东.有关1:25万航磁系列图的几个问题.物探与化探,2007,31(5):459~464
[119]王万银.位场总水平导数极值位置空间变化规律研究.地球物理学报,2010,53(9):2257~2270
[120]王万银,邱之云,刘金兰,等.位场数据处理中的最小曲率扩边和补空方法研究.地球物理学进展,2009,24(4):1327~1338
[121]王万银,邱之云,杨永,等.位场边缘识别方法研究进展.地球物理学进展,2010,25(1):196~210
[122]王想,李桐林. Tilt梯度及其水平导数提取重磁源边界位置.地球物理学进展,2004,19(3):625～630
[123]王绪本,毛立峰,高永才.电磁导数场概率成像方法研究.成都理工大学学报(自然科学版),2004,31(6):679～683
[124]韦昌山,蔡锦辉,李林,等.华南成矿区成矿规律和找矿方向综合研究成果报告.北京:全国地质资料馆内部资料,2003
[125]吴宣志,刘光海,薛光琦,等.富里叶变换和位场谱分析方法及其应用.北京:测绘出版社,1987
[126]吴中海,叶澄清,潘云鹤．一个改进的简单多边形凸包算法.计算机辅助设计与图形学学报,1997,9(1):9～13
[127]夏玲燕,吴汉宁,柏冠军,等.柴达木盆地航磁资料做弱信息增强技术研究及在线性构造识别中的应用.地球物型学进展,2008,23(4):1058～1062
[128]肖锋.重力数据处理方法的研究及其在钾盐矿勘探中的应用:[博士学位论文].吉林:吉林大学,2009
[129]谢勇,张虎生,张立为,等.江西省全省重磁综合编图与找矿靶区优选成果报告.北京:全国地质资料馆内部资料,2007
[130]许令周,关继腾,房文静.高次导数的概率成像原理.青岛大学学报,2003,16(4):32～36
[131]许令周,关继腾,房文静.减小概率成像分辨宽度的一种方法.物探化探计算技术,2005,27(2):138～140
[132]严加永,吕庆田,孟贵祥,等．铜陵矿集区中酸性岩体三维成像及对深部找矿方向的指示.矿床地质,2009,28(6):838～849
[133]严良俊,胡文宝,姚长利．重磁资料面积处理中的滤波增强技术与应用．勘探地球物理进展，2006,29(2):102~103
[134]杨高印．位场数据处理的一项新技术—小子域滤波法．石油地球物理勘探,1995,30(2):240~244
[135]姚长利，管志宁，吴其斌，等．欧拉反演方法分析及实用技术改进．物探与化探，2004，28(2)：150～155
[136]姚长利,郝天珧等.重磁反演约束条件及三维物性反演技术策略.物探与化探,2002,26(4):253～257
[137]姚长利,郝天珧等.重磁遗传算法三维反演中高速及有效存储方法技术.地球物理学报,2003,46(2):252～258
[138]姚长利,黎益仕,管志宁.重磁异常正反演可视化实时方法技术改进.现代地质,1998,12(1):115～122
[139]姚长利,郑元满等.重力遗传算法三维反演中动态数组优化方法.物探化探计算技术,2002.24(3):240～245
[140]余钦范,楼海.水平梯度法提取重磁源边界位置.物探化探计算技术,1994,16(4):363～367
[141]袁奎荣.隐伏花岗岩预测及深部找矿.北京:科学出版,1990
[142]翟丽娜.重磁数据处理方法在东北地区构造识别中的应用:[博士学位论文].吉林:吉林大学,2012
[143]张贵宾.重力异常的正则化等维反演研究.重磁数据处理解释应用研讨会论文集,2008:138～141
[144]张季生.用三维解析信号技术处理华南航磁异常.物探与化探,2000,24(3):190～196
[145]张季生.增强解析信号技术的应用条件.物探与化探,1999,23(4):296～300
[146]张明华,乔计花,黄金明,等.重磁电数据处理解释软件RGIS.北京:地质出版社,2011
[147]张维宸,刘建芬,谢连文.利用航磁数据推断隐伏(半隐伏)岩体.东华理工大学学报(自然科学版),2008,41(4):349～356
[148]赵希刚,吴汉宁,柏冠军,等.重磁异常解释断裂构造的处理方法及图示技术．地球物理学进展,2008,23(2):414~421
[149]郑广如,乔春贵,刘英会.高精度航磁资料圈定隐伏岩体的效果.物探与化探,2003,27(1):18~22
[150]郑元满.重磁三维物性实时可视化反演关键技术研究:[博士学位论文].北京:中国地质大学（北京）,2011
[151]钟清.区域重力资料在地质填图中的边界定位问题研究:[博士学位论文].北京:中国地质大学（北京）,2006
[152]钟清,孟小红,刘士毅．重力资料定位地质体边界问题的探讨.物探化探计算技术,2007,29(增刊):35~38
[153]周风桐,赵永贵,曹洛华,等.勘查地球物理勘查地球化学文集第24集隐伏矿地球物理研究.北京:地质出版社，2004
[154]朱国器,黎海龙,黎旭勤,等.广西桂西地区1:20万区域重力综合报告,北京:全国地质资料馆内部资料,2005