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.

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.

Geochronological and geochemical constraints on the granitic gneiss in the Huozhou Complex: implications for the tectonic evolution of the Trans-North China Orogen

The Trans-North China Orogen is a major Neoarchean Paleoproterozoic collisional orogenic belt above the North China Craton, formed due to prolonged and complex processes. Even though many NeoarcheanPaleoproterozoic magmatic and metamorphic activities have been reported, due to the Huozhou Complex’s small outcropping range, little attention has been paid to the origin of various igneous rocks of the Huozhou Complex in the center of the Trans-North China Orogen. The Huozhou Complex, located south of the Luè liang, Wutai, and Hengshan complexes, is an important window into the Early Precambrian structure and evolution of the North China Craton. Its magma and metamorphism are crucial to understanding the development of the structural evolution of the Trans-North China Orogen. The Huozhou metamorphic complex area exposes a range of Precambrian metamorphic rocks, among which the most extensively dispersed is felsic biotite plagioclase gneiss. In this study comprehensive geological field survey, micropetrology,chronology, geochemistry, and Hf isotope analysis were carried out for the Qinggangping and Anziping gneiss in the north of the Huozhou Complex. The results show that the magmatic zircon age of the Qinggangping gneiss is2196 ± 14 Ma, and its protolith is I-type granite, formed by partial melting of igneous rocks in the absence of weathering. Its source is mainly the juvenile crust from depleted mantle dating 2431–2719 Ma, with a small amount of mantle-derived material. The Anziping gneiss has a metamorphic zircon age of 1931 ± 13 Ma with an S-type granite protolith belonging to peraluminous granite.The Anziping gneiss is formed by recycling pre-existing crustal components at 2613–2848 Ma. A minor quantity of mantle-derived magma is also introduced to the crust simultaneously. The samples of Qinggangping gneiss and Anziping gneiss show the characteristics of obvious negative Nb, Ti, and P elements in the spider diagram of primitive mantle standardization. This implies that the rocks have the characteristics of magmatic rocks in an island arc or subduction environment, which could have formed in the tectonic environment of the continental margin arc.

Evaluation and Driving Force Analysis of Cultivated Land Quality in Black Soil Region of Northeast China

Cultivated land is an important natural resource to ensure food,ecological and economic security.The cultivated land quality evaluation(CQE)is greatly significant for protecting and managing cultivated land.In this study,320 counties in the black soil region of Northeast China(BSRNC)represent the research units used to construct the CQE system measuring the soil properties(SP),cultivated land productivity(CLP),ecological environment(EE)and social economy(SE).The total of 19 factors were selected to calculate the integrated fertility index(IFI)and divided into grades.Simultaneously,we used the coupling coordination degree model to comprehensively analyze the spatial pattern of the cultivated land quality(CLQ)in the BSRNC,and use the structural equation model(SEM)to analyze the driving mechanism.The results show the following:1)The CLQ of 262 counties in the BSRNC is in a state of coupling and coordination,and the coupling and coordination degree presents a spatial distribution pattern of‘high in the southwest and northeast,low in the northwest and southeast’.The coordinated development degree of 271 counties is between 0.4 and 0.6,which is in a transitional state between coordination and disorder.2)The CLQ in the BSRNC is generally good,with an average grade of 3.High-quality cultivated land accounts for 58.45%of all counties,middle-and upper-quality cultivated land accounts for 27.05%,and poor-quality cultivated land accounts for 14.49%.3)The SEM analysis shows that the SP,CLP,EE,and SE all influence the CLQ.Among them,the SP has the largest driving force on the CLQ,while the SE has the smallest driving force on it.The results confirm that the main factors affecting the evaluation results are crop productivity level,normalized difference vegetation index,ratio vegetation index,difference vegetation index,and organic carbon content.When implementing protection measures in counties with a low CLQ,considering a balanced coordination of multiple systems and reasonably controlling the quality degradation are important.This study provides the current situation and driving factors of the CLQ in the BSRNC and will play an important role in black soil governance and utilization.

Parameters of new three-water model based on nuclear magnetic experiment and optimization algorithm

Clastic rock reservoir is the main reservoir type in the oil and gas field.Archie formula or various conductive models developed on the basis of Archie’s formula are usually used to interpret this kind of reservoir,and the three-water model is widely used as well.However,there are many parameters in the threewater model,and some of them are difficult to determine.Most of the determination methods are based on the statistics of large amount of experimental data.In this study,the authors determine the value of the parameters of the new three-water model based on the nuclear magnetic data and the genetic optimization algorithm.The relative error between the resistivity calculated based on these parameters and the resistivity measured experimentally at 100%water content is 0.9024.The method studied in this paper can be easily applied without much experimental data.It can provide reference for other regions to determine the parameters of the new three-water model.

Analysis on advanced transient EM detectability of coal mine roadway

The hidden water-bearing structures near the roadway tunnelling face are very likely to cause water seepage accidents in coal mines.Currently,transient electromagnetic(EM)technology has be-come an important method to detect water damage in advance of roadway excavation.In this paper,the time-domain finite element algorithm based on unstructured tetrahedron grids is used to accurate-ly simulate the geological body in front of the roadway excavation face and analyze its response.The authors detect the distance between the roadway excavation face and the low-resistivity water-bearing body,the resistivity difference between the low-resistivity body and surrounding rock,and the influence of the size of the low-resistivity body on the transient EM response.Furthermore,the common types of low-resistivity bodies in the roadway drivage process are used for modeling to analyze the attenuation of the detected EM response when there are low-resistivity bodies in front of the roadway.The research in this paper can help effectively detecting the water-bearing low-resistivity body in front of the roadway drivage and lay a foundation for reducing the risk of water seepage accidents.

Determination of multi-component content and construction of digital cores based on CT grey thresholds of altered igneous rocks

Altered igneous reservoirs have low porosity and permeability,compact structure and certain heterogeneity.A simple digital core with certain generality and multi-parameter constraints can be con-structed to characterize the microscopic pore structure and mineral composition.In this paper,based on core X-ray,CT images and whole-rock mineral analysis,threshold segmentation of mass content and grayscale distribution of various minerals in different lithologies of igneous rocks in the buried hill of Huizhou depression is carried out to construct digital core of altered igneous rocks.The results show that after converting the mineral mass content into volume content,the minerals of altered igneous rocks in Huizhou depression can be classified into components.According to the range of grayscale value,components can be divided into six parts.Due to the difference of the content of components in different lithologies of igneous rocks,differentiated grayscale threshold segmentation is needed to obtain the digital core for a single lithology.The final digital core generation process includes two steps:building a single component digital core,and stacking and combining.This kind of universal digital core model can support the subsequent pore scale numerical simulation and comprehensive rock physics research.

Classification of vegetative types in Changbai Mountain based on optical and microwave remote sensing data

Highly accurate vegetative type distribution information is of great significance for forestry resource monitoring and management.In order to improve the classification accuracy of forest types,Sentinel-1 and 2 data of Changbai Mountain protection development zone were selected,and combined with DEM to construct a multi-featured random forest type classification model incorporating fusing intensity,texture,spectral,vegetation index and topography information and using random forest Gini index(GI)for optimization.The overall accuracy of classification was 94.60%and the Kappa coefficient was 0.933.Comparing the classification results before and after feature optimization,it shows that feature optimization has a greater impact on the classification accuracy.Comparing the classification results of random forest,maximum likelihood method and CART decision tree under the same conditions,it shows that the random forest has a higher performance and can be applied to forestry research work such as forest resource survey and monitoring.

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.