频率域波形反演中与频率相关的影响因素分析

波动方程深度偏移是解决复杂地质体成像的关键技术,基于波动方程的速度建模为其提供更为精确的速度模型.频率域波形反演是目前研究最为广泛的波动方程速度建模方法之一,它推动了波形反演在勘探尺度下的应用.本文通过对频率域波形反演的实现,分析对比了其有效执行过程中与频率相关的影响因素.介绍了时间域的多尺度反演方法在频率域的一种实现方式,对比分析了输入数据的频点带宽和应用的子波频带范围不同时对反演结果的影响.本文通过设计的山地地质模型对频率域波形反演进行了测试和对比,得到的结论为频率域波形反演的有效计算提供了依据和参考.     

Blocky regularization schemes for Full‐Waveform Inversion★

With ill‐posed inverse problems such as Full‐Waveform Inversion, regularization schemes are needed to constrain the solution. Whereas many regularization schemes end up smoothing the model, an undesirable effect with FWI where high‐resolution maps are sought, blocky regularization does not: it identifies and preserves strong velocity contrasts leading to step‐like functions. These models might be needed for imaging with wave‐equation based techniques such as Reverse Time Migration or for reservoir characterization. Enforcing blockiness in the model space amounts to enforcing a sparse representation of discontinuities in the model. Sparseness can be obtained using the ?¹ norm or Cauchy function which are related to long‐tailed probability density functions. Detecting these discontinuities with vertical and horizontal gradient operators helps constraining the model in both directions. Blocky regularization can also help recovering higher wavenumbers that the data used for inversion would allow, thus helping controlling the cost of FWI. While the Cauchy function yields blockier models, both ?¹ and Cauchy attenuate illumination and inversion artifacts.     

基于聚类分析的多子波优化重定算法

提出了基于聚类分析的多子波优化重定算法,即利用聚类分析对波形进行分类,然后在分类的基础上对每一聚类的波形分别采用β参数拟合算法和门槛算法对其进行重定,最后以误差平方和最小为依据选择最优的作为最后的结果。经过基于聚类分析的多子波优化算法重定后的卫星测高数据与验潮站数据进行了比较,发现该算法能在一定程度上改善近海卫星测高数据的精度,提高重定成功率,降低重定次数。     

薄互层地震波形特征研究——以博兴洼陷沙四段滩坝砂为例

陆相断陷湖盆中薄互层非常普遍,如滩坝砂岩,具有横向变化快、单层厚度薄等特点,加之受地震资料分辨率低等因素的影响,致使其地震预测难度极大。文中在分析博兴洼陷沙四段滩坝砂岩发育区各种岩性地球物理特征的基础上,从正演模型入手,结合实际验证,分析了不同岩性及其组合的地震反射特征,找出了研究区砂体组合反射在地震波形结构、极大值两侧面积比等地震属性上的表现形式,总结了薄互层及其组合模式与地震反射各属性的对应关系,以此进行地震属性预测,取得良好效果,在东营凹陷西部沙四段滩坝砂岩储层的勘探中较准确地预测出储层的发育带。     

基于Born敏感核函数的速度、密度双参数全波形反演

速度、密度之间的相互耦合使得密度在多参数全波形反演中较难获得.本文将截断高斯-牛顿法用于声介质速度、密度双参数全波形反演,通过考虑近似Hessian矩阵中反映速度、密度相互作用的非主对角块元素,有效解决了多参数全波形反演中速度、密度之间的耦合问题,在不采用反演策略的情况下,仍能够获得精度较高的速度、密度反演结果.常规的截断牛顿类全波形反演通常利用一阶伴随状态法求取目标函数对模型参数的梯度,利用二阶伴随状态法或有限差分法求解Hessian-向量乘,在每一步内循环迭代过程中需要额外求解两次正演问题,计算量较大.本文基于Born近似,将梯度计算中的核函数-向量乘表示为具有明确物理意义的向量-标量乘的累加运算,同时将Hessian-向量乘转化为两次核函数-向量乘,无需额外求解正演问题,有效降低了计算量.数值实验证明了本文提出的方法的有效性.     

墨西哥湾GC955H井天然气水合物储层声波衰减特征

墨西哥湾GC955H井钻遇两种不同的水合物储层,浅层裂隙充填型水合物和深层砂岩型水合物.浅层水合物充填在细粒泥质沉积物的裂隙中,由电阻率测井计算的饱和度平均值为25%.深层水合物充填在砂岩孔隙中,由电阻率计算的饱和度平均为65%.基于声波全波形数据,本文计算了GC955H井储层的声波衰减大小,结果发现两种水合物层对声波衰减的影响不同.泥岩水合物层中的声波衰减与上覆背景泥岩沉积层基本相当.砂岩水合物层的声波衰减大于0.1,最大0.15,远高于上覆和下伏的背景砂岩层.对声波频率的分析发现,GC955H井泥岩层中水合物对声波频率的影响不大,砂岩水合物层的声波频谱与上下背景砂岩层相比发生变化,主频增大.通过对比两种水合物储层的特征,本文初步定性分析了造成水合物对储层声波传播影响不同的原因,包括岩性、水合物饱和度、水合物赋存方式等;但对此的定量描述需要未来更详尽地讨论和研究水合物的声波衰减机制.     

Modelling of rutting of two flexible pavements with the shakedown theory and the finite element method

This paper presents a finite element program, for the modelling of rutting of flexible pavements. In its present version, the program incorporates a permanent deformation model for unbound granular materials based on the concept of the shakedown theory developed by Zarka for metallic structures under cyclic loadings and has been used to estimate the permanent deformations of unbound granular materials (UGM) subjected to traffic loading. The calculation is performed in two steps: the first step consists in modelling the resilient behaviour of the pavement in 3D, using non-linear elastic models, to determine the stress field in the pavement. Then stress paths are derived and used to calculate the permanent deformations and the displacements, using a Drucker–Prager yield surface. An application to the prediction of the permanent deformations of experimental pavements with an unbound granular base, tested on the LCPC pavement testing facility is presented.     

井间断裂构造的计算机分析方法

井间断裂构造分析是根据测井资料,以地层分层数据为基础,利用计算机技术实现井间断裂分析的一种计算机自动分析方法.动态波形匹配算法能够很好地建立井间地层之间的对比关系,可以通过曲线拟合、特征提取、匹配代价计算等步骤,自动绘制井间地层对比的路径图解.通过对各种构造的路径图分析,总结出正断层、逆断层、不整合、同沉积断层、犁式断层和地层尖灭等34种情况下的路径图模式,提出路径图模式及相关概念.     

Retracking of SARAL/AltiKa Radar Altimetry Waveforms for Optimal Gravity Field Recovery

The accuracy of the marine gravity field derived from satellite altimetry depends on dense track spacing as well as high range precision. Here, we investigate the range precision that can be achieved using a new shorter wavelength Ka-band altimeter AltiKa aboard the SARAL spacecraft. We agree with a previous study that found that the range precision given in the SARAL/AltiKa Geophysical Data Records is more precise than that of Ku-band altimeter by a factor of two. Moreover, we show that two-pass retracking can further improve the range precision by a factor of 1.7 with respect to the 40 Hz-retracked data (item of range_40 hz) provided in the Geophysical Data Records. The important conclusion is that a dedicated Ka-band altimeter-mapping mission could substantially improve the global accuracy of the marine gravity field with complete coverage and a track spacing of <6 km achievable in ～1.3 years. This would reveal thousands of uncharted seamounts on the ocean floor as well as important tectonic features such as microplates and abyssal hill fabric.     

Pore structure characteristics of tight sandstones in the northern Songliao Basin,China

Understanding the pore structure characteristics of tight gas sandstones is the primary purpose of reservoir evaluation and efforts to characterize tight gas transport and storage mechanisms and their controls. Due to the various pore types and multi-scale pore sizes in tight reservoirs, it is essential to combine several techniques to characterize pore structure. Scanning electron microscopy (SEM), nitrogen gas adsorption (N₂GA), mercury intrusion porosimetry (MIP) and nuclear magnetic resonance (NMR) were conducted on tight sandstones from the Lower Cretaceous Shahezi Formation in the northern Songliao Basin to investigate pore structure characteristics systematically (e.g., type and size distribution of pores) and to establish how significant porosity and permeability are for different pore types. The studied tight sandstones are composed of intergranular pores, dissolution pores and intercrystalline pores. The integration of N₂GA and NMR can be used as an efficient method to uncover full pore size distribution (PSD) of tight sandstones, with pore sizes ranging from 2 nm to dozens of microns. The full PSDs indicate that the pore sizes of tight sandstones are primarily distributed within 1.0 μm. With an increase in porosity and permeability, pores with larger sizes contribute more to porosity. Intercrystalline pores and intergranular/dissolution pores can be clearly distinguished on the basis of mercury intrusion and surface fractal. The relative contribution of intercrystalline pores to porosity ranges from 58.43% to 91.74% with an average of 79.74%. The intercrystalline pores are the primary contributor to pore space, whereas intergranular/dissolution pores make a considerably greater contribution to permeability. A specific quantity of intergranular/dissolution pores is the key to producing high porosity and permeability in tight sandstone reservoirs. The new two permeability estimation models show an applicable estimation of permeability with R² values of 0.955 and 0.962 for models using D_max (pore diameter corresponding to displacement pressure) and D_f (pore diameter at inflection point), respectively. These results indicate that both D_max and D_f are key factors in determining permeability.