北极冰-水耦合系统中声波传播特性

利用固体和流体介质中波传播理论,导出了冰-水两层复合结构中导波频散方程。进一步,利用二分法对频散方程进行了数值求解,得到了ω-k频散曲线(ω与k分别为圆频率和波数),以及相速度和群速度频散曲线。结果表明:冰-水两层复合结构中导波由具有相同厚度水层和冰层中导波耦合而成,但与水层和冰层中导波频散曲线相比,复合结构中导波频散曲线除第1阶模式外,其余高阶模式均发生了很大变化。从原水层第1阶模式的截止频率开始,复合结构第2阶模式的相速度曲线被压低,各高阶(大于2阶)模式的相速度曲线出现一个跃变点,群速度曲线出现一个极大和一个极小值。水层越厚,复合结构各高阶模式的截止频率越低,相同频带内导波模式越丰富。水层厚度保持不变时,复合结构各阶模式的相速度和群速度曲线均随冰层厚度的增加而向低频方向移动。另外,还进一步分析了冰-水复合结构的导波波结构,发现第1阶导波模式的能量主要集中在冰层内和海表面附近,而2阶以上高阶导波模式的振动位移幅度随深度方向呈现周期性特征,并且模式阶数越高,振动越复杂。      

充水有限长圆柱薄壳声散射:Ⅰ.理论

研究水中充水有限长圆柱薄壳的声散射。采用弹性薄壳理论结合边界条件导出了有限长圆柱薄壳倾斜入射时散射声场的解析解,内部充水的作用反映在系统总阻抗中增加了内部流体阻抗。数值计算表明,在反向散射声场中,刚性散射部分仅在正横附近有比较大的贡献,而斜入射时起主要作用的是弹性散射部分。散射波相速度频散曲线的计算表明,与内部真空的情况相比内部充水后除正横附近的纵波和剪切波以外,增加了一组内部流体引起的附加波,其对散射声场贡献非常重要。在计算的频率-角度谱图中,内部流体附加波呈现"碗"形共振曲线。      

无限长粘弹性圆柱管中轴对称波的传播模式和衰减

在Kelvin-Voigt线性粘弹性模型框架内研究了无限长粘弹性圆柱管中轴对称波的传播和衰减。通过一个复根搜索程序对频散方程的求解得到无量纲化相速度频散曲线和衰减曲线。比较了引入损耗因子后粘弹性圆柱管与弹性圆柱管中轴对称波传播的特性。给出了圆柱管的外内径之比、损耗因子和材料参数对相速度和衰减的影响。分析结果表明与弹性圆柱管相比粘弹性圆柱管中轴对称波的各阶传播模式均存在衰减,高阶传播模式并无严格意义的截止频率。损耗因子对第一阶传播模式的相速度影响很小,而对衰减的影响则比较大,穿孔率对波传播的相速度和衰减有相当大的影响。     

液态金属快堆主管道中超声导波传播特性研究

讨论了钠冷快堆(Sodium-cooled Fast Reactor,SFR)主管道的整体温度和内部液态金属钠流动速度的变化对管道导波传播特性的影响。推导了充液管道中导波频散方程的一般形式,并给出了管道内液态金属钠处于流动状态下的导波频散方程。采用数值计算方法获得了管内液态金属钠处于不同温度和不同流速时的导波纵向模式频散曲线和导波时域波形。结果表明,温度变化对基阶纵向模式的影响较小,但对高阶纵向模式的影响较大;液态钠流速增大会使导波频散曲线向高频轻微移动,但在实际检测中可以忽路管内液体流动速度的影响。通过对时域接收波形的模拟计算,进一步考察了液态金属钠的温度及流动速度变化对导波传播的影响,并通过对比不同模态的激发特点和不同频段的导波时域波形特点,结合导波频散曲线,给出了适用于SFR管道超声无损检测的导波模态和声源激发频段选择方案。      

偶极子声波测井在套管井中的声场特征*

偶极子声波测井在套管井中的应用越来越多,本文重点研究了套管井不同胶结状况下偶极子声源在套管井中激发的声场特征,通过解析方法分析了水泥环第一界面和第二界面不同胶结状况下井孔模式波的频散和衰减。计算结果表明,在水泥环第一界面或第二界面胶结差时井孔声场中可观测地层弯曲波和套管弯曲波;与胶结良好时相比,地层弯曲波的频散曲线向低频偏移,水泥环第二界面胶结差时偏移较大;套管弯曲波对水泥界面存在的流体环厚度很敏感,随着流体环厚度的增加套管弯曲波的相速度逐渐增大,激发幅度也明显增强,这一特征为利用套管弯曲波的相速度估算流体环厚度提供了可能。本文的计算结果为利用套管井中越来越多的偶极子声波测井资料评价水泥胶结质量奠定了理论基础。     

弯曲模式声波在充液井孔中的传播

本文由波动方程讨论了充液井孔中,弯曲模式声波在弹性和非弹性、无限大井壁介质以及有限大井壁介质情况下沿井轴方向的传播。数值计算了它们的相速度频散曲线和井孔的频率响应特性,计算了衰减正弦信号在井孔中激励的波形。利用PZT压电晶片的弯曲振动模式,制作了偶极子测井换能器,在混凝土浇注的人工模拟井中进行了测试,试验结果与理论计算结果作了对比。     

含孔隙分层半空间中瑞利波的传播特性研究

为了研究含孔隙介质分层半空间中瑞利波的传播规律,分析孔隙介质参数对瑞利波频散曲线的影响,本文进行了数值模拟研究。采用传递矩阵算法,计算了含孔隙分层半空间中一定频率范围内瑞利波所有模式的频散及激发强度曲线,并与均匀弹性固体分层半空间情况作了类比分析,在含孔隙覆盖层的两层模型和含低速孔隙夹层的三层模型下,详细研究了孔隙度、渗透率、层厚度等参数对瑞利波各模式的影响,发现孔隙度及层厚度的变化对频散曲线影响较大,而渗透率的变化对频散曲线影响较小。      

充粘液管材中超声纵向轴对称导波的频散特性分析

模拟分析充粘液管内径-壁厚比和超声纵向轴对称导波波包宽度的关系。考虑了衰减的情况下,在各传播距离、频厚积和激发脉冲周数上,采用脉冲回波法,分析了导波的频散特性。结果表明:对于特定的材料,当频厚积一定时,充粘液管材中传播的导波特性只与内径-壁厚比有关;各高阶模式的截止频厚积随内径-壁厚比的增大而减小。用L(0,1)模式检测内径-壁厚比在4附近的充粘液管时灵敏度较低;而用L(0,2)模式检测内径-壁厚比在9.2附近的管时灵敏度较低。对于不同内径-壁厚比的充粘液管,检测中应当用不同的纵向轴对称导波模式,所用的激发脉冲周数和频厚积也应当不同。     

孔隙介质包裹的充液管道结构中导波传播特性

研究孔隙介质包裹的充液管道中纵向导波传播特性,分析孔隙介质参数对频散曲线的影响。建立了孔隙介质包裹充液管道的结构模型,利用孔隙介质弹性波动理论,建立对应的频散方程,数值模拟计算得到该模型的频散曲线和时域波形,并分析了孔隙介质参数以及管壁厚度对频散曲线的影响。结果表明孔隙介质的渗透率对于导波频散的影响较小,孔隙度的改变对时域波形的位移幅度影响较大。同时,导波存在衰减,且衰减随孔隙度增大而增大。所得结论为埋地管道无损检测方面提供一定理论参考。      

二维矩形波导管内和频差频声场的非线性传播

本文采用声场部分波理论和二阶微扰方法，在二维矩形波导管内，讨论了两个具有不同频率的基波模式相互作用产生的和频差频声场的非线性传播问题。当两个模式的相速相等时，和频差频声场才具有非线性性质；和频差频声场的对称性取决于两个基波模式的对称性。由二阶非线性波的初始激发条件，本文还得到了和频差频声场的解析表达式。     

Borehole guided waves in a non-Newtonian (Maxwell) fluid-saturated porous medium

<正>The property of acoustic guided waves generated in a fluid-filled borehole surrounded by a non-Newtonian (Maxwell) fluid-saturated porous formation with a permeable wall is investigated.The influence of non-Newtonian effects on acoustic guided waves such as Stoneley waves,pseudo-Rayleigh waves,flexural waves,and screw waves propagations in a fluid-filled borehole is demonstrated based on the generalized Biot-Tsiklauri model by calculating their velocity dispersion and attenuation coefficients.The corresponding acoustic waveforms illustrate their properties in time domain.The results are also compared with those based on generalized Biot's theory.The results show that the influence of non-Newtonian effect on acoustic guided wave,especially on the attenuation coefficient of guided wave propagation in borehole is noticeable.     

Numerical study on scanning radiation acoustic field in formations generated from a borehole

The directivity of acoustic transducers used in conventional acoustic logging tools is uncontrollable[1,2], which inevitably affects investigation depth and resolution. At present, deep and wide range of investigation in petroleum exploration is urgently re- quired. It is important to improve the exploration capability to find more complex and fine reservoirs[3], for which the direction of the radiated acoustic energy is a direct factor. Acoustic field in the formations generated by the source…     

Leaky modes and the first arrivals in cased boreholes with poorly bonded conditions

The generation mechanism of the first arrivals in the cased boreholes for the poorly bonded conditions is investigated. Based on the analyses of the Riemann surface structure of the characteristic function, the dispersion features, excitation spectra and contributions of modes excited in the cased boreholes with different cementing types are studied. The phase velocity dispersion studies of leaky modes show that high-order modes form "plateau" regions with one approximate velocity denoted by v separated by their cutoff frequencies, in which the phase velocity changes little with a considerable frequency range, while the group velocity keeps a relatively constant high value. Usually, the operation frequency range of a specific cementing evaluation acoustic logging tool is covered by such a "plateau" region. Mode excitation and contribution analyses show that the first arrivals in the cased boreholes for the poorly bonded conditions are the contributions from leaky modes, where the traveling velocity of the first arrivals processed by slowness time coherence(STC) method is equal to the approximated velocity v. Analyses on generation of leaky modes in the cased boreholes supplement the understanding of the generation mechanism of the first arrivals.     

Inversion of guided-wave dispersion data with application to borehole acoustics

The problem of inferring unknown geometry and material parameters of a waveguide model from noisy samples of the associated modal dispersion curves is considered. In a significant reduction of the complexity of a common inversion methodology, the inner of two nested iterations is eliminated: The approach described does not employ explicit fitting of the data to computed dispersion curves. Instead, the unknown parameters are adjusted to minimize a cost function derived directly from the determinant of the boundary condition system matrix. This results in an efficient inversion scheme that, in the case of noise-free data, yields exact results. Multimode data can be simultaneously processed without extra complications. Furthermore, the inversion scheme can accommodate an arbitrary number of unknown parameters, provided that the data have sufficient sensitivity to these parameters. As an important application, we consider the sonic guidance condition for a fluid-filled borehole in an elastic, homogeneous, and isotropic rock formation for numerical forward and inverse dispersion analysis. We investigate numerically the parametric inversion with errors in the model parameters and the influence of bandwidth and noise, and examine the cases of multifrequency and multimode data, using simulated flexural and Stoneley dispersion data.     

The group velocity variation of Lamb wave in fiber reinforced composite plate

Experimentally measured Lamb wave group velocities in composite materials with anisotropic characteristics are not the same as the theoretical group velocities which is calculated with the Lamb wave dispersion equation. This discrepancy arises from the fact that the angle between the group velocity direction and the phase velocity direction in anisotropic materials exists. Wave propagation in a composite material with anisotropic characteristics should be considered with respect to magnitude correction in addition to direction correction. In this study, S0 mode phase velocity dispersion curves are depicted with the variation of degree with respect to the fiber direction using a Lamb wave dispersion relation in the unidirectional, bidirectional, and quasi-isotropic composite plates. Slowness surface is sketched by the reciprocal value of the phase velocity curves. The magnitude and direction of the group velocity could be calculated from the slowness surface. The recalculated group velocities with consideration of the magnitude and direction from the slowness surface are compared with experimentally measured group velocities. The proposed method shows good agreements with theoretical and experimental results.     

倾斜地层下的井孔声场分析及声速校正图版制作

理解声波在倾斜井孔下的传播规律,对于准确获得目的层真速度有着至关重要的作用。通过有限差分和几何声学分析方法对比了倾斜地层下纵横波首波速度的变化规律,旨在建立倾斜地层下的声速校正图版。研究结果表明,当倾斜地层速度"上低下高"时,从全波中得到的声波速度不再是地层的真实速度,均大于目的层真速度。尤其是当地层界面倾角大于45°时,即便是源距超过了4 m,仍对首波速度有影响,必须对其校正。相反,当倾斜地层速度"上高下低"时,任何源距和地层倾角下的视速度值都接近目的层真速度,可以不予校正。根据以上分析,建立了倾斜地层下的声速校正图版,并指出下层速度的变化对于校正图版有着较大影响,源距的变化对于图版影响很小,在不同源距下可以使用同一个图版。      

Tube-wave propagation in a fluid-filled borehole generated by a single point force applied to the surrounding formation

Propagation of tube waves in an infinite fluid-filled borehole, generated by a single-force point source placed in the elastic surrounding formation, is analyzed in the long-wave approximation. Integral representations of the precise solution are obtained both for fast and slow formations. An asymptotic analysis of tube-wave propagation in the fluid-filled borehole is performed on the basis of these two integral representations. The complete asymptotic wave field in the borehole fluid for a fast formation consists of P and SV phases and the lowest eigenmode of the Stoneley wave (tube wave). For a slow formation the conical Stoneley wave (Mach wave) is generated. It appears only behind the critical angle defined by the ratio of the S wave velocity in the formation to the low-frequency Stoneley wave velocity and decays weakly with an offset. Asymptotic wave forms are in good agreement with wave forms obtained by straightforward calculations.     

Theoretical simulation of electroacoustic borehole logging in a fluid-saturated porous formation

Electroacoustic (E-A) logging describes the acoustic response to an electromagnetic (EM) source in a fluid-filled borehole surrounded by a porous medium. The E-A response is simulated by two different methods in this paper. In the coupled method, the EM field and the acoustic field are modeled using Pride's model, which couples Maxwell's equations and Biot's equations. In the uncoupled method, the EM field is uninfluenced by the converted acoustic field, resulting in separate acoustic formulation with an electrokinetic source term derived from the primary EM field. The difference of the transient full waveforms between the above two methods is remarkably small for all examples, thus confirming the validity of using the computationally simpler uncoupled method. It is shown from the simulated waveforms that an EM-accompanying acoustic field is coupled to the EM field and appears with an apparent phase velocity of the EM wave in the formation. Acoustic waves with the conventional acoustic velocities are also seen in the converted full waveforms. For the sandstone models used in this paper, when permeability is less than 1 Darcy, the E-A Stoneley wave amplitude increases with porosity, which is different from that in conventional acoustic-to-acoustic logging.     

利用气枪声源数据的地声参数反演*

本文通过对东海采集的气枪声源数据进行处理,基于地声参数对不同声场观测量的敏感性差异,利用分步反演的策略进行分析。首先选择距离接收阵七个不同的距离点数据,通过warping变换准确提取第三阶至第八阶简正波的频散曲线;利用海底声速对频散曲线敏感的特性来反演声速;由Hamilton经验公式求得海底密度;通过传播损失拟合获得海底衰减。不同距离点数据反演的海底声速与密度一致性较高。实验提取的频散曲线和和反演参数仿真结果、实验获取的传播损失与反演参数获得传播损失均拟合较好。