各向异性黏弹性介质伪谱法波场模拟

当地震信号通过复杂地球介质时，地层除了表现为各向异性，还表现为内在的黏弹性特征 因此，为准确描述地震波在地球介质中的传播特征，理想的地球介质模型应该能够模拟岩石 的各向异性特征和衰减特征.给出了各向异性黏弹性介质模型的波动方程，推导了伪谱法波场正演模拟的递推公式，并利用伪谱法实现了地震波波场数值模拟.表明了该介质模型中地 震波场特征与各向异性主轴方位和介质的黏滞性参数之间的关系. 关键词： 黏弹性 各向异性 伪谱法     

方位各向异性黏弹性介质波场有限元模拟

随着裂缝性储集层的地位越来越重要，黏弹性各向异性介质中地震波传播理论的研究也越来越受到重视.在弱黏滞性条件下，针对标准线性黏弹性体模型，采用时间增量法，在每一个时间步长确定黏弹性变形及载荷，并在计算下一个时间步长的黏弹性变形时把它加载进去，建立类似于解弹性问题的黏弹性方位各向异性介质波动方程的有限元递推公式.数值模拟给出了黏弹性波在地球介质中的波场传播特征. 关键词： 黏弹性 方位各向异性 有限元 时间增量法     

线性黏弹性各向异性介质速度频散和衰减特征研究

地层岩石既非各向同性的,也非完全弹性的,正确地描述波在地层中的传播需要搞清楚岩石的各向异性和黏弹性特征.针对裂缝性地层的裂缝发育方位问题,建立具有任意方位角的方位各向异性黏弹性本构关系.使用均匀平面简谐波分析方法研究其频散关系,得到Christoffel方程,进而获得均匀平面波的复速度、相速度、衰减系数和品质因子的表达式.通过黏弹性方位各向异性页岩和砂岩进行模拟,研究了波场在地层中的传播特征,如相速度、衰减系数和品质因子等随频率、方位和入射角的变化关系 关键词： 黏弹性 各向异性 相速度 衰减系数 品质因子     

各向异性粘弹性孔隙介质地震波场伪谱法正演模拟

地球介质各向异性、粘弹性以及孔隙特征是地震资料描述油气储层时应综合考虑的. 建立了各向异性粘弹性孔隙介质模型，导出了各向异性粘弹性孔隙介质的弹性波波动方程，采用伪谱法正演模拟了各向异性粘弹性孔隙介质地震波，进行了其波场特征分析. 该研究有益于加深对地震波在实际地球介质中传播规律的认识.     

双相各向异性随机介质伪谱法地震波场特征分析

各向异性、双相孔隙以及非均匀性是描述油气储层时应综合考虑的. 结合随机介质理论和双相介质模型建立了双相各向异性随机介质模型,采用伪谱法模拟了双相各向异性随机介质地震波场,结果表明：双相各向异性随机介质地震波场存在散射波、旅行时扰动等复杂的波场特征,这些特征强烈依赖于随机介质模型参数. 在大非均匀空间尺度下,非均匀幅度主要影响波的旅行时扰动;在小非均匀空间尺度下,非均匀幅度主要影响波的散射. 该研究使人们有可能在统计意义下反演油气储层的非均匀特征,有益于加深对地震波在油气储层中传播规律的认识. 关键词： 双相各向异性 随机介质 伪谱法 地震波     

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

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

各向异性海底地层海洋可控源电磁响应三维积分方程法数值模拟

采用积分方程法对各向异性海底地层海洋可控源电磁（MCSEM）响应进行三维数值模拟.首先利用算子理论给出各向异性海底地层中的压缩积分方程,由于满足压缩映射条件,该积分方程在任意参数下总是迭代收敛的.然后为了提高计算效率,引入分区域多重网格准线性近似技术.通过具体算例验证了所述算法在计算精度与效率方面的有效性.最后利用该算法考察并分析了海底地层的各向异性对MCSEM三维响应特征的影响. 关键词： 海洋可控源电磁法 各向异性 三维模拟 积分方程法     

相控圆弧阵声源在TI地层井孔中产生的声场

研究了相控圆弧阵声源在各向异性地层井孔中产生的声场的特征,重点分析了该声源对于地层横波各向异性的方位敏感性,旨在探索一种新的地层声学各向异性测量方法。结果显示:相控圆弧阵声源在地层中可激发多种方位阶数的声场,主要成分为单极子波场和偶极子波场;在HTI地层竖直井中,方位角不同的圆弧阵声源激发的波动的速度不同,相控圆弧阵声源对于地层横波速度的方位各向异性的敏感性主要是它所激发的弯曲波的贡献。相比于传统的正交偶极子声波测井而言,采用相控圆弧阵声源可以评价井附近地层的周向非均质性,可以提高井周测量覆盖次数,提高反演结果的可靠性。另外,相控圆弧阵声源兼具备多极子声源的功能,优化的相位控制方式可以实现方位各向异性的最佳测量。      

利用积分方程法的各向异性地层频率测深三维模拟

应用一种新的积分方程算法对各向异性地层频率域电磁测深三维问题进行正演模拟.利用范数小于或等于1的修正Green算子得到各向异性地层中的新积分方程,由于满足压缩映射条件,该积分方程在任意参数条件下总是迭代收敛的.提出一种应用迭代法求解新积分方程的迭代初值优化选择方法,并通过具体算例对比说明该方法可使频率域电磁测深三维模拟的计算效率得到有效提高.最后应用所述算法对层状各向异性大地轴向频率测深视电阻率响应进行三维数值模拟,考察分析地层的各向异性对视电阻率响应特征的影响,得到一些重要结论.     

各向异性超常媒质矩形波导的导波特性研究

研究中通过建立和求解填充各向异性超常媒质的矩形波导TE模式和TM模式的波动方程和场分布方程,讨论了各向异性超常媒质矩形波导的有效介电常数、有效磁导率、相位常数以及传输系数与所填充超常媒质本构参数之间的变化关系.研究发现,填充了各向异性超常媒质的矩形波导的导波特性在截止频率之上的频率区域,其导波特性由所填充超常媒质的横向本构参数决定,而与超常媒质的纵向本构参数无关.但是,在截止频率之下的频率区域,该矩形波导的导波特性则由所填充超常媒质的横向本构参数和纵向本构参数共同决定.研究中通过典型的计算实例对所得到的结 关键词： 各向异性 超常媒质 矩形波导 导波特性     

激光激发瑞利波检测表面倾斜缺陷的研究

研究激光激发瑞利波检测样品的表面倾斜缺陷。基于频域热弹耦合方程,采用有限元方法建立激光激发瑞利波检测倾斜缺陷的数值模型,研究倾斜缺陷的检测机理。数值计算含不同的长度及倾斜角度缺陷的样品中瞬态位移波形,分析瑞利波在倾斜缺陷处模式转换的过程,研究各种瑞利波的传播路径。在此基础上,研究缺陷宽度与材料黏性对瑞利波传播及缺陷检测的影响。结果表明:瑞利波在缺陷处产生的反射及透射瑞利波的到达观测点的时间可以检测缺陷位置和长度,瑞利波在缺陷的底部发生模式转换产生的切变波可以检测缺陷倾斜角度。数值结果和已有的实验结果一致,从而为表面倾斜缺陷的检测提供有效的理论依据。      

Asymptotic analysis of Bloch-Floquet waves in a thin bi-material strip with a periodic array of finite-length cracks

We present an asymptotic algorithm to solve a problem of wave propagation in a thin bi-material strip with an array of cracks situated at the interface between two materials. For small frequencies we construct an asymptotic solution which takes into account the singular behavior near the crack tips and the smooth nature of the oscillation far away from them. We construct the boundary layer solutions near the crack tips. The boundary layers are harmonic solutions in scaled domains. Dispersion equations are derived and solved within the frame of the asymptotic model.     

Asymptotic analysis of Bloch–Floquet waves in a thin bi-material strip with a periodic array of finite-length cracks

We present an asymptotic algorithm to solve a problem of wave propagation in a thin bi-material strip with an array of cracks situated at the interface between two materials. For small frequencies we construct an asymptotic solution which takes into account the singular behavior near the crack tips and the smooth nature of the oscillation far away from them. We construct the boundary layer solutions near the crack tips. The boundary layers are harmonic solutions in scaled domains. Dispersion equations are derived and solved within the frame of the asymptotic model.     

Detection of fatigue-induced micro-cracks in a pipe by using time-reversed nonlinear guided waves: a three-dimensional model study

Localization of fatigue-related micro-cracks in pipelines is of increasing importance in industrial applications. A three-dimensional (3D) fatigue-crack imaging technique combining nonlinear guided waves with time reversal is proposed in this paper for potential applications in pipeline inspections. By using this method, the non-classical nonlinear guided waves generated from micro-cracks with hysteretic behavior are recorded, and the third harmonic waves are used to reconstruct the fatigue-crack images in a pipe by using a time reversal (TR) process. The feasibility of this method is examined by the imaging simulations for a steel pipe with varied defect areas. A finite-difference time-domain (FDTD) code is programmed to solve the wave equations under cylindrical coordinates, and simulate the experimental process of wave propagation. The results show that: (1) the proposed technique has excellent spatial retrofocusing capability; (2) the accuracy of defect localization and sizing depends on the crack orientation and the adopted guided wave mode; and (3) different displacement/stress components have varied sensitivities to the crack orientation.     

On propagation in continuous random media

The analysis of wave propagation in continuous random media typically proceeds from the parabolic wave equation with back scatter neglected. A closed hierarchy of moment equations can be obtained by using the Novikov-Furutsu theorem. When the same procedure is applied in the spatial Fourier domain, one obtains a closed hierarchy of coupled moment equations for the forward- and back-scattered wavefields that is not restricted to narrow scattering angles nor to small local perturbations. The general equations are difficult to solve, but a Markov-like approximation is suggested by the form of the scattering terms. Simple algebraic solutions can be obtained if a narrow-angle-scatter approximation is then invoked. Thus, three distinct approximations are explicit in this analysis, namely closure, Markov and narrow-angle scatter.

The results show that the extinction of the coherent wavefield has a distinctly different form from the corresponding result for propagation in a sparse distribution of discrete scatteres. Furthermore, when the scatter is constrained to narrow forwardand back-scattered cones, there is no back-scatter enhancement. These results are discussed within the context of the extension of the spectral-domain formalism to discrete random media. The general continuous-media moment equations are developed but not solved. The results correct and extend an earlier analysis that used a perturbation approach to compute the scattering functions rather than the Novikov-Furutsu theorem.     

各向异性渗流条件下弹性波的传播特征

研究了弹性波在非均匀裂纹孔隙介质中的传播特性,建立了各向异性喷射流模型.当弹性波通过裂纹孔隙介质时,由于波的扰动及裂纹和孔隙几何结构的不一致,导致在裂纹内部及裂纹与周边孔隙之间同时存在着流体压力梯度.此时的弹性波波动响应中包含着裂纹内连通性特征和背景孔隙渗透率信息.流体的动态流动过程使得介质的等效弹性参数为复数(非完全弹性),并且具有频率依赖性.当弹性波为低频和高频极限时,介质为完全弹性;当处于中间频段时,波有衰减和频率依赖.裂纹孔隙介质的各向异性连通性(渗透率)对应着各向异性特征频率(当渗流长度等于非均匀尺度时的弹性波频率),波的传播受到裂纹内连通性的影响.在一定频段内,随着裂纹厚度的增加,将出现第二峰值,峰值大小同时受到裂纹厚度和半径的影响.     

Temporal Solitons in Nonlinear Media Modeled by?Modified Complex Ginzburg Landau Equation Under?Collective Variable Approach

In the present paper, we have investigated the possibility of the existence of soliton solution in media modeled by the modified complex Ginzburg Landau equation. We have employed the technique of collective variables (CVs) to obtain a set of six coupled ordinary differential equations, one each for all the CVs included in the ansatz for the pulse. The coupled differential equations for the collective variables have been numerically solved to reveal the pulse dynamics which show stable soliton propagation.     

飞行器外流场对光传输的影响

针对细光束在飞行器外流场中传输受流场影响的问题，发展了光波标量衍射方程和气体运动的完全可压缩Navier-Stokes方程耦合的方法，并结合长度加权插值技术，对光束在飞行器定常外流场中传输的情形进行了数值仿真分析.结果表明,通过飞行器外流场的光束会发生偏转和离焦，为保证飞行器上光学装置的正常工作，须考虑必要的调校或补偿.