基于均方差比归一化垂向梯度法的位场边界检测

针对位场数据处理中边界定位精度和识别能力的问题,提出一种新的边界检测方法——均方差比归一化垂向梯度法。均方差比(MSER)可进行边界识别,是针对边界点异常方向性和均方差衡量数据波动性提出的,对全区数据点4个方向的均方差进行归一化后选择各个数据点均方差比的最大值作为滤波输出实现的;均方差比归一化垂向梯度(NVD-MSER)可进行边界增强,是通过均方差比的垂向梯度及其总梯度的比值实现的。模型试验对比分析结果表明,NVD-MSER方法具有计算稳定性强、反映的边界位置连续性好、与实际模型边界对比偏差小的优点,这说明NVD-MSER法有较强的边界检测能力。用NVD-MSER法可以检测出黑龙江虎林盆地19条断裂,而欧拉反褶积只能识别出11条断裂,说明NVD-MSER法增强了对断裂平面位置的识别能力。

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.

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.

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.

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.