核磁共振测井的正则化-启发式阈值降噪研究

核磁共振(NMR)测井是在高温、高压的极端条件下进行,样品体积小,回波幅度微弱.受到来源复杂的噪声影响,NMR信号通常被淹没在噪声中,测量数据的信噪比较低.提出正则化-启发式阈值算法(Regularization Heursure, R-Heursure)对小波分解后的细节系数阈值降噪,正则化因子的选取与地层孔隙结构和测量数据的原始信噪比相关.通过最大相关系数能量准则选取最优化的母小波函数、消失矩和分解层次,采用正则化因子约束估计的阈值,使选取的阈值恰好能大于噪声水平而不损失小孔(或微孔)的响应信息.数值模拟和实际测井资料处理验证了R-Heursure算法的降噪性能,NMR测井数据的信噪比得以有效改善,为储层评价提供更准确的信息.     

水平层状介质中双侧向测井资料的迭代Tikhonov正则化反演

根据非线性反演理论与Morozov偏差原理研究建立从双侧向测井(DLL)资料中同时重构地层原状电阻率、侵入带电阻率、侵入半径、层界面位置以及井眼泥浆电阻率的迭代正则化算法.首先利用Tikhonov正则化反演理论将双侧向测井资料的反演问题转化为含有稳定泛函的非线性目标函数的极小化问题,并利用Gauss-Newton算法确定极小化解.为得到稳定的反演结果并有效实现测井资料的最佳拟合,在迭代过程中将Morozov偏差原理和Cholesky分解技术相结合,建立了一套后验选择正则化因子的方法.最后通过理论模型和大庆油田实际测井资料的处理结果,验证了该算法能够取得更为满意的反演效果.     

低信噪比核磁共振T2谱反演算法研究（英文）

正则化因子的选择方法决定了正则化方法的稳定性和准确度。在分析改进的奇异值分解法与正则化方法的关系的基础上,给出一个正则因子计算公式。数值模拟试验表明改进的奇异值分解法和正则化方法适应低信噪比。当信噪比低于30时正则化方法比奇异值分解法结果更好,当信噪比高于30时奇异值分解法更优。采用本文提出的正则因子的正则化方法可以适用于实际核磁测井的低信噪比(5相似文献   

感应-横向组合测井联合反演及多解性研究

电阻率测井反演问题解的局部收敛性和非唯一性增加了测井解释的难度，影响了实际的应用，把测量原理迥异的感应型和电流型测井结果联合起来反演，是减小这些影响的有效方法，利用经测井理论和反演技术，建立了适用感应－横向组合测井响应的联合演算法，研究了单一的横向测井反演的感应－横向联合反演的局部收敛性、非唯一性及解的偏差，结果表明，联合反演更有利于同时确定地层电阻率、侵入半径和侵入带电阻率，当初值地较大范围内变化时，反演结果均能在真值或真值附近收敛，改善了反演结果的局部收敛性，增加了反演结果的可信度，避免了差异较大乃至相互矛盾的解，现场实例证实了这一结论。     

垂直井眼中三维阵列感应测井资料自适应井眼校正

三维感应测井受井眼环境影响很大,如何消除这些不利影响是有效获取各向异性地层电导率真值并提高资料处理和解释质量的一项重要工作.本文以一维柱状各向异性地层模型为基础,通过最小平方反演技术从实际测量资料中快速反演井径、偏心距、泥浆电导率、地层水平电导率和各向异性系数等模型参数,设法减小井眼环境对测量结果的影响.首先,利用仪器偏心情况下不同旋转角电导率张量的关系方程,建立旋转角提取与校正方法,获取无旋转角三维感应测井数据.在此基础上,根据事先计算好的井眼校正库并结合多维非线性有限元逼近技术,建立快速计算三维感应测井响应以及Jacobi导数矩阵的插值公式,并利用自适应正则化迭代反演技术不断修改模型参数,实现理论合成资料与输入资料的最佳拟合.最后,通过反演出的模型参数计算三维感应测井资料所有分量的校正量,实现三维感应资料的井眼校正.理论模型和实际资料的处理结果均证实,在低阻井眼泥浆情况下,一维柱状模型中的反演方法仍然可以提取出较可靠的地层水平和垂直电导率,且井眼环境校正对于受井眼环境影响相对较小的测井曲线具有非常好的校正效果.     

估计辐射参数的井间电磁波层析成像技术

为了避免使用不合理初始辐射场强和方向性因子带来的误差,研究了估计辐射参数的井间电磁波层析成像技术.通过时域有限差分法模拟表明,天线长度与波长的比值、钻孔充填情况、钻孔周围介质的物性均会影响偶极天线的初始辐射场强或方向性因子;为此结合已知的分层资料,将它们设为未知参数,并设定初始辐射场强与发射点位置相关,方向性因子随射线角度而变化;采用正则化反演方法,由钻孔资料建立了模型方差目标函数,使得反演结果与钻遇的地质特征保持一致.通过理论模型试验和实例应用分析表明,相对于传统射线层析成像方法,估计辐射参数的正则化层析成像技术有助于提高反演的准确性.     

基于水库气枪震源数据的频域水准比例因子反褶积方法分析

基于2016年福建南一水库气枪实验资料,利用频率域水准因子反褶积方法和互相关时延检测技术,研究水准比例因子及台站的背景噪声水平对反褶积计算结果的影响,并对反褶积法消除不同枪压引起的震源影响效果进行分析。结果表明:(1)水准比例因子越小,反褶积计算后的波形信噪比越小,当水准比例因子取值过小时,走时识别误差较大;(2)台站气枪记录的信噪比越大,走时识别精度越高,参考台记录的信噪比对结果的影响远大于远场台,当远场台记录的信噪比大于10时,走时误差一般小于6ms,而当参考台记录的信噪比为30左右时,走时误差可能达到20ms;(3)气枪震源差异较大时,频域水准反褶积方法去除震源效应的效果较好,而在震源差异较小时,可能引入方法误差。     

基于Backus-Gilbert理论的孔隙介质核磁共振弛豫反演

孔隙介质核磁共振(NMR)弛豫数据的多指数反演在NMR测井和岩心分析中均非常重要.本文基于Backus-Gilbert(BG)理论,提出一种NMR弛豫多指数反演的新方法.从解的非唯一性出发,不仅构造出一种解估计,更重要的是评价各种可接受的解估计,通过引入解估计分辨率和解估计方差对解进行评价,找出最佳折中解.通过算例,比较了新方法与传统正则化方法的效果,结果表明新方法具有明显优越性,在低信噪比条件下解稳定.最后分析了新方法的影响因素.     

沃尔什函数法测井曲线高分辨率反演

沃尔什函数法首先将地层物理量展开成沃尔什级数，其次利用线性模型参数估计法由测井曲线和响应方程估计出展开系数，进而求出地层物理量的真值。本文介绍了该方法的原理、参数估计和误差控制的方法，并对深感应测井的模型数据和实际资料进行了试算，获得了满意的效果。     

水平层状介质中双侧向资料的全参数正则化迭代反演与应用

水平层状介质中的模型向量通常由地层原状电阻率、侵人带电阻率、侵入半径以及地层界面深度等参数组成,整个测井记录实质上是模型向量的非线性函数.本文将研究双侧向仪器Born近似表达式并建立用双侧向测井资料同时反演所有模型的全参数正则化迭代算法(RWPII).双侧向电极系上的电流和电位必须满足特定的聚焦条件,其正演模拟与Born近似表达式与非聚集型仪器完全不同.应用摄动原理,首先研究双侧向仪器响应的Born近似,给出电极系上的电流密度、电流和电位微小变化与电导率摄动之间的积分关系.然后,根据层状地层中电阻率分布特征以及半解析正演模拟技术,建立Fréchet导数矩阵的快速算法,该Fréchet导数矩阵包含了视电阻率对所有模型参数的偏导数.在此基础上,研究同时反演所有模型参数的正则化迭代算法,并给出正则化因子的具体计算方法.最后,应用理论模型和实际资料反演结果验证该反演算法的有效性.     

多分辨率地震信号反褶积

基于二进小波变换提出了一种新的反褶积方法─-多分辨率地震信号反褶积．在地震信号二进小波变换域中的各尺度上分别进行其分辨率随小波尺度变化的反褶积,利用不同分辨率反褶积结果之间的相关性,以及测量噪声随尺度的衰减特性,从低分辨率反褶积结果逼近高分辨率反褶积结果．理论分析和实验表明,该方法有较高的精度,并且在较低信噪比情况下有好的效果．     

Analysis of the Frequency Domain Water-level Deconvolution Method Based on Reservoir Airgun Source Data

Based on the 2016 airgun experimental data of the Fujian Nanyi reservoir, we adopted the frequency domain water-level deconvolution method and cross-correlation time delay detection technique to study the influence of level scaling factor and the background noise level of the station on deconvolution calculation results, and analyze the effect of deconvolution on eliminating the influence of the source caused by different air-gun pressures. The results show that:(1) When the level scaling factor is smaller, the signal to noise ratio of the waveform after the deconvolution is smaller, and when the level scaling factor is over smaller, the identification error of travel time is greater. (2) When the SNR of the station record is higher, the recognition accuracy of travel time is higher, the influence of SNR on the reference station record is far greater than the far station, when the SNR of the far station record is more than 10, the error of travel time is within 6ms, but when the SNR of the reference station record is 30, the travel time error may reach to 20ms. (3) When the airgun source difference is big, the frequency domain water-level deconvolution method has better effect on eliminating the source influence, but the method error may be introduced when the source difference is small.     

Vibroseis deconvolution: A comparison of pre and post correlation vibroseis deconvolution data in real noisy data

Vibroseis is a source used commonly for inland seismic exploration. This non-destructive source is often used in urban areas with strong environmental noise. The main goal of seismic data processing is to increase the signal/noise ratio where a determinant step is deconvolution. Vibroseis seismic data do not meet the basic minimum-phase assumption for the application of spiking and predictive deconvolution, therefore various techniques, such as phase shift, are applied to the data, to be able to successfully perform deconvolution of vibroseis data.This work analyzes the application of deconvolution techniques before and after cross-correlation on a real data set acquired for high resolution prospection of deep aquifers. In particular, we compare pre-correlation spiking and predictive deconvolution with Wiener filtering and with post-correlation time variant spectral whitening deconvolution. The main result is that at small offsets, post cross-correlation spectral whitening deconvolution and pre-correlation spiking deconvolution yield comparable results, while for large offsets the best result is obtained by applying a pre-cross-correlation predictive deconvolution.     

基于超级混合虚折射干涉法的阵列声波处理方法及其在套管井中的应用

在声学测井勘探时,复杂多变的测井条件(如胶结不好的套管井和随钻测井等)下的强直达干扰信号会导致地层信号的信噪比降低,使得地层有用信号的提取十分困难.本文在超级虚折射法的基础上,利用声源与接收器的互易原理,对原始波场信号进行互相关和褶积重建,并对所有可能的位形产生的信号进行叠加,将该叠加波场与超级虚折射波场进行叠加,重构出超级混合虚折射波场,使不同接收器上的折射波信噪比在超级虚折射的基础上进一步得到提高,显著提高了有用信号的信噪比,同时解决了超级虚折射近场能量不足的问题.套管井测井的理论模拟及实际处理表明,该方法可以很好地压制套管波信号、提高地层信号的信噪比.     

浅层地震勘探在沿海地区隐伏断层探测中的应用

中国沿海地区的第四系厚度变化较大,在大部分地区存在淤泥、黏土和砂层的多重互层情况,这种地层结构对地震波的传播十分不利,容易造成对高频信号能量的强烈耗散。沿海地区又具有人口密度大、工业干扰强烈的特点,地震波的激发、接收条件较差。文中给出了在粤东、天津和唐山覆盖层厚度变化较大的地区,利用浅层地震勘探开展隐伏断层探测的应用实例。较系统地介绍了野外工作、资料处理和解释方法。阐述了在强干扰环境下提高信噪比和分辨率的有效技术途径。所给实例中的反射波组能量强,隐伏断层的反应明显,经跨断点钻探验证,结果证实地震勘探确定的界面深度、断点位置和基岩错距等参数都是可靠的。所述技术方法对类似地区的城市活断层探测和工程物探工作都有较重要的参考价值     

Widely linear denoising of multicomponent seismic data

Seismic data processing is a challenging task, especially when dealing with vector-valued datasets. These data are characterized by correlated components, where different levels of uncorrelated random noise corrupt each one of the components. Mitigating such noise while preserving the signal of interest is a primary goal in the seismic-processing workflow. The frequency-space deconvolution is a well-known linear prediction technique, which is commonly used for random noise suppression. This paper represents vector-field seismic data through quaternion arrays and shows how to mitigate random noise by proposing the extension of the frequency-space deconvolution to its hypercomplex version, the quaternion frequency-space deconvolution. It also shows how a widely linear prediction model exploits the correlation between data components of improper signals. The widely linear scheme, named widely-linear quaternion frequency-space deconvolution, produces longer prediction filters, which have enhanced signal preservation capabilities shown through synthetic and field vector-valued data examples.     

THE PHASE PROPERTY OF A FINITE REALIZATION OF RANDOM NOISE AND ITS INFLUENCE ON DECONVOLUTION *

A finite realization of a discrete random noise process may be considered as a one-sided energy signal. Its phase property can then be described by means of the center position. The samples of such a realization are the components of a random signal vector and the center position is therefore a random variable. A statistical analysis shows that the expected value of the center position equals half the time duration of the realization. This implies that the Z-transform of the realization may be expected to have an equal number of poles and zeros inside and outside the unit circle. The standard deviation from the expected value of the center position is shown to depend on the time duration of the realization and on the autocorrelation of the process. It follows that, for processes that can be described by the convolution of a white series and a disturbance wavelet, the center position is independent of the phase property of the wavelet. A conclusion based on these results is that the homomorphic technique of wavelet estimation through cepstrum stacking must give questionable outcomes. Another conclusion is that the super-position of a realization of random noise on a minimum phase wavelet will in general give a mixed phase resulting signal. It is pointed out that schemes for the derivation of deconvolution filters do not take account of this phenomenon.     

用遗传算法实现地震信号反褶积

遗传算法作为寻优手段具有全局优化和很好的稳定性．本文将遗传算法用于地震信号反褶积处理，与已往方法相比它具有更好的分辨率和稳定性我们采用Bernoulli－Gaussian模型和ARMA模型分别描述地震反射系数序列和地震子波，用最大似然和最小预测误差准则分别构造用于估计反射系数序列和地震子波的目标函数，用遗传算法优化目标函数，以实现地震信号反褶积．     

小波变换的过零点特性与地震勘探信号的信噪比和分辨率

利用小波变换研究地震勘探信号小波变换的过零点特性，本文提出了用小波变换的过零点特性和地震勘探信号相邻道的横向相关性提高信号分辨率和信噪比的新方法．该方法包括两个主要步骤：①利用相邻地震道信号具有很好相关性，而噪音相关性差的特点以及小波变换的过零点特性得到有效反射波同相轴随空间坐标的变化信息．②利用奇异值分解和最小二乘（SVD-TLS）方法沿同相轴对振幅进行多项式拟合去噪并增加信号高频提高信号分辨率．     

Phase inversion deconvolution for long and short period multiples attenuationl1

A new approach to deconvolution has been developed to improve the attenuation of multiple energy. This approach to deconvolution is unique in that it not only eliminates the usual assumptions of a minimum phase lag wavelet and a random distribution of impulses, but also overcomes the noise limitation of the homomorphic deconvolution and its inherent instability to phase computation. We attempt to analyse the continuous alteration of the acoustic waveform during the propagation through a linear system. Based on the results of this analysis, the surface-related measurements are described as a convolution of the impulse response of the system with the non-stationary forward wavelet which includes all multiple terms generated within the system. The amplitude spectrum of the forward wavelet is recovered from the amplitude spectrum of the recorded signal, using the difference between the rate of decay of the source wavelet and the duration of the measurement. The phase spectrum of the forward wavelet is estimated using the Hilbert transform and the fact that the mixed phase lag wavelet can be presented as a convolution of the minimum and maximum phase lag wavelets. The multiples are discriminated from primaries by comparison of the phase spectrum of the seismic signal and the inverse of the forward wavelet. Therefore, the technique is called phase inversion deconvolution (PID). This approach requires no velocity information in order to recognize and attenuate multiple energy. Therefore, primary energy is recovered in the near-offset region where the velocity differential between primary and multiple energies is very small.