A method for gravity anomaly separation based on preferential continuation and its application

基于 Pawlowski (1995 ) 建议的优先的继续方法，我们在各种各样的来源互相是 uncorrelated 与，继续高度足够大的情况中与优先的向上的继续操作员为严肃异例分离建议一个方法和过程。我们也为估计最佳介绍一个方法向上继续的高度，与不同继续高度基于分析一个综合严肃异例的优先的向上的继续操作员的特征变化。方法在铁存款上在未加工的 Bouguer 严肃数据上被测试。结果证明方法高效地并且清楚地把数据分开成地区性的异例和剩余异例。

Numerical simulation of the dual laterolog for carbonate cave reservoirs and response characteristics

洞碳酸盐形成被异质描绘，它使电的木头预言困难。知道怎么使用双 laterolog 精确地识别当前是重要的并且份量上评估洞。使用数字模拟计算电的木头回答能为洞鉴定和评估提供一个理论基础。在这份报纸，基于双 laterolog 原则，我们首先学习球形的洞用有限元素方法建模的不同尺寸(女性)，决定在在包括地与不同充满的材料学习洞模型的木头回答以后充满的抵抗力和洞之间的一种关系，并且最后学习洞上的回答用页岩和页岩填满的双 laterolog 。数字结果为碳酸盐洞水库的鉴定和评估提供一个理论基础。

Three-dimensional conjugate gradient inversion of magnetotelluric full information data

基于阻抗张肌数据， tipper 数据，和conjugate 坡度算法的分析，我们发展一三维( 3D )为转换结合坡度算法完整的信息数据决定了从的 magnetotelluric 五电并且磁场部件并且讨论方法为 3D 倒置结果的量的解释使用完整的信息数据。从合成数据的 3D 倒置的结果显示从转换的结果更好联合张肌和 tipper 数据是的阻抗比的完整的信息数据源于转换仅仅阻抗张肌数据(或 tipper 数据) 在改进分辨率和可靠性。合成例子也表明这个 3D 倒置算法的有效性和稳定性。

A joint high-resolution processing method and its application for thin inter-beds

Seismic processing characterizing thickness and borders of thin inter-beds has gradually evolved fi＇om post-stack migration to pre-stack migration, and the latter considers both vertical and lateral resolutions. As the key processing methods for improving vertical and lateral resolution, conventional deconvolution and pre-stack time migration （PSTM） are not simply dominated by the estimation and compression of the wavelet because of its instability. Therefbre, considering the variations of wavelet frequency belbre, during and alter PSTM can obtain good common reflection point （CRP） gathers and imaging profiles of thin inter-beds. Based on the frequency characteristics of the wavelet before, during and after PSTM, a joint high-resolution processing method for thin inter-beds is proposed in this paper, including inverse Q filtering for high-frequency compensation before PSTM, optimum weighting Kirchhoff PSTM for preserving high-frequencies during PSTM, and wavelet harmonizer deconvolution tier consistent processing and frequency-band broadening after PSTM. An application to real data characterized by mudstone beds in the Oriente Basin proved that the joint high-resolution processing method is effective for determining the thickness and borders of thin inter-beds and is favorable for subsequent reservoir prediction and seismic inversions.

Gaussian fitting in gamma-ray spectral decomposition

In order to extract the information of 662-keV 137Cs spectra from the overlapping peaks with 583-keV 208Tl,609-keV 214Bi,the overlapping peaks are measured by in-situ γ-ray spectrometer using a NaI(Tl) detector.The spectral model is optimized by the Gaussian fitting algorithm,and the optimized fitting indexes for fitting/original value are from 0.96 to 0.99.Gaussian fitting verified by experiment is feasible for γ-ray spectrum analysis.The full energy peak of 137Cs is extracted correctly from the overlapping peaks,it is important for in-situ γ-ray spectrometer to estimate contamination of 137Cs in radiated environment and nuclear accident.

Experimental Verification and Research for the Distortion in the Integrated Frequency Responses of the High-Pressure Sealed Cabin and Magnetic Field Sensor

Although magnetotelluric sounding method applied to the land is advanced, there are many difficulties when it is applied to marine environment, one of which is how to lay magnetic field sensors down to the seafloor to complete measurements. To protect the magnetic field sensors from intense erosion and high pressure, suitable high-pressure sealed cabins must be designed to load them. For the consideration of magnetic measurement and marine operation, the sealed pressure cabin should be nonmagnetic and transportable. Among all optional materials, LC4 super.hard aluminum alloy has the highest performance of price/quality ratio to make the sealed pressure cabin. However, it does not mean that the high-pressure sealed cabin made using LC4 will be perfect in performance. In fact, because of its weak magnetism, the pressure cabin made using LC4 has distorting effect on frequency responses of the magnetic field sensors sealed in it. This distorting effect does not affect the use of the magnetic field sensor, but if we want to eliminate its effect, we should study it by experimental measurements. In our experiment tests, frequency sweep magnetic field as excitation signal was used, and then responses of the magnetic field sensor before and after being loaded into the high-pressure sealed cabin were measured. Finally, normalized abnormal curves for the frequency responses were obtained, through which we could show how the high-pressure sealed cabin produces effects on the responses of the magnetic field sensor. Experimental results suggest that the response distortion induced by the sealed pressure cabin appears on mid- and high-frequency areas. Using experimental results as standardization data, the frequency responses collected from seafloor magnetotelluric measurements can be corrected to restore real information about the seafloor field source.

Study of sedimentary sequence cycles by well-seismic calibration

In order to solve the problems of the fine division of sedimentary sequence cycles and their change in two-dimensional space as well as lateral extension contrast, we developed a method of wavelet depth-frequency analysis. The single signal and composite signal of different Milankovitch cycles are obtained by numerical simulation. The simulated composite signal can be separated into single signals of a single frequency cycle. We also develop a well-seismic calibration insertion technology which helps to realize the calibration from the spectrum characteristics of a single well to the seismic profile. And then we determine the change and distribution characteristics of spectrum cycles in the two-dimensional space. It points out the direction in determining the variations of the regional sedimentary sequence cycles, underground strata structure and the contact relationship.

Love wave tomography in Italy from seismic ambient noise

We estimate Love wave empirical Green＇s functions from cross-correlations of ambient seismic noise to study the crust and uppermost mantle structure in Italy. Transverse-component ambient noise data from October 2005 through March 2007 recorded at 114 seismic stations from the Istituto Nazionale di Geofisica e Vulcanologia （INGV） national broadband network, the Mediterranean Very Broadband Seismographic Network （MedNet） and the Austrian Central Institute for Meteorology and Geodynamics （ZAMG） yield more than 2 000 Love wave group velocity measurements using the multiple-filter analysis technique. In the short period band （5-20 s）, the cross-correlations show clearly one-sided asymmetric feature due to non-tmiform noise distribution and high local activities, and in the long period band （〉20 s） this feature becomes weak owing to more diffusive noise distribution. Based on these measurements, Love wave group velocity dispersion maps in the 8-34 s period band are constructed, then the SH wave velocity structures from the Love wave dispersions are inverted. The final results obtained from Love wave data are overall in good agreement with those from Rayleigh waves. Both Love and Rayleigh wave inversions all reveal that the Po plain basin is resolved with low velocity at shallow depth, and the Tyrrhenian sea is characterized with higher velocity below 8 km due to its thin oceanic crust.

Electrical Structure and Fault Features of Crust and Upper Mantle beneath the Western Margin of the Qinghai-Tibet Plateau:Evidence from the Magnetotelluric Survey along Zhada-Quanshui Lake Profile

The magnetotelluric （MT） survey along the Zhada （札达）-Quanshui （泉水） Lake profile on the western margin of the Qinghai （青海）-Tibet plateau shows that the study area is divided into three tectonic provinces by the Yalung Tsangpo and Bangong （班公）-Nujiang （怒江） sutures. From south to north these are the Himalayan terrane, Gangdise terrane, and Qiangtang （羌塘） terrane. For the study area, there are widespread high-conductivity layers in the mid and lower crust, the top layers of which fluctuate intensively. The high-conductivity layer within the Gangdise terrane is deeper than those within the Qiangtang terrane and the Himalaya terrane, and the deepest high-conductivity layer is to the south of the Bangong-Nujiang suture. The top surface of the high-conductivity layer in the south of the Bangong-Nujiang suture is about 20 km lower than that in the north of it. The high-conductivity layer within the Gangdise terrane dips toward north and there are two high-conductivity layers within the crust of the southern Qiangtang terrane. In the upper crust along the profile, there are groups of lateral electrical gradient zones or distortion zones of different scales and occurrence indicating the distribution of faults and sutures along the profile. According to the electrical structure, the structural characteristics and space distribution of the Yalung Tsangpo suture, Bangong.Nujiang suture, and the major faults of Longmucuo （龙木错） and Geerzangbu are inferred.

A Study on the Classification and Well-Logging Identification of Eclogite in the Main Hole of Chinese Continental Scientific Drilling Project

Eclogite, one of the important lithologies in the main hole of the Chinese Continental Scientific Drilling （CCSD） Project, exists above the depth of 3 245 m and has distinctive responses of gamma-ray, compensating density and neutron well-logging, and so on. In this study, according to the diversities of minerals and chemical components and well-logging responses, edogites are classified from three aspects of origin, content of oxygen, and sub-mineral. We studied the logging identification method for eclogite sub-classes based on multi-element statistics and reconstructed 11 kinds of eclogite. As a result, eclogites can be divided into 6 types using well logs. In the light of this recognition, the eclogite in the main hole is divided into 20 sections, and the distribution characters of all sub-classes of eclogite are analyzed, which will provide important data for geological research of CCSD.

Nuclear magnetic resonance(NMR) microscopic simulation based on random-walk: Theory and parameters analysis

The microscopic response characteristics of nuclear magnetic resonance(NMR) are widely used for characterizing complex pore structures of rocks. Due to the prohibitive NMR experiment cost, numerical simulation was employed as an alternative approach to verify some theoretical aspects of NMR responses. Firstly, the basic principles of pore-scale NMR simulation based on random-walk method(RWM) were introduced. The RWM-simulated results were benchmarked with the analytical results for an ideal spherical pore model. Then, the effects of two numerical parameters, namely diffusion radius and walk numbers, were studied on the simulation accuracy. The simulation method is then applied to various pore models with different pore sizes and pore shapes filled with different fluids to study the microscopic NMR response characteristics. The numerical experiments are useful for understanding and interpreting NMR measurements and the simulation code provides a numerical tool to perform pixel-based digital rock analysis.

Lithospheric Electrical Structure across the Eastern Segment of the Altyn Tagh Fault on the Northern Margin of the Tibetan Plateau

Project INDEPTH （InterNational DEep Profiling of Tibet and the Himalaya） is an interdisciplinary program designed to develop a better understanding of deep structures and mechanics of the Tibetan Plateau. As a component of magnetoteUuric （MT） work in the 4th phase of the project, MT data were collected along a profile that crosses the eastern segment of the Altyn Tagh fault on the northern margin of the plateau. Time series data processing used robust algorithms to give high quality responses. Dimensionality analysis showed that 2D approach is only valid for the northern section of the profile. Consequently, 2D inversions were only conducted for the northern section, and 3D inversions were conducted on MT data from the whole profile. From the 2D inversion model, the eastern segment of the Altyn Tagh fault only appears as a crustal structure, which suggests accommodation of strike slip motion along the Altyn Tagh fault by thrusting within the Qilian block. A large-scale off-proffie conductor within the mid-lower crust of the Qilian block was revealed from the 3D inversion model, which is probably correlated with the North Qaidam thrust belt. Furthermore, the unconnected conductors from the 3D inversion model indicate that deformations in the study area are generally localized.

Research of the Conductive Structure of Crust and the Upper Mantle beneath the South-Central Tibetan Plateau

With the super-wide band magnetoteiluric sounding data of the JUong （吉隆）-Cuoqin （措勤） profile （named line 800） which was completed in 2001 and the Dingri （定日）-Cuomai （措迈） profile （named line 900） which was completed in 2004, we obtained the strike direction of each MT station by strike analysis, then traced profiles that were perpendicular to the main strike direction, and finally obtained the resistivity model of each profile by nonlinear conjugate gradients （NLCG） inversion. With these two models, we described the resistivity structure features of the crust and the upper mantle of the center-southern Tibetan plateau and its relationship with Yalung Tsangpo suture： the upper crust of the research area is a resistive layer with resistivity value range of 200-3 000 Ω.m. The depth of its bottom surface is about 15-20 km generally, but the bottom surface of resistive layer is deeper in the middle of these two profiles. At llne 900, it is about 30 km deep, and even at line 800, it is about 38 km deep. There is a gradient belt of resistivity at the depth of 15-45 km, and a conductive layer is beneath it with resistivity even less than 5 Ω.m. This conductive layer is composed of individual conductive bodies, and at the south of the Yalung Tsangpo suture, the conductive bodies are smaller with thickness about 10 km and lean to the north slightly. However, at the north of the Yalung Tsangpo suture, the conductive bodies are larger with thickness about 30 km and also lean to the north slightly. Relatively, the conductive bodies of line 900 are thinner than those of line 800, and the depth of the bottom surface of line 900 is also shallower. At last, after analyzing the effect factors to the resistivity of rocks, it was concluded that the very conductive layer was caused by partial melt or connective water in rocks. It suggests that the middle and lower crust of the center-southern Tibetan plateau is very thick, hot, flabby, and waxy.

An Explicit Fractal Interpolation Algorithm for Reconstruction of Seismic Data

Based on our previous investigation and the pioneering work of other researchers, a novel explicit fractal interpolation method based on affine transform is proposed, in which we approximate the vertical scaling factors by the locally explicit expression. Numerical experiments indicate that the explicit fractal interpolation method shows great accuracy of reconstruction of the seismic profile and yields significant improvement over wave-equation based trace interpolation methods （unified approach）.

Experimental study of the relationship between fluid density and saturation and sonic wave velocity of rock samples from the WXS Depression,South China Sea

The relationship between fluid density and saturation and sonic wave velocity of rock samples taken from the WXS Depression in the South China Sea was studied by an oil-water replacement experiment under simulated in-situ temperature and pressure conditions.Two kinds of low-density oils（0.691 and 0.749 g/cm^3） and two kinds of high-density oils（0.834 and 0.873 g/cm^3） were used to saturate the rock samples at different oil-saturation states,and the saturated P- and S-wave velocities were measured.Through Gassmann＇s equation,the theoretical P- and S-wave velocities were also calculated by the fluid replacement method.With the comparison of the measured values and the theoretical values, this study comes to the following conclusions.（1） With the increase of oil saturation and the decrease of water saturation,the P-wave velocity of rock samples saturated by low-density oil increases and the changing rule is in accord with the effective fluid theory;the P-wave velocity of rock samples saturated by high-density oil decreases and the changing rule goes against the theory.（2） With the increase of oil density（namely 0.691→0.749→0.834→0.873 g/cm^3） when oil saturation is unchanged,P-wave velocity increases gradually.（3） The S-wave velocity is always stable and is not affected by the change of oil density and saturation.The results can be used to constrain pre-stack seismic inversion,and the variation rule of sonic wave velocity is valuable for hydrocarbon identification in the study area.

A Preliminary Study on Paleomagnetism and Rock Magnetism of Eclogite from the Maobei Area

A preliminary study of paleomagnetism and rock magnetism has been done on 55 eclogite samples collected from the Chinese Continental Scientific Drilling （CCSD） site at the Maobei （毛北） area, Donghai （东海） County, Jiangsu （江苏） Province. Also the isothermal remanence, hysteresis loop, magnetic fabric, thermal susceptibility were measured, and analyses were made by electron-probe and scanning electric microscope on some samples synchronously. The result indicates that there are two groups of stable remanence, the normal and reversed magnetization. The remanence orientations are： D=94.3°, I=-29.1° and D=273.7°, I=15.4°, respectively. The magnetization intensity and the density of the samples which carry the normal magnetization are very different from those bearing reversed magnetization. The magnetic anisotropy is weak, and the minimum axis is hardly determined. The isothermal remanence and the hysteresis loop show that the magnetic carriers of the eclogite are likely SD （single domain） and PSD （pseudo.single domain） magnetite. According to the magnetic property, the cause of formation of magnetic carriers, the mechanism of the remanence, and the significance for the tectonics are discussed.

The Offset-Domain Prestack Depth Migration with Optimal Separable Approximation

The offset-domain prestack depth migration with optimal separable approximation, based on the double square root equation, is used to image complex media with large and rapid velocity variations. The method downward continues the source and the receiver wavefields simultaneously. The mixed domain algorithm with forward Fourier and inverse Fourier transform is used to construct the double square root equation wavefield extrapolation operator. This operator separates variables in the wave number domain and variables in the space domain. The phase operation is implemented in the wave number domain, whereas the time delay for lateral velocity variation is corrected in the space domain. The migration algorithm is efficient since the seismic data are not computed shot by shot. The data set test of the Marmousi model indicates that the offset-domain migration provides a satisfied seismic migration section on which complex geologic structures are imaged in media with large and rapid lateral velocity variations.

Application of Reservoir Seismic Inversion to the Damintun Sag in the Liaohe Oilfield

According to the special geologic conditions of the Damintun （大民屯） sag in the Liaohe （辽河） basin, with a complex structure and rapid lateral change of thin interbeds, the technique of logging-constraint seismic inversion based on prestack high-resolution seismic data was used in the description of oil-gas reservoirs. Reservoir seismic inversion can effectively identify underground complex geologic structures and seismic anomalous reflection volumes and quantitatively predict the distribution of sandstones in space and their variant law in combination with lithologic interpretation. This work studies the wave impedance inversion of high-resolution prestack seismic data, and logging multi-attribute data inversion, and applies these methods to the Damintun sag. As a result, the vertical resolution of reservoir prediction is raised, ability of identifying thin-interbed sand bodies is enhanced, reliability of reservoir prediction is improved, and favorable lithologic traps of this area are further confirmed. These effects are of significance in the exploration of hidden hydrocarbons in this oilfield.

High Precision Time Domain Forward Modeling for Crosshole Electromagnetic Tomography

To improve the resolution of crosshole electromagnetic tomography, high precision of forward modeling is necessary. A pseudo-spectral time domain （PSTD） forward modeling was used to simulate electromagnetic wave propagation between two boreholes. The PSTD algorithm is based on the finite difference time domain （FDTD） method and uses the fast Fourier transform （FFT） algorithm for spatial derivatives in Maxwell＇s equations. Besides having the strongpoint of the FDTD method, the calculation precision of the PSTD algorithm is higher than that of the FDTD method under the same calculation condition. The forward modeling using the PSTD method will play an important role in enhancing the resolution of crosshole electromagnetic tomography.

Simultaneous Structure-Coupled Joint Inversion of Gravity and Magnetic Data Based on a Damped Least-Squares Technique

The structure-coupled joint inversion method of gravity and magnetic data is a powerful tool for?developing improved physical property models with high resolution and compatible features;?however, the conventional procedure is inefficient due to the truncated singular values decomposition?(SVD) process at each iteration. To improve the algorithm, a technique using damped leastsquares?is adopted to calculate the structural term of model updates, instead of the truncated SVD. This?produces structural coupled density and magnetization images with high efficiency. A so-called?coupling factor is introduced to regulate the tuning of the desired final structural similarity level.?Synthetic examples show that the joint inversion results are internally consistent and achieve?higher?resolution than separated. The acceptable runtime performance of the damped least squares?technique used in joint inversion indicates that it is more suitable for practical use than the truncated SVD method.