三维波动方程正演及模型应用研究
三维波动方程
三维波动方程
近年来，随着莫比乌斯旅行器技术的发展，三维波动方程在科学界和互联网领域中越来越受到重视。
三维波动方程（3D Wave Equation）又称三维Kirchhoff波动方程，是一种基于原点的数学模型，一般用于研究被称为波动的物理现象，比如声音、光等。
这种方程是用来描述一维、二维或三维电磁学波在介质中传播的高级模型，被应用到声学、电磁学、地震学和热力学等多个学科领域中。
三维波动方程在互联网行业发挥着越来越重要的作用，如在图像传输方面，3D Wave Equation可以把原本静态的图片转化为动态的响应帧，使图片显示更加生动活泼，增强用户体验。
此外，三维波动方程也被应用到音频行业，帮助实现更加生动的立体声。
而且，三维波动方程技术可以在全息图像、游戏开发、空间导航等方面实现进一步扩展。
三维波动方程在保证内容质量的同时还保证了高精度度和高稳定性，以更加精确准确的方式模拟物理世界，因此不仅仅在互联网领域具有重要意义，而且在更多领域有着广泛的应用前景。
根据技术发展趋势，三维波动方程将在互联网行业越来越受人关注，并开始发挥更大的作用。
基于三维克希霍夫正演的照明度分析
量分析 。基于照明度分析 的观测系统设计， 其判断观测 系统是否合理的依据为 ： 各种照明能量是否均匀 。
２ 实验 结果及 分 析
实验一 ： 理论倾斜界面模型 ， 模型分为 ３ ： 层 第一层 为倾斜界面， 速度为 １９ｍ／ ￣ ４２ ｓ第二层为倾斜界面 ， 速度 为 ２４ ｍ／ ； 三层 为水 平 界 面 ， 度 为 ２７ｍ／。采 ０３ ｓ第 速 ６Ｏ ｓ
用 ¨。照 明度模 拟是观 测 系统优化设 计所 采用 的一种 正演 方法 。它是观 测 系统优化设 计所 采 用的一
种 新 的分析评 价方 法 。把 三维 波动方程 克希霍 夫积 分 解和 射 线理论 结合起 来 ， 实现 正演模 拟 的 照 明
度模 拟 ， 用于指导和 优化观 测 系统设计 。
５ Ｏ
西 部探 矿工程
２ １ 年第 ｌ ０１ ２期
基 于 三维克 希霍 夫 正 演 的照 明度 分 析
未 睨
（ 西南石 油大学 ，四川 成 都 ６ ００ ） １５０ 摘 要 ： 今 ， 着勘探 对 象的 复 杂 ， 统 的 基 于共 中心 点和 共反 射 点重 合 的 观测 系统 设计 已不 适 现 随 传
来 指 导观 测 系统设 计 。 １５ 炮 点 一面 元 一检 波点对 应 能量分 析 ．
追踪出所有通过 目的层反射 的能量。从射线角度 上讲 ， 一个炮点 和一个 检波点就可 以确定一个地下 面 元， 记录每个炮 一面元一检波点对应的能量 。只有记录 炮 一面元一检波点对应 能量才可以进行真共反射点 能
２１ ０ １年第 １ 期 ２
移。
西部 探矿 工程
５ ｌ
１４ 炮 点入 射能 量分 析 ．
追踪 出所有通过 目的层反射的能量， 叠加当前炮点
含裂隙储层模型数值模拟对比分析
0引言地层介质中，孔隙作为一种储存空间，是油气资源运移或者扩散的重要通道，控制油气资源的产能。
深部介质受到上覆介质的压力，水平或者接近水平的裂隙空间被压实，高角度或者垂直的裂隙空间得到保留，HTI 模型是与高角度裂隙特别是垂直裂隙相适应的等效地质模型。
含裂隙岩石模型的发展大致经历了建立含裂隙储层模型、建立含裂隙储层模型的理论基础、含裂隙储层数值模拟3个阶段。
第一阶段，韩媛媛[1]对Eshelby 等效介质模型进行了综合论述和模型建立（如图1所示）。
图1 Eshelby等效介质模型第二阶段，在各向同性介质的背景下，加入单一结构裂隙（如定向排列的垂直或平行裂隙）和多元结构裂隙（如同时存在多个单一结构的裂隙系统），但各向同性的背景在实际地质应用中不占优势，对各向异性背景的研究更有价值[2]。
第三阶段，用数值模拟方法研究各种含裂隙模型的响应特征和适用范围[3]，常见的数值模拟方法有波动方程法和有限差分法[4]。
熊晓军等[5]用三维波动方程对缝洞地质模型进行了数值模拟。
本文采用波动方程方法对Hudson 模型和Cheng 模型进行数值模拟，比较、分析模型的优、缺点及适用范围，为研究含裂隙储层AVO 响应特征提供模型参考。
1裂隙储层等效模型分析1.1Hudson 模型Hudson [6-7]用椭球状裂隙近似模拟岩石介质中的扁平状裂隙，发展了含裂隙介质中的弹性波场理论。
Hudson 模型的假设有4点。
第一，弹性波的波长远大于定向排列的裂隙尺度。
第二，裂隙分布稀疏、均匀，裂隙所占比重小。
第三，裂隙分布不连续，每一个椭球状的裂隙彼此独立。
第四，裂隙厚度超过裂隙长度。
对于具有三轴对称的横向各向同性材料，弹性刚度张量由c 11、c 13、c 33、c 44、c 66 5个独立的弹性常数表示：=66444433131313111213121100000000000000000000000c c c c c c c c cc c c c（1）其中，c 12=c 11-2c 66。
三维波动方程的解法
三维波动方程的解法随着科技的发展，数字化软件的使用越来越广泛，计算机模拟已经成为了许多领域中不可或缺的重要工具，如气象预报、油气勘探、地震预测等。
在这些领域中，三维波动方程的解法是一项至关重要的任务，本文将介绍一些常见的解法。
一、有限差分法有限差分法是一个经典的数值解法，其基本思路是对微分方程进行离散化处理，然后求解离散化方程组。
在三维波动方程中，有限差分法的思路是将连续的空间进行网格化，将时间轴分成若干个时间步长。
这样，在每个时间步长内，将每个空间点一一计算，从而得到下一个时间步长的解。
有限差分法的优点是简单易懂，方便实现，但是它的精度可能会受到差分步长的影响，因此需要注意选择适当的差分步长。
二、有限元方法有限元方法是一类广泛应用于各个领域的数值分析方法，它的基本思路是将求解区域分割成若干小单元，然后通过求解每个小单元内的问题，从而得到整个求解区域的解。
在三维波动方程中，有限元方法可以将求解区域分割成若干四面体单元或者六面体单元。
通过计算每个单元的刚度矩阵和质量矩阵，可以得到离散化的三维波动方程。
然后通过求解离散化方程组，就可以得到整个求解区域的解。
有限元方法的优点是可以适应复杂的求解区域，精度高、收敛速度快。
当然，它的复杂度也比有限差分法高一些，需要更多的计算资源。
三、边界元法边界元法是一种利用边界条件而不是体系方程求解问题的方法，它的基本思路是将求解区域的边界分解成若干离散的小元素，然后通过求解每个小元素之间的关系，从而得到整个求解区域内部的解。
在三维波动方程中，边界元法可以将求解区域的边界分解成若干小面单元，然后通过求解相邻面积之间的关系，从而得到三维波动方程的解。
边界元法的优点是可以减少求解区域的离散化，大大降低了计算量。
然而，边界元法需要求解每个小元素之间的关系，该关系矩阵中存在一个对角主元，会导致矩阵的条件数很大，因此需要特殊的算法来求解。
结语：三维波动方程的解法有许多种，但是每种方法都有自己的优缺点。
三维波动方程有限差分正演方法
三维波动方程有限差分正演方法三维波动方程是描述地震波传播和地震波场的方程，是地震勘探、地震监测等领域中常用的数值模拟方法。
三维波动方程的有限差分正演方法是一种常用的数值解法，通过将连续的偏微分方程离散化，转化为差分方程进行求解，可以得到地震波场的时空分布情况。
∂²P/∂t²=c²(∂²P/∂x²+∂²P/∂y²+∂²P/∂z²)其中，P是地震波场的压力波动，t是时间，c是地震波传播速度，x、y、z是空间坐标。
有限差分正演方法通过离散化空间和时间进行数值求解。
对空间进行离散化，将地震波场的x、y、z坐标分别划分为Nx、Ny、Nz个网格点，其中Nx、Ny、Nz分别表示x、y、z方向的网格数。
对时间进行离散化，将t划分为Nt个时间步长，其中Nt表示时间步数。
将地震波场的压力P表示为P(i,j,k,n)，其中i、j、k表示网格点坐标，n表示时间步长。
在有限差分正演方法中，采用中央差分格式对空间导数进行离散化，采用二阶精度的差分格式对时间导数进行离散化，得到如下差分方程：P(i,j,k,n+1)=2P(i,j,k,n)-P(i,j,k,n-1)+c²Δt²(∂²P/∂x²+∂²P/∂y²+∂²P/∂z²)其中，Δt是时间步长。
上述差分方程可以通过迭代求解，从而得到地震波场的时空分布。
有限差分正演方法的优点是简单、直观，容易编程实现。
但是，由于空间和时间的离散化会引入数值误差，因此需要合理选择网格大小和时间步长，以保证数值解的精度和稳定性。
此外，由于三维方程的求解计算量较大，所以需要进行高效的并行计算。
总结来说，三维波动方程的有限差分正演方法是一种常用的地震波场数值模拟方法，通过离散化地震波场的空间和时间，将偏微分方程转化为差分方程进行求解。
波动方程数学模型
有关“波动方程”的数学模型
有关“波动方程”的数学模型如下：
波动方程是揭示波动现象奥秘的数学工具，描述了波在时间和空间中的传播行为，并在物理学、工程学等领域得到广泛应用。
无论是声波、电磁波还是水波，波动方程都能揭示它们的传播速度和振幅变化规律。
波动方程的一般形式可以表示为：∂²u / ∂t² = v² ∂²u / ∂x²，其中u表示波的位移或振幅关于时间t和空间位置x（以及可能的其他坐标）的函数，v表示波的传播速度。
这个方程在一维情况下描述了波沿着x轴正向传播且保持形状不变，而在二维或三维情况下，则需要考虑更多维度上的变化。
此外，波动方程还可以通过差商代替微商的方式进行表述，例如向前差商和二阶差商等。
总之，波动方程是理解和研究波动现象的重要工具，其数学模型为我们提供了深入探索波动行为本质的基础。
三维波动方程有限差分正演方法
三维波动方程有限差分正演方法三维波动方程有限差分正演方法是地球物理勘探领域中常用的数值计算方法之一，其主要应用于地震波传播与反演等领域。
一、三维波动方程有限差分正演方法原理三维波动方程的一般形式可以表示为：\[ \frac{\partial^2 p}{\partial t^2} = \nabla^2 p +f(x,y,z,t) \]其中$p$表示波场，$f(x,y,z,t)$表示源项函数，$\nabla^2$表示拉普拉斯算子。
对于三维波动方程的有限差分正演方法，其基本的数值离散形式如下：\[ \frac{p_{i,j,k}^{n+1} - 2p_{i,j,k}^n + p_{i,j,k}^{n-1}}{\Delta t^2} = c_x^2 \frac{p_{i+1,j,k}^n - 2p_{i,j,k}^n + p_{i-1,j,k}^n}{\Delta x^2} + c_y^2 \frac{p_{i,j+1,k}^n -2p_{i,j,k}^n + p_{i,j-1,k}^n}{\Delta y^2} + c_z^2\frac{p_{i,j,k+1}^n - 2p_{i,j,k}^n + p_{i,j,k-1}^n}{\Delta z^2} + f_{i,j,k}^n \]其中$p_{i,j,k}^n$表示波场在离散网格点$(i,j,k)$处的值，$\Delta t,\Delta x,\Delta y,\Delta z$分别表示时间和空间的离散步长，$c_x,c_y,c_z$分别表示波速在$x,y,z$方向上的离散形式，$f_{i,j,k}^n$表示源项在离散网格点$(i,j,k)$处的值。
该有限差分正演方法可以通过迭代求解，即根据当前时刻$t^n$的波场值$p_{i,j,k}^n$，计算当前时刻$t^{n+1}$的波场值$p_{i,j,k}^{n+1}$。
在迭代过程中，需要进行边界条件处理和源项的更新等操作，以确保该方法的数值计算精度和稳定性。
基于GPU的三维波动方程有限差分正演
基于GPU的三维波动方程有限差分正演[摘要]三维地震资料的处理和解释都需要有效的三维正演模型予以验证，实际工作中三维地震资料模型较大且结构复杂，在普通的桌面级计算机上难以完成正演运算，必须借助工作站甚至大型机，使其计算成本增加，成为实际生产应用中的瓶颈。
这里将利用GPU技术，研究波动方程正演模拟的高性能计算方法，使之能够满足生产中海量数据处理的需要。
[关键字] 三维声波方程有限差分正演GPU并行计算CUDA0前言三维地震勘探技术是一种信息量大、精度高的石油地球物理勘探方法，并已经在实际生产中大规模应用。
在复杂构造地区，三维资料的处理与解释需要与三维正演理论模型进行对比验证，国内外学者做了大量的研究工作并提出了许多高精度的三维正演算法。
一般来说，正演方法在网格规模快速增长的时候，计算量急速增长，使得普通的桌面级计算机无法承受运算耗时，必须通过工作站和大型机完成运算。
GPU并行运算可以有效降低运算耗时，使大数据量三维正演可以在桌面级计算机上完成。
1 三维有限差分算法的GPU算法三维GPU并行算法的访存模式比较复杂，因为要使用二维线程块逐一计算4个场分量，同一Block的不同线程不能重复利用处于不同面（即不同）网格内的波场值，所以这里的数据不能载入共享存储器。
在TFL模式中，在沿i方向逐步运算的过程中，当前线程会两次利用y和z场分量，因此只能使用寄存器存储当前和上一层的y和z场分量，当前层的波场值在计算过程中逐步转为上一层的波场值。
用方向场值举例递推TFL模式算法步骤如下：（1）存储当前层y、z方向场值，设置其为上一层场值；（2）从全局存储器读取当前层y、z方向场值，将其存储到寄存器；（3）对当前Block中所有线程进行同步；（4）将y、z场值从寄存器读出，将其载入共享存储器；（5）同步当前Block内的所有线程；（6）更新x并写入global存储器；（7）将当前层号加1，返回步骤（l）重复执行上述步骤，直到最后一层；2 数值实验所选测试模型为三维均匀速度模型，网格模型大小为109*109*109，xyz三个方向的网格步长均为1m，时间步长为0.2ms，介质速度为250m/s，有限差分精度为8阶。
砂泥岩互层地质模型中砂岩层偏移成像的特征分析
0前言随着松辽盆地南部长岭断陷勘探的深入，大面积的薄砂泥岩互层岩性油气藏成为主要勘探目标，该类油气藏的主要特点是储层厚度小，地层岩性和厚度横向变化较大，对预测精度要求高。
长岭断陷砂岩层的厚度2-10m，薄层内的多次波、转换波与有效反射波发生干涉现象。
薄互层反射的干涉效应和储层横向变化引起地震记录振幅、相位和频率也发生变化。
变化规律受薄互层内各岩层速度、单层厚度等因素的影响。
薄互层地震反射波场特征也会给后续资料处理及解释带来困难。
本文从正演出发，分析了简单砂泥岩模型中薄砂岩层的成像特征，继而通过建立长岭断陷的3D 地质模型，设计强化的观测系统，正演出三维单炮进行偏移处理，分析不同砂岩层厚度的成像特征。
1模型建立3D 模型建立分为两部分，构造建模和速度建模。
构造建模是根据解释成果，采用断裂恢复法由解释的断层和地层数据直接网格化生成地层面，通过简单定义断面边界、断面间主辅关系即可快速建立错综复杂的断层网络系统；速度建模是利用井资料、VSP 以及速度谱等资料，在构造、岩相和沉积层序模型的约束下，把井的纵向分辨率高和速度谱的横向分辨率高的特点结合起来，在三维空间沿层插值生成高精度的速度模型。
采用两步建模法，特别是利用井资料，显著提高了模型精度，速度互层明显，符合长岭断陷地区薄砂泥岩互层的特征，为准确反映目标地区的地质结构提供了保证。
可以刻画出砂泥岩互层的地质特征（如图1）。
2基于三维模型的波动方程正演三维波动方程正演采用的是最常用的有限差分法作为正演计算的数值模拟方法，它是将波动方程中波场函数的空间导数和时间导数用相应的空间、时间的差分来代替。
它的基本原理就是是把求解区域划分为差分网格，然后用有限的网格节点代替连续的求解域，利用微商与差商的近似关系将描述介质传播的微分方程转化为差分方程进行求解。
有限差分法具有计算速度快、占用内存小等优点，该方法对于近远场及复杂边界都有广泛的适用性，能够准确地模拟波在各种介质及复杂结构地层中的传播规律。
波动方程正演模型的研究与应用
波动方程正演模型的研究与应用郑鸿明* 娄 兵 蒋 立（新疆油田公司勘探开发研究院地物所）摘要野外采集的地震数据是经过大地滤波后的畸变信号，处理的地震剖面只是间接地反映了地下构造和地质体的特征，虽然目前有很多方法和手段可以分析并提取相关的地质信息，但由于处理对波场的改造和噪声的存在以及方法本身的多解性问题降低了识别地质信息的可靠性。
处理中每一步对有效信息的影响有多大，对地震属性解释的影响有多大，没有一个定量的标准，只能凭经验和认识来定性地判断。
正演模型在弹性波理论指导下，遵循严格的数学公式，可以最佳模拟地下各种情况。
各种处理方法和不同的处理流程所得到的结果能否符合或最佳逼近波动方程建立的数学模型，正演模型是判断处理工作合理性的良好准则。
主题词地质模型波动方程正演模型地震响应模块测试1 引 言随着地震勘探的不断深入，地震勘探也由构造型油气藏勘探进入精细的岩性勘探阶段，要求地震勘探能够反映地下地质体岩性变化，以及识别含油、气、水的地震响应特征，分辨薄互层、低幅度构造的能力。
地球物理学家们在长期的实践中已经研究开发了很多相关的技术，虽然理论上这些方法都能够成立，这些技术应用成功的实例也很多，但也不乏有失败的教训，往往产生多解性，或与钻探的结论不符。
这里除了复杂地表和复杂地下构造形成的复杂地震波场而不满足建立在简单地质模型处理理论的因素外，与处理过程对地震波场的改造也有很大关系。
从地震数据的采集到最终处理的地震剖面，整个过程是一个系统工程，地下地质结构、地质体的岩性变化以及含流体的性质，对处理人员来说是看不见、摸不着的“黑匣子”，我们所看到的只是经过大地滤波后产生畸变的地震波场，如何从这个畸变的地震波场中去伪存真、恢复真实的构造形态、提取储层的相关地震属性信息，这是岩性处理的最终目标。
处理中的每一步环环相扣、相互影响、相互制约，而我们对处理中的每一步产生的中间结果所应达到的标准只是凭经验、感觉进行定性判定，加入了很多人为因素，这些因素或多或少影响着我们对解释成果的正确认识。
三维波动方程有限差分正演方法
三维波动方程有限差分正演方法摘要：本文主要介绍了三维波动方程有限差分正演方法的基本原理和实现过程，并对其优缺点进行了分析。
通过数值实验验证了该方法的可行性和准确性，为地震勘探、地下水文学等领域的研究提供了参考。
关键词：三维波动方程、有限差分、正演方法、数值实验一、引言三维波动方程是地震勘探、地下水文学等领域的基础理论，其求解方法对于相关领域的研究具有重要意义。
有限差分正演方法是求解三维波动方程常用的数值方法之一，其具有计算速度快、精度高等优点，在实际应用中得到了广泛的应用。
本文将介绍三维波动方程有限差分正演方法的基本原理和实现过程，并对其优缺点进行分析，同时通过数值实验验证该方法的可行性和准确性。
二、三维波动方程有限差分正演方法的基本原理三维波动方程可以表示为：其中，u(x,y,z,t)为波动场，c(x,y,z)为介质速度，ρ(x,y,z)为介质密度，t为时间。
有限差分正演方法通过将空间和时间离散化，将三维波动方程转化为差分方程，进而求解波动场在不同时刻的数值解。
具体而言，有限差分正演方法将空间和时间分别离散化，将空间网格点和时间网格点相结合，构成一个四维网格空间，其中每个网格点对应一个波动场的数值解。
有限差分正演方法的求解过程可以分为以下几个步骤：1. 空间离散化对于三维空间，可以将其分为x、y、z三个方向，分别进行离散化。
假设在x方向上，空间步长为Δx，在y方向上，空间步长为Δy，在z方向上，空间步长为Δz。
则可以将空间网格点表示为(xi,yj,zk)，其中i=1,2,...,nx，j=1,2,...,ny，k=1,2,...,nz，nx、ny、nz分别为x、y、z方向的网格数。
2. 时间离散化假设时间步长为Δt，则在t时刻波动场的数值解可以表示为ui,j,k^n，其中n表示时间步数，i、j、k表示空间网格点的索引。
3. 有限差分近似将三维波动方程中的导数项用有限差分近似表示，例如：其中，ui,j,k^n表示在t时刻，在(xi,yj,zk)处的波动场数值解，c(xi,yj,zk)表示在(xi,yj,zk)处的介质速度，ρ(xi,yj,zk)表示在(xi,yj,zk)处的介质密度，Δx、Δy、Δz、Δt分别为空间和时间步长。
三维洛仑兹变换式用于三维波动方程
三维洛仑兹变换式用于三维波动方程三维洛伦兹变换式用于三维波动方程引言：在物理学中，波动方程是描述波动现象的一种基本方程。
在经典力学和电动力学中，波动方程可以描述振动的机械波和电磁波的传播。
为了研究波动现象在不同参考系下的行为，我们需要引入三维洛伦兹变换式，以便在不同参考系下分析波动方程的解。
一、波动方程的基本概念波动方程是一种描述波动现象的偏微分方程，通常用来描述波动的传播和干涉。
在三维空间中，波动方程的一般形式可以表示为：□^2φ - (1/v^2) ∂^2φ/∂t^2 = 0其中，φ表示波函数，□^2表示拉普拉斯算符，v表示波速。
这个方程描述了波函数在空间和时间上的变化关系。
二、三维洛伦兹变换式的引入为了研究波动方程在不同参考系下的行为，我们需要引入三维洛伦兹变换式。
三维洛伦兹变换式是狭义相对论中的基本变换式，用于描述时间和空间坐标在不同参考系间的变换关系。
在三维波动方程中，我们使用三维洛伦兹变换式将波函数的空间和时间坐标从一个参考系转换到另一个参考系。
三、三维洛伦兹变换式的具体形式三维洛伦兹变换式可以表示为：x' = γ(x - vt)y' = yz' = zt' = γ(t - vx/c^2)其中，x、y、z表示原参考系中的空间坐标，t表示时间坐标；x'、y'、z'表示新参考系中的空间坐标，t'表示新参考系中的时间坐标；γ表示洛伦兹因子，v表示参考系之间的相对速度，c表示光速。
四、三维洛伦兹变换式在波动方程中的应用将波函数的空间坐标和时间坐标代入三维洛伦兹变换式，我们可以得到波函数在不同参考系中的表达式。
通过对波函数在不同参考系中的变换，我们可以研究波动现象在不同参考系中的特性和行为。
五、结论三维洛伦兹变换式在三维波动方程中起到了重要的作用。
它描述了波函数在不同参考系中的变换关系，使我们能够更全面地理解和分析波动现象的特性。
三维波动方程的解
ｌ ｆ ， ＝０ （ ｚ ＋王 ， ¨＋ ） ， ｔ ＞ ０
方程 （ ４ ） 是亥 姆 霍兹 程 ，（ ４ ） 中 的 ｕ可 以代表 电场强
； ， Ｙ ， ｚ ， ０ ） ＝０ ， 一 ∞＜ ， ｙ ， ＜ ∞
｛ ， ， Ｙ ， 。 ， ０ ） ＝ （ — ， ｙ 一 叼 ，
Ｃ
Ｚ Ｃ 一，
ｎ
收 稿 日期 ： ２ ０ １ ２ — １ １ — １ ５
作者简介： 张子珍０９ ６ ５ 一 ） ， 女 ，Ｉ Ｊ Ｊ 西大 同人 ，敦授 ， 研究力‘ ： 论物理 。
“ 西大 Ｍ大学学报（ 自然科学 版 ） Ｅ一： Ｃｓ ｉ ｎ
ａ
Ｘ Ｓ ｉ １ ３
【 “ ， Ｙ ， ， ０ ） ＝ ， Ｙ ， ）
（ ３ ）
ｆ ４ １
先 求 出点 源 引 起 的基 本 解 ，再 利 用 叠 加 原 理
求 出连续 源的解 。点源 引起 的基 本解１ ｌ 所 满足 的
方 程 是 ｌ
方程 ｆ ３ １ 的解是
Ｔ （ ｔ ） Ｏ Ｃ ． ｅ － ｉ ｋ ＝ｅ ，
Ｍ ， Ｙ ， ， ｆ ） ＝ Ｍ ， Ｙ ， ｚ ） Ｔ （ ｔ ） ，
得
ｆ ｆ ） ＋ｋ ２ ａ ７ ＝０ ，
Ｍ＋七 “ ＝ ０。
ｉ “ ｆ ， ＝ａ ２ （ ｔ ｉ ＋ｕ ”＋“ ） ， ｔ ＞ Ｏ
｛ ｔ Ｉ （ ｘ ， Ｙ ， ， ０ ） ＝ ， Ｙ ， ｚ ） ， － ｏ ￣ ＜ ｘ ， Ｙ ， ｚ ＜ ｏ ｏ （ １ ）
标系下， 分别求解波动方程 ，并将各种情形赋予相 应的物 理 意义 ，以求 对 该 问题有 一个 整体 的认 识 。
鲁Ｊ 一 Ｊ 一 Ｊ 一 （ 州） ， 叼 ， ） ｄ 叼 ＝ 鲁
三维波动方程有限差分正演方法
三维波动方程有限差分正演方法
首先，将三维波动方程进行离散化处理，将连续的空间和时间域离散
为网格点的集合。
在空间域中，将地下介质划分为多个均匀网格，每个网
格点上对应一个地震波场的值。
在时间轴上，将时间域离散为多个时间步长，每个时间步长表示一个时刻。
然后，利用有限差分算子将三维波动方程转化为离散形式。
常用的有
限差分格式有常系数差分格式和变系数差分格式。
常系数差分格式适用于
各向同性介质，而变系数差分格式适用于各向异性介质。
在三维波动方程有限差分正演方法中，需要考虑各向异性介质的模型。
各向异性介质与各向同性介质不同，其波速和波阻抗在不同方向上具有不
同的值。
因此，在离散化过程中，需要引入特殊的差分算子来考虑各向异
性的影响。
最后，利用迭代求解的方法，按时间步长依次求解离散化的波动方程。
利用差分算子和初始条件，根据时空的变化逐步更新波场的数值，并计算
出每个网格点的波动值。
通过不断迭代求解，最终可以得到地震波在地下
的传播和反射结果。
三维波动方程有限差分正演方法在地震勘探中具有重要的应用价值。
它可以模拟地震波在地下介质中的传播和反射过程，帮助研究人员了解地
下结构、油气储层等。
同时，它也可以用来研究地震波在各向异性介质中
的传播规律，提高地震勘探的准确性和效率。
总之，三维波动方程有限差分正演方法是一种常用的数值模拟方法，
可以用来模拟地震波在地下介质中的传播和反射。
它具有广泛的应用价值，对地震勘探、岩石物理、地质测井等领域具有重要的意义。
三维波动方程基本解的一个求法
三维波动方程基本解的一个求法随着数学技术的发展，三维波动方程也得到了重视，已经成为许多科学研究的重要工具。
三维波动方程是一个模拟物理现象的数学模型，其基本解十分重要，它可以用来研究物理问题的准确性和评价模型的有效性。
本文将分析三维波动方程的基本解以及一种求解该基本解的方法。
首先，我们来了解一下三维波动方程的定义。
三维波动方程是一个常微分方程，包括了三个空间维度，空间变量可以由x, y, z三个坐标表示，时间变量可以由t表示。
空间和时间变量可以组合一个四元组(x, y, z, t)，这个四元组用来定义三维波动方程的解，这个解可以用下面的形式表示：u(x, y, z, t)=F(x, y, z, t)其中F(x, y, z, t)是一个函数，可以用来计算三维波动方程的基本解。
三维波动方程的基本解是一个空间动力盛行的系统，它产生的结果是分布在空间上的，因而其基本解有“定常解”和“波动解”之分。
定常解是指当t=0时，F(x, y, z, t)的值仅依赖于(x,y,z)三个空间变量，而随着时间变量t的变化，F(x, y, z, t)的值不会发生变化。
波动解是指F(x, y, z, t)的值随时间变化不断发生改变，但是空间变量(x, y, z)的变化不会影响函数F(x, y, z, t)的值。
两种解形式都十分重要，分析这两种解形式涉及到函数F(x, y, z, t)的研究。
求解三维波动方程的基本解的一种方法是分析函数F(x, y, z, t)的定义，并以此为基础构建一组数学模型。
根据空间变量(x, y, z)的定义，将其表示为由模型参数构成的向量，再根据时间变量t对F(x, y, z, t)进行求导，可以构成一组联立方程系统。
根据系统联立方程的具体形式，可以用几何技术、迭代法、梯度下降法等求解出一组解。
一旦求得解，就可以判断F(x, y, z, t)的变化，从而得出三维波动方程的基本解。
本文简要地介绍了三维波动方程的基本解和一种求解它的方法，但是求解三维波动方程的基本解是一个极其复杂的过程，它涉及到多种数学技巧以及数值分析技术。
三维波动方程正演技术在地震采集施工方法设计优化中的应用
２ 三 维 波 动 方 程 正 演 模 拟 对 采 集 施 工 方 法 优 化 设 计 及 验 证
针对 苏北 盆地 勘探 部署 几 上几 下 的 困难 地 区 ，即面积 超过 ｌ ｋ 的 Ｊ 城 区提 出基于地 质 目标 的三 Ｏ ｍ。 Ｈ 维波 动方 程 正演模 拟 及采 集方 法设 计 优化 研究 ，希望 通过 三 维地震 波 场模 拟来 指导 和 验证地 震采 集施 工
果 不符 合地震 波 的动 力学 特 征 ，针对 复 杂地质 构 造会 出现 盲 区 。为此 ，应 用 地震 波 的运 动学 和动力 学原 理 ，计算 给 定地质 模 型 的地震 响应 ，以准 确分 析 油气 藏地 震 波 的波场 特征 即三维 波动 方程 正演 ，是 地球
物理 学 家一 直努力 的 目标
・１９・ ９
三 维 波 动 方 程 正 演 技 术 在 地 震 采 集 施 工 方 法 设 计 优 化 中 的 应 用
成云 要苎魏 一，
，
燕 ，徐 继 伟
， 蔫 江石勘局球理探，苏州２７ 璧 曼 （ ～ 地物勘处江扬２０ 苏油探……一 … … ～ … ５） ０ 一～ … …
按 照反射 和 透射规 律 将波 场分 解 为上 行 波和 下行 波 ，即反 射 波和 透射 波 。该方 法有 两 个优 点 ，一是计 算 速度 较其 他数 值算 法 如有 限差 分法 、有 限元法 等 要快 ；二 是节 省 内存 ，对 于 二维 问题 ，只需 要若 干个 一
维数 组存 储 波场 ；对 于 三维 问题 ，只需要 若 干个 二维 数组 存储 波 场 ，即可 完成 正演 计算 ，正确模 拟 一次
第三章 三维波动方程的定解问题-2
t 0
2 F (r a t )
t 0
(ru ) r u r a t
t 0
( 3.35)
考虑到先前有一个动作 ：令 r 0 ，得到了（ 3.33 ）：
u (0, t ) 2F (a t )
另一方面，令 r 0 ，正是
u( M , t ) lim u ( r , t ) u (0, t )
深圳大学电子科学与技术学院
物理意义
1 1 1 u(M , t ) (M ' ) dS (M ' ) dS 2 2 4 a t t S M 4a t S M at at
T0 d
M'
D
M
at
S
M at
如果D< at，u(M,t) = 0 (扰动阵尾已过)
( 3.32)
上式的结果代入（ 3.30 ）式，得
u (0, t ) 2F (a t )
( 3.33)
u (r u ) u (r , t ) r F ( r a t ) G( r a t ) r r u r aF ( r a t ) aG( r a t ) t
d
, , 0
,
0 2 .
dS si n dd 球立体角元 r2
x
dS r 2 sin d d 球面上的面积元 dV dSdr r 2dr sin d d 球的体积元
dV dSdr r 2dr sin d d 球的体积元 r 2drd 球的体积元
因此，以下我们将先求 出 u (r , t ) 。 （这比求解 u( M , t ) 方便得多了） 。
波动方程正演模型及应用
波动方程正演模型及应用吴清岭 张 平 施泽龙3(大庆石油管理局勘探开发研究院)摘 要 地震资料解释经常用到正演模型。
常规的褶积模型不能模拟地震波的动力学特征。
本文采用声波方程,通过四阶有限差分近似,实现了复杂地质构造零炮检距的数值模拟。
文中同时展示了实际应用效果。
主题词 正演模型 有限差分 零炮检距剖面作者简介 吴清岭,男,1962年生,1983年毕业于华东石油学院勘探系,硕士,高级工程师,现从事地震方法研究工作。
地址:(163712)黑龙江省大庆市让胡路区勘探开发研究院。
3 参加本工作的还有杨有林同志。
在地震资料解释中,人们力图得到能够保持地震波的运动学与动力学特征的波动方程正演模型,以达到精确模拟地震波传播特性的目的。
在求解波动方程的2种数值解法(有限差分法和有限元法)中,有限差分法是一种快速有效的方法,并且地质模型的复杂程度不影响运算速度。
本文介绍了对声波方程采用四阶有限差分近似制作零炮检距剖面的基本过程及应用效果。
一、基本原理1,计算公式在二维空间域内,二维声波方程为1C 292u 9t 2=92u 9x 2+92u9z 2式中 C ———声学介质下地震波的纵波速度;u ———声压。
设Δh 为空间采样步长;Δt 为时间采样步长;m 、n 、l 分别为正整数;则有x =m ・Δh z =n ・Δh t =l ・Δt 对时间域采用二阶有限差分;对空间域采用四阶有限差分(推导过程略),其数值计算公式为u (m ,n ,l +1)=(A 2/12){16[u (m +1,n ,l )+u (m -1,n ,l )+u (m ,n +1,l )+u (m ,n -1,l )]-[u (m +2,n ,l )+u (m -2,n ,l )+u (m ,n +2,l )+u (m ,n -2,l )]}+(2-5A 2)[u (m ,n ,l )-u (m ,n ,l -1)]其中 A 2=C 2(m ,n )Δt 2/Δh 2式中 C (m ,n )———介质速度的空间离散值;Δt ———时间离散步长;Δh ———空间离散步长。
基于傅氏变换有限差分的三维波动方程正演模拟并行算法
Ａ一 ＋ ）Ａ一 一 ） （ 芝＋ 。 （一 毒＋ 专 ｚ 詈 ；
极大地提高了运算效率 ， 且易于在不同的设备上实
现， 为其 实用化 奠定 了 良好 的基础 。
薏嚣 ） １—． ＋ ６＋ （ ’
基金 项 目 ：国 家 ９ ３ 题 （０ ７ Ｂ ０ ６ ８ ； 球 探 测 与信 息技 术教 育 部 重 点 实验 室 开放 基金 （４ 。９ ８ ０ ３ ７课 ２０ Ｃ ２ ９０ ）地 ０１０ＯＯ５） 收稿 日期 ：２ １ 一Ｏ —２ ０１ ９ ６ 改 回 日期 ：２ １ ～Ｏ —２ Ｏ２ ４ ４
（． １地球 探 测与信 息技 术教 育部 重点 实验 室 成都理 工 大学 ， 都 成
８ １０ ） ４ ０ ０
６０ ５ ； １０ ９
２ ｑ 国石 油塔 里木 油 田公 司研 究院 ， 尔勒 ．－ 库
摘
要 ：傅 氏变换有 限差 分 （ Ｆ Ｄ） 即 Ｆ 混合 波 场延 拓 法 波动 方 程 正 演模 拟 精 度 高，＿￣ 维 Ｆ Ｄ ４ ｖ － Ｆ
因其 计算 复杂 、 据存 储 巨大而难 以推广 。这里提 出三 维 Ｆ Ｄ 一 步 法正 演模 拟 并 行 算 法 ， 一 数 Ｆ 把 个 串行 的大规模 频 率循环 分割 成 并行 的 多个小规 模 频 率循 环 ， 并发 执 行 多 个波 场延 拓 来提 高运
算速 度 ， 并通过 由给 定 的资 源设 定 并行数 ， 以及 适 时 重排 三 维数 据 ， 大大 增 强 了算 法 的适 应性 和 时效 性 。通 过 对 比并行和 非并 行 正演 实验 ， 证 了算 法的有 效性 。 验 关键 词 ：Ｆ Ｄ 混合延拓 ；三 雏数据 重排 ；三维 正 演模 拟 ；频 率循 环 ；并行 算法 Ｆ
波动方程正演在地震勘探设计中的应用_丁大伟
为 入 射 角 ，θ２为 透 射 角 ，ｖ１为 入 射 波 波 速 ，ｖ２为 透 射 波 波
速，Ｖ反为反射波单位向量，Ｖ透为透射波单位向量。其中，透
射角由Ｓｎｅｌｌ定律直接计算，公式如下：
θ２＝ｓｉｎ－（１
ｖ２ｓｉｎθ１） ｖ１
三、应用效果对比
（一）小断块发育的复杂构造区岩性勘探实例分析
某区块小断裂发育，地层破碎，储层砂体薄，地表过大 型障碍区，如何能够得到好的小断层，识别薄砂体，同时，
坠τｘｚ ＝μ（ 坠ｖｘ ＋ 坠ｖｚ ）
坠ｔ
坠ｚ 坠ｘ
其中：τｘｘ＝τｘ（ｘ ｘ，ｚ，ｔ），τｚｚ＝τｚ（ｚ ｘ，ｚ，ｔ），τｘｚ＝τｘ（ｚ ｘ，ｚ， ｔ）是应力张量；ρ＝ρ（ｘ，ｚ）是密度；ｖｘ＝ｖ（ｘ ｘ，ｚ，ｔ），ｖｚ＝ｖ（ｚ ｘ，ｚ， ｔ）是速度向量；λ＝λ（ｘ，ｚ），μ＝μ（ｘ，ｚ）是拉梅系数。
【教法研究】
波动方程正演在地震勘探设计中的应用
丁大伟１，柳 溪２，李 宏２
（１．河北涿州石油物探学校，河北 涿州 ０７２７５０； ２．华北油田地球物理勘探研究院，河北 任丘 ０６２５５２）
摘要：随着地震勘探程度的提高，勘探难度逐渐增加，勘探精度要求越来越高，针对地质目标的采集方法研究越来
越重要，利用现有的地震资料，结合地震反演建立合理的地质模型，运用波动方程正演技术设计地震勘探采集方法，研
坠 坠ｚ
（ ρ１
坠ρ 坠ｚ
）＝
１ Ｋ
坠２ρ 坠ｔ２
＋ｆ（ｘ，ｚ，ｔ）δ
（ｘ－ｘ０，ｚ－ｚ０）
在三维空间域内，声波方程为：
坠坠ｘ（ρ１ 坠坠ρｘ ）＋ 坠坠ｙ（ρ１ 坠坠ρｙ ）＋ 坠坠ｚ（ρ１
坠坠ρｚ ）＝
１ Ｋ
坠２ρ 坠ｔ２
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
三维波动方程正演及模型应用研究1熊晓军，贺振华，黄德济成都理工大学油气藏地质与开发工程国家重点实验室（610059）E-mail:xiongxiaojun@摘 要：为了真实准确地反映三维地质体的波场特征，在频率-波数域将二维波场延拓算子推广到三维空间，采用三维波动方程延拓方法实现了三维地质模型的快速叠后正演。
该方法采用傅立叶变换进行计算不仅计算迅速，算法稳定，而且可以采用相位移加插值方法处理一定的横向变速的情况，为更加灵活方便地模拟地下复杂的三维地质体提供了一种有效的工具。
作为实例，首先进行了三维French 模型的数值模拟，得到了和实际物理模型实验结果相一致的正演记录，并对比分析了三维偏移剖面和二维偏移剖面的偏移效果，验证了该方法的正确性和适用性。
然后进行了三维缝洞地质模型的正演计算，得到了高信噪比的正演记录，对认识和解释三维缝洞地质体的地震波场特征很有帮助。
关键词：三维地震正演，波动方程，波场延拓，French 模型，缝洞模拟1. 引 言在地震资料采集、处理和解释中通常要进行地震波场的数值模拟，即假设已知地下的地质情况，应用地震波的运动学和动力学的基本原理，计算给定地质模型的地震响应，其对正确认识地震波的运动学和动力学特征，以及准确地分析油气藏的反射波场特征有着重要的指导意义。
常规的正演方法主要有波动方程法和几何射线法两大类[1]。
几何射线法[2]将地震波波动理论简化为射线理论，主要考虑地震波传播的运动学特征，所得的地震波的传播时间比较准确，但是计算结果很难保持地震波的动力学特征，而且对复杂的地质构造会出现盲区。
波动方程法通过求解地震波波动方程，建立地下地震波运动的波场，它包含了地震波传播的所有信息，在地震模拟中占有重要地位。
目前，常规的波动方程正演方法，如有限差分法[3]、积分法[4]以及F-K 域[5]等方法，主要进行二维地质体的数值模拟，而实际的地下构造往往是三维的，因此，十分有必要研究基于三维地质体的数值模拟方法。
本文在二维正演方法的基础上，将二维波场延拓算子推广到三维空间，在频率—波数域进行三维叠后正演模拟。
该方法不仅算法稳定，计算速度快，具有处理一定横向变速的能力，而且可以得到高信噪比的零偏移距正演记录。
2. 方法原理利用波动方程进行地震波场数值模拟的核心是波场延拓，对于垂向变速介质，利用二维标量波动方程，在频率—波数域可以得到各个深度间隔内的相位移延拓的正演和偏移公式[6]， i zi Z ik i x i x e z k P z k P ∆+=),,(),,(1ωω (1)i zi Z ik i x i x e z k P z k P ∆−+=),,(),,(1ωω (2)1 本课题得到高等学校博士学科点专项科研基金资助课题（SRFDP ）资助，编号：20040616001。
- 1 -)1(22222ωωv k v k x zi −= (3)其中式(1)为二维正演延拓公式，其延拓方向为由下至上，式(2)为二维偏移延拓公式，其延拓方向为由地面向下延拓。
式中表示x k x 方向的波数，表示深度间隔i 内在方向的波数，zi k z ω表示频率，表示当前深度延拓的深度间隔。
z ∆ 本文在上面的二维波场延拓公式的基础上，采用三维波动方程进行计算，将二维波场延拓算子推广到三维空间，得到了以下的垂向变速的三维地质体的正演和偏移成像的波场延拓公式，i zi Z ik i y x i y x e z k k P z k k P ∆+='),,,(),,,(1ωω (4)i zi Z ik i y x i y x e z k k P z k k P ∆−+='),,,(),,,(1ωω (5)))(1(222222'ωωv k k v k y x zi +−= (6)其中式(4)为三维正演延拓公式，(5)为三维偏移延拓公式。
式中表示方向的波数，表示深度间隔i 内在方向的波数。
y k y 'zi k z 当速度横向变换时，本文采用了Gazdag 等人提出的相位移加插值方法[7]的思想，在三维空间选取三个参考速度进行拉格朗日插值计算，这样使该方法具有了更大的适用性。
3. 理论模型计算3.1三维French 模型为了检验本文提出的方法的正确性，我们进行了三维French 模型的数值模拟计算。
为了与原物理模型（参见附图1， French ，1974）有所区别，在原物理模型的上部施加了一个水平层。
图1（a ）是本文采用的French 模型和自激自收的观测系统。
该模型共包含6个地质体，“1”代表自上往下的第一个水平界面，“2”代表自上往下的第二个水平界面，“3”代表自上往下的第三个水平界面，“4”代表斜面，“5”和“6”分别代表两个大小和形态都相同的穹隆。
本次采用的自激自收的观测系统共有128条测线，每条测线包含128个检波点。
图1（b ）是去除第一个水平界面后的平面图，在该图上添加了观测系统的线道号。
图1（c ）是相应的速度模型，为一沿x 方向的垂直切片。
对于穹隆5和6，它们的速度都是3000m/s. 图2是第39号测线的时间剖面图，图中仅包含有一次反射波和绕射波。
从图2（a ）可以看到，似乎L39通过了两个穹隆的近顶部，但是从图1（b ）可以看到，L39仅通过代号为“6”的穹隆。
图2（b ）是该条测线的三维偏移剖面，其实现了正确的偏移归位，这也证明了本文提出的三维波动方程正演方法的正确性。
此外，为了说明三维偏移的效果，我们还对L39的正演记录进行了二维偏移，如图2（c ）所示，图中显示了两个穹隆，与实际不符合，这就造成了解释陷阱。
- 2 -（a）French模型和观测系统（b）去除水平界面1的平面图（c）速度模型，其中穹隆5和6的速度都为3000m/s图1 用于计算的三维模型（a）L39的三维正演记录- 3 -（b）L39的频率—波数域三维波动方程偏移记录（c）L39的频率—波数域二维波动方程偏移记录图2 Line 39的时间剖面图3.2三维French模型（a）三维洞穴模型（b）经过洞穴的水平切片（c）速度模型，其中洞穴体的速度为1000m/s图3 简单的三维洞穴模型图3是一个简化的三维缝洞模型，它采用和图1（a）相同的自激自收观测系统，共有128条测线，每条测线包含128个检波点。
图3（a）由三个地质体组成，“1”代表至上往下的- 4 -第一个水平层界面，“2”代表自上往下的第二个水平层界面，“3”代表位于两个水平层界面之间的洞穴体。
图3（b）是一个经过洞穴体的水平切片，并在该图上标识了观测系统的线道号。
图3（c）是与该三维缝洞模型对应的速度模型，其中洞穴体的速度为1000m/s。
图4（a）是第64号测线的三维正演记录，从图中可以清楚的观察洞穴引起的绕射波和反射波的特征及其影响范围，如CDP45~CDP85之间的圆弧所示。
此外，图中CDP64附近出现了多个强的绕射同相轴，可以解释为洞穴体的其它绕射点综合影响的结果。
图4（b）是第64号测线的三维偏移记录，图中的绕射全部正确归位，缝洞体的形态十分清晰。
图4（c）是第64号测线的反射系数剖面。
通过对比图4（b）和图4(c)，我们验证了本文采用的三维正演和偏移方法的正确性，同时也说明采用该方法进行缝洞地质体识别和研究的可行性。
（a）L64的三维正演记录（b）L64的频率—波数域三维波动方程偏移记录（c）L64的反射系数剖面图4 Line 64的时间剖面图4. 结束语通过理论模型计算的结果表明，本文采用的三维波动方程正演方法是一种切实可行的数- 5 -值模拟方法。
该方法不仅计算速度快，算法稳定，而且可以处理一定的横向变速的情况，具有很大的适用性。
其得到的正演记录仅包含了地震波的反射和绕射特征，信噪比高，而且能够通过常规的偏移方法实现其正确的偏移归位，为今后进一步研究缝洞储层的地震波场响应以及检测和识别缝洞地质体提供了一种新的思路和有效的工具。
参考文献[1] 张永刚. 地震波场数值模拟方法.石油物探，42： 143~148. 2003[2] 张钋，刘洪，李幼铭.射线追踪方法的发展现状.地球物理学进展， 15：36~45. 2000[3] 陈伟. 起伏地表条件下二维地震波场的数值模拟.勘探地球物理进展，42： 26~31. 2005[4] 贺振华，《反射地震资料偏移处理与反演方法》，重庆：重庆大学出版社,1989.5[5] 熊小兵，贺振华.相位移加有限差分法波动方程正演模拟.石油物探，37：22~28. 1998[6] Lai J, Gardner G H. Forward modeling using the PSPI method. Seismic acoustic laboratory, Annual progressreview, 14: 271~300, 1984[7] Gazdag J F, Sguazzero P. Migration of seismic data by phase shift plus interpolation. Geophysics, 49:124~131,19823-D wave equation forward modelingand its application of modelsXiaoJun XIONG ZhenHua HE DeJi HUANG State Key Lab. of Oil and Gas Reservoir Geology and Exploration, Chengdu University ofTechnology, Chengdu 610059, ChinaAbstractTo correctly study the characters of 3D seismic wave-field, the paper extends 2D wave-field extrapolating operator to 3D space in F-K domain, and applies 3D wave equation to realize fast post-stack forward modeling of 3D geologic body. The method has fast calculating efficiency by FFT, stabile algorithm and can handle variable velocity in landscape orientation by PSPI, so it provides a new effective tool for conveniently simulating complex 3D underground geologic body. As the examples, the paper firstly simulates 3D French model, which can get the same results as practical physical model, and tests its correctness and applicability by contrastively analyzing the migrated profiles of 3D and 2D. Then the paper simulates the 3-D geological fractured model and gets its seismic profiles with high S/N, which is helpful for us to recognize and explain the characters of the model’s seismic wave-field.Keywords: 3-D seismic forward modeling, wave equation, wave-field extrapolation, fractured simulation, French mode作者简介：熊晓军，男。