非均匀饱含黏性流体孔隙介质中声波传播及井孔声场分析

声波在饱含流体孔隙介质中的传播特性与流体的黏滞性及孔隙介质的非均匀性密切相关.本文在Biot理论基础上,考虑了孔隙流体的剪切应力及孔隙结构的非均匀性,采用含黏性流体孔隙介质中的波动理论,研究了孔隙介质中四种体波的频散和衰减特性,分析了慢横波对快纵波转换散射的影响,进一步推导了孔隙地层井孔中的模式波及其声场的解析解,研究了非均匀孔隙介质中井孔模式波和波列的特征.研究结果表明,含黏性流体孔隙介质中存在慢横波,慢横波的频散很强,其传播特征受到介质孔隙度、渗透率及孔隙流体黏度的影响.在非均匀孔隙介质中,与慢横波相关的剪切应力平衡过程不仅导致快纵波的频散和衰减,还会影响井孔伪瑞利波及斯通利波的传播特征.本文的工作完善了孔隙介质中声波传播的物理机制,为孔隙地层井孔声波的解释与应用提供了理论指导.     

孔隙地层井壁上的声波首波及其诱导电磁场的原因

声波测井时孔隙地层中的声波首波平行于井轴沿井壁传播，它既有轴向位移分量，又有垂直于井壁的位移分量.这种以快纵波速度传播的波, 不仅含有由快、慢纵波势给出的梯度场,而且还含有由横波势给出的旋度场.慢纵波势的梯度是渗流位移首波的主要构成成分, 也是声电效应测井响应中存在伴随声波首波电场的主要原因.首波包含有旋度位移场，是存在伴随声波首波磁场的原因. 关键词： 孔隙介质 声波首波 诱导电磁场 测井     

基于动电效应的震电测井方法研究进展*

与双电层和孔隙流体渗流有关的声波-电磁场耦合效应(动电耦合波)在油气储层勘探、地震电磁场等领域有着潜在应用价值。该文简述动电耦合波在理论模型以及井孔动电耦合波的实验测量、模拟计算方面的研究进展,并对油气储层测井领域进一步的动电耦合波研究进行了展望。基于动电信号的探测方法同时接收声波和电磁场两类信号,可避免单一类别信号方法的不足。实验获得了岩心的动电耦合系数,但其表征的孔隙岩石物理性质还有待进一步认识。基于耦合控制方程,开展了震电波场的计算和分析。结果表明,声源激发后,可接收到伴随声波的电磁信号和早于声波的电磁首波。这两类信号都既对影响声波特性的岩石模量、孔隙度和渗透率等参数敏感,又与电导率等岩石电学性质密切相关。非饱和岩石动电耦合波理论和基于动电耦合波的参数反演方法等问题有待进一步研究。     

软地层随钻测井仪器偏心情况下低频四极子源 激发波场的数值分析*

为了探索测井仪器偏心对随钻低频四极子波场的影响,利用有限元法(多物理场耦合有限元软件包)对软地层井孔中随钻仪器偏心情况下低频四极子声源激发的波场进行了数值模拟。通过对随钻四极子测井波形的分析可知,当采用低频随钻四极子偏心声源激发时,测井波列中除传播速度与低速地层横波速度相当的地层四极子波模式外,还存在一个幅度非常低、传播速度略高于井孔流体中声速的模式波,该模式被对应钻铤弯曲波模式,且该模式波的幅度随着随钻四极子源偏离井孔中心距离的增大呈二次方规律增大,地层四极子模式波幅度呈三次方规律增加。     

流体饱和砂岩中弹性波速度的一种动力学计算模型

本文在怀利(Wyllie)平均时间公式的基础上,结合声速测井的条件,引入固体骨架和孔隙流体的平均等效长度等三个假设,采用普通弹性波动力学的方法,导出流体饱和砂岩地层中弹性波传播方程,得到纵、横波速度的计算公式,从而将怀利的计算纵波速度的运动学模型,发展成计算纵、横波速度的动力学模型。本文结果与Raymer等人1980年提出的经验公式十分接近,也与Wyllie等人测定的空气饱和砂岩的声速数据相吻合。     

含少量气泡流体饱和孔隙介质中的弹性波

考虑孔隙流体中含有少量气泡,且气泡在声波作用下线性振动,研究声波在这种孔隙介质中的传播特性.本文先由流体质量守恒方程和孔隙度微分与流体压力微分的关系推导出了含有气泡形式的渗流连续性方程;在处理渗流连续性方程中的气体体积分数时间导数时,应用Commander气泡线性振动理论导出气体体积分数时间导数与流体压强时间导数的关系,进而得到了修正的Biot形式的渗流连续性方程;最后结合Biot动力学方程求得了含气泡形式的位移场方程,便可得到两类纵波及一类横波的声学特性.通过对快、慢纵波的频散、衰减及两类波引起的流体位移与固体位移关系的考察,发现少量气泡的存在对快纵波和慢纵波的传播特性影响较大.     

等离子体加载耦合腔慢波结构色散分析

利用磁化等离子体介电张量和纵向场分量法对任意大小轴向磁场中等离子体填充耦合腔慢波结构进行场分析，得到磁化等离子体填充电子注通道电磁波场分量的轴对称精确解.在此基础上，建立耦合腔分区模型、采用场匹配方法建立色散方程，并数值计算得出不同等离子体密度及磁场下的耦合腔色散曲线.对不同密度等离子体填充情况下的耦合腔色散特性、混合模式的形成机理以及等离子体空间电荷波进行了分析讨论. 关键词： 耦合腔慢波结构 等离子体 混合模式 色散特性     

基于Biot-喷射流统一模型Maxwell流体饱和孔隙介质中的弹性波

推广Biot-Tsiklauri声学模型的同时借鉴Dvorkin和Nur的工作,建立了具有任意孔径分布并顾及喷射流动机制的非牛顿流体饱和孔隙介质声学模型,研究了非牛顿流体(Maxwell流体)饱和孔隙介质中的弹性波的衰减和频散特性.着重讨论充孔隙Maxwell流体的非牛顿流效应对弹性波的频散和衰减的影响.研究表明,饱和流体的非牛顿流效应和喷射流动机制均是引起弹性波波频散和衰减的重要因素.依据非牛顿流体(Maxwell流体)饱和各向同性孔隙介质的Biot-喷射流声学模型,喷射流动只影响纵波的频散和衰减,而饱和流体的非牛顿流效应不仅影响纵波,而且还影响横波的频散和衰减.     

裂缝诱导TTI双孔隙介质波场传播特征

本文在考虑油藏内流体影响的基础上, 进一步讨论了裂缝诱导各向异性的极化角和方位角的影响, 对裂缝诱导TTI (tilted transverse isotropy)双孔隙介质模型进行了研究. 在裂缝诱导HTI (horizontal transverse isotropy)双孔隙介质理论的基础上, 用Bond变换推导了裂缝诱导TTI双孔隙介质的柔度系数矩阵和耗散系数矩阵, 从而建立了介质的一阶速度应力方程. 采用交错网格高阶有限差分法及PML边界条件, 对xoz平面内的2.5维矢量波动方程进行了数值模拟. 结果表明,裂缝的极化角和方位角的存在都会导致横波分裂, 而在双层裂缝诱导TTI双孔隙介质模型的分界面上,又会产生转换波的分裂和横波的再分裂现象, 这就增加了波场的复杂性, 从而为进一步研究实际地球介质的地震波场传播特征奠定了基础. 关键词： 裂缝诱导TTI 双孔隙 裂缝极化角 裂缝方位角     

弹性波斜入射声子晶体的缺陷模及其转型

推导出弹性波斜入射固-固掺杂结构声子晶体的转移矩阵和透射系数公式.计算固-固掺杂结构声子晶体中弹性波的透射系数.得到当横波斜入射时,透射波中横波的缺陷模随着入射角的增大而减弱,横波向纵波的转型随着入射角的增大而增强.当纵波斜入射时,透射波中纵波的缺陷模随着入射角的增大而减弱,纵波向横波的转型随着入射角的增大而增强.     

Acoustic mode waves and individual arrivals excited by a dipole source in fluid-filled boreholes

An approach of separating individual wave arrivals for a dipole logging is presented. The branch points are treated as a type of logarithm and the Riemann surface structure of the multivalued function is studied, that is, the displacement potential within the borehole. Based on the theoretical analysis, the complex poles contributing to the wave field on various Riemann sheets are investigated in detail for the case of a fast formation. It is shown that poles on Riemann sheet (0,0) are real and form branches of modes with dispersion. Mathematically, it is demonstrated that the flexural mode has no cutoff frequency, which is different from the traditional point of view. Poles on other relevant Riemann sheets are complex and form many branches on the complex frequency-wavenumber plane. Further investigation on the pole and branch cut contributions indicates that the vertical branch cut integration method has limitations in separating wave arrivals. By properly taking into account the complex poles on various Riemann sheets together with branch cut integrations, wave arrivals are separated from the full waveforms effectively for both the fast and slow formation models. Specially, there are complex poles on Riemann sheet (0,−1) in the vicinity of the compressional branch cut for a slow formation with a large Poisson’s ratio, which have small imaginary parts and contribute a lot to the p-wave arrival. Supported by the National Natural Science Foundation of China (Grant No. 10534040)     

Pseudo interface waves observed at the fluid/porous-medium interface. A comparison of two methods

At the fluid/porous-medium interface the pseudo-Rayleigh (pR) and pseudo-Stoneley (pSt) waves exist. The relation with the corresponding poles in the slowness plane is not unambiguous, depending on the choice of branch cuts. For a point-force excitation, the far-field Green's functions are computed using vertical branch cuts (method I) implying that the pR- and pSt-poles obey the radiation condition. Then, a separate pseudo interface wave is entirely captured by the corresponding pole residue because the loop integral along a branch cut contributes to a body wave only. When hyperbolic branch cuts are used (method II) the poles lie on the "principal" Riemann sheet. Then, also the loop integrals necessarily contribute to the pR-wave because the pR-pole is different from that in method I. They do not contribute to the pSt-wave when the pSt-pole lies on the principal Riemann sheet because the pole is identical to that in method I. When the pSt-pole has migrated to another Riemann sheet, however, the pSt-wave is fully captured by the loop integrals. In conclusion, the phase velocity and attenuation of a separate pseudo interface wave can be computed from the pole location in method I, but should be extracted from the full response in method II.     

Numerical Simulation of Shear-Horizontal-Wave-Induced Electromagnetic Field in Borehole Surrounded by Porous Formation

Seismoelectric fieM excited by purely torsional loading applied directly to the borehole wall is considered. A brief formulation and some computed waveforms show the advantage of using shear-horizontal （SH） transverseelectric （TE） seismoelectric waves logging to measure shear velocity in a fluid-saturated porous formation. By assuming that the acoustic field is not influenced by its induced electromagnetic field due to seismoeleetric effect, the coupling governing equations for electromagnetic field are reduced to Maxwell equations with a propagation current source. It is shown that this simplification is valid and the borehole seismoelectric conversion efficient is mainly dependent on the electrokinetic coupling coefficient. The receivers to detect the conversion electromagnetic field and to obtain shear velocity can be set in the borehole fluid in the SH-TE seismoelectric wave log.     

Resonance blocking and passing effects in two-dimensional elastic waveguides with obstacles

Resonance localization of wave energy in two-dimensional (2D) waveguides with obstacles, known as a trapped mode effect, results in blocking of wave propagation. This effect is closely connected with the allocation of natural resonance poles in the complex frequency plane, which are in fact the spectral points of the related boundary value problem. With several obstacles the number of poles increases in parallel with the number of defects. The location of the poles in the complex frequency plane depends on the defect's relative position, but the gaps of transmission coefficient plots generally remain in the same frequency ranges as for every single obstacle separately. This property gives a possibility to extend gap bands by a properly selected combination of various scatterers. On the other hand, a resonance wave passing in narrow bands associated with the poles is also observed. Thus, while a resonance response of a single obstacle works as a blocker, the waveguide with several obstacles becomes opened in narrow vicinities of nearly real spectral poles, just as it is known for one-dimensional (1D) waveguides with a finite number of periodic scatterers. In the present paper the blocking and passing effects are analyzed based on a semi-analytical model for wave propagation in a 2D elastic layer with cracks or rigid inclusions.     

Leaky modes and their contributions to the compressional head wave in a borehole excited by a dipole source

A method of studying the contributions of leaky modes to the wave field is presented based on the analysis of the Riemann surface structure of the characteristic function, and the sensitivities of con- tributions to various factors of interest are examimed. Numerical results show that their contributions to the compressional head wave are related to the distributions of complex poles on (-1, -1) and (0, -1) Riemann sheets on the frequency-wavenumber (ω - k) plane. For fast formations, their contributions ar...     

Stability problem in nonlinear wave propagation

An explicit expression for the excitation spectrum of the stationary solutions of a nonlinear wave equation is obtained. It is found that all branches of many-valued solutions of a nonlinear wave equation between the (2K+1,2K+2) turning points (branch points in the complex plane of the nonlinearity parameter) are unstable. Some parts of branches between the (2K,2K+1) turning points are also unstable. The instability of the latter is related to the possibility that pairs of complex conjugate eigenvalues cross the real axis in the κ plane. Zh. éksp. Teor. Fiz. 114, 1487–1499 (October 1998) Published in English in the original Russian journal. Reproduced here with stylistic changes by the Translation Editor.     

Excitation of a tube wave in a borehole by an external isotropic point source

The excitation of a tube wave in an infinite fluid-filled borehole by an external isotropic point source is considered. The solution to the problem is obtained in the form of a double integral with respect to the ray parameter (slowness) and frequency. The integral with respect to the slowness is transformed to a contour integral in the complex slowness plane and then reduced to the integral over the edges of the cut of the vertical slowness function and the semiresidues at the poles. An asymptotic expression for the wave field in the borehole is obtained with allowance for the radiation condition at infinity. It is shown that, when a longitudinal spherical wave is incident on the borehole, only one low-frequency Stoneley wave is excited and not two, as was assumed earlier [1].     

Electroseismic waves excited by vertical magnetic dipole in borenole

Acoustic and electromagnetic fields are coupled in a fluid saturated porous medium due to seismoelectric effect. Seismoelectric well logging method has been proposed to detect deep target formation utilizing such effect. Because of uncoupling of SH waves with P-SV waves, a simple and forthright way to get shear waves information is possible, especially for soft or slow formation whose shear wave velocity is lower than the velocity of borehole fluid. We consider the wave fields excited by a vertical magnetic dipole (VMD) source. Two methods are used to simulate, one is the coupled method based on Pride model and the other is the uncoupled method. For two methods, the frequency wavenumber domain representations of the acoustic field and associated seismoelectric field are formulated. The full waveforms of acoustic waves and electromagnetic wave induced SH waves excited by VMD source in the time domain propagation in borehole are simulated and analyzed.     

Adaptive optics wave function reconstruction and phase unwrapping when branch points are present

An algorithm is presented for reconstruction of adaptive optics wavefront sensor data which produces results that contain the part of the wavefront distortion associated with the hidden phase as well as the scalar phase part that is captured by an ordinary least-mean-square-error reconstructor. The algorithm uses a multigrid formulation and multiplication of complex phasors representing the measured phase differences to reconstruct the distorted wave function. The algorithm is formulated to work with Hudgin-geometry data, but in the appendix a variant of this algorithm is described that allows operation with Fried-geometry data. Also an algorithm is presented that generates a phase function from the reconstructed wave function, a phase function that has its branch cuts placed so that the 2π discontinuities of the branch cut occur where the optical intensity is much lower than the average intensity. The reconstructor algorithm is formulated as a noise-variance-weighted reconstructor. It is found that the algorithm's noise gain is only slightly greater than that of a noise-variance-weighted least-mean-square-error reconstructor so long as the noise variance for the input phase difference measurement data is less than about 0.25 rad².     

超短脉冲复宗量辛格高斯光束

利用理论解析推导的方法，在傍轴近似条件下，给出了一组新的超短脉冲光束的解析解，称为超短脉冲复宗量辛格高斯光束.此脉冲光束解的每个频率分量都是复宗量高斯光束，时间脉冲的形状为辛格函数.对这种超短脉冲光束及其在自由空间中的传输过程进行了较为细致的研究，讨论了超短脉冲复宗量辛格高斯光束的轴上光强、光强的横向分布、脉冲极性反转、脉冲延迟等性质. 关键词： 脉冲光束 缓变包络近似 脉冲传输