Estimating Field Source Parameters of Gravity Change in North China Based on the Euler Deconvolution Method

Based on the absolute and relative gravity observations in North China from 2009 to 2014, spatial dynamic variations of the regional gravity field are obtained. We employed the Euler deconvolution method and the theoretical model to get the best estimates of parameters. Gravity field change caused by the depth and distribution in North China is calculated by back analysis. The results show the structural index that equals 1 is suitable for inversion of the gravity variation data. The inversion results indicate that the depths of anomaly field sources are spread over the Hetao fault. The research method of this paper can be used in the quantitative study on the field source and may shed new light on the interpretations of gravity change, and also provide quantitative basis for earthquake prediction index criterions based on the gravity change.     

双层界面重磁联合反演方法研究

Based on the synchronous joint gravity and magnetic inversion of single interface by Pilkington and the need of revealing Cenozoic and crystalline basement thickness in the new round of oil-gas exploration, we propose a joint gravity and magnetic inversion methodfor two-layer models by concentrating on the relationship between the change of thicknessI and position of the middle layer and anomaly and discuss the effects of the key parameters. Model tests and application to field data show the validity of this method.     

利用大地电磁资料对玄武岩成像的可行性

The magnetotelluric （MT） method has been among the favorite supporting tools for seismic imaging of sub-salt and sub-basalt targets. In this paper we present an example from Kachchh, India （where basaltic rocks overlie Mesozoic sedimentary rocks）, and discuss the feasibility of using MT method as an exploration tool in this geological setting. Our results highlight the difference in magnetotelluric response caused by the thin intrabasalt layering. The key issue addressed in this paper is what MT can and cannot provide in such geological settings. First, we compute apparent resistivity and phase response curves using representative resistivity-depth models and borehole data from the study area. Later, we compare these results to assess the plausibility of using MT to image the sub-volcanic sediments at Kachchh. Finally, we substantiate our discussion through one-dimensional inversion of the field observed MT data from this region that exhibits poor sensitivity of MT for thin basalt layers.     

依据密度分布频率特性的重力反演方法研究

Three-dimensional gravity inversion based on the mass property model is very popular in recent years. The time and efficiency of inversion algorithms is relative to the magnitude of the target mesh. One approach is to search over the entire solution space for a more refined result. However, the inversion will be difficult with the increased parameters in the large search space and the number of computations increases exponentially. ｜n this paper, we propose a novel approach based on the frequency characteristics of the density distribution over the mesh. The purposes of our study are to reduce the parameters of three- dimensional gravity inversion and to lighten the image quality of the inversion result. The results show that the new method can expedite the inversion processing and get a better geological interpretation than tradition methods.     

Geological structure of the yellow sea area from regional gravity and magnetic interpretation

Regional gravity and aeromagnetic data covering the area of 32°- 38° N, 118°-127° E at the scale of 1：1,000,000 are coordinated and integrated in a synthetic study of the South China Yellow Sea and adjacent areas. Depth to magnetic crystalline basement and its structure are determined by magnetic anomaly inversion. Depth to and thickness of the Paleozoic rock are also revealed by gravity anomaly inversion with constrains of the basement and known seismic information from several profiles. Structure units, main faults, basin boundaries, and sub-suppressions are outlined on the basis of gravity data interpretation.     

Mantle plume,large igneous province and continental breakup——Additionally discussing the Cenozoic and Mesozoic mantle plume problems in East China

Based on the former workers study results such as numerical simulation of fluid mechanics,seismic tomography of the whole earth and igneous rocks,the basic characteristics of mantle plumes are summarized in detail,namely the mantle plume,from the D″layer near the core-mantle bouundary(CMB)of 2900 km deep,is characterized by the spape of large head and thin narrow conduit,by the physical property of high temperature and low viscosity.The LIP(large igneous province)is the best exhibition when the mantle plume ascends to the surface.According to the basic characteristics of the mantle plumes and the LIP,as well as the temporal-spatial relationships between the mantle plume and continental breakup,the detailed research on petrology,geochemistry,temporal-spatial distribution,tectonic background of the Cenozoic-Mesozoic igneous rocks and gravity anomaly distribution in East China has been done.As a result,the Mesozoic igneous rocks in Southeast China should not be regarded as an example of typical LIP related to mantle plumes.for their related characteristics are not consitent with those of the typical LIPs related to mantle plumes.The Cenozoic igneous rocks in Northeast China have no the typical characteristics of mantle plumes and hotspots,so the Cenozoic volcanism in Northeast China might have no the direct relationships with the activity of mantle plumes.     

Gravity effect calculation of three-dimensional linear density distribution and its application

An algorithm for calculating gravity effect of three-dimensional (3D) linear density distribution is presented in this paper. The linear continuous density distribution is represented with 3D grid model, which has a resemblance to the velocity model used in some seismic tomography codes. The consensus in representation method of density model and velocity model facilitates the seismic-gravity-integrated interpretation or simultaneous inversion. The numerical test of synthetic data shows that although the analytical gravity formula for linear density distribution is more complex than that for piecewise constant density distribution, it takes less time to calculate the gravity effect with linear density model than that with piecewise constant density model. In addition, this method is used in the integrated interpretation of 3D seismological tomography and gravity data in Dabie Mountain area.     

The application of the high-frequency electromagnetic sounding method in the exploration of underground thermal water

High-frequency electromagnetic sounding is an electromagnetic exploration method using the natural high-frequency electromagnetic field as a field source. It has higher resolution and greater depth penetration than the direct current method and is especially fit for geothermal energy exploration and low- and mid-level groundwater detection. We introduce a successful application of high-frequency electromagnetic sounding for evaluating geothermal water resources. The high frequency electromagnetic system （MT-USA with a frequency range from 10 KHz to 1 Hz） is first applied to sample field data from China. A remote reference station is used to assure sampled data quality. We then perform 2D inversion image processing with the electrical method data. The results basically indicate the spatial distribution of underground geothermal water and provide favorable clues to finding the sources of the subsurface geothermal water in this area.     

Regional gravity field modeling based on rectangular harmonic analysis

Regional gravity field modeling with high-precision and high-resolution is one of the most important scientific objectives in geodesy,and can provide fundamental information for geophysics,geodynamics,seismology,and mineral exploration.Rectangular harmonic analysis(RHA)is proposed for regional gravity field modeling in this paper.By solving the Laplace’s equation of gravitational potential in local Cartesian coordinate system,the rectangular harmonic expansions of disturbing potential,gravity anomaly,gravity disturbance,geoid undulation and deflection of the vertical are derived,and so are the formula for signal degree variance and error degree variance of the rectangular harmonic coefficients(RHC).We also present the mathematical model and detailed algorithm for the solution of RHC using RHA from gravity observations.In order to reduce the edge effects caused by periodic continuation in RHA,we propose the strategy of extending the size of computation domain.The RHA-based modeling method is validated by conducting numerical experiments based on simulated ground and airborne gravity data that are generated from geopotential model EGM2008 and contaminated by Gauss white noise with standard deviation of 2 mGal.The accuracy of the 2.5′×2.5′geoid undulations computed from ground and airborne gravity data is 1 and 1.4cm,respectively.The standard error of the gravity disturbances that downward continued from the flight height of 4 km to the geoid is only 3.1 mGal.Numerical results confirm that RHA is able to provide a reliable and accurate regional gravity field model,which may be a new option for the representation of the fine structure of regional gravity field.     

Seismic Signatures of Rock Pore Structure

Rock pore structure is one of the important parameters in controlling both seismic wave velocity and permeability in sandstones and carbonate rocks. For a given porosity of two similar rocks with different pore structures, their acoustic wave speeds can differ 2 km/s, and permeability can span nearly six orders of magnitude from 0.01 mD to 20 D in both sandstone and limestone. In this paper, we summarize a two-parameter elastic velocity model reduced from a general poroelastic theory, to characterize the effect of pore structures on seismic wave propagation. For a given mineralogy and fluid type of a reservoir, this velocity model is defined by porosity and a frame flexibility factor, which can be used in seismic inversion and reservoir characterization to improve estimation of porosity and reserves. The frame flexibility factor can be used for quantitative classification of rock pore structure types (PST) and may be related to pore connectivity and permeability, using both poststack and prestack seismic data. This study also helps explain why amplitude versus offset analysis (AVO) in some cases fails for the purpose of fluid detection: pore structure effect on seismic waves can mask all the fluid effects, especially in carbonate rocks.     

The inversion of density structure by graphic processing unit (GPU) and identification of igneous rocks in Xisha area

Organic reefs, the targets of deep-water petroleum exploration, developed widely in Xisha area. However, there are concealed igneous rocks undersea, to which organic rocks have nearly equal wave impedance. So the igneous rocks have become interference for future exploration by having similar seismic reflection characteristics. Yet, the density and magnetism of organic reefs are very different from igneous rocks. It has obvious advantages to identify organic reefs and igneous rocks by gravity and magnetic data. At first, frequency decomposition was applied to the free-air gravity anomaly in Xisha area to obtain the 2D subdivision of the gravity anomaly and magnetic anomaly in the vertical direction. Thus, the distribution of igneous rocks in the horizontal direction can be acquired according to high-frequency field, low-frequency field, and its physical properties. Then, 3D forward modeling of gravitational field was carried out to establish the density model of this area by reference to physical properties of rocks based on former researches. Furthermore, 3D inversion of gravity anomaly by genetic algorithm method of the graphic processing unit (GPU) parallel processing in Xisha target area was applied, and 3D density structure of this area was obtained. By this way, we can confine the igneous rocks to the certain depth according to the density of the igneous rocks. The frequency decomposition and 3D inversion of gravity anomaly by genetic algorithm method of the GPU parallel processing proved to be a useful method for recognizing igneous rocks to its 3D geological position. So organic reefs and igneous rocks can be identified, which provide a prescient information for further exploration.     

APPLICATION OF 3D HIGH PRECISION GRAVITY METHOD TO THE EXPLORATION OF DEEP BURIED HILL

Due to the influence of volcanic rock or weathering crust coverage, the quality of deep seismic data is often not reliable and it is difficult to identify small deep buried hill structure beneath the cover. However, the gravity exploration technology can use the remarkable density differences between the object and the overlying strata to identify this special geological reservoir. Although recently several rounds of seismic exploration have been carried out in Raoyang depression, Huabei Oilfield, North China to determine the existence of Hu 8 north deep buried hill, whether the buried hill really exists or not as well as its scale is still in controversial. In this paper, based on the information of seismic data and formation density, deep processing of 3D high precision ground gravity data has been carried out for 3D forward and inversion computation. The dimensional gravity forward calculation results show that the density of the deep anomaly body forms a relatively low gravity anomaly in the earth's surface. By using the potential processing method of vertical second derivative, and sliding filtering, the residual anomaly is separated from regional or background field, which illustrates the existence of a local high gravity anomaly at Hu 8 north area. According to the amplitude of residual gravity anomaly and formation density difference modeling, through a number of 3D forward calculations and 3D inversion of gravity data, the existence of Hu 8 north buried hill and its possible scale are basically determined. The results prove that the 3D high precision gravity method is effective in determining the deep buried hill structure in case that seismic data is not reliable.     

基于二阶导数的磁源边界与顶部深度快速反演

为实现磁异常反演,首先提出了磁异常非参数快速反演的概念,即无需提供先验信息.在Nabighian提出的磁场通用梯度公式的基础上,推导实现了基于磁异常垂向二阶导数的非参数快速反演(V2D_depth).它不仅可以获取场源的边界信息,同时可以反演场源的埋深.通过与Tilt_depth方法对比,本文方法计算的场源边界更清晰,反演的深度也更接近真实深度,同时较大程度上克服了Tilt_depth方法受叠加异常的影响.理论模型验证了方法的有效性,并通过准噶尔盆地某区块磁异常数据的处理,提取了受强大的区域背景场掩盖的石炭系火成岩产生的弱异常,突出了构造分区、断裂分布等信息,获取了火成岩的位置及埋深参数,为该区的火成岩油气藏勘探提供了有效的处理途径.     

综合物探技术识别地震反射异常体属性

埋藏在6000m以下的地震反射异常体是礁体还是火成岩体？或者是海底滑塌体？在高质量重磁电资采集的基础上，通过对研究区域内的地层岩石物性细致研究，利用CEMP资料界面约束二维反演、重磁力异常剥离技术和LCT重磁力地震联合反演解释技术细致研究了异常体的物性参数，从密度、磁化率、电阻率、速度等多种特征推测判断了深层异常体地质属性，为钻井部署提供了间接依据。     

重磁三维反演伊犁盆地中部密度和磁性结构

伊犁盆地中部区域地处多个构造单元的交汇区,构造复杂,缺乏高精度地质资料,制约了对盆地构造演化的进一步认识.本文联合使用多种先验信息作为约束条件,利用高精度重磁数据反演获得了研究区地下海拔-10 km以浅的三维密度和磁性结构,增强了反演结果的可靠性.结果揭示以白石墩次凸为中心发育了一个"北断南超"的南西-北东向不对称型凹陷,北部沉积厚度大于南部;区内主要断裂具有高角度特征,生烃中心受断裂控制,阿吾拉勒山前凹陷浅部可能发育逆掩断裂;火成岩以侵入的基性火成岩为主,伴有部分喷出火成岩和中酸性火成岩.     

Improved preconditioned conjugate gradient algorithm and application in 3D inversion of gravity-gradiometry data

With the continuous development of full tensor gradiometer (FTG) measurement techniques, three-dimensional (3D) inversion of FTG data is becoming increasingly used in oil and gas exploration. In the fast processing and interpretation of large-scale high-precision data, the use of the graphics processing unit process unit (GPU) and preconditioning methods are very important in the data inversion. In this paper, an improved preconditioned conjugate gradient algorithm is proposed by combining the symmetric successive over-relaxation (SSOR) technique and the incomplete Choleksy decomposition conjugate gradient algorithm (ICCG). Since preparing the preconditioner requires extra time, a parallel implement based on GPU is proposed. The improved method is then applied in the inversion of noisecontaminated synthetic data to prove its adaptability in the inversion of 3D FTG data. Results show that the parallel SSOR-ICCG algorithm based on NVIDIA Tesla C2050 GPU achieves a speedup of approximately 25 times that of a serial program using a 2.0 GHz Central Processing Unit (CPU). Real airborne gravity-gradiometry data from Vinton salt dome (southwest Louisiana, USA) are also considered. Good results are obtained, which verifies the efficiency and feasibility of the proposed parallel method in fast inversion of 3D FTG data.     

南海重磁异常特征及火成岩分布(英文)

南海火成岩油气藏具有广阔的勘探前景,综合利用地球物理方法圈划与识别火成岩体、研究火成岩分布是火成岩油气藏研究的基础。针对南海重磁场特征,采用低纬度、变倾角化极技术进行了磁异常化极处理,利用优选延拓方法实现重磁异常分离并提取南海海域浅部火成岩重磁异常信息,利用磁异常三维相关成像给出南海火成岩的三维空间等效分布,在重磁梯度突出局部异常边界信息的基础上,通过梯度加权的重磁相关分析勾画不同类型火成岩的平面展布,火成岩的分布特征显示出受地壳深部结构及断裂构造的控制与影响。     

四川地区流动重力资料的位场 分离与异常特征提取

用于监测四川地区重力场随时间变化的流动重力测网于2008年5月12日汶川M_S8.0地震后进行了全面改造, 形成了新的闭合环线并对整网点位每年实施2次重力观测工作. 该文以2010—2012年获取的5期流动重力资料为基础, 分别采用小波多尺度分解和三维视密度反演等位场分离方法, 对测区内重力场的动态变化进行多层次分解及异常特征提取, 以探求四川测区范围内重力异常动态变化特征以及地震前后不同深度范围内密度变化情况与重力场响应之间的关系. 结果表明, 对同一期资料进行小波分解得到不同阶次的细节图像, 揭示了不同深度处重力异常动态变化特征; 分解结果同时也表明了四川测区范围内重力场的变化形态与川内特有的断裂构造格架和深部结构环境及过渡带特征密切相关. 通过对不同时期的差分重力异常资料进行三维视密度反演对比研究发现, 在地震前后四川测区范围内由浅到深密度变化呈现出不同程度的差异, 震区附近出现较为明显的横向密度变化特征且浅层物质密度局部性的异常变化差异较为突出, 测区范围内中上地壳深度层密度反演切片显示出的趋势性差异基本上反映了地震前后测区内地下物质的动态变化趋势, 因此, 深刻认识不同深度范围内重力场的动态演化特征, 对研究壳内物质密度变化以及与地震的孕育和发展有关的深部介质环境具有积极的意义.      

基于GPU并行的重力、重力梯度三维正演快速计算及反演策略

利用NVIDIA CUDA编程平台,实现了基于GPU并行的重力、重力梯度三维快速正演计算方法.采用当前在重力数据约束反演或联合反演中流行的物性模型(密度大小不同、规则排列的长方体单元)作为地下剖分单元,对任意三维复杂模型体均可用很多物性模型进行组合近似,利用解析方法计算出所有物性模型在计算点的异常值并累加求和,得到整个模型体在某一计算点引起的重力(或重力梯度)值.针对精细的复杂模型体产生的问题,采用GPU并行计算技术,主要包括线程有效索引与优化的并行归约技术进行高效计算.在显卡型号为NVIDIA Quadro 2000相对于单线程CPU程序,重力和重力梯度U_xx、U_xy正演计算可以分别达到60与50倍的加速.本文还讨论了GPU并行计算在两种反演方法中的策略,为快速三维反演技术提供了借鉴.     

基于地球物理方法预测二连盆地及邻区晚古生代潜在烃源岩分布

二连盆地及邻区是我国油气勘探的重点区域之一,普遍认为烃源岩层主要分布在中生代地层,对其晚古生代地层的烃源岩潜力研究较少,本文对研究区内晚古生代碳酸盐岩及泥岩进行了有机地球化学实验并分析其烃源岩潜力,评价结果表明二连盆地及邻区内晚古生代地层具有一定的生烃潜力,受高成熟度的影响,剩余潜力以生气能力为主.依据研究区内电法剖面、布格重力异常特征及磁性资料认为晚古生代潜在烃源岩层具有"低电阻、高密度、弱磁性"的地球物理性质,并初步预测了二连盆地及邻区晚古生代的深度大致分布在3～5km的范围内,进一步利用优化后的小子域滤波重力异常分离技术对晚古生代地层进行了厚度预测,刻画其平面分布特征,认为研究区内的晚古生代地层具有两个分布中心,其中心区厚度范围为1000～2000m之间,与野外地质剖面的厚度基本一致,晚古生代潜在烃源岩的平面预测可以为二连盆地油气勘探工作提供新方向.