0 引言

声-弹耦合介质在现实世界中普遍存在，包括海洋环境、河流及坝体和船舶等.地震波在这种声-弹耦合介质中的数值模拟问题，一直是人们的研究热点和难点.Hou 等(2012)将地震波在该情况下的数值模拟方法分为两种：(1)单一控制方程法，即在流相和固相中采用相同的地震波波动方程，利用一阶速度-应力方程针对声-弹耦合介质进行地震波数值模拟(Virieux, 1986; Burns and Stephen 1990; Dougherty and Stephen 1991; Stephen 1996; Van Vossen et al. 2002，Basabe and Sen, 2015)，这种方法不必考虑界面的几何特征，实现方法简单，但是由于方程中的变量较多(各向同性条件下二维情况为5个，三维情况9个)，计算效率相对较低.(2)双方程耦合方法，即分别用声波方程和弹性波方程来描述地震波在液相和固相介质中的传播并利用连续性条件将两者介质联系起来.该方法需要精确考虑声-弹耦合介质的界面信息.Stephen(1983，1988)和Sochacki等(1991)利用声-弹界面的应力、应变连续性条件采用有限差分方法对声学和弹性介质的纯位移方程耦合正演模拟研究，虽然短时间内波场传播十分稳定，一方面，随着传播时间增加，界面耦合近似误差积累逐渐增加导致地震波传播不稳定，另一方面，界面处耦合差分方法很难适用于高阶格式，极大的限制了空间差分步长.Lee等(2009)结合声-弹界面连续性条件利用单元网格有限差分数值模拟了水平界面情况下的地震波传播，同样该方法很难推广到高阶差分格式，限制了实际应用.曲英明等(2015)和Qu等(2017)利用一阶速度-应力方程将该方法推广到海底起伏模型情况，克服了起伏界面模拟时的阶梯散射波场干扰.Komatitsch等(2000)利用谱元法数值模拟了水平及起伏海底条件的地震波传播，模拟结果相对稳定，Carcione和Helle(2004)利用伪谱法分析研究了海底不同形式波的传播特征.Zhang(2004)结合有限差分和有限元方法进行了单元格子法声-弹介质地震波数值模拟方法研究，该方法对起伏不规则海底界面适用性相对较强.宣领宽等(2013)考虑了四边形单元的双线性特性并结合变步长思想进一步提高了该方法的数值模拟精度.Choi等(2008)在二阶声压波动理论和弹性位移波动理论框架下采用频率域有限元方法正演模拟了声-弹耦合介质地震波传播，并将该方法应用到不同形式下的全波形反演研究中，进一步验证了该方法实现的稳定性和有效性(Kim et al，2009；Bae et al，2010，2011).Yu (2016a，2016b)通过声压与应力参数的定量关系将传统的一阶速度与应力方程变换为一阶速度-应力-压力方程并开展了海上各向同性和VTI介质的逆时偏移研究.上述方法中采用一阶速度-应力波动方程模拟精度相对较高且声-弹界面模拟稳定，但是在现有计算环境下进行实际应用时，势必会面临计算量大、计算效率相对较低的问题(Choi et al.，2008).基于有限元、谱元和单元格子法采用双方程耦合方法可以有效的模拟声-弹边界处的地震波传播，但是在该方法需要进行精确的网格剖分，尤其针对复杂三维情况下的网格剖分工作需要特别留意，以免引入不必要的正演模拟误差(Xu and McMechan, 2014).

频率域数值模拟方法不仅可以提取模拟结果的单频信息和无时间频散以及对于介质的泊松比适用范围广，还便于进行并行计算和黏弹及双相介质的地震波传播研究，在实际应用模拟中逐渐受到人们的重视.Lysmer 和Drake (1972)首先提出了频率域有限元地震波正演模拟方法，随后，Marfurt(1984)与Marfurt和Shin(1989)对该方法的实施过程以及差分格式、精度和稳定性等特点进行了系统的阐述和分析.为了提高空间模拟精度，Jo(1996)、Shin和Sohn(1998)以及Min(2000)分别提出了加权平均9点以及25点差分格式，有效的提高了数值模拟精度.刘璐等(2013)、Tang(2015)进一步提出了基于优化15点、17点和27点差分算子的频率域正演模拟方法.殷文等(200)、张衡等(2014)、刘财等(2014)和唐祥德等(2015)针对频率域数值模拟方法做了大量的分析研究工作.随着计算机数值计算能力的提高，频率域数值模拟方法逐渐在向三维问题以及实际资料应用等(Pratt, 1999；and Operto, 2009；高凤霞，2014；Li et al. 2015，2016；Amestoy et al.，2016)方向转变.

2018年培训内容和形式已逐渐趋向成熟。考虑以后将手机、电脑联机，操作通过扫描大屏二维码注册使用中、外文数据库，演示检索过程，让新员工有更直观的认识和更深入地参与。

为尽可能的减少频率域正演模拟的计算内存，提高计算效率，本文针对声-弹耦合介质，借助等效交错网格思想并充分考虑密度参数空间变化对地震波传播的影响(Bartolo et al., 2012，2015)，推导了声-弹耦合地震波波动方程.在流相介质和固相介质中分别采用非均质情况频率域二阶声压标量波、二阶纯位移控制方程，为保证流、固相介质间地震波能量的稳定传输和有效交换，提出了声-弹耦合界面处转换过渡层方法，并详细阐述了过渡层与上下介质空间差分具体耦合方法.在与非均质纯位移波动方程正演结果对比分析的基础上，首先采用各向同性单层流相介质模型进行正演模拟验证了声-弹耦合方程数值模拟中过渡层策略的有效性和准确性，随后又数值模拟了地震波在声-弹耦合介质简单模型和复杂Marmousi2模型中的传播，进一步验证了本文方法稳定性和准确性，同时该方法可以简单的推广到三维情况.

1 声-弹耦合介质正演模拟方法

1.1 声-弹耦合介质波动方程

一阶位移(速度)-应力方程是定量描述地震波在三维各向同性完全弹性介质中传播的基本理论，由运动平衡微分方程和应力-应变方程组成.其中，时间域运动平衡方程为

(1)

应力-应变方程为

T=C LTU,

(2)

其中，U=(U,V,W)T为质点位移向量，T=(τxx,τyy,τzz,τxy,τxz,τyz)T为正应力和剪切应力组成的应力分量，t为时间，ρ为密度，L和C分别为偏微分算子矩阵和各向同性弹性介质中的刚度张量矩阵，其具体表达式如下：

(3)

(4)

当地震波在流体(声学)介质中传播时，其剪切模量为零，即剪切应力消失，且声压P为正应力如下定量关系：

方程(11)和(12)中的其他偏微分项的差分格式和方程(13)与(14)类似，不再另行表述，此外，我们将采用空间4阶有限差分格式开展下文声-弹耦合介质频率域正演模拟研究工作.

(5)

其中，λ、μ为拉梅常数.

P=-τxx=-τyy=-τzz.

(6)

把公式(6)代入方程(1—2)，则得流体介质中一阶位移-压力方程：

(7)

其中，P为声压参数，Tf=(-P,-P,-P)T为声压向量，流相介质中的偏微分算子Lf和刚度矩阵Cf的具体表达为

(8)

+ρ(x)fx(x,ω)=0,

(9a)

(9b)

(9c)

其中，ω为圆频率， S(x,ω)为震源函数，x为(x,z).

扭转冲击工具顺时针冲击结束状态如图3b所示。此时，工具液流分为3部分：第部分直接通过节流喷嘴到达钻头；第2部分通过分流器到达上冲击腔，此时冲击锤已和下冲击面接触，顺时针冲击完毕；第3部分是启动器内的压力通过卸荷槽联通得到平衡。此时，因为启动器周围的环境压力处于近似相等的状态，启动器借助着顺时针转动的惯性继续转动。

(2) 设置的辅助点区域须在液相介质内.由于辅助点是用一阶位移-压力方程进行表征的，该方程组只在液相中成立，故过渡层中的辅助点须设置在液相介质中.若过渡层在固相介质中，则辅助点则用一阶位移-应力方程进行表征，本文暂不讨论该种情况.

(10)

含震源项其二维形式的表达式为

(11)

同样，对方程(5)进行傅里叶变换，其纯位移波动方程的二维形式为

对方程(7)进行傅里叶变换，表征流相介质的二维频率域一阶位移-压力方程组可写为

在全国组织工作会议上，习近平总书记发表重要讲话，首次提出新时代的组织路线：全面贯彻新时代中国特色社会主义思想，以组织体系建设为重点，着力培养忠诚干净担当的高素质干部，着力集聚爱国奉献的各方面优秀人才，坚持德才兼备、以德为先、任人唯贤，为坚持和加强党的全面领导、坚持和发展中国特色社会主义提供坚强组织保证。

(12a)

+ρ(x)fz(x,ω)=0,

统计学软件为SPSS 13.0，组间均数比较采用χ2检验，检验标准为α=0.05，P＜0.05为具有显著性差异。

(12b)

图1 交错网格差分格式示意图 (a) 传统交错网格; (b) 声压标量波等效交错; (c) 弹性波方程等效交错. Fig.1 The diagram of staggered-grid method (a)conventional staggered-grid method; (b、c) equivalent staggered-grid method for acoustic and elastic equation

其中，F=(fx,fz)为震源函数.

近红外光谱仪（Near Infrared Spectrum Instrument，NIRS）是介于可见光（Vis）和中红外（MIR）之间的电磁辐射波，美国材料检测协会（ASTM）将近红外光谱区定为780～2 526nm，是人们在吸收光谱中发现的第一个非可见光区。

1.2 等效交错网格思想

以二维为例，传统的交错网格差分格式是在空间和时间上应力(τxx,τzz,τxz)和速度(vx,vz)交错更新，每次都需要保存至少两个时刻的所有参数值，而等效交错网格在利用其周围应力或应变信息时，又可以采用其本身的差分格式来代替，从而实现仅含有纯压力(应力)或纯位移交错差分思想(如图1所示)，下面给出非均质情况二阶纯位移弹性波方程(12)x方向空间四阶的差分格式：

(13)

(14)

将公式(2)代入(1)中可求得三维非均质情况下二阶纯位移方程表达式：

1.3 声-弹介质频率域声-弹耦合波动方程

图2 频率域基于交错网格过渡层与上、 下等效交错网格耦合过程示意图 Fig.2 Diagram of staggered-grid method for transition layer in frequency domain

为尽可能的减少频率域数值计算中求解的线性方程组时所占计算内存，流相介质和固相介质分别主要采用二阶纯声压方程(公式(11))、纯位移方程(公式(12))进行定量表征，并提出了如图2所示的过渡层思想来保证两相之间的地震波能量的传输与交换，二阶纯声压和纯位移方程控制区域的离散点分别用图2中的红色矩形和红色三角形表示.在两种弹性波动方程控制区域的交界处，为保证流相介质中的地震波能量准确传输到固相介质中，在二阶纯声压方程控制区域边界处设置三层用一阶位移-压力方程(9)中的(9b)与(9c)公式表征的离散点(绿色三角形)作为辅助，将地震波能量传输到其下面二阶纯位移方程控制的区域介质中.反之，为保证固相介质中的地震波能量准确传输到液相介质中，在二阶纯位移方程控制的边界处同样设置三层用一阶位移-压力方程中的(9a)公式表征的离散点(黑色矩形)作为辅助，将地震波能量传输到其上面二阶纯声压方程控制的区域介质中，从而完成地震波能量在声-弹耦合介质中的数值模拟.则声-弹耦合频率域波动方程可表达为如下形式：

随着丝绸之路的兴起，各种商贸活动在这条古道上蓬勃开展，中国古代货币也随之被带到境外流通使用，起到了类似于国际货币的重要作用，进一步推动了商贸活动的顺利进行，促进了沿线各国的商贸和经济的发展以及社会文化的交流，并在此过程中，因材质、形制、文字、铸造技术等的选择不同而逐渐形成了具有时代性、继承性和传播性的钱币文化。这种钱币文化随同丝绸之路两端的东方中华文化、西方的希腊和罗马文化，以及丝路沿线两侧的印度文化、伊斯兰文化而彼此传播、互相影响、彼此交融，对东西方的社会经济文化发展都起到了重要的推动作用，其影响非常广泛和深远，中国的丝绸之路也因此闻名于世。

(15)

其中，S为流相介质中的震源项，为单位矩阵.

在实际处理边界层时，以下两点需要特别注意：

(1) 辅助点的层数由方程空间离散的阶数决定.为保证地震波能量在两相之间准确交换时，二阶纯声压和二阶纯位移方程控制区域边界点处的空间离散需要利用到辅助点的值，当采用的空间离散阶数逐渐升高，则辅助点的层数也要相应的增加.其中，当采用4、6、8、10阶空间离散时，则所需辅助点的层数依次为3、5、7、9.

将(7)式整理化简可得无震源项非均质情况下的二阶声压方程

在流、固相介质和过渡层内部同样采用交错网格差分格式，实现了整个数值模拟区域都采用交错网格差分思想，一方面保证了地震波在流-固相之间传播的稳定性，同时也保证了数值模拟的精度.本文空间4阶差分格式来模拟研究，该方法同样可以推广到更高的差分阶数.为模拟无限远区域的假设，采用完全匹配层(PML)吸收边界条件.

2 数值模拟结果

2.1 流相介质单层模型

图3 流相介质单层模型 Fig.3 The model containing fluid

另外，为进一步直观说明声-弹耦合频率域方程与纯位移波动方程计算优势，表1给出了声-弹耦合介质单层模型两种方法的待求参数个数及组建的阻抗矩阵大小结果，本文提出的声-弹耦合频率域波动方程数值模拟方法阻抗矩阵有所减小，计算效率相对提高，而且当模型的水层的厚度占比越大时，本文的计算优势体现的更加明显.

谈到研究会的发展离不开航天710所（以下简称710所），因为710所是研究会的发起者、建设者和研究会重要研究活动组织者、参与者及资源的提供者，研究会是710所面向社会经济领域重要窗口和科研平台，710所是研究会依托单位。

通过本文频率域数值模拟方法，首先计算求取不同参数不同频率的空间分布，图4为流相介质单层模型在30 Hz、60 Hz时声压(P)，图5为基于本文声-弹耦合方程和纯位移方程求取的水平位移(U)和垂直位移(W)不同频率下的频谱幅值.图6为抽取模型中第20道频率域频谱幅值的比较结果，对模型中各点进行短时傅里叶反变换(IFFT)后可获得地震波不同时间波场快照(图7、图8).通过图6可知，本文所描述方法(ESG)与传统交错网格方法(SSG)求取的正演模拟结果基本相同，并且在过程层界面没有任何人工噪声的产生，进一步通过对比两种模拟结果的时间域波场快照，验证了本文方法的准确性.

2.2 声-弹耦合介质模型

为进一步验证和分析本文频率域声-弹耦合波动方程对地震波在声-弹界面处传播特征刻画准确性，故设计一简单的声-弹耦合介质单层模型，模型如图9所示，海底界面为水平界面.模型大小为201×201，x、z方向的网格间距为4 m，水层的厚度为200 m，水下4 m激发主频为25 Hz的雷克子波，检波器的设置在水下16 m，PML吸收边界的厚度为120 m.基于本文方法数值计算不同频率下声压参数和位移参数的频率信息，同时为验证数值模拟结果的准确与否，同样进行了非均质情况纯位移波动正演模拟计算.图10给出了不同频率下的声压频谱快照，两种方法不同频率位移分量频率快照模拟结果如图11所示，并进一步抽取了图11中第20道的模拟结果进行对比分析，不同频率计算的数值模拟结果完全一致，如图12所示.对数值模拟结果经过短时傅里叶反变换，求得不同参数不同时间的波场快照(图13、14).由图可见，地震波以声压方式在液相介质中传播时，海底反射清晰可见.地震波以位移形式在固相介质中传播时，地震波转换为纵波和横波，与传统的交错网格正演模拟结果相比，两者的模拟效果一致，验证了本文方法的准确性.

再次，就贴牌生产的问题而言。重庆自贸区大约占地120平方公里。其中包括两江新区、沙坪坝区、重庆西部物流园、西永微电子产业园区。今年新建的重庆自由贸易区拥有全面开放的特点，这使得通关和出境条件便利，利于自由贸易区内定牌加工行为的增加，这也导致知识产权侵权现象会更加普遍。

图4 流相介质单层模型声压参数不同频率的分布情况 (a) 30 Hz； (b) 60 Hz. Fig.4 Snap shots of acoustic model in frequency domain

图5 流相介质单层模型频率域耦合方程(a—d)及纯位移方程(e—h)不同参数不同频率(a、c、e、g为30 Hz, b、d、f、h为60 Hz)的分布情况，从左到右依次为水平位移(a、b、e、f)及垂向位移(c、d、g、h) Fig.5 Snap shots of acoustic model using the coupled (a—d) and displacement (e—h) in frequency domain (a,c,e,g for 30 Hz, b,d,f,h for 60 Hz) equations of the horizontal (a,b,e,f) and vertical (c,d,g,h) displacement

首先，通过流相介质单层模型数值模拟验证本文数值模拟方法的准确性.图3为流相介质单层模型示意图，五角星代表炮点位置，距离水面8 m，三角形代表接收点位置，距离水面40 m，x与z方向的空间间隔为4 m，选用主频为30 Hz的雷克子波激发，时间采样间隔为0.00195 s.图3中虚线为过渡层的位置，距离水面280 m.

表1 声-弹耦合单层模型两种方法求解方程组大小比较 Table 1 The comparison for different methodsin the elastic media overlying the fluid phase

声-弹耦合波动方程纯位移波动方程待求参数个数P=80×261U=W=186×261U=W=261×261阻抗矩阵大小117972×117972136242×136242

图6 流相介质单层模型第20道频率域耦合方程及纯位移方程数值模拟比较 (a—b) 水平位移； (c—d) 垂向位移. Fig.6 The 20th trace for the coupled and displacement equations (a—b) The horizontal displacement； (c—d) The vertical displacement.

图7 流相介质单层模型声压参数波场快照 (a) 0.225 s； (b) 0.275 s. Fig.7 Snap shots of acoustic model in time domain (a) 0.225 s； (b)0.275 s.

为进一步验证和说明本文方法对复杂模型的适用性和稳定性，采用复杂Marmousi2模型进行测试，如图15所示，密度参数和纵波速度之间的变化趋势是一致的，故未在图中展示.需要特别指出，为满足模型的模拟精度，本文未考虑原始模型中海底松软低速层(高泊松比).模型大小为201×141，x,z方向的网格间距为5 m，水层的厚度为200 m，水下5 m激发主频为25 Hz的雷克子波，检波器的设置在水下10 m和海底，分别接收地震波传播时的声压和位移信息.通过声-弹耦合方程和纯位移方程频率域正演模拟计算，获取声压和位移参数的频率信息，图16和17给出了两种方法30 Hz时声压和位移分量的空间频谱快照，并抽取了图17中两种模拟方法第20道的位移分量进行比较，如图18所示，两者模拟结果完全相同.对频率域的正演信息进行短时傅里叶反变换，可求得时间域地震波声压和位移参数的传播快照(图19、20).在Marmousi2模型中，海底与弹性介质中的地层呈不整合接触，使得海底附近的构造相对复杂，通过与传统的交错网格正演模拟方法比较，本文的声-弹耦合数值模拟方法同样可对复杂海底地质条件环境进行有效的模拟，进一步说明了本文方法对复杂地下声-弹耦合环境具有较强的适应性.

图8 基于频率域耦合方程(a、b、e、f)及纯位移方程(c、d、g、h)流相介质单层模型波场快照(a—d) 0.225 s; (e—h) 0.275 s Fig.8 Snap shots using the coupled (a,b,e,f)and displacement (c,d,g,h)equations for different time (a—d) 0.225 s; (e—h) 0.275 s

图9 声-弹耦合介质单层模型示意图 Fig.9 The elastic model overlying the fluid phase

图10 声压参数不同频率的分布情况 (a) 30 Hz； (b) 60 Hz. Fig.10 Snap shots of acoustic model in frequency domain

图11 频率域耦合方程(a—d)及纯位移方程(e—h)不同参数不同频率(a、c、e、g为30 Hz, b、d、f、h为60 Hz)的分布情况， 从左到右依次为水平位移(a、b、e、f)及垂向位移(c、d、g、h) Fig.11 Snap shots of acoustic model using the coupled (a—d) and displacement (e—h) in frequency domain (a,c,e,g for 30 Hz, b,d,f,h for 60 Hz) equations of the horizontal (a,b,e,f) and vertical (c,d,g,h) displacement

图12 第20道频率域耦合方程及纯位移方程数值模拟比较 (a—b) 水平位移； (c—d) 垂向位移. Fig.12 The 20th trace for the coupled and displacement equations (a—b) The horizontal displacement； (c—d) The vertical displacement.

图13 流相介质单层模型0.3 s (a)、 0.35 s (b) 声压参数波场快照 Fig.13 Snap shots of acoustic model in time domain (a) 0.3 s； (b) 0.35 s

3 分析与结论

本文针对声-弹耦合介质，借助等效交错网格思想并充分考虑密度参数空间变化对地震波传播的影响，为尽可能的减少频率域正演模拟的计算内存，提高计算效率，推导了频率域声-弹耦合地震波波动方程，在流相介质和固相介质中分别采用非均质情况频率域二阶声压标量波、二阶纯位移控制方程，为保证声-弹相介质间地震波能量的稳定传输和有效交换，提出了声-弹耦合界面处转换过渡层策略，并详细说明了过渡层与上下介质空间差分具体耦合方法.基于不同模型通过对该方法的试算结果分析研究，取得了一定的认识如下：

(1) 与传统的交错网格正演模拟，本文提出的声-弹耦合介质频率正演模拟在保持与传统交错网格方法相同空间精度的前提下，尽可能的压缩了方程中的参数变量，降低了计算运行所需的内存，提高了计算效率.

(2) 通过简单声-弹单层模型和复杂Marmousi2模型试算并与非均质纯位移波动方程正演结果进行对比以及和解析解进行对比分析，不仅验证了基于等效交错网格频率域声-弹耦合介质数值模型的可行性，而且说明了本文所述方法的适用性和稳定性.

综上所述，孕产妇教育普及率仍待提高，年龄、参加孕妇教育及月收入是孕产妇健康素养的重要影响因素；母亲选择喂养方式与其对母乳喂养的态度及行为密切相关，同时健康素养合格的母亲更倾向于选择母乳喂养。

(3) 该方法同样可以简单的推广到三维情况下，并且相对于传统的交错网格有限差分方法，本文的声-弹耦合地震波数值模拟方法对于海上的地震偏移和全波形反演等方法和实际研究有很好的适用性，尤其适用于当前海上地震勘探OBC地震采集资料(同时采集的海水中的声压资料和海底的位移或速度资料).

致谢 十分感谢两位匿名专家对稿件提高的宝贵意见以及SPICE无私分享基于弹性波波动方程的地震波传播解析解程序.

References

Amestoy P, Brossier R, Buttari A, et al. 2016. Fast 3D frequency-domain full-waveform inversion with a parallel block low-rank multifrontal direct solver: Application to OBC data from the North Sea. Geophysics, 81(6): R363-R383.

Bae H S, Shin C, Cha Y H, et al. 2010. 2D acoustic-elastic coupled waveform inversion in the Laplace domain. Geophysical Prospecting, 58(6): 997-1010.

图14 基于频率域耦合方程(a、b、e、f)及纯位移方程(c、d、g、h)流相介质单层模型波场快照(a—d) 0.3 s; (e—h) 0.35 s Fig.14 Snap shots using the coupled (a,b,e,f) and displacement (c,d,g,h) equations for different time (a—d) 0.3 s; (e—h) 0.35 s

图15 Marmousi2模型 (a) 纵波速度； (b) 横波速度. Fig.15 Matmousi2 model (a) P-velocity； (b) S-velocity.

图16 Maimousi2模型声压参数30 Hz频率域波场快照 Fig.16 Maimousi2 model in frequency domain of 30 Hz

图19 Marmousi2模型声压参数波场快照 Fig.19 Snap shots of Marmousi2 model in time domain

图17 Maimousi2模型频率域耦合方程(a—b)及纯位移方程(c—d)不同参数30 Hz分布情况， 从左到右依次为水平位移(a、c)及垂向位移(b、d) Fig.17 Snap shots of acoustic model using the coupled (a—b)and displacement(c—d)equations of 30Hz of horizontal (a,c) and vertical (b,d) displacements

图18 Marmousi2模型第20道数值模拟比较 (a) 水平位移； (c) 垂向位移. Fig.18 The 20th trace for the coupled and displacement equations (a) The horizontal displacement； (b) The vertical displacement.

图20 基于频率域耦合方程(a、b)及纯位移方程(c、d)Marmousi2模型波场快照 Fig.20 Snap shots using the coupled (a,b) and displacement (c,d) equtions

Bae H S, Pyun S, Chung W, et al. 2011. Frequency-domain acoustic-elastic coupled waveform inversion using the Gauss-Newton conjugate gradient method. Geophysical Prospecting, 60(3): 413-432.

Burns D R, Stephen R A. 1990. Three-dimensional numerical modeling of geoacoustic scattering from seafloor topography. The Journal of the Acoustical Society of America, 88(5): 2338-2345.

Carcione J M, Helle H B. 2004. The physics and simulation of wave propagation at the ocean bottom. Geophysics, 69(3): 825-839.

Choi Y, Min D J, Shin C. 2008. Two-dimensional waveform inversion of multi-component data in acoustic-elastic coupled media. Geophysical Prospecting, 56(6): 863-881.

De Basabe J D, Sen M K. 2015. A comparison of finite-difference and spectral-element methods for elastic wave propagation in media with a fluid-solid interface. Geophysical Journal International, 200(1): 278-298.

Di Bartolo L, Dors C, Mansur W J. 2012. A new family of finite-difference schemes to solve the heterogeneous acoustic wave equation. Geophysics, 77(5): T187-T199.

Di Bartolo L, Dors C, Mansur W J. 2015. Theory of equivalent staggered-grid schemes: application to rotated and standard grids in anisotropic media. Geophysical Prospecting, 63(5): 1097-1125.

Dougherty M E, Stephen R A. 1991. Seismo/acoustic propagation through rough seafloors. The Journal of the Acoustical Society of America, 90(5): 2637-2651.

Gao F X. 2014. Multiparameter inversion from wave equations by frequency domain full waveform inversion [Ph. D. thesis] (in Chinese). Changchun: Jilin University.

Hou G N, Wang J, Layton A. 2012. Numerical methods for fluid-structure interaction-a review. Communications in Computational Physics, 12(2): 337-377.

Jo C H, Shin C, Suh J H. 1996. An optimal 9-point, finite-difference, frequency-space, 2-D scalar wave extrapolator. Geophysics, 61(2): 529-537.

Kim M H, Choi Y, Cha Y H, et al. 2009. 2-D frequency-domain waveform inversion of coupled acoustic-elastic media with an irregular interface. Pure and Applied Geophysics, 166(12): 1967-1985.

Komatitsch D, Barnes C, Tromp J. 2000. Wave propagation near a fluid-solid interface: A spectral-element approach. Geophysics, 65(2): 623-631.

Lee H Y, Lim S C, Min D J, et al. 2009. 2D time-domain acoustic-elastic coupled modeling: a cell-based finite-difference method. Geosciences Journal, 13(4): 407-414.

Li Y, Han B, Métivier L, et al. 2016. Optimal fourth-order staggered-grid finite-difference scheme for 3D frequency-domain viscoelastic wave modeling. Journal of Computational Physics, 321: 1055-1078.

Li Y, Métivier L, Brossier R, et al. 2015. 2D and 3D frequency-domain elastic wave modeling in complex media with a parallel iterative solver. Geophysics, 80(3): T101-T118.

Liu C, Yang Q J, Lu Q, et al. 2014. Simulation of wave propagation in two-phase porous media using a frequency-space domain method. Chinese Journal of Geophysics (in Chinese), 57(9): 2885-2899, doi: 10.6038/cjg20140914.

Liu L, Liu H, Liu H W. 2013. Optimal 15-point finite difference forward modeling in frequency-space domain. Chinese Journal of Geophysics (in Chinese), 56(2): 644-652, doi: 10.6038/cjg20130228.

Lysmer J, Drake L. 1972. A finite element method for seismology. ∥Bolt B A ed. Methods in computational physics II Seismology: Surface Waves and Earth Oscillations. New York: Academic Press Inc.

Marfurt K J. 1984. Accuracy of finite-difference and finite-element modeling of the scalar and elastic wave equations. Geophysics, 49(5): 533-549.

Marfurt K J, Shin C S. 1989. The future of iterative modeling in geophysical exploration. ∥Eisner E ed. Handbook of Geophysical Exploration: Seismic Exploration. Oxford: Pergamum Press.

Min D J, Shin C, Kwon B, et al. 2000. Improved frequency-domain elastic wave modeling using weighted-averaging difference operators. Geophysics, 65(3): 884-895.

Pratt R G. 1999. Seismic waveform inversion in the frequency domain, Part 1: Theory and verification in a physical scale model. Geophysics, 64(3): 888-901.

Qu Y M, Huang J P, Li Z C, et al. 2015. The application of variable coordinate system forward modeling using cell-based staggered grids in acoustic-elastic coupled medium. Geophysical Prospecting for Petroleum (in Chinese), 54(5): 582-591.

Qu Y M, Huang J P, Li Z C, et al. 2017. A hybrid grid method in an auxiliary coordinate system for irregular fluid-solid interface modelling. Geophysical Journal International, 208(3): 1540-1556.

Shin C, Sohn H. 1998. A frequency-space 2-D scalar wave extrapolator using extended 25-point finite-difference operator. Geophysics, 63(1): 289-296.

Sochacki J S, George J H, Ewing R E, et al. 1991. Interface conditions for acoustic and elastic wave propagation. Geophysics, 56(2): 168-181.

Stephen R A. 1983. A comparison of finite difference and reflectivity seismograms for marine models. Geophysical Journal International, 72(1): 39-57.

Stephen R A. 1988. A review of finite difference methods for seismo-acoustics problems at the seafloor. Reviews of Geophysics, 26(3): 445-458.

Stephen R A. 1996. Modeling sea surface scattering by the time-domain finite-difference method. The Journal of the Acoustical Society of America, 100(4): 2070-2078.

Tang X D, Liu H, Zhang H, et al. 2015. An adaptable 17-point scheme for high-accuracy frequency-domain acoustic wave modeling in 2D constant density media. Geophysics, 80(6): T211-T221.

Tang X D, Liu H, Zhang H. 2015. Frequency-space domain acoustic modeling based on average-derivative and GPU implementation. Chinese Journal of Geophysics (in Chinese), 58(4): 1341-1354, doi: 10.6038/cjg20150421.

Van Vossen R, Robertsson J O A, Chapman C H. 2002. Finite-difference modeling of wave propagation in a fluid-solid configuration. Geophysics, 67(2): 618-624.

Virieux J. 1986. P-SV wave propagation in heterogeneous media: Velocity-stress finite-difference method. Geophysics, 51(4): 889-901.

Virieux J, Operto S. 2009. An overview of full-waveform inversion in exploration geophysics. Geophysics, 74(6): WCC1-WCC26.

Xu K, McMechan G A. 2014. 2D frequency-domain elastic full-waveform inversion using time-domain modeling and a multistep-length gradient approach. Geophysics, 79(2): R41-R53.

Xuan L K, Zhang W P, Ming P J, et al. 2013. Numerical simulation of acoustic-elastic interaction based on multi-time-step unstructured finite volume time domain method. Chinese Journal of Solid Mechanics (in Chinese), 34(2): 158-168.

Yin W. 2008. Forward modeling and parallel Algorithm Based on high-order finite-difference method in frequency domain. Journal of Jilin University (Earth Science Edition) (in Chinese), 38(1): 144-151.

Yu P F, Geng J H, Li X B, et al. 2016a. Acoustic-elastic coupled equation for ocean bottom seismic data elastic reverse time migration. Geophysics, 81(5): S333-S345.

Yu P F, Geng J H, et al. 2016b. Separating quasi-P-wave in transversely isotropic media with a vertical symmetry axis by synthesized pressure applied to ocean-bottom seismic data elastic reverse time migration. Geophysics, 81(6): C295-C307.

Zhang J F. 2004. Wave propagation across fluid-solid interfaces: a grid method approach. Geophysical Journal International, 159(1): 240-252.

Zhang H, Liu H, Liu L, et al. 2014. Frequency domain acoustic equation high-order modeling based on an average-derivative method. Chinese Journal of Geophysics (in Chinese), 57(5): 1599-1611, doi: 10.6038/cjg20140523.

附中文参考文献

高凤霞. 2014. 频率域波动方程多参数全波形反演方法研究[博士论文]. 长春: 吉林大学.

刘财, 杨庆节, 鹿琪等. 2014. 双相介质中地震波的频率-空间域数值模拟. 地球物理学报, 57(9): 2885-2899, doi: 10.6038/cjg20140914.

刘璐, 刘洪, 刘红伟. 2013. 优化15点频率-空间域有限差分正演模拟. 地球物理学报, 56(2): 644-652, doi: 10.6038/cjg20130228.

曲英铭, 黄建平, 李振春等. 2015. 基于单元交错网格的变坐标系正演模拟方法在声-弹介质中的应用. 石油物探, 54(5): 582-591.

唐祥德, 刘洪, 张衡. 2015. 基于加权平均导数的频率-空间域正演模拟及GPU实现. 地球物理学报, 58(4): 1341-1354, doi: 10.6038/cjg20150421.

宣领宽, 张文平, 明平剑等. 2013, 基于不同时间步长时域非结构有限体积法模拟声-弹性耦合问题. 固体力学学报, 34(2): 158-168.

殷文. 2008. 基于频率域高阶有限差分法的正演模拟及并行算法. 吉林大学学报(地球科学版), 38(1): 144-151.

张衡, 刘洪, 刘璐等. 2014. 基于平均导数方法的声波方程频率域高阶正演. 地球物理学报, 57(5): 1599-1611, doi: 10.6038/cjg20140523.