Study of formation boundary and dip attribute extraction based on edge detection technology

In the seismic profile interpretation process,as the seismic data are big and the small geological features are difficult to identify,improvement of the efficiency is needed. In this study,structure tensor method in computer image edge detection processing is applied into the 2D seismic profile. Coherent attribute is used to extract formation edge. At the same time,extracting the eigenvalues and eigenvectors to calculate the seismic geometric properties which include dip and apparent dip,automatic identification is achieved. Testing the Gaussian kernel function with synthetic models and comparing the coherent attribute and dip attribute extraction results before and after,the conclusion that Gaussian filter can remove the random noise is obtained.

Implementation technology of Kirchhoff demigration

Demigration refers to directly applying a specific imaging technique to a migrated section. It is applied primarily to seismic data mapping. In a previous research study, a time-efficient implementation technology of demigration was expounded. In the present study, the Fast Marching Method （FMM） used for traveltime computation in the isochrone-staek demigration, is developed. Furthermore, other key techniques （ such as selection of aperture and antialiasing filtering factor） are analyzed in detail. Besides, the detail implementation method and program flow are given, which is shown their good computational efficiency and high-quality demi- gration effect. This implementation technique is illustrated with both the V（z） model and Marmousi model. It provides a basic method for implementing demigration in the application of seismic data mapping.

Application of spectral decomposition technology in shallow gas detection based on Wigner-Ville distribution

Shallow gas is composed of all kinds of shallow buried natural gas resources( < 1500 m) with relatively small reserve for each gas resource. It has some advantages such as shallow burial depth,good physical properties and the huge accumulations. Based on the Wigner-Vill distribution,a general spectral decomposition method is applied in the shallow gas detection. Cone-shaped kernel function filtering method is used to suppress cross-terms of the Wigner-Ville distribution,which is tested on field seismic data. Because of shallow gas reservoir has a characteristic that low frequency energy is stronger and high frequency energy is weaker,it indicates the presence of shallow gas successfully.

Research on compound plugging removal technology and its application in Xinmu oilfield of Jilin

Most fault-block reservoirs in Xinmu oilfield belong to heterogeneous sandstone in characters which has low permeability, and reservoir pollution is a common phenomenon in this area. Acidizing deplugging in oil wells has become one of the major measures to improve production efficiency in the field. A compound deplug- ging technology in high efficiency low corrosion is developed for this kind of low permeability sandstone reser- voir. On the basis of profoundly understanding of the reservoir＇s physical properties and sensitivity, along with comprehensive analysis of the cause for jams in oil wells, a series of experiments are carried out in order to in- vestigate the dissolution reaction among samples and deplugging inhibitor, sample dissolving speed, formation fluid compatibility, reduction of secondary pollution, etc. Considered reservoir condition in nearby wells the op- timized compositional deplugging liquid formula is selected for this reservoir. It is featured by reducing the de- plugging reaction speed, extending solution for processing radius, preventing secondary damage in dissolution processing, and removing plug pollution effectively. To implement this high efficiency low corrosion deplugging measure based on reservoir condition in the borehole and nearby wells, a relatively better result of deplugging and production increasing is achieved, which enriches the measures to increase production in Xinmu oilfield, and this method can be applied to other similar oiffields for the purpose of maintaining the crude oil production and providing assistant for increasing the production significantly.

A method to interpret fracture aperture of rock slope using adaptive shape and unmanned aerial vehicle multi-angle nap-of-the-object photogrammetry

The aperture of natural rock fractures significantly affects the deformation and strength properties of rock masses,as well as the hydrodynamic properties of fractured rock masses.The conventional measurement methods are inadequate for collecting data on high-steep rock slopes in complex mountainous regions.This study establishes a high-resolution three-dimensional model of a rock slope using unmanned aerial vehicle(UAV)multi-angle nap-of-the-object photogrammetry to obtain edge feature points of fractures.Fracture opening morphology is characterized using coordinate projection and transformation.Fracture central axis is determined using vertical measuring lines,allowing for the interpretation of aperture of adaptive fracture shape.The feasibility and reliability of the new method are verified at a construction site of a railway in southeast Tibet,China.The study shows that the fracture aperture has a significant interval effect and size effect.The optimal sampling length for fractures is approximately 0.5e1 m,and the optimal aperture interpretation results can be achieved when the measuring line spacing is 1%of the sampling length.Tensile fractures in the study area generally have larger apertures than shear fractures,and their tendency to increase with slope height is also greater than that of shear fractures.The aperture of tensile fractures is generally positively correlated with their trace length,while the correlation between the aperture of shear fractures and their trace length appears to be weak.Fractures of different orientations exhibit certain differences in their distribution of aperture,but generally follow the forms of normal,log-normal,and gamma distributions.This study provides essential data support for rock and slope stability evaluation,which is of significant practical importance.

Multi-scale data joint inversion of minerals and porosity in altered igneous reservoirs—A case study in the South China Sea

There are abundant igneous gas reservoirs in the South China Sea with significant value of research,and lithology classification,mineral analysis and porosity inversion are important links in reservoir evaluation.However,affected by the diverse lithology,complicated mineral and widespread alteration,conventional logging lithology classification and mineral inversion become considerably difficult.At the same time,owing to the limitation of the wireline log response equation,the quantity and accuracy of minerals can hardly meet the exploration requirements of igneous formations.To overcome those issues,this study takes the South China Sea as an example,and combines multi-scale data such as micro rock slices,petrophysical experiments,wireline log and element cutting log to establish a set of joint inversion methods for minerals and porosity of altered igneous rocks.Specifically,we define the lithology and mineral characteristics through core slices and mineral data,and establish an igneous multi-mineral volumetric model.Then we determine element cutting log correction method based on core element data,and combine wireline log and corrected element cutting log to perform the lithology classification and joint inversion of minerals and porosity.However,it is always difficult to determine the elemental eigenvalues of different minerals in inversion.This paper uses multiple linear regression methods to solve this problem.Finally,an integrated inversion technique for altered igneous formations was developed.The results show that the corrected element cutting log are in good agreement with the core element data,and the mineral and porosity results obtained from the joint inversion based on the wireline log and corrected element cutting log are also in good agreement with the core data from X-ray diffraction.The results demonstrate that the inversion technique is applicable and this study provides a new direction for the mineral inversion research of altered igneous formations.

Effect of gas saturation on P-wave velocity in tight sandstone

By measuring the variation of the P-and S-wave velocities of tight sandstone samples under water saturation,it was confirmed that with the decrease in water saturation,the P-wave velocity first decreased and then increased.The variation in velocity was influenced by the sandstone’s porosity.The commonly used Gassmann equation based on fluid substitution theory was studied.Comparing the calculated results with the measured data,it was found that the Gassmann equation agreed well with the measured data at high water saturation,but it could not explain the bending phenomenon of P-wave velocity at low saturation.This indicated that these equations could not accurately describe the relationship between fluid content and rock acoustic velocity.The reasons for this phenomenon were discussed through Taylor’s expansion.The coefficients of the fitting formula were calculated and verified by fitting the measured acoustic velocity changes of the cores.The relationship between P-wave velocity and saturation was discussed,which provides experimental support for calculating saturation using seismic and acoustic logging data.

3D time-domain forward modeling of airborne transient electromagnetism considering superparamagnetic effect

The superparamagnetic effect arises from the superparamagnetism exhibited by a multitude of nano-sized magnetic mineral particles under an external electric field.This phenomenon manifests in transient electromagnetic data primarily as a deceleration in the attenuation rate of late-stage signals,a characteristic difficult to discern directly from airborne transient electromagnetic signals,consequently leading to significant misinterpretations of subterranean electrical structures.This study embarks on 3D forward modeling of airborne electromagnetic responses in the frequency domain,accounting for the superparamagnetic effect,utilizing an unstructured finite element method.Superparamagnetic responses in the time domain were obtained through frequency-time conversion.This investigation explores the influence of various parameters-such as magnetic susceptibility,time constants,and flight altitude-on the superparamagnetic effect by examining the response characteristics of typical targets.Findings indicate that in its late stages,the superparamagnetic effect can induce a relative anomaly of up to 300%.There is a positive correlation between magnetic susceptibility and the strength of the superparamagnetic effect.The influence of the time constant's upper and lower limits on the superparamagnetic effect is minimal;however,the range between these limits significantly affects the effect,showing a negative correlation with its intensity.Higher flight altitudes weaken the superparamagnetic signal.The impact is most pronounced when superparamagnetic minerals are shallowly buried,effectively shielding the underlying geology with the characteristics of a good conductivity anomaly,but this effect diminishes with greater depth.The insights from this study provide a theoretical framework for a deeper understanding of the superparamagnetic effect in transient electromagnetic signals and for more accurate interpretations of subterranean geological and electrical structures.

UAV-based transient electromagnetic 3D forward modeling and inversion and analysis of exploration capability

Unmanned aerial vehicle transient electromagnetic(UAV-TEM)is a novel airborne exploration method that offers advantages such as low cost,simple operation,high exploration efficiency and suitability for near-surface exploration in complex terrain areas.To improve the accuracy of data interpretation in this method,the authors conducted a systematic three-dimensional(3D)forward modeling and inversion of the UAV-TEM.This study utilized the finite element method based on unstructured tetrahedral elements and employed the second-order backward Euler method for time discretization.This allowed for accurate 3D modeling and accounted for the effects of complex terrain.Based on these,the influence characteristics of flight altitudes and the sizes,burial depths,and resistivities of anomalies are compared and analyzed to explore the UAV-TEM systems’exploration capability.Lastly,four typical geoelectrical models of landslides are designed,and the inversion method based on the Gauss-Newton optimization method is used to image the landslide models and analyze the imaging effect of the UAV-TEM method on landslide geohazards.Numerical results showed that UAV-TEM could have better exploration resolution and fine imaging of nearsurface structures,providing important technical support for monitoring,early warning,and preventing landslides and other geological hazards.

Horizontal-to-vertical spectral ratio inversion method based on multimodal forest optimization algorithm

The exploration of urban underground spaces is of great significance to urban planning,geological disaster prevention,resource exploration and environmental monitoring.However,due to the existing of severe interferences,conventional seismic methods cannot adapt to the complex urban environment well.Since adopting the single-node data acquisition method and taking the seismic ambient noise as the signal,the microtremor horizontal-to-vertical spectral ratio(HVSR)method can effectively avoid the strong interference problems caused by the complex urban environment,which could obtain information such as S-wave velocity and thickness of underground formations by fitting the microtremor HVSR curve.Nevertheless,HVSR curve inversion is a multi-parameter curve fitting process.And conventional inversion methods can easily converge to the local minimum,which will directly affect the reliability of the inversion results.Thus,the authors propose a HVSR inversion method based on the multimodal forest optimization algorithm,which uses the efficient clustering technique and locates the global optimum quickly.Tests on synthetic data show that the inversion results of the proposed method are consistent with the forward model.Both the adaption and stability to the abnormal layer velocity model are demonstrated.The results of the real field data are also verified by the drilling information.

Geoid interpolation in Gongzhuling City and the 1cm quasi geoid refinement

The refinement of the 1cm geoid holds significant importance in engineering applications,where the accuracy of the geoid is frequently constrained by its resolution.However,there has been limited exploration into the specific relationship between geoid resolution and accuracy.This article aims to address this gap by thoroughly examining said relationship.This study employs the CapRCR modification to calculate the geoid of Gongzhuling City.The findings indicate that the accuracy can be enhanced by 2%to 9%through encryption of the geoid resolution from 2.5′×2.5′to 1.5′×1.5′.Furthermore,this improvement can be augmented by 15%to 21%through encryption of the gravity anomaly to the same resolution.The accuracy of the geoid exhibits a linear relationship with the resolution of gravity anomalies.The theoretical accuracy of the geoid,excluding integration errors,is determined to be 1.21 cm.In engineering applications,opting for a resolution of 30''×30''can result in the geoid containing only an integration error of 2 mm.However,to attain an accuracy of 1 cm,it becomes imperative to further mitigate data errors.Consequently,the final refined quasi-geoid accuracy is established at 0.56 mm.

3D magnetotelluric inversions with unstructured finite-element and limited-memory quasi-Newton methods

Traditional 3D Magnetotelluric(MT) forward modeling and inversions are mostly based on structured meshes that have limited accuracy when modeling undulating surfaces and arbitrary structures. By contrast, unstructured-grid-based methods can model complex underground structures with high accuracy and overcome the defects of traditional methods, such as the high computational cost for improving model accuracy and the difficulty of inverting with topography. In this paper, we used the limited-memory quasi-Newton(L-BFGS) method with an unstructured finite-element grid to perform 3D MT inversions. This method avoids explicitly calculating Hessian matrices, which greatly reduces the memory requirements. After the first iteration, the approximate inverse Hessian matrix well approximates the true one, and the Newton step(set to 1) can meet the sufficient descent condition. Only one calculation of the objective function and its gradient are needed for each iteration, which greatly improves its computational efficiency. This approach is well-suited for large-scale 3D MT inversions. We have tested our algorithm on data with and without topography, and the results matched the real models well. We can recommend performing inversions based on an unstructured finite-element method and the L-BFGS method for situations with topography and complex underground structures.

Three-dimensional inversion of controlled-source audio-frequency magnetotelluric data based on unstructured fi nite-element method

We propose a new 3D inversion scheme to invert the near-and transition-zone data of CSAMT with topography accurately.In this new method,the earth was discretized into unstructured tetrahedra to fit the ragged topography and the vector fi nite-element method was adopted to obtain precise responses and good sensitivity.To simulate the attitude and shape of the transmitter,we divided a long-grounded transmitter into dipoles and integrated these dipoles to obtain good responses in the near-and transition-fi eld zones.Next,we designed an L2 norm-based objective functional and applied a standard quasi-Newton method as the optimization method to solve the inverse problem and guarantee steady convergence.We tested our 3D inversion method first on synthetic data and then on a field dataset acquired from select sites near Changbai Mountain,China.In both tests,the new inversion algorithm achieved excellent fitting between the predicted and observed data,even in near-and transition-fi eld zones,and the inversion results agreed well with the true model.These fi ndings reveal that the proposed algorithm is eff ective for 3D inversion of CSAMT data.

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.

Phase velocity in 2D TTI media

We derive an expression for phase velocity in 2D tilted transverse isotropy （TTI） media. Snapshots of phase velocity in TTI and transverse isotropy （TI） model media are simulated and analyzed using the derived expression. In addition, the x-component character differences between the modeled phase velocities of the two media models are compared and analyzed.

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 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.

A rock physics model for the characterization of organic-rich shale from elastic properties

Kerogen content and kerogen porosity play a significant role in elastic properties of organic-rich shales. We construct a rock physics model for organic-rich shales to quantify the effect of kerogen content and kerogen porosity using the Kuster and Toks6z theory and the selfconsistent approximation method. Rock physics modeling results show that with the increase of kerogen content and kerogen-related porosity, the velocity and density of shales decrease, and the effect of kerogen porosity becomes more obvious only for higher kerogen content. We also find that the Poisson＇s ratio of the shale is not sensitive to kerogen porosity for the case of gas saturation. Finally, for the seismic reflection responses of an organic-rich shale layer, forward modeling results indicate the fifth type AVO re- sponses which correspond to a negative intercept and a positive gradient. The absolute values of intercept and gradient increase with kerogen content and kerogen porosity, and present predictable variations associated with velocities and density.

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.