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.

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.

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.

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.

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.

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.

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.

Lithological Characterization and Its Application Based on Three-Dimensional Structure-Coupled Joint Inversion of Gravity and Magnetic Data

Incorporating structural-coupling constraint, known as the cross-gradients criterion, helps to improve the focussing trend in cross-plot of multiple physical properties. Based on this feature, a?post-processing technique is studied to characterize the lithological types of subsurface geological materials after joint inversion. A simple domain transform, which converts two kinds of participant physical properties into an artificial complex array, is adopted to extract anomalies manually from homogenous host rock. A synthetic example shows that structure-coupled joint inverted results tend to concentrate on the feature trends in the cross-plot, and the main geological targets are recovered well by a radius-azimuth plot. In a field data example, the lithological characterization?reveals that the main rock types interpreted in the study area agree with the geological information, thus demonstrating the feasibility of this technique.

A MATLAB-Based Numerical and GUI Implementation of Cross-Gradients Joint Inversion of Gravity and Magnetic Data

The cross-gradients joint inversion technique has been applied to multiple geophysical data with a significant improvement on compatibility, but its numerical implementation for practical use is rarely discussed in the literature. We present a MATLAB-based three-dimensional cross-gradients joint inversion program with application to gravity and magnetic data. The input and output information was examined with care to create a rational, independent design of a graphical user interface (GUI) and computing kernel. For 3D visualization and data file operations, UBC-GIF tools are invoked using a series of I/O functions. Some key issues regarding the iterative joint inversion algorithm are also discussed: for instance, the forward difference of cross gradients, and matrix pseudo inverse computation. A synthetic example is employed to illustrate the whole process. Joint and separate inversions can be performed flexibly by switching the inversion mode. The resulting density model and susceptibility model demonstrate the correctness of the proposed program.

Logging Evaluation Method of Low Resistivity Reservoir——A Case Study of Well Block DX12 in Junggar Basin

The Hutubi （呼图壁） River reservoir of well block DX12 is a lithologic hydrocarbon reservoir that is under tectonic settings. The main oil-bearing sand body in this area is thin and has a poor transverse connectivity. Because of the complexity of the oil-water relationship, the oil reservoir presents a low resistivity feature, which brings great difficulties to hydrocarbon reservoir identification. This article develops an effective method of well log interpretation that can meet the requirement of low resistivity reservoir well logging evaluation. The authors combine the oil reservoir geology feature, the oil well logging curve characteristics and chemical analytical data to analyze the reasons for low resistivity, then establish the appropriate reservoir parameter explanation model, which uses different saturation computational methods according to different generations. When the clay content is more than 5%, we select W-S dual water model; when the shale content is more than 13%, we use the Schlumberger formula; when the shale content is less then 13%, we use Archie＇s formula. The well logging evaluation method of low resistivity reservoir has been improved by the irreducible water saturation formula which is established by the permeability, the porosity, the coefficient of pore structure and the shale content, hydrocarbon reservoir recognition charts, and the non-resistivity logging methods （repeat formation test （RFT）; modular dynamic test （MDT）, etc.）. The coincidence rate for this arrangement of the well logging integrated interpretation is 82.6% in the well block DX12. It is a powerful direction for low resistivity well log interpretation.

Rayleigh Wave Group Velocity Distribution in Ningxia

In this article, seven months ambient noise data and 10 events recorded at seven digital stations from the Ningxia （宁夏） regional seismic network and 5 500-t controlled source explosion data recorded by 15 temporary and 7 permanent seismic stations are used to measure dispersion curves of fundamental mode Rayleigh waves. The study region was divided into grids with 0.1°×0.1°; group velocity distributions of Rayleigh waves from 6-22 s were determined with the Occam＇s inversion technique. These velocity distribution maps show the lateral velocity variations in the study area, and the velocity structures are correlated with surface geology and tectonic units. The Yinchuan （银川） basin is clearly featured with low velocities, and the Helan （贺兰） Mountain and southern mountain areas are revealed with high velocities.

Modeling of Borehole Radar for Well Logging Using Pseudo-spectral Time Domain Algorithm

In this article, numerical modeling of borehole radar for well logging in time domain is developed using pseudo-spectral time domain algorithm in axisymmetric cylindrical coordinate for proximate true formation model. The conductivity and relative permittivity logging curves are obtained from the data of borehole radar for well logging. Since the relative permittivity logging curve is not affected by salinity of formation water, borehole radar for well logging has obvious advantages as compared with conventional electrical logging. The borehole radar for well logging is a one-transmitter and two-receiver logging tool. The conductivity and relative permittivity logging curves are obtained successfully by measuring the amplitude radio and the time difference of pulse waveform from two receivers. The calculated conductivity and relative permittivity logging curves are close to the true value of surrounding formation, which tests the usability and reliability of borehole radar for well logging. The numerical modeling of borehole radar for well logging laid the important foundation for researching its logging tool.

PML Absorbing Boundary Condition for Seismic Numerical Modeling by Convolutional Differentiator in Fluid-Saturated Porous Media

The perfectly matched layer（PML） was first introduced by Berenger as an absorbing boundary condition for electromagnetic wave propagation.In this article,a method is developed to ex-tend the PML to simulating seismic wave propagation in fluid-saturated porous medium.This non-physical boundary is used at the computational edge of a Forsyte polynomial convolutional differenti-ator（FPCD） algorithm as an absorbing boundary condition to truncate unbounded media.The incor-poration of PML in Biot＇s equations is given.Numerical results show that the PML absorbing bound-ary condition attenuates the outgoing waves effectively and eliminates the reflections adequately.

Computation of Magnetic Anomalies and Gradients for Spatial Arbitrary Posture Regular Body

In the interaction computation for 3D gravity and magnetic anomalies due to arbitrarily shaped homogenous magnetized polyhedron model composed of triangular facets, there are many difficult points, such as mass computing, absence of a mature computer technique in 3D geological body modeling, inconvenient human-computer interaction, hard program coding, etc.. Based on the formulae of the magnetic field due to horizontal regular bodies, and by applying forward theory with the three-dimensional Cartesian coordinate system transformation, the forward problems of magnetic anomalies and gradient tensors for arbitrary slantwise regular bodies were solved. It is shown that the magnetic calculating expressions of the arbitrary posture regular body are corrected by comparing results with the homogeneous polyhedral body model outcome data. Furthermore, in the same condition, the former significantly reduced forward time. Applying a new forward method of regular body expressions in arbitrary posture, developed software for interaction computation between the 3D geological body model and magnetic field has advantages of fast calculation speed, easy manipulation, etc..

Characterization of Fault Zones by Analysis of Aftershock Waveform Data

Large property contrasts between materials in a fault zone and the surrounding rock are often produced by repeating earthquakes. Fault zones are usually characterized by fluid concentration, clay-rich fault gouge, increased porosity, and dilatant cracks. Thus, fault zones are thought to have reduced seismic velocities than the surrounding rocks. In this article, we first investigated the synthetic waveforms at a linear array across a vertical fault zone by using 3D finite difference simulation. Synthetic waveforms show that when sources are close to, inside, or below the fault zone, both arrival times and waveforms of P-and S-waves vary systematically across the fault zone due to reflections and transmissions from boundaries of the low-velocity fault zone. The arrival-time patterns and waveform characteristics can be used to determine the fault zone structure. Then, we applied this method to the aftershock waveform data of the 1992 Landers M7.4 and the 2008 Wenchuan （汶川） M8.0 earthquakes. Landers waveform data reveal a low-velocity zone with a width of approximately 270-370 m, and P-and S-wave velocity reductions relative to the host rock of approximately 35%-60%; Wenchuan waveform data suggest a low-velocity zone with a width of approximately 220-300 m, and P-and S-wave velocities drop relative to the host rock of approximately 55%.