基于波动方程有限差分算法的接收函数正演与偏移
波动方程时空域有限差分数值解及吸收边界条件研究进展
２ 波动方程时空域有限差分数值解法
在波动方程时空域 有 限 差 分 数 值 求 解 中 ， 通常 低阶或高阶差分计算 采用二阶差分计算时间 导 数 、 空间导数 。 常 规 有 限 差 分 方 法 通 常 基 于 空 间 域 频 散关系设计空间差分系数 ；基于时空域频散关系的 有限差分方 法 则 是 基 于 时 空 域 频 散 关 系 设 计 空 间 差分系数 ， 因而能 获 得 更 高 的 精 度 ；优 化 有 限 差 分 方法通过优化方法设计 有 限 差 分 系 数 ， 能够进一步 提高精度 。 下面分别介 绍 常 规 有 限 差 分 方 法 、 基于 时空域频散 关 系 的 有 限 差 分 方 法 和 优 化 有 限 差 分 方法 。
３ ７
度为二阶 ， 而采用新 的 空 间 差 分 系 数 求 解 一 维 声 波 方程精度为 ２ 求解二维 、 三维波动 方程 可分别 犕 阶， 在８ 个、 ４ ８个方向达到２ 犕 阶精度
［ ］ ３ ８
频散速度与真实速度的比 值 ， 比值 越 接 近 １， 表示数 值求解频散越小 ， 精度越高 。 由图 １ 可见 ， 常规基于 空间域频散关系的差分方法精度低于基于时空域频 前者的频散曲线随速度变化很 散关系的差分方法 ， 大（ 图１ ） ， 而后者的频散曲线随速度变化则很小 ａ （ 图１ ） 。 ｂ 图２ 为 基 于 空 间 域 频 散 关 系 的 差 分 方 法、 基于 时空域频散关系的差分方法模拟得到的多层水平介
［ ］ １ ８， １ ９ ［ １ ５， １ ６］ ［ １ ７］ ［ ］ １ ２， １ ３
量大大增加 。 ２． ２ 基于时空域频散关系的有限差分方法 目前大部分差分方法在空间域确定空间差分算 子， 但 是 地 震 波 传 播 计 算 通 常 是 在 时 间—空 间 域 进 行的 。 如果将这些空间差分算子直接应用于求解波 动方程 ， 频散总是存在而且有时较大 。 例如 ， 当常规 二阶时间差分和 ２ 犕 阶空间差分算子应 用于求解波
基于FFT的波动方程VOFFLC控制
基于FFT的波动方程VOFFLC控制汪洋;江厚顺;汪杰;许冬进;尹彪【期刊名称】《工程数学学报》【年(卷),期】2024(41)1【摘要】针对复杂波动方程的无穷维特性,基于Simulink平台利用FFT(Fast Fouri-er Transform)方法将其从时域PDE(Partial Differential Equations)模型转化为频域ODE(Ordinary Differential Equation)模型,并在频域上搭建类似于集中参数的控制系统。
通过FFT和成熟FDM(Finite Difference Method)模拟实验结果的对比,证明采用FFT原理模拟PDE波动方程的思路正确;在频域ODE模型上施加自适应VOFFLC闭环控制,并设计了两种控制反馈规则。
其中,采用乘法法则的VOFFLC控制时,波动呈现和原有形态一致、而周期缩短和振幅减小的现象;采用减法法则的VOFFLC控制时,可以实现类似边界控制的结果,然而在空间维度上可以实现向量级控制,即实现对该维度上任意函数形状、插值函数或者散点的向量级别控制,而这是边界控制做不到的。
因而,基于FFT的波动方程VOFFLC控制有进一步的研究意义和广阔的实用价值。
【总页数】13页(P186-198)【作者】汪洋;江厚顺;汪杰;许冬进;尹彪【作者单位】长江大学石油工程学院;扬州工业职业技术学院化学工程学院【正文语种】中文【中图分类】TP13【相关文献】1.基于波动方程的聚焦点控制照明叠前深度偏移2.基于FFT的电压波动与闪变测量算法3.基于FFT的电压波动和闪变算法应用研究4.资源性商品价格短期波动预测理论实证分析——基于动态控制方程5.变系数波动方程所决定的控制系统的最小能量控制和快速控制问题因版权原因，仅展示原文概要，查看原文内容请购买。
基于频散介质标量波方程有限差分的正演模拟方法及装置[发明专利]
![基于频散介质标量波方程有限差分的正演模拟方法及装置[发明专利]](https://img.taocdn.com/s3/m/b6be470a524de518974b7d49.png)
专利名称：基于频散介质标量波方程有限差分的正演模拟方法及装置
专利类型：发明专利
发明人：吴国忱,杨凌云,印兴耀,宗兆云,曹丹平,张佳佳
申请号：CN202011497848.1
申请日：20201217
公开号：CN112578450A
公开日：
20210330
专利内容由知识产权出版社提供
摘要：本发明公开了一种基于频散介质标量波方程有限差分的正演模拟方法及装置。
该方法包括：设定炮点的震源位置参数和角频率参数；根据震源位置参数和角频率参数，构建引入衰减因子的频率域频散介质标量波方程；对频率域频散介质标量波方程进行离散，以得到频率域频散介质标量波方程的差分格式；根据频率域频散介质标量波方程的差分格式计算角频率参数对应的波场值。
本发明通过频率域进行波动方程的正演模拟方法，稳定性较高，能够更加精确地描述频散介质中地震波场的频率变化。
申请人：中国石油大学(华东)
地址：266580 山东省东营市北二路271号
国籍：CN
代理机构：北京润平知识产权代理有限公司
更多信息请下载全文后查看。
一种利用有限差分来正演模拟声波波形的方法
一种利用有限差分来正演模拟声波波形的方法
王晓飞;刘海涅
【期刊名称】《科技风》
【年(卷),期】2012(000)022
【摘要】声波波形对于研究并旁地层的情况有着非常重要的意义.有限差分是常用的正演方法.本文利用有限差分的方法对软地层和硬地层的不同模型进行计算,并对结果进行了分析.
【总页数】3页(P42-44)
【作者】王晓飞;刘海涅
【作者单位】中海油田服务股份有限公司油田技术事业部,北京市101149;中海油田服务股份有限公司油田技术事业部,北京市101149
【正文语种】中文
【相关文献】
1.利用哈特莱变换进行井间声波波场正演模拟 [J], 刘迎曦;张霖斌
2.流固边界耦合介质高阶有限差分地震正演模拟方法 [J], 吴国忱;李青阳;吴建鲁;梁展源
3.一种新型有限差分网格剖分方法在大地电磁一维正演中的应用 [J], 张辉;唐新功
4.利用远震波形反演和宽频带地震波正演模拟推断2008年汶川地震的破裂过程[J], Takeshi Nakamur;Seiji Tsuboi;Yoshiyuki Kaneda;Yoshiko Yamanaka;付萍杰;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
5.基于三维有限差分方法的三分量\r感应测井正演模拟 [J], 郭晨;陈晓亮;卢圣鹏
因版权原因，仅展示原文概要，查看原文内容请购买。
波动方程偏移
中国石油大学(华东)
基本原理
利用一系列脉冲测试得出结论：归一化的成像条件可以衰减浅层的成像噪声，并能补偿深层能量， 最终改善成像效果。 对比陡倾角处的能量可以看出检波归一化对陡倾角的能量稍微强一点，但是震源归一化补偿的角 度更大一些。
图3-22 不同成像条件的脉冲测试结果 (a) 互相关成像条件；(b) 震源归一化互相关成像条件；(c) 检波归一化互相关成像条件 中国石油大学(华东)
记录波场和震源波场的振幅比成像条件：I(x,y,z)=U(x,y,z,t)/D(x,y,z,t)
tmax
I ( x, y, z)
tmax
0
S ( x, y, z,t ) R( x, y, z,t )dt
振幅比成像条件提供了最
好的分辨率，确定上述成像条件的方程仅是成像时间的形式，它的定义与震源波场最大值一致，但无法保
基本原理
对该方案做脉冲测试，选取的参数为：介质速度为2000m/s，速度模型横纵向采样 率都是4m，脉冲放在x=400m，t=360ms处。理想情况下的脉冲响应应为一个半径 为360m的半圆，如图中的白色虚线所示，实验得到的实际的RTM脉冲响应与理论值 吻合比较好，特别是在大倾角的地方，也有能量存在。 CDP CDP
引 言
近年来，随着向量并行计算机及高性能集群的快速发展 和高精度速度分析与建模技术的实现，使得计算存储需求和 速度场要求对RTM的制约逐步减小，RTM重新进入人们的视 线。
相对其他方法而言，RTM用全波方程对波场延拓，避免对
波动方程的近似，无倾角限制，可以实现复杂构造高角度成 像、可以成像回转波、棱柱波等各种体波，可以进行多次反 射波成像，使得多次波收敛和聚焦；因此，RTM也凭借其对 复杂记录波场和地质构造的适应能力，为多次波的积极利用-
光纤Bragg 光栅地震检波器是利用光栅的波长调制原理
光纤Bragg 光栅地震检波器是利用光栅的波长调制原理,即利用外界的微扰振动来改变光栅的栅距,再转化为对应的波长变化量,通过检测波长的变化来测量加速度的大小. 由图(2) 可以看出,该检波器的结构是将光纤布拉格光栅固定在弹簧片上构成的. 当地面振动时,固定于地面的检波器随之运动,检波器中的光栅与检波器外壳产生相对运动,在光电接收器件的接收面上产生运动的莫尔条纹,经处理后,以数字信号形式输出[9～13 ] .研究表明,光纤Bragg 光栅以其特有的特性(高灵敏度和大动态范围,较强的抗干扰能力和具有一定的智能作用,与大地耦合作用好,谐波失真小等) 比压电检波器(此种检波器是根据某些物质的压电效应制成的. 当沿一定方向对某些电介质施力而使其变形时,介质内部就产生极化现象.压电检波器正是利用这种压电效应,将地震波引起的压电效应转变为电信号的一种机电转换装置. ) 不仅具有高的灵敏度系数,并且动态范围也足以实现地震勘探领域的不同频率段的需求,在地质石油勘探领域具有广阔的应用前景. [17～29 ]5 结论地震检波器是地震勘探中的重要环节,地震采集数据的品质基本上取决于检波器本身的品质,埋置环境与记录数据系统的性能. 目前,数据采集的动态范围已经达到了120 dB ,所以检波器本身的动态范围越来越成为地球物理勘探技术中的首要技术要求[30～40 ] .光纤布拉格光栅以其高分辨力,高信噪比,高精确度,高可靠性,为振动的测量提供了理想的技术手段. 利用光纤布拉格光栅作为检波器的敏感元件,可具有动态范围宽,抗电磁干扰,信号频带宽等特点,易于满足高精度,高分辨率的地震勘探要求. 光纤Bragg 光栅地震检波器在地震石油勘探领域,桥梁结构健康检测领域和海洋开发领域等具有广泛的应用价值,是值得研究并推广的一项新技术.[ 1 ] 陶果,多雪峰. 我国地球物理测井技术的发展与战略初探[J ] .地球物理学进,2001 ,16 (3) :98～101.[ 2 ] 阮驰,崔崧,高应俊. 光纤Bragg 光栅与石油仪器[ J ] . 石油仪器,2001 ,15 (6) :1～4.[ 3 ] 陆文凯,丁文龙,张善文,肖焕钦,赵铭海. 基于信号子空间分解的三维地震资料高分辨率处理方法[J ] . 地球物理学报,2005 ,48 (4) :896～901.[ 4 ] 陈祖传. 地球物理勘探技术的进展[ J ] . 地球物物理学进展,2000 ,10 (3) :1～19.[5 ] 张向林,郭果,刘新茹. 油气地球物理勘探技术进展[J ] . 地球物理学进展,2006 ,21 (1) :143～151.[ 6 ] Kersey A D. Fiber grating sensors [ J ] . Lightwave Technol ,1997 ,15 (8) :1442～1463.[ 7 ] Hill K O , Meltz G. Fiber Bragg Grating Technology Funda2mentals andOverviews[J ] . Journal of Lighwave Tet hnology ,1997 ,15 :1263～1274.[ 8 ] 苑立波. 光纤光栅原理和应用[J ] . 光通信技术,1998 ,22 (1) :70～72.[9 ] 梁磊,周雪芳. 新型光纤Bragg 光栅地震检波器的研究[J ] . 承德石油高等专科学校学报,2003 ,5 (1) :4～7.[ 10 ] 高华,李淑清,南忠良,陶知非,蒋诚志. 光栅地震检波器的研究[J ] . 航空精密制造技术,2003 ,39 (1) :40～42.[ 11 ] 王红落,常旭,陈传仁. 基于波动方程有限差分算法的接收函数正演与偏移[J ] . 地球物理学报,2005 ,48 (2) :415～422.[ 12 ] 袁子龙,韦丹宁,李婷婷. 高分辨率地震勘探智能程控型前置放大器的设计[J ] . 地球物理学进展,2006 ,21 (1) :300～303.[ 13 ] 姚陈. 地震三维矢量反射波场[J ] . 地球物理学进展,2006 ,21(2) :430～439 . [ 14 ] 于常青,祁晓明,朴永红,于文芹. 联合高分辨率地震和精细测井资料的剩余油气分析[ J ] . 地球物理学进展,2005 ,20 (2) :张向林,陶果. 油气井生产测井中的光纤传感技术[J ] . 地球物理学进展,2005 ,20(3) :796～800.[16 ] 周振安,刘爱英. 光纤光栅传感器用于高精度应变测量研究[J ] . 地球物理学进展,2005 ,20 (3) :864～866.[ 17 ] 陈海峰,肖立志,张元中,付建伟,赵小亮. 光纤Bragg 光栅在油气工业中的若干应用及进展[J ] . 地球物理学进展,2006 ,21(2) :572～577.[ 18 ] 陈相府,安西峰,王高伟. 浅层高分辨地震勘探在采空区勘测中的应用[J ] . 地球物理学进展,2005 ,20 (2) :437～439.[ 19 ] 徐朝繁,张先康,刘宝金,等. 高分辨折射地震资料处理方法及其应用[J ] . 地球物理学进展,2005 ,20 (4) :1052～1058.[ 20 ] 张军华,吕宁,田连玉,陆文志,钟磊. 地震资料去噪方法、技术综合评述[J ] . 地球物理学进展,2005 ,20 (4) :1083～1091.[ 21 ] 崔若飞,武旭仁,陈同俊. 煤矿地震数据管理系统的开发[J ] .地球物理学进,2005 ,20 (2) :374～376.[22 ] 周灿灿,王昌学. 水平井测井解释技术综述[J ] . 地球物理学进展,2006 ,21 1) :152～160.[ 23 ] 石建新,王延光,毕丽飞,等. 多分量地震资料处理解释技术研究[J ] . 地球物理学进展,2006 ,21 (2) :505～511.[ 24 ] 于景邨,刘志新,岳建华,刘树才. 煤矿深部开采中的地球物理技术现状及展望[ J ] .地球物理学进展,2007 ,22 ( 2) : 586 ～591.[ 25 ] 彭富清,霍立业. 海洋地球物理导航[ J ] . 地球物理学进[ 26 ] 孟庆生,楚贤峰,郭秀军,樊玉清,贾永刚. 高分辨率数据处理技术在近海工程地震勘探中的应用[ J ] . 地球物理学进展,2007 ,22 (3) :1006～1010.[27 ] 杨文采,于常青. 深层油气地球物理勘探基础研究[J ] . 地球物理学进展,2007 ,22 (4) :1238～1242.[ 28 ] 金翔龙. 海洋地球物理研究与海底探测声学技术的发展[J ] .地球物理学进展,2007 ,22 (4) :1243～1249.[29 ] 张向林,刘新茹,李健,卢涛. 我国油气开发监测技术进展[J ] .地球物理学进展,2007 ,22 (4) :1360～1363.[ 30 ] 李淑清,陶知非. 未来地震检波器理论分析[ J ] . 物探装备,2003 ,13 (3) :152～156.[ 31 ] 刘光林,刘泰生,高中录,李刚,姚光凯. 地震检波器的发展方向[J ] . 勘探地球物理进展,2003 ,26 (3) :178～185.[ 32 ] 柏冠军,吴汉宁,赵希刚,王靖华. 地震资料预测薄层厚度方法研究与应用[J ] . 地球物理学进展,2006 ,21 (2) :554～558.[33 ] 张进铎. 地震解释技术现状及发展趋势[J ] . 地球物理学进展,2006 ,21 (2) :578～587.[34 ] 原宏壮,陆大卫,张辛耘,孙建孟. 测井技术新进展综述[J ] . 地球物理学进展,2005 ,20 (3) :786～795.[ 35 ] 冷元宝,朱萍玉,周杨,王送来. 基于分布式光纤传感的堤坝安全监测技术及展望[J ] . 地球物理学进展,2007 ,22 (3) :1001～1005.[ 36 ] 陈会忠. 地震信息系统发展综述[J ] . 地球物理学进展,2007 ,22 (4) :1142～1146. [ 37 ] 池顺良. 深井宽频钻孔应变地震仪与高频地震学———地震预测观测技术的发展方向,实现地震预报的希望[J ] . 地球物理学进展,2007 ,22 (4) :1164～1170.[38 ] 陆其鹄,彭克中,易碧金. 我国地球物理仪器的发展[J ] . 地球物理学进展,2007 ,22 (4) :1332～1337.[ 39 ] 赖晓玲,刁桂苓,孙译. 用近场余震观测资料研究张北地震的发震构造[J ] . 地球物理学进展,2007 ,22 (1) :63～67.[ 40 ] 张福明,查明,邵才瑞,印兴耀. 天然气的测井勘探与评价技术[J ] . 地球物理学进展,2007 ,22 (1) :179～185.。
频率波数域波动方程偏移
（二）频率波数域波动方程偏移序：有限差分法是在时空域进行偏移，利用付氏变换可在频率波数域实现偏移。
1．偏移公式① 速度减半后的波动方程：042222222=∂∂-∂∂+∂∂tuV z u x u (6.4-67)② 对上式进行关于x 和t 的二维付氏变换，速度用常数，得0)4(22222=-+U k Vdz U d xω (6.4-77) 式中),,(ωz k U U x =是波场函数u(x,z,t)的二维付氏变换。
③ 求解（6.4-77），有两个解，分别对应着上行波和下行波。
偏移研究的是上行波的向下延拓问题，所以只取上行波解为：])4(exp[),0,(),,(21222z k Vj k U z k U x x x -=ωωω （6.4-78）物理意义：用地面波场的付氏变换),0,(ωx k U ,可求出地下任何深度处的波场的付氏变换),,(ωz k U x ，是频率波数域内的常速波场延拓公式。
④ 求地下任意深度处的波场u(x,z,t) 对（6.4-78）进行反付氏变换，得x x k t j x dk d e z k U t z x u x ωωπω)(),,(21),,(++∞+∞-∞+∞-⎰⎰=（6.4-79）⑤ 成像取t=0时刻的波场，由（6.4-79）得x x jk x dk d e z k U z x u x ωωπ+∞+∞-∞+∞-⎰⎰=),,(21)0,,(⎰⎰∞+∞-∞+∞-=π21x x x x dk d z k Vx k j k U ωωω]})4([exp{),0,(21222-+（6.4-80） 2．频波域波动方程偏移的特点优点：①利用快速付氏变换，偏移效率高。
②适合于大倾角的地区。
缺点：①速度横向变化大的地区不能用。
②必须注意采样间隔，以免出现假频。
（三）克希霍夫积分偏移 1.用克希霍夫积分解求解波动方程 2．维波动方程克希霍夫积分解(x,y,z,t)P13图6.1-12 克希霍夫积分示意图如果围绕着震源的封闭曲面Q ，已知 Q 面上波动的位移位φ(x 1,y 1,z 1,t)及其对时间对空间的导数,且这些值是连续的没有奇点。
接收函数波动方程叠后偏移方法
接收函数波动方程叠后偏移方法介绍陈凌，张耀阳中国科学院地质与地球物理研究所接收函数波动方程叠后偏移方法是我们借鉴勘探地震学中发展成熟的反射波偏移成像技术，2005年以来发展的一种新的天然地震接收函数偏移成像方法（Chen et al., 2005a; 2005b）。
该方法包括两个基本步骤：时间域CCP 叠加和深度域波场反向延拓。
与反射地震学中处理反射波资料的CMP 叠加类似，CCP 叠加用以提高资料的信噪比。
波场反向延拓则是一个将产生Ps 转换波（图1d）的转换波源偏移至其真实位置的过程。
在CCP 叠加中，用一个1D 参考模型将来自单个地震事件的接收函数按垂直入射情况做时间校正后，根据其转换点位置进行组合叠加。
所得到的CCP 叠加道集可以近似看作是零偏移距资料，即转换波从速度间断面的转换点处垂直向上传播至地表而被记录到的波场。
对于复杂结构，只采用CCP 叠加将会给成像结果带来假像，而偏移对于转换界面的正确成像是必不可少的（Ryberg and Weber, 2000; Poppeliersand Pavlis, 2003a）。
我们采用的波场偏移原理与反射地震学中普遍采用的爆炸反射面原理（Claerbout, 1985; Sheriff and Geldart, 1995）类似。
在爆炸反射面模型中：每一个反射界面都被当作是布满了密度正比于界面反射强度的爆炸源。
将这些爆炸源同时在时间t = 0 激发，地震波以实际速度的二分之一向上传播以满足双向传播走时（图1a, b）。
这样在地表上接收的波场可以模拟由CMP 叠加得到的零偏移距（零炮点-接收点距离）道集。
倒转上述过程，即将地表观测的零偏移距波场反向传播至t = 0 时刻，就可以得到一个使所有反射体正确归位的深度偏移成像结果。
在接收函数偏移中，通过CCP 叠加产生的零偏移距接收函数道集可以用来模拟地表记录到的将所有地下转换波源同时在t = 0 时刻激发得到的波场，而源的强度正比于Ps 的转换系数（图1c）。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
矿产资源开发利用方案编写内容要求及审查大纲
矿产资源开发利用方案编写内容要求及《矿产资源开发利用方案》审查大纲一、概述
㈠矿区位置、隶属关系和企业性质。
如为改扩建矿山, 应说明矿山现状、
特点及存在的主要问题。
㈡编制依据
(1简述项目前期工作进展情况及与有关方面对项目的意向性协议情况。
(2 列出开发利用方案编制所依据的主要基础性资料的名称。
如经储量管理部门认定的矿区地质勘探报告、选矿试验报告、加工利用试验报告、工程地质初评资料、矿区水文资料和供水资料等。
对改、扩建矿山应有生产实际资料, 如矿山总平面现状图、矿床开拓系统图、采场现状图和主要采选设备清单等。
二、矿产品需求现状和预测
㈠该矿产在国内需求情况和市场供应情况
1、矿产品现状及加工利用趋向。
2、国内近、远期的需求量及主要销向预测。
㈡产品价格分析
1、国内矿产品价格现状。
2、矿产品价格稳定性及变化趋势。
三、矿产资源概况
㈠矿区总体概况
1、矿区总体规划情况。
2、矿区矿产资源概况。
3、该设计与矿区总体开发的关系。
㈡该设计项目的资源概况
1、矿床地质及构造特征。
2、矿床开采技术条件及水文地质条件。