2013年吉林前郭M_S5.5震群的双差法重新定位及震源机制

基于区域数字地震台网的震相资料、松辽盆地中部地壳速度结构,使用双差定位方法对2013年10月31日开始的吉林省前郭爆发性震群序列进行重新定位。根据重新定位后的地震群空间分布图像结合震中区附近的物探资料探讨发震构造的特征。重定位后EW、NS和UD 3个方向大致的定位精度分别为0.9km、0.7km和1.2km,地震序列震呈NW向密集条带状分布,长度约12km,正交NE向,宽约6km。震群序列中地震震源深度较浅,80%的地震分布在6~8km深度范围内,在深度剖面上震群序列表现为明显的密集现象。根据重新定位后震群序列的空间分布特征、序列中等震连发特征以及较大级别地震的震源机制解,推断发震构造为克山-大安断裂带内的推覆逆冲构造。冲断面向E倾斜,近地表倾角较陡,深部几乎平直,表现为典型的铲形逆冲断层的特征,断层面纵向破裂长度大于横向破裂长度。根据发震构造的特征可以推断3个5级地震发生在下层冲断面上,余震主要在冲断面上部触发。

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 noise- contaminated 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 airbome gravity-gradiometry data from Vinton salt dome （south- west 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.

Seismic imaging of the middle and upper crust by double-difference tomography:the Haicheng earthquake (Ms 7.3) in Liaoning Province

The Haicheng earthquake （Ms 7.3） occurred in Liaoning Province （39°N-43°N, 120°E-126°E ）, China on February 4, 1975. The mortality rate was only 0.02% owing to the first timely and accurate prediction, although the area affected by the earthquake was 9200 km^2 and covered cities with a population density of 1000 p/km^2. In this study, the double- difference （DD） tomography method was used to obtain high-resolution three-dimensional （3D） P- and S-wave velocity （Vp and Vs） structures and Vp/Vs as well as the earthquake locations. Tomography results suggest that velocity structure at shallow depth coincides well with topography and sediment thickness. The earthquake locations form a northwest-striking zone associated with the Jinzhou（JZ） Fault and a northeast-striking zone associated with the Haichenghe-Dayanghe （HD） Fault, and suggest that the JZ Fault consists of three faults and the Ms 7.3 Haicheng earthquake originated at the intersection of the JZ and the Faults. Low- velocity zones （LVZs） with low Vp/Vs are observed at 15-20 km depth beneath the Haicheng （HC） region. We interpret the LVZs in the middle crust as regions of fluids, suggesting rock dehydration at high temperatures. The LVZs and low Vp/Vs in the upper crust are attributed to groundwater-filled cracks and pores. We believe that large crustal earthquakes in this area are caused by the combination of faulting and fluid movement in the middle crust.

Efficient finite-volume simulation of the LWD orthogonal azimuth electromagnetic response in a three-dimensional anisotropic formation using potentials on cylindrical meshes

In this study,the cylindrical finite-volume method(FVM)is advanced for the efficient and high-precision simulation of the logging while drilling(LWD)orthogonal azimuth electromagnetic tool(OAEMT)response in a three-dimensional(3 D)anisotropic formation.To overcome the ill-condition and convergence problems arising from the low induction number,Maxwell’s equations are reformulated into a mixed Helmholtz equation for the coupled potentials in a cylindrical coordinate system.The electrical fi eld continuation method is applied to approximate the perfectly electrical conducting(PEC)boundary condition,to improve the discretization accuracy of the Helmholtz equation on the surface of metal mandrels.On the base,the 3 D FVM on Lebedev’s staggered grids in the cylindrical coordinates is employed to discretize the mixed equations to ensure good conformity with typical well-logging tool geometries.The equivalent conductivity in a non-uniform element is determined by a standardization technique.The direct solver,PARDISO,is applied to efficiently solve the sparse linear equation systems for the multi-transmitter problem.To reduce the number of calls to PARDISO,the whole computational domain is divided into small windows that contain multiple measuring points.The electromagnetic(EM)solutions produced by all the transmitters per window are simultaneously solved because the discrete matrix,relevant to all the transmitters in the same window,is changed.Finally,the 3 D FVM is validated against the numerical mode matching method(NMM),and the characteristics of both the coaxial and coplanar responses of the EM field tool are investigated using the numerical results.

Airborne electromagnetic data denoising based on dictionary learning

Time-domain airborne electromagnetic(AEM)data are frequently subject to interference from various types of noise,which can reduce the data quality and affect data inversion and interpretation.Traditional denoising methods primarily deal with data directly,without analyzing the data in detail;thus,the results are not always satisfactory.In this paper,we propose a method based on dictionary learning for EM data denoising.This method uses dictionary learning to perform feature analysis and to extract and reconstruct the true signal.In the process of dictionary learning,the random noise is fi ltered out as residuals.To verify the eff ectiveness of this dictionary learning approach for denoising,we use a fi xed overcomplete discrete cosine transform(ODCT)dictionary algorithm,the method-of-optimal-directions(MOD)dictionary learning algorithm,and the K-singular value decomposition(K-SVD)dictionary learning algorithm to denoise decay curves at single points and to denoise profi le data for diff erent time channels in time-domain AEM.The results show obvious diff erences among the three dictionaries for denoising AEM data,with the K-SVD dictionary achieving the best performance.

3D density inversion of gravity gradiometry data with a multilevel hybrid parallel algorithm

The density inversion of gravity gradiometry data has attracted considerable attention;however,in large datasets,the multiplicity and low depth resolution as well as efficiency are constrained by time and computer memory requirements.To solve these problems,we improve the reweighting focusing inversion and probability tomography inversion with joint multiple tensors and prior information constraints,and assess the inversion results,computing efficiency,and dataset size.A Message Passing Interface(MPI)-Open Multi-Processing(OpenMP)-Computed Unified Device Architecture(CUDA)multilevel hybrid parallel inversion,named Hybrinv for short,is proposed.Using model and real data from the Vinton Dome,we confirm that Hybrinv can be used to compute the density distribution.For data size of 100×100×20,the hybrid parallel algorithm is fast and based on the run time and scalability we infer that it can be used to process the large-scale data.

3D forward modeling and response analysis for marine CSEMs towed by two ships

A dual-ship-towed marine electromagnetic （EM） system is a new marine exploration technology recently being developed in China. Compared with traditional marine EM systems, the new system tows the transmitters and receivers using two ships, rendering it unnecessary to position EM receivers at the seafloor in advance. This makes the system more flexible, allowing for different configurations （e.g., in-line, broadside, and azimuthal and concentric scanning） that can produce more detailed underwater structural information. We develop a three-dimensional goal-oriented adaptive forward modeling method for the new marine EM system and analyze the responses for four survey configurations. Ocean-ottom topography has a strong effect on the marine EM responses; thus, we develop a forward modeling algorithm based on the finite-element method and unstructured grids. To satisfy the requirements for modeling the moving transmitters of a dual-ship-towed EM system, we use a single mesh for each of the transmitter locations. This mitigates the mesh complexity by refining the grids near the transmitters and minimizes the computational cost. To generate a rational mesh while maintaining the accuracy for single transmitter, we develop a goal-oriented adaptive method with separate mesh refinements for areas around the transmitting source and those far away. To test the modeling algorithm and accuracy, we compare the EM responses calculated by the proposed algorithm and semi-analytical results and from published sources. Furthermore, by analyzing the EM responses for four survey configurations, we are confirm that compared with traditional marine EM systems with only in-line array, a dual-ship-towed marine system can collect more data.

Translation-invariant wavelet denoising of full-tensor gravity-gradiometer data

Denoising of full-tensor gravity-gradiometer data involves detailed information from field sources, especially the data mixed with high-frequency random noise. We present a denoising method based on the translation-invariant wavelet with mixed thresholding and adaptive threshold to remove the random noise and retain the data details. The novel mixed thresholding approach is devised to filter the random noise based on the energy distribution of the wavelet coefficients corresponding to the signal and random noise. The translation- invariant wavelet suppresses pseudo-Gibbs phenomena, and the mixed thresholding better separates the wavelet coefficients than traditional thresholding. Adaptive Bayesian threshold is used to process the wavelet coefficients according to the specific characteristics of the wavelet coefficients at each decomposition scale. A two-dimensional discrete wavelet transform is used to denoise gridded data for better computational efficiency. The results of denoising model and real data suggest that compared with Gaussian regional filter, the proposed method suppresses the white Gaussian noise and preserves the high-frequency information in gravity-gradiometer data. Satisfactory denoising is achieved with the translation-invariant wavelet.

3D inversion modeling of joint gravity and magnetic data based on a sinusoidal correlation constraint

Joint inversion based on a correlation constraint utilizes a linear correlation function as a structural constraint.The linear correlation function contains a denominator,which may result in a singularity as the objective function is optimized,leading to an unstable inversion calculation.To improve the robustness of this calculation,this paper proposes a new method in which a sinusoidal correlation function is employed as the structural constraint for joint inversion instead of the conventional linear correlation function.This structural constraint does not contain a denominator,thereby preventing a singularity.Compared with the joint inversion method based on a cross-gradient constraint,the joint inversion method based on a sinusoidal correlation constraint exhibits good performance.An application to actual data demonstrates that this method can process real data.

3D parallel inversion of time-domain airborne EM data

To improve the inversion accuracy of time-domain airborne electromagnetic data, we propose a parallel 3D inversion algorithm for airborne EM data based on the direct Gauss-Newton optimization. Forward modeling is performed in the frequency domain based on the scattered secondary electrical field. Then, the inverse Fourier transform and convolution of the transmitting waveform are used to calculate the EM responses and the sensitivity matrix in the time domain for arbitrary transmitting waves. To optimize the computational time and memory requirements, we use the EM ＂footprint＂ concept to reduce the model size and obtain the sparse sensitivity matrix. To improve the 3D inversion, we use the OpenMP library and parallel computing. We test the proposed 3D parallel inversion code using two synthetic datasets and a field dataset. The time-domain airborne EM inversion results suggest that the proposed algorithm is effective, efficient, and practical.

Passive multiple reverse time migration imaging based on wave decomposition and normalized imaging conditions

With the development of seismic exploration,passive-source seismic data has attracted increasing attention.Ambient noise passive seismic sources exists widely in nature and industrial production.Passive seismic data is important in logging while drilling(LWD),large-scale structural exploration,etc.In this paper,we proposed a passive multiple reverse time migration imaging(PMRTMI)method based on wavefield decomposition and normalized imaging conditions method.This method differs from seismic interferometry in that it can use raw passive seismic data directly in RTM imaging without reconstruction of virtual active gather,and we use the wavefield decomposition method to eliminate the low frequency noise in RTM.Further,the energy normalized imaging condition is used in full wavefield decomposition,which can not only enhance the image quality of both edge and deep information but also overcome the wrong energy problem caused by uneven distribution of passive sources;furthermore,this method exhibits high efficiency.Finally,numerical examples with the Marmousi model show the effectiveness of the method.

Three-dimensional forward modeling for the SBTEM method using an unstructured fi nite-element method

In this study,we propose a three-dimensional(3D)forward modeling algorithm of surface-to-borehole transient electromagnetic(SBTEM)fields based on an unstructured vector fi nite-element method to analyze the characteristics of SBTEM responses for complex geoelectrical models.To solve the double-curl diff usion equation for the electric fi eld,we use an unstructured tetrahedral mesh to discretize the model domain and select the unconditionally stable backward Euler scheme to discretize the time derivative.In our numerical experiments,we use a grounded wire as a transmitting source.After validating the algorithm’s eff ectiveness,we first analyze the diffusion characteristics and detectability of the electromagnetic field.After that,we focus our attention on the distribution and the cause of zero bands for Ex and dBy/dt components with the hope of guiding future field surveys.Finally,by simulating diff erent models,we analyze the capability of the SBTEM method in detecting typical mineral veins so that we can provide a reference for mineral resource exploration in the deep earth.

Seismic AVAZ inversion for orthorhombic shale reservoirs in the Longmaxi area, Sichuan

Seismic AVAZ inversion method based on an orthorhombic model can be used to invert anisotropy parameters of the Longmaxi shale gas reservoir in the Sichuan Basin..As traditional seismic inversion workfl ow does not suffi ciently consider the infl uence of fracture orientation,we predict fracture orientation using the method based on the Fourier series to correct pre-stacked azimuth gathers to guarantee the accuracy of input data,and then conduct seismic AVAZ inversion based on the VTI constraints and Bayesian framework to predict anisotropy parameters of the shale gas reservoir in the study area.We further analyze the rock physical relation between anisotropy parameters and fracture compliance and mineral content for quantitative interpretation of seismic inversion results.Research results reveal that the inverted anisotropy parameters are related to P-and S-wave respectively,and thus can be used to distinguish the effect of fracture and fl uids by the joint interpretation.Meanwhile high values of anisotropy parameters correspond to high values of fracture compliance,so the anisotropy parameters can refl ect the development of fractures in reservoir.There is two sets of data from different sources,including the content of brittle mineral quartz obtained from well data and the anisotropy parameters inverted from seismic data,also show the positive correlation.This further indicates high content of brittle mineral makes fractures developing in shale reservoir and enhances seismic anisotropy of the shale reservoir.The inversion results demonstrate the characterization of fractures and brittleness for the Longmaxi shale gas reservoir in the Sichuan Basin.

A fast-imaging method for airborne gravity gradient data based on tensor invariants

Airborne gravity gradient data contain additional short-wavelength information about the buried geological bodies.This study develops a fast interpretation method based on the gravity gradient data for the sources’spatial location and physical property parameters.This study analyzes the advantages of the source parameter inversion method based on tensor invariants.It proposes a normalized fast-imaging method based on tensor invariants to quickly estimate the spatial location parameters of sources through the local maximum value position of the imaging results.First,the tensor invariant characteristics and the imaging method’s effect in a simple model are analyzed using a theoretical model.Second,to analyze the imaging method’s application effect in complex model conditions,the method’s applicability is quantitatively analyzed using the data added with noise,superimposed anomalies of adjacent sources,and anomalies of deep and shallow geological bodies.The theoretical model’s simulation results show that the model’s imaging results in this study have satisfactory performance on the spatial position estimation of the sources.Finally,the method is applied to the gravity anomaly data corresponding to the Humble salt dome.The imaging results can effectively estimate the distribution of the salt dome’s horizontal and depths,verifying the practicability of the method.

Continuous time-varying Q-factor estimation method in the time-frequency domain

The Q-factor is an important physical parameter for characterizing the absorption and attenuation of seismic waves propagating in underground media,which is of great signifi cance for improving the resolution of seismic data,oil and gas detection,and reservoir description.In this paper,the local centroid frequency is defi ned using shaping regularization and used to estimate the Q values of the formation.We propose a continuous time-varying Q-estimation method in the time-frequency domain according to the local centroid frequency,namely,the local centroid frequency shift(LCFS)method.This method can reasonably reduce the calculation error caused by the low accuracy of the time picking of the target formation in the traditional methods.The theoretical and real seismic data processing results show that the time-varying Q values can be accurately estimated using the LCFS method.Compared with the traditional Q-estimation methods,this method does not need to extract the top and bottom interfaces of the target formation;it can also obtain relatively reasonable Q values when there is no eff ective frequency spectrum information.Simultaneously,a reasonable inverse Q fi ltering result can be obtained using the continuous time-varying Q values.

Magnetic interface forward and inversion method based on Padé approximation

The magnetic interface forward and inversion method is realized using the Taylor series expansion to linearize the Fourier transform of the exponential function. With a large expansion step and unbounded neighborhood, the Taylor series is not convergent, and therefore, this paper presents the magnetic interface forward and inversion method based on Pade approximation instead of the Taylor series expansion. Compared with the Taylor series, Pade＇s expansion＇s convergence is more stable and its approximation more accurate. Model tests show the validity of the magnetic forward modeling and inversion of Pade approximation proposed in the paper, and when this inversion method is applied to the measured data of the Matagami area in Canada, a stable and reasonable distribution of underground interface is obtained.

Three-dimensional gravity and magnetic data acquisition and study on its joint gradient Euler deconvolution method

Three-dimensional(3 D)gravitational and magnetic exploration is performed using aerial measurement tools,however,this has difficulties with measuring-height design and the construction of a joint-interpretation scheme.At present,the height in such experiments is set according to the measurement scale,and the distribution characteristics of anomalies are not fully considered.Here,we present the idea of using the attenuation characteristics of a singular-value spectrum to evaluate the contributions of various measurement heights and multi-height combinations for inversion to correctly and reasonably design appropriate measuring heights and the number of various measurement heights to be set.The jointgradient Euler-deconvolution method can accurately obtain the distribution of geological bodies from 3 D gravitational and magnetic data at an improved resolution,and experimental tests confirm these findings.Therefore,an actual 3 D aeromagnetic-data-acquisition and inversion test were carried out in the vicinity of the Zhurihe Iron Mine in Inner Mongolia.The fl ight-height diff erence was set to 60 m,and the specifi c distribution of lodes was obtained by the joint-gradient Euler-deconvolution method.This provides a reliable basis for future detailed exploration and proves that the methods presented in this paper have good practicalapplication eff ects and prospects.