2008年12
月第43卷 第6期
・正演技术・
一阶弹性波方程的变网格高阶有限差分数值模拟
李振春3 张 慧 张 华
(①中国石油大学地球资源与信息学院, 山东东营257061)
李振春, 张慧, 张华. 一阶弹性波方程的变网格高阶有限差分数值模拟. 石油地球物理勘探, 2008, 43(6) :
711~716
摘要 使用可变网格的有限差分法进行地震模拟有许多独特的优点, 主要表现为对地质模型的离散化更为合理, 在低速带和复杂构造区域, 可将局部网格划分得相对精细些, 不仅提高了模拟精度, 消除了因采样不足导致的频散现象, 而且可以减少计算机内存需求, 保持模型计算的灵活性。本文提出一种新的基于高阶交错网格技术的弹性波数值模拟方法, 通过改变网格的空间步长实现了局部网格加密技术, 弥补了常规网格的缺陷和不足。试算结果表明, 本文提出的算法稳定性较好, 且能够提高模拟精度, 减少计算时间, 提高计算效率。关键词 变网格 交错网格 高阶有限差分 弹性波
1 引言
法技术(地震采集。为适应复杂勘探条件的需要, 、高效率的地震数值模拟技术变得越来越重要。
自Virieux [1]率先采用交错网格有限差分法模拟P 2SV 波在非均匀介质中的传播以来, 交错网格成为了有限差分法正演模拟中普遍使用的网格剖分技术。同二阶位移方程相比, 使用交错网格的一阶速度—应力方程, 无论是模拟精度还是模拟的稳定性都有了较大提高。但是如果局部地区构造较为复杂, 或者存在低速层时, 为了提高复杂构造的正演精度, 通常对整个模型使用较小的采样网格, 以保证对频散现象的抑制和对稳定性的保障。这对保持模型计算的灵活性非常不利, 还会对模型的其他部分造成过采样, 增加惊人的计算量。所以, 设计一种更合理的离散化方式进行正演模拟是非常必要的。
从空间采样的角度考虑, 既要提高模拟精度、又要降低计算机内存需求的最有效方法就是在模型的不同区域采用不同步长的网格, 即变网格(也称为不规则网格) 。为此J ast ram 等提出了基于二维声波方程针对某一深度变网格步长的算法[2];J ast ram 等
[3]; [4Pitarka 实现了低阶精度的不[5]。
本文在前人研究的基础上提出一种高阶变网格数值模拟方法, 实现了高精度的局部网格加密技术, 另外还改进了前人对变化网格过渡带采用的波场插值算法, 避免了在过渡带区域进行大量的插值计算, 也不需要平滑函数进行平滑, 同时还保持了算法的稳定性。
2 方法原理
文中所提出的算法是基于Virieux [1]和Graves [6]
提出的一阶速度—应力方程, 对于各向同性介质, 有
τxx +5τ5t v x =b (5x z xz ) +f x
τxz +5τ5t v z =b (5x z zz ) +f z λ+2μ) 5x v x +λ5τ5z v z t xx =(
λ+2μ) 5z v z +λ5τ5x v x t zz =(
(5z v x +5x v z ) 5τt xz =μ
(1)
ττ式中:v x , v z 分别表示速度; τxx , zz , xz 分别表示应力向量; 5x 、5z 、5t 分别表示一阶偏微分算子5/5x 、5/5z 、5/5t ; b 表示密度的倒数; f x , f z 为体力向量;
λ及μ为拉梅常数。
3山东省东营市中国石油大学地球资源与信息学院,257061
本文于2008年1月8日收到, 修改稿于同年7月30日收到。本课题由国家973专题(2007CB209605) 、国家863专题(2006AA06Z206) 和CNPC 物探重点实验室中国石油大学(华东) 研究室资助。
712 石油地球物理勘探2008年
通过对时间导数采用二阶近似, 得到弹性波一阶速度—应力方程的离散形式为
v x v
n+2
i , j
=v x
n-2
i , j
τxx +D τ+Δtb x (D x z xz ) |n
i , j
n
Δ 1=d x i /2, Δ2=d x i /2
n-1
n-1k =1
n+2
z i+2, j+2
=v
n-2
z i+2, j+2
τ+Δtb z (D x xz +D τz zz ) |i+, j+22
n+2i+2, j n+2i+2, j
Δn-1=
k =1
d x i+k+d x i /2, Δn =∑d x i-k +d x i /2∑
ττ
n+1
xx i +2, j n+1
zz i +2, j
=τ
n
xx i +2, j n
zz i +2, j
λ+2μ) D x v x +λ+Δt[(D z v z ]
图1 计算对称点为i +
的差分系数所需的变网格节点2
n-1
=τ) D z v z +λ+Δt[(λ+2μD x v x ] 同理, 以i 为对称点的差分算子为Δn-1=
n-1k =1
n+1n
τμxz (D z v x +D x v z ) ]xz i , j +2=τxz i , j +2+Δt[(
n+2
i , j +2
d x ∑
i+(k-1)
+d x i+n-1/2, Δn =
k =1
d x ∑
i-k
+d x i-n /2。
(2)
式中:D x , D z 分别表示对x , z 的一阶微分算子;
b x =(b i , j +b i +1, j ) /2, b z =(b i , j +b i , j +1) /2; μxz =(1/μi , j +1/μi +1, j +1/μi , j +1+1/μi +1, j +1) /4。
用变量g 来表示速度v x 、v z 或应力向量τxx 、ττzz 、xz , 则在常步长的交错网格技术中, 用下式计算式
(2) 中的g 关于x 的2n 阶空间偏导数[7]
D x g (x , z ) =
C [g (x +Δx
Δx ∑
n i
i =1
下面给出变网格算法的任意偶数阶精度差分近似式及差分系数的求取方法, 即
n
D x g (x , z ) =
i -1
[c ∑
2i-1
g (x +Δ2i-1, z ) +
(4)
+c 2i -Δ2i , z ) ]
n
, ) -
i ; i 是空间差分算子, 它
d g 关于x 。令g (x , z ) =g z exp (i k x ) , 则式(4) 可写成
n
(3) -g (x -Δx i-1z x 2n 式中, 系数C i n 由x (x 2i -1, z ) 和g (x -Δx 2i -1, z ) 。
与常网格不同, 变网格的差分算子是空间变化
i k =
i -1
[c ∑
2i-1
Δ2i-1) +c 2i exp (-i k Δ2i ) ]exp (+i k
(5)
的。给定地质模型, 在确定网格剖分后就可求出各
网格节点的差分算子供递推计算时调用。在交错网
格技术中差分算子有两个对称点, 即i +或i , 不
2
同的对称点对应的Δi 也不同。下面以对称点i +
2
为例求取Δi , 如图1所示。
对式(5) 中的指数项进行泰勒展开, 它的2n 阶泰勒展开式为
2
Δi ) ≈1+i k Δi +(i k ) 2Δexp (i k i +…+
2
2n-12n
(k Δi ) 2n-1+o (Δ(6) +i )
(2n -1) ! 把式(6) 代入式(5) 并整理, 可以得到
2222
ΔΔ) () (c 1ΔΔ2 i k =(c 1+c 2+…+c 2n ) +i k (c 1Δ-c +…-c +i 1222n 2n 1+c 2Δ2+…+c 2n 2n ) +
2
+(i ) 3
3
6
332n-1
(c 1ΔΔ31+c 2Δ2+…-c 2n 2n ) +…+(i )
2
n-12n-12n-1
(c 1Δ-c 2Δ+…-c 2
n Δ2n ) 12
(2n -1) !
(7)
2n-1
写成矩阵的形式
1i Δ1
1
-i Δ2
…
…ω
1i Δ2n-1
1
-i Δ2n
c 1 c 2 …=
2n-22n-1
0i
…
(i ) 2n-2Δ1
2n-2
2n-1
(i ) 2n-1Δ1
…
2n-2
(i ) 2n-2Δ22n-12n-1
Δ-(i ) 2
……
(i ) 2n-2Δ2n
2n-1
Δ-(i ) 2n
2n-2
…(i ) 2n-2Δ2n-12n-1…(i ) 2n-1Δ2n-1
c 2n-1c 2n
00
(8)
第43卷 第6期李振春等:一阶弹性波方程的变网格高阶有限差分数值模拟 713
通过解方程组(8) , 可以得到有限差分算子D x 的系数c i 。同样地, 差分算子D z 的差分系数也是通过这种方法求得。对于交错的网格, 由于变量所定义的网格点的位置不同, 对于每一个差分算子D x 或D z 都有两组不同的差分系数。
在使用变网格时, 网格步长的变化可能会导致数值反射的出现。究其原因是因为波场离散后, 相速度是网格步长的函数。当相速度梯度较大时, 即使速度和密度都没有变化, 入射波的能量也会部分反射回来, 导致了数值反射现象。因为频散程度与频率有关, 所以当入射波的频率较低时, 产生的数值反射较高频时弱。在网格变化的区域对波场进行插值计算虽可以在某种程度上压制数值反射, 但是插值算法计算量大, 且效果并不令人满意, 所以本文没有采用插值计算。本文提出的算法是根据各点所对应的步长的变化计算其差分系数, 在一个循环内完成所有点的波场计算, 不会产生数值反射, 并且减小了计算量。
式中:h =max (d x i , d z i ) 表示最大的空间网格步长;
v min 是地震波速的最小值。
本文所提算法需要满足的稳定性条件为
≤
min (d x i , d z i ) n i =1
∑
(-1) i-1C i
式中:V max 为最大速度值; ρ为密度; C i 为式(3) 中做2n 阶空间展开时的系数。
4 数值模拟实例
为了验证提出的变网格差分算法的稳定性及有效性, 设计两个模型正演实例来说明文中所述方法的实际效果。4. 1 。模×3000m , 纵波速度为3500m/s , 横0125来计算, 密度为常数。P 波震源激发, 震源主频为30Hz , 位于(1500m ,1500m ) 处。变网格步长为10m ,5m 。与其进行对比的常规交错网格步长为5m 。数值模拟中, 差分解法的精度为o (Δt 2, Δx 10
) 。弹性波波场快照和单炮记录分别记录于图2和图3中。
从上面的弹性波波场快照和单炮记录中可以看出, 变网格的过渡区没有产生数值边界反射, 说明本文提出的算法成功地消除了空间网格变化造成的数值反射, 进而说明了该算法的合理性和良好的稳定性。
3 法较好地解决了边界反射问题。使用透射边界能够吸收绝大部分边界反射的能量, 剩余能量通过衰减边界再次减小, 这样计算量增加不大, 吸收也较彻底, 可以较好地消除人为边界反射。
稳定性条件是有限差分数值模拟中十分重要的问题[9]。经过验证, 利用本文所提出的算法能够精确模拟的最大的地震波频率为
f max
图2a 变网格波场(左) 和常规网格波场(右) 水平分量快照(t =360ms )
714 石油地球物理勘探2008年
第43卷 第6期李振春等:一阶弹性波方程的变网格高阶有限差分数值模拟 715
412 背斜模型
模型大小为3000m ×2000m , 模型中横波速度取泊松比为0. 25来计算; P 波震源激发, 震源子波主频为20Hz 。采用精度为o (Δt 2, Δx 10) 的常规交错网格有限差分法和变
网格有限差分法进行模拟对比。模型参数及单炮记录如图4和图5所示。
从上面的模拟效果对比分析得知:在用常规交错网格有限差分法模拟的炮记录中, 由于背斜的影响, 自背斜顶界面反射以下出现了严重的数值频散; 而用变网格算法进行模拟时, 由于在背斜低速区域采用小网格步长计算, 故由背斜影响导致的频散得到了较好地压制, 充分显示了此算法较高的模拟精度
。
图4
背斜模型
716 石油地球物理勘探2008年
5 结论
(1) 本文在网格过渡区域, 没有使用插值计算, 而
是根据各点所对应的步长的变化计算其差分系数, 在一个循环内完成对所有点的波场计算, 成功地消除了数值反射, 且同波场插值计算相比减少了计算量。
(2) 对于常规网格算法难以解决的地下复杂构造的问题, 采用变网格算法可较好地解决, 且能较灵活地适应地层模型的变化, 对低速体的适应性较常规网格算法有了较大改善。另外, 相对常规网格算法而言, 变网格算法节省了计算机内存, 减少了计算时间, 且模拟的区域越大, 节省的时间就会越多, 避免了对整个模型的精细采样造成的资源浪费。
(3) 本文提出的算法精度可达任意偶数阶精度, 高阶变网格算法的提出, 尤其是与交错网格有限差分法相结合, 可以有效地提高模拟精度和计算效率。
参考文献
[1] Virieux J. P 2media :velocity 2stress finite 2difference method. Geo 2
physics , 1986, 51, 88~901[2] J astram C and Behle A. Acoustic modeling on a verti 2
cally verying grid. Geop hys Pros p , 1992, 40, 157~169[3] J astram C and Tessmer E. Elastic modeling on a grid of
vertically varying spacing. Geop hys Pros p , 1994, 42, 357~370
[4] Falk J , Tessmer E and Gajewski D. Tube wave model 2
ing by the finite 2difference method with varying grid spacing. Pros p , 1996, 148, 77~93
[5] Pitarka A. 3D elastic finite 2difference modeling of seis 2
mic motion using staggered grids with nonuniform spacing. B ulletin of Seismological of A merica , 1999, 89, 54~68
[6] Graves R W. Simulating seismic wave propagation
in 3D elastic using staggered 2grid finite 2differences. B ull Seism S oc A m , 1996, 86, 1091~1107[7] 董良国, 马在田, 曹景忠. 一阶弹性波方程交错网格高
阶差分解法. 地球物理学报,2000, 43(3) :411~419[8] 廖振鹏, 周正华, 张艳红. ,2002,45(4) :534~545[ , 地球物理学报,2000, 43(6) :(本文编辑:冯杏芝)
(上接第688页)
虚线框在横向上的分布与图8中虚线框在横向上的分布比较可以得知, 阻抗低值处并非都富集天然气。通过时频分析技术和阻抗反演技术相结合, 最终确定储层的分布范围为图8中白色实线框内。
从对实际资料的应用和分析中发现, 将时频分析技术和波阻抗反演技术相结合可以更加精细地刻画礁滩相储集层的空间分布、油气的富集范围, 计算油气储层的厚度, 估算油气的储量, 从而有效地对礁滩相储层进行预测, 降低钻探风险。
只对垂直剖面进行了探讨, 建议在实际的生产研究中也应该考虑对沿层切片、水平时间切片及波阻抗剖面和单频剖面叠合图进行分析, 从而更加准确地对礁滩相油气藏进行预测。
参考文献
[1] Thomas Englehart ; 石天译. 墨西哥湾中部气云四分量
地震数据的解释. 石油勘探开发情报,2001(6) :33~43[2] 李等华, 唐跃, 殷积峰等. 川东黄龙场构造上二叠统长
兴组生物礁特征与潜伏礁预测. 中国地质, 2006, 33(2) :427~434[3] Castagna J P ,Sun Shengjie ,Siegf ried R W. Instantane 2
ous spectral analysis :Detection of low 2f requency shad 2ows associated with hydrocarbons. T he L eading Ed ge , 2003, 120~127[4] 董臣强, 王军, 张金伟. 时频分析技术在三角洲层序分
析中的应用. 断块油气,2002,9(2) :18~20[5] 张宏彬, 何樵登等. 宽带约束反演. 石油物探, 1995,
34(1) :1~10[6] 张霖斌, 姚振兴等. 快速模拟退火算法及应用. 石油地
球物理勘探,1997,32(5) :654~659[7] 王赟, 王妙月, 彭苏萍. 地球物理随机联合反演. 地球物
理学报,1999,42(Supp. 1) :142~151[8] 寇晓丽, 刘三阳. 基于模拟退火的粒子群算法求解约
束优化问题. 吉林大学学报(工学版) , 2007, 37(1) :136~140
5 结束语
通过对礁滩相油气储集层的特征分析和研究, 提出了将时频分析技术和波阻抗反演技术相结合对礁滩相油气藏进行预测的方法。文中分析了单独应用两种技术对礁滩相油气藏预测的优缺点。综合两种技术的优点, 可以更加有效地对礁滩相油气藏进行预测。通过对实际资料的应用发现, 采用文中提出的礁滩相预测研究方法和思路, 可以有效识别储集层的位置, 精细刻画储集层的范围, 计算储集层的厚度, 减少钻井的风险, 更好地指导油气开发。文中
(本文编辑:冯小球)
Vol. 43 No. 6Abstracts Ⅶ
V ariable 2grid high 2order f inite 2difference numeric simulation of f irst 2order elastic w ave equation. Li Zhen 2chun 1, Zhang H ui 1and Zhang H u a 1. OGP , 2008,43(6) :711~716
U sing variable 2grid finite 2difference algorit hm to carry out t he seismic simulation has many u 2nique advantages t hat mainly are more reasonable for discretization of geologic model , and relatively fine division of local grid can be implemented in low 2velocity zone and complex st ruct ural area , which not only improves t he simulated p recision and eliminates t he dispersion resulted f rom insuffi 2cient sampling , but also can reduce t he require 2ment of comp uter memory , keeping t he flexibility of model comp utatio n. The paper presented a new elastic wave numeric simulatio n met hod based on high 2order staggered 2gird technique , which real 2ized t he technique densifying local grid by changing t he spatial step of grid , remedying t shortcoming of ordinary showed t hat t he t he paper is characterized by , improved p reci 2sion of simulation , reducing co mp utational time and improving comp utational efficiency.
K ey w ords :variable grid , staggered grid , high 2or 2der finite difference , elastic wave
1. College of Eart h Resource and Informatio n , Chi 2na University of Pet roleum , Dongying City , Shan 2do ng Province , 257061, China
Elastic w ave simulation of P 2and S 2w aves separation by 22D staggered grid and application. Zhang Jian 2lei 1, Tian Zhen 2ping 1and W ang Cheng 2xiang 1. OGP , 2008,43(6) :717~722
Through systematically st udying on t he met h 2od separating P 2wave from S 2wave by second 2order elastic wave equation and first 2order staggered grid algorit hm , t he paper presented t he formula rep re 2sentation simulating P 2and S 2waves separation in first 2order staggered grid. The P 2and S 2waves are nat urally separated but do not be coupled toget her in a simulating p rocess of t he met hod , p roducing f ull P 2wave and f ull S 2wave simulated records. The result s of numeric simulation showed t he correct 2ness of t he met hod. There is important p ractical
meaning for knowing p ropagation rule of elastic wave and guiding numeric simulation of reservoir by st udy of f ull P 2wave and f ull S 2wave simulation. K ey w ords :elastic wave simulation , P 2and S 2wave separatio n , staggered grid
1, Technical Research Center of Geop hysical Pros 2pecting , B GP , Zhuozhou City , Hebei Province , 072751, China
Staggered 2grid pseudo 2spectrum simulation of frac 2ture and cave reservoir. Liu Jiong , Ma Jian 2w ei 1, Yang H ui 2zhu 1and B a Jing 1. OGP , 2008, 43(6) :723~727
The paper carried out t he analysis and correla 2tion of t raditional p seudo 2spect rum and improved staggered 2grid p seudo 2spect rum f rom t heoretical met hod cases , t he superiority of 2met hod. The more of differential localization of p seudo 2spect rum is showed by first 2order differential comp utational case of p ulse se 2quence ; t he correct ness of staggered 2grid p seudo 2spect rum met hod is proved by 22D acoustic wave comp utational case ; t he precision of t he staggered 2grid p seudo 2spect rum met hod being higher t han t hat of traditional p seudo 2spect rum is proved by t he co mp utational case of elastic wave in non 2hom 2ogeneo us medium. Finally , t he propagation issue of elastic wave in non 2homogeneous f ract ure and cave reservoir is simulated by staggered 2grid p seu 2do 2spect rum met hod.
K ey w ords :numeric simulation of wave equation , t raditional p seudo 2spect rum , staggered 2grid p seu 2do 2spect rum , fract ure and cave reservoir
1. Instit ute of Seismic Wave Exploration and De 2velop ment , Tsinghua University , Beijing City , 100084, China
Method for identifying fluid property in volcanite reservoir. Zhang Li 2hu a 1, P an B ao 2zhi 1, Shan G ang 2yi 1and H e Li 2. OGP , 2008,43(6) :728~730
The vulcanite reservoir is characterized by complex lit hology and variety of pore st ruct ure , which brought about difficulty in identifying and e 2valuating t he lit hology of volcanite reservoir. Tak 2ing deep volcanite reservoir in Daqing Oilfield as
2008年12
月第43卷 第6期
作者介绍
吕公河 教授级高工,1963年生; 1984年毕业于长春地质
学院物探专业,1996年获得青岛海洋大学海洋地质专业硕士学位,2005年毕业于西安交通大学工程力学专业, 获得博士学位; 现任胜利油田地球物理勘探高级专家, 长期从事滩海以及高分辨率地震勘探技术研究工作。陈浩林 高级工程师,1969年生;1991年毕业于中国地质
大学(武汉) 地球物理专业, 获学士学位;2001年毕业于中国石油大学(北京) 地球物理专业, 获工程硕士学位; 2005年毕业于中国矿业大学(北京) 地球物理专业, 获博士学位。长期从事地震资料采集、处理方法及理论研究。王振华 硕士,1958年生;1982年毕业于华东石油学院勘探
地球物理专业,1985年毕业于华东石油学院北京研究生部应用地球物理专业。主要从事地球物理勘探技术的开发与应用。现在北京华阳捷泰技术有限公司工作。单联瑜 高级工程师, 1959年生;1982年毕业于华东石油
学院地球物理勘探专业, 究和地震资料处理研究所在读博士王玉贵 高级工程师; 年毕业于石油大学
(华东) ,2006年获中国石油大学应用地球物理专业硕士学位, 现为中国地质大学地球物理与信息技术学院在读博士。一直从事VSP 资料处理、解释及方法研究工作。裴正林 博士, 副研究员,1962年生;1984年毕业于西安地
质学院物探专业,2000年毕业于中国地质大学(北京) , 获博士学位,2000~2003年在中国石油大学(北京) 从事博士后研究。主要从事地震波传播理论及数值模拟, 井间地震层析成像, 小波域地震属性提取和油气识别, 以及小波理论应用研究。获省部级科技进步一等奖1项, 二等奖1项, 三等奖2项。合著1部, 发表论文50余篇。段云卿 副教授,1958年生;1982年毕业于华东石油学院物
探专业, 现在中国地质大学(北京) 从事教学和地震资料处理及解释方法研究。已发表论文多篇。
何兵寿 副教授,1973年生;1996年毕业于中国矿业大学
(徐州) 物探专业,2002年获石油大学(北京) 应用地球物理专业博士学位, 现在中国海洋大学海洋地球科学学院从事多波多分量地震技术、纵波高分辨率处理技术研究和相关教学工作。樊德华 高级工程师,1967年生;1991年毕业于成都地质学
院, 现在中石化胜利油田从事石油勘探工作。刘振武 教授, 英国Heriot 2Watt 大学石油工程专业博士,
1949年生; 中国石油天然气集团公司副总工程师, 享受国务院政府专家特殊津贴。编著有《二十一世纪初中国
油气关键技术》《、油气应用基础研究展望》和《企业技术
创新管理》《、中国石油“十五”科技进展丛书》等著作。陈 林 博士研究生,1980年生;2004年毕业于江汉石油学
院, 获学士学位; 现为中国科学院地质与地球物理研究所博士研究生; 主要从事海洋地球物理和盆地模拟等方面的研究。赵丽平 工程师,1962年生;1988年毕业于胜利油田职工大
学石油地质专业,2004年毕业于石油大学资源勘察专业, 长期从事地震资料综合研究。蔡涵鹏 1983生; 现为成都理工大学信息工程学院硕士研
究生, 就读专业为地球探测与信息技术; 研究方向为地震资料处理与解释。彭文绪 高级地质师,1972年生;1995年毕业于江汉石油学
院, 获学士学位。1998年毕业于石油大学(北京) , 获硕。。 年生;1988年毕业于石油
(, 获学士学位,2005年在中国科学院广州地球化学所获构造地质博士学位。现在主要从事油气勘探综合研究工作。张海燕 副教授,1968年生;1990年毕业于清华大学精密
仪器专业, 获工学学士学位;1993年毕业于中国海洋大学海洋物理专业, 获理学硕士学位;2008年1月毕业于中国海洋大学海洋地质专业, 获理学博士学位。现在中国海洋大学信息科学与工程学院从事AVO 叠前反演、计算智能研究。已经在国内外核心学术期刊上发表论文十余篇, 其中有多篇被EI 收录。李振春 理学博士, 教授, 博士生导师; 中国石油大学(华东)
地球物理系主任。主要从事地震波传播与正演模拟、地震成像与偏移速度分析、多尺度地震资料联合反演与CFP 2AV P 分析理论与方法的教学与研究工作。张建磊 硕士,1978年生;2001年毕业于石油大学(华东) 地
球探测与信息技术学院,2004年获得石油大学(华东) 地球探测与信息技术专业硕士学位, 现在东方地球物理公司物探技术研究中心从事地震正演和偏移方法研究。刘 炯 1980年生, 清华大学在读博士, 主要从事复杂介质
中地震波传播机理的研究。张丽华 工程师,1974年生;1996年毕业于长春地质学院,
现为吉林大学地球探测科学与技术专业博士研究生, 研究方向为石油地球物理测井数据处理与解释研究。柳永坡 博士生, 副教授,1971年生;1995年毕业于大庆石
油学院, 主要研究领域为:软件工程、测试技术、性能故障分析。王 月 1985年生; 现在中国石油大学(北京) 攻读硕士学
位, 主要从事石油物探市场分析及国际石油经济等方面的研究, 已发表论文4篇。
2008年12
月第43卷 第6期
・正演技术・
一阶弹性波方程的变网格高阶有限差分数值模拟
李振春3 张 慧 张 华
(①中国石油大学地球资源与信息学院, 山东东营257061)
李振春, 张慧, 张华. 一阶弹性波方程的变网格高阶有限差分数值模拟. 石油地球物理勘探, 2008, 43(6) :
711~716
摘要 使用可变网格的有限差分法进行地震模拟有许多独特的优点, 主要表现为对地质模型的离散化更为合理, 在低速带和复杂构造区域, 可将局部网格划分得相对精细些, 不仅提高了模拟精度, 消除了因采样不足导致的频散现象, 而且可以减少计算机内存需求, 保持模型计算的灵活性。本文提出一种新的基于高阶交错网格技术的弹性波数值模拟方法, 通过改变网格的空间步长实现了局部网格加密技术, 弥补了常规网格的缺陷和不足。试算结果表明, 本文提出的算法稳定性较好, 且能够提高模拟精度, 减少计算时间, 提高计算效率。关键词 变网格 交错网格 高阶有限差分 弹性波
1 引言
法技术(地震采集。为适应复杂勘探条件的需要, 、高效率的地震数值模拟技术变得越来越重要。
自Virieux [1]率先采用交错网格有限差分法模拟P 2SV 波在非均匀介质中的传播以来, 交错网格成为了有限差分法正演模拟中普遍使用的网格剖分技术。同二阶位移方程相比, 使用交错网格的一阶速度—应力方程, 无论是模拟精度还是模拟的稳定性都有了较大提高。但是如果局部地区构造较为复杂, 或者存在低速层时, 为了提高复杂构造的正演精度, 通常对整个模型使用较小的采样网格, 以保证对频散现象的抑制和对稳定性的保障。这对保持模型计算的灵活性非常不利, 还会对模型的其他部分造成过采样, 增加惊人的计算量。所以, 设计一种更合理的离散化方式进行正演模拟是非常必要的。
从空间采样的角度考虑, 既要提高模拟精度、又要降低计算机内存需求的最有效方法就是在模型的不同区域采用不同步长的网格, 即变网格(也称为不规则网格) 。为此J ast ram 等提出了基于二维声波方程针对某一深度变网格步长的算法[2];J ast ram 等
[3]; [4Pitarka 实现了低阶精度的不[5]。
本文在前人研究的基础上提出一种高阶变网格数值模拟方法, 实现了高精度的局部网格加密技术, 另外还改进了前人对变化网格过渡带采用的波场插值算法, 避免了在过渡带区域进行大量的插值计算, 也不需要平滑函数进行平滑, 同时还保持了算法的稳定性。
2 方法原理
文中所提出的算法是基于Virieux [1]和Graves [6]
提出的一阶速度—应力方程, 对于各向同性介质, 有
τxx +5τ5t v x =b (5x z xz ) +f x
τxz +5τ5t v z =b (5x z zz ) +f z λ+2μ) 5x v x +λ5τ5z v z t xx =(
λ+2μ) 5z v z +λ5τ5x v x t zz =(
(5z v x +5x v z ) 5τt xz =μ
(1)
ττ式中:v x , v z 分别表示速度; τxx , zz , xz 分别表示应力向量; 5x 、5z 、5t 分别表示一阶偏微分算子5/5x 、5/5z 、5/5t ; b 表示密度的倒数; f x , f z 为体力向量;
λ及μ为拉梅常数。
3山东省东营市中国石油大学地球资源与信息学院,257061
本文于2008年1月8日收到, 修改稿于同年7月30日收到。本课题由国家973专题(2007CB209605) 、国家863专题(2006AA06Z206) 和CNPC 物探重点实验室中国石油大学(华东) 研究室资助。
712 石油地球物理勘探2008年
通过对时间导数采用二阶近似, 得到弹性波一阶速度—应力方程的离散形式为
v x v
n+2
i , j
=v x
n-2
i , j
τxx +D τ+Δtb x (D x z xz ) |n
i , j
n
Δ 1=d x i /2, Δ2=d x i /2
n-1
n-1k =1
n+2
z i+2, j+2
=v
n-2
z i+2, j+2
τ+Δtb z (D x xz +D τz zz ) |i+, j+22
n+2i+2, j n+2i+2, j
Δn-1=
k =1
d x i+k+d x i /2, Δn =∑d x i-k +d x i /2∑
ττ
n+1
xx i +2, j n+1
zz i +2, j
=τ
n
xx i +2, j n
zz i +2, j
λ+2μ) D x v x +λ+Δt[(D z v z ]
图1 计算对称点为i +
的差分系数所需的变网格节点2
n-1
=τ) D z v z +λ+Δt[(λ+2μD x v x ] 同理, 以i 为对称点的差分算子为Δn-1=
n-1k =1
n+1n
τμxz (D z v x +D x v z ) ]xz i , j +2=τxz i , j +2+Δt[(
n+2
i , j +2
d x ∑
i+(k-1)
+d x i+n-1/2, Δn =
k =1
d x ∑
i-k
+d x i-n /2。
(2)
式中:D x , D z 分别表示对x , z 的一阶微分算子;
b x =(b i , j +b i +1, j ) /2, b z =(b i , j +b i , j +1) /2; μxz =(1/μi , j +1/μi +1, j +1/μi , j +1+1/μi +1, j +1) /4。
用变量g 来表示速度v x 、v z 或应力向量τxx 、ττzz 、xz , 则在常步长的交错网格技术中, 用下式计算式
(2) 中的g 关于x 的2n 阶空间偏导数[7]
D x g (x , z ) =
C [g (x +Δx
Δx ∑
n i
i =1
下面给出变网格算法的任意偶数阶精度差分近似式及差分系数的求取方法, 即
n
D x g (x , z ) =
i -1
[c ∑
2i-1
g (x +Δ2i-1, z ) +
(4)
+c 2i -Δ2i , z ) ]
n
, ) -
i ; i 是空间差分算子, 它
d g 关于x 。令g (x , z ) =g z exp (i k x ) , 则式(4) 可写成
n
(3) -g (x -Δx i-1z x 2n 式中, 系数C i n 由x (x 2i -1, z ) 和g (x -Δx 2i -1, z ) 。
与常网格不同, 变网格的差分算子是空间变化
i k =
i -1
[c ∑
2i-1
Δ2i-1) +c 2i exp (-i k Δ2i ) ]exp (+i k
(5)
的。给定地质模型, 在确定网格剖分后就可求出各
网格节点的差分算子供递推计算时调用。在交错网
格技术中差分算子有两个对称点, 即i +或i , 不
2
同的对称点对应的Δi 也不同。下面以对称点i +
2
为例求取Δi , 如图1所示。
对式(5) 中的指数项进行泰勒展开, 它的2n 阶泰勒展开式为
2
Δi ) ≈1+i k Δi +(i k ) 2Δexp (i k i +…+
2
2n-12n
(k Δi ) 2n-1+o (Δ(6) +i )
(2n -1) ! 把式(6) 代入式(5) 并整理, 可以得到
2222
ΔΔ) () (c 1ΔΔ2 i k =(c 1+c 2+…+c 2n ) +i k (c 1Δ-c +…-c +i 1222n 2n 1+c 2Δ2+…+c 2n 2n ) +
2
+(i ) 3
3
6
332n-1
(c 1ΔΔ31+c 2Δ2+…-c 2n 2n ) +…+(i )
2
n-12n-12n-1
(c 1Δ-c 2Δ+…-c 2
n Δ2n ) 12
(2n -1) !
(7)
2n-1
写成矩阵的形式
1i Δ1
1
-i Δ2
…
…ω
1i Δ2n-1
1
-i Δ2n
c 1 c 2 …=
2n-22n-1
0i
…
(i ) 2n-2Δ1
2n-2
2n-1
(i ) 2n-1Δ1
…
2n-2
(i ) 2n-2Δ22n-12n-1
Δ-(i ) 2
……
(i ) 2n-2Δ2n
2n-1
Δ-(i ) 2n
2n-2
…(i ) 2n-2Δ2n-12n-1…(i ) 2n-1Δ2n-1
c 2n-1c 2n
00
(8)
第43卷 第6期李振春等:一阶弹性波方程的变网格高阶有限差分数值模拟 713
通过解方程组(8) , 可以得到有限差分算子D x 的系数c i 。同样地, 差分算子D z 的差分系数也是通过这种方法求得。对于交错的网格, 由于变量所定义的网格点的位置不同, 对于每一个差分算子D x 或D z 都有两组不同的差分系数。
在使用变网格时, 网格步长的变化可能会导致数值反射的出现。究其原因是因为波场离散后, 相速度是网格步长的函数。当相速度梯度较大时, 即使速度和密度都没有变化, 入射波的能量也会部分反射回来, 导致了数值反射现象。因为频散程度与频率有关, 所以当入射波的频率较低时, 产生的数值反射较高频时弱。在网格变化的区域对波场进行插值计算虽可以在某种程度上压制数值反射, 但是插值算法计算量大, 且效果并不令人满意, 所以本文没有采用插值计算。本文提出的算法是根据各点所对应的步长的变化计算其差分系数, 在一个循环内完成所有点的波场计算, 不会产生数值反射, 并且减小了计算量。
式中:h =max (d x i , d z i ) 表示最大的空间网格步长;
v min 是地震波速的最小值。
本文所提算法需要满足的稳定性条件为
≤
min (d x i , d z i ) n i =1
∑
(-1) i-1C i
式中:V max 为最大速度值; ρ为密度; C i 为式(3) 中做2n 阶空间展开时的系数。
4 数值模拟实例
为了验证提出的变网格差分算法的稳定性及有效性, 设计两个模型正演实例来说明文中所述方法的实际效果。4. 1 。模×3000m , 纵波速度为3500m/s , 横0125来计算, 密度为常数。P 波震源激发, 震源主频为30Hz , 位于(1500m ,1500m ) 处。变网格步长为10m ,5m 。与其进行对比的常规交错网格步长为5m 。数值模拟中, 差分解法的精度为o (Δt 2, Δx 10
) 。弹性波波场快照和单炮记录分别记录于图2和图3中。
从上面的弹性波波场快照和单炮记录中可以看出, 变网格的过渡区没有产生数值边界反射, 说明本文提出的算法成功地消除了空间网格变化造成的数值反射, 进而说明了该算法的合理性和良好的稳定性。
3 法较好地解决了边界反射问题。使用透射边界能够吸收绝大部分边界反射的能量, 剩余能量通过衰减边界再次减小, 这样计算量增加不大, 吸收也较彻底, 可以较好地消除人为边界反射。
稳定性条件是有限差分数值模拟中十分重要的问题[9]。经过验证, 利用本文所提出的算法能够精确模拟的最大的地震波频率为
f max
图2a 变网格波场(左) 和常规网格波场(右) 水平分量快照(t =360ms )
714 石油地球物理勘探2008年
第43卷 第6期李振春等:一阶弹性波方程的变网格高阶有限差分数值模拟 715
412 背斜模型
模型大小为3000m ×2000m , 模型中横波速度取泊松比为0. 25来计算; P 波震源激发, 震源子波主频为20Hz 。采用精度为o (Δt 2, Δx 10) 的常规交错网格有限差分法和变
网格有限差分法进行模拟对比。模型参数及单炮记录如图4和图5所示。
从上面的模拟效果对比分析得知:在用常规交错网格有限差分法模拟的炮记录中, 由于背斜的影响, 自背斜顶界面反射以下出现了严重的数值频散; 而用变网格算法进行模拟时, 由于在背斜低速区域采用小网格步长计算, 故由背斜影响导致的频散得到了较好地压制, 充分显示了此算法较高的模拟精度
。
图4
背斜模型
716 石油地球物理勘探2008年
5 结论
(1) 本文在网格过渡区域, 没有使用插值计算, 而
是根据各点所对应的步长的变化计算其差分系数, 在一个循环内完成对所有点的波场计算, 成功地消除了数值反射, 且同波场插值计算相比减少了计算量。
(2) 对于常规网格算法难以解决的地下复杂构造的问题, 采用变网格算法可较好地解决, 且能较灵活地适应地层模型的变化, 对低速体的适应性较常规网格算法有了较大改善。另外, 相对常规网格算法而言, 变网格算法节省了计算机内存, 减少了计算时间, 且模拟的区域越大, 节省的时间就会越多, 避免了对整个模型的精细采样造成的资源浪费。
(3) 本文提出的算法精度可达任意偶数阶精度, 高阶变网格算法的提出, 尤其是与交错网格有限差分法相结合, 可以有效地提高模拟精度和计算效率。
参考文献
[1] Virieux J. P 2media :velocity 2stress finite 2difference method. Geo 2
physics , 1986, 51, 88~901[2] J astram C and Behle A. Acoustic modeling on a verti 2
cally verying grid. Geop hys Pros p , 1992, 40, 157~169[3] J astram C and Tessmer E. Elastic modeling on a grid of
vertically varying spacing. Geop hys Pros p , 1994, 42, 357~370
[4] Falk J , Tessmer E and Gajewski D. Tube wave model 2
ing by the finite 2difference method with varying grid spacing. Pros p , 1996, 148, 77~93
[5] Pitarka A. 3D elastic finite 2difference modeling of seis 2
mic motion using staggered grids with nonuniform spacing. B ulletin of Seismological of A merica , 1999, 89, 54~68
[6] Graves R W. Simulating seismic wave propagation
in 3D elastic using staggered 2grid finite 2differences. B ull Seism S oc A m , 1996, 86, 1091~1107[7] 董良国, 马在田, 曹景忠. 一阶弹性波方程交错网格高
阶差分解法. 地球物理学报,2000, 43(3) :411~419[8] 廖振鹏, 周正华, 张艳红. ,2002,45(4) :534~545[ , 地球物理学报,2000, 43(6) :(本文编辑:冯杏芝)
(上接第688页)
虚线框在横向上的分布与图8中虚线框在横向上的分布比较可以得知, 阻抗低值处并非都富集天然气。通过时频分析技术和阻抗反演技术相结合, 最终确定储层的分布范围为图8中白色实线框内。
从对实际资料的应用和分析中发现, 将时频分析技术和波阻抗反演技术相结合可以更加精细地刻画礁滩相储集层的空间分布、油气的富集范围, 计算油气储层的厚度, 估算油气的储量, 从而有效地对礁滩相储层进行预测, 降低钻探风险。
只对垂直剖面进行了探讨, 建议在实际的生产研究中也应该考虑对沿层切片、水平时间切片及波阻抗剖面和单频剖面叠合图进行分析, 从而更加准确地对礁滩相油气藏进行预测。
参考文献
[1] Thomas Englehart ; 石天译. 墨西哥湾中部气云四分量
地震数据的解释. 石油勘探开发情报,2001(6) :33~43[2] 李等华, 唐跃, 殷积峰等. 川东黄龙场构造上二叠统长
兴组生物礁特征与潜伏礁预测. 中国地质, 2006, 33(2) :427~434[3] Castagna J P ,Sun Shengjie ,Siegf ried R W. Instantane 2
ous spectral analysis :Detection of low 2f requency shad 2ows associated with hydrocarbons. T he L eading Ed ge , 2003, 120~127[4] 董臣强, 王军, 张金伟. 时频分析技术在三角洲层序分
析中的应用. 断块油气,2002,9(2) :18~20[5] 张宏彬, 何樵登等. 宽带约束反演. 石油物探, 1995,
34(1) :1~10[6] 张霖斌, 姚振兴等. 快速模拟退火算法及应用. 石油地
球物理勘探,1997,32(5) :654~659[7] 王赟, 王妙月, 彭苏萍. 地球物理随机联合反演. 地球物
理学报,1999,42(Supp. 1) :142~151[8] 寇晓丽, 刘三阳. 基于模拟退火的粒子群算法求解约
束优化问题. 吉林大学学报(工学版) , 2007, 37(1) :136~140
5 结束语
通过对礁滩相油气储集层的特征分析和研究, 提出了将时频分析技术和波阻抗反演技术相结合对礁滩相油气藏进行预测的方法。文中分析了单独应用两种技术对礁滩相油气藏预测的优缺点。综合两种技术的优点, 可以更加有效地对礁滩相油气藏进行预测。通过对实际资料的应用发现, 采用文中提出的礁滩相预测研究方法和思路, 可以有效识别储集层的位置, 精细刻画储集层的范围, 计算储集层的厚度, 减少钻井的风险, 更好地指导油气开发。文中
(本文编辑:冯小球)
Vol. 43 No. 6Abstracts Ⅶ
V ariable 2grid high 2order f inite 2difference numeric simulation of f irst 2order elastic w ave equation. Li Zhen 2chun 1, Zhang H ui 1and Zhang H u a 1. OGP , 2008,43(6) :711~716
U sing variable 2grid finite 2difference algorit hm to carry out t he seismic simulation has many u 2nique advantages t hat mainly are more reasonable for discretization of geologic model , and relatively fine division of local grid can be implemented in low 2velocity zone and complex st ruct ural area , which not only improves t he simulated p recision and eliminates t he dispersion resulted f rom insuffi 2cient sampling , but also can reduce t he require 2ment of comp uter memory , keeping t he flexibility of model comp utatio n. The paper presented a new elastic wave numeric simulatio n met hod based on high 2order staggered 2gird technique , which real 2ized t he technique densifying local grid by changing t he spatial step of grid , remedying t shortcoming of ordinary showed t hat t he t he paper is characterized by , improved p reci 2sion of simulation , reducing co mp utational time and improving comp utational efficiency.
K ey w ords :variable grid , staggered grid , high 2or 2der finite difference , elastic wave
1. College of Eart h Resource and Informatio n , Chi 2na University of Pet roleum , Dongying City , Shan 2do ng Province , 257061, China
Elastic w ave simulation of P 2and S 2w aves separation by 22D staggered grid and application. Zhang Jian 2lei 1, Tian Zhen 2ping 1and W ang Cheng 2xiang 1. OGP , 2008,43(6) :717~722
Through systematically st udying on t he met h 2od separating P 2wave from S 2wave by second 2order elastic wave equation and first 2order staggered grid algorit hm , t he paper presented t he formula rep re 2sentation simulating P 2and S 2waves separation in first 2order staggered grid. The P 2and S 2waves are nat urally separated but do not be coupled toget her in a simulating p rocess of t he met hod , p roducing f ull P 2wave and f ull S 2wave simulated records. The result s of numeric simulation showed t he correct 2ness of t he met hod. There is important p ractical
meaning for knowing p ropagation rule of elastic wave and guiding numeric simulation of reservoir by st udy of f ull P 2wave and f ull S 2wave simulation. K ey w ords :elastic wave simulation , P 2and S 2wave separatio n , staggered grid
1, Technical Research Center of Geop hysical Pros 2pecting , B GP , Zhuozhou City , Hebei Province , 072751, China
Staggered 2grid pseudo 2spectrum simulation of frac 2ture and cave reservoir. Liu Jiong , Ma Jian 2w ei 1, Yang H ui 2zhu 1and B a Jing 1. OGP , 2008, 43(6) :723~727
The paper carried out t he analysis and correla 2tion of t raditional p seudo 2spect rum and improved staggered 2grid p seudo 2spect rum f rom t heoretical met hod cases , t he superiority of 2met hod. The more of differential localization of p seudo 2spect rum is showed by first 2order differential comp utational case of p ulse se 2quence ; t he correct ness of staggered 2grid p seudo 2spect rum met hod is proved by 22D acoustic wave comp utational case ; t he precision of t he staggered 2grid p seudo 2spect rum met hod being higher t han t hat of traditional p seudo 2spect rum is proved by t he co mp utational case of elastic wave in non 2hom 2ogeneo us medium. Finally , t he propagation issue of elastic wave in non 2homogeneous f ract ure and cave reservoir is simulated by staggered 2grid p seu 2do 2spect rum met hod.
K ey w ords :numeric simulation of wave equation , t raditional p seudo 2spect rum , staggered 2grid p seu 2do 2spect rum , fract ure and cave reservoir
1. Instit ute of Seismic Wave Exploration and De 2velop ment , Tsinghua University , Beijing City , 100084, China
Method for identifying fluid property in volcanite reservoir. Zhang Li 2hu a 1, P an B ao 2zhi 1, Shan G ang 2yi 1and H e Li 2. OGP , 2008,43(6) :728~730
The vulcanite reservoir is characterized by complex lit hology and variety of pore st ruct ure , which brought about difficulty in identifying and e 2valuating t he lit hology of volcanite reservoir. Tak 2ing deep volcanite reservoir in Daqing Oilfield as
2008年12
月第43卷 第6期
作者介绍
吕公河 教授级高工,1963年生; 1984年毕业于长春地质
学院物探专业,1996年获得青岛海洋大学海洋地质专业硕士学位,2005年毕业于西安交通大学工程力学专业, 获得博士学位; 现任胜利油田地球物理勘探高级专家, 长期从事滩海以及高分辨率地震勘探技术研究工作。陈浩林 高级工程师,1969年生;1991年毕业于中国地质
大学(武汉) 地球物理专业, 获学士学位;2001年毕业于中国石油大学(北京) 地球物理专业, 获工程硕士学位; 2005年毕业于中国矿业大学(北京) 地球物理专业, 获博士学位。长期从事地震资料采集、处理方法及理论研究。王振华 硕士,1958年生;1982年毕业于华东石油学院勘探
地球物理专业,1985年毕业于华东石油学院北京研究生部应用地球物理专业。主要从事地球物理勘探技术的开发与应用。现在北京华阳捷泰技术有限公司工作。单联瑜 高级工程师, 1959年生;1982年毕业于华东石油
学院地球物理勘探专业, 究和地震资料处理研究所在读博士王玉贵 高级工程师; 年毕业于石油大学
(华东) ,2006年获中国石油大学应用地球物理专业硕士学位, 现为中国地质大学地球物理与信息技术学院在读博士。一直从事VSP 资料处理、解释及方法研究工作。裴正林 博士, 副研究员,1962年生;1984年毕业于西安地
质学院物探专业,2000年毕业于中国地质大学(北京) , 获博士学位,2000~2003年在中国石油大学(北京) 从事博士后研究。主要从事地震波传播理论及数值模拟, 井间地震层析成像, 小波域地震属性提取和油气识别, 以及小波理论应用研究。获省部级科技进步一等奖1项, 二等奖1项, 三等奖2项。合著1部, 发表论文50余篇。段云卿 副教授,1958年生;1982年毕业于华东石油学院物
探专业, 现在中国地质大学(北京) 从事教学和地震资料处理及解释方法研究。已发表论文多篇。
何兵寿 副教授,1973年生;1996年毕业于中国矿业大学
(徐州) 物探专业,2002年获石油大学(北京) 应用地球物理专业博士学位, 现在中国海洋大学海洋地球科学学院从事多波多分量地震技术、纵波高分辨率处理技术研究和相关教学工作。樊德华 高级工程师,1967年生;1991年毕业于成都地质学
院, 现在中石化胜利油田从事石油勘探工作。刘振武 教授, 英国Heriot 2Watt 大学石油工程专业博士,
1949年生; 中国石油天然气集团公司副总工程师, 享受国务院政府专家特殊津贴。编著有《二十一世纪初中国
油气关键技术》《、油气应用基础研究展望》和《企业技术
创新管理》《、中国石油“十五”科技进展丛书》等著作。陈 林 博士研究生,1980年生;2004年毕业于江汉石油学
院, 获学士学位; 现为中国科学院地质与地球物理研究所博士研究生; 主要从事海洋地球物理和盆地模拟等方面的研究。赵丽平 工程师,1962年生;1988年毕业于胜利油田职工大
学石油地质专业,2004年毕业于石油大学资源勘察专业, 长期从事地震资料综合研究。蔡涵鹏 1983生; 现为成都理工大学信息工程学院硕士研
究生, 就读专业为地球探测与信息技术; 研究方向为地震资料处理与解释。彭文绪 高级地质师,1972年生;1995年毕业于江汉石油学
院, 获学士学位。1998年毕业于石油大学(北京) , 获硕。。 年生;1988年毕业于石油
(, 获学士学位,2005年在中国科学院广州地球化学所获构造地质博士学位。现在主要从事油气勘探综合研究工作。张海燕 副教授,1968年生;1990年毕业于清华大学精密
仪器专业, 获工学学士学位;1993年毕业于中国海洋大学海洋物理专业, 获理学硕士学位;2008年1月毕业于中国海洋大学海洋地质专业, 获理学博士学位。现在中国海洋大学信息科学与工程学院从事AVO 叠前反演、计算智能研究。已经在国内外核心学术期刊上发表论文十余篇, 其中有多篇被EI 收录。李振春 理学博士, 教授, 博士生导师; 中国石油大学(华东)
地球物理系主任。主要从事地震波传播与正演模拟、地震成像与偏移速度分析、多尺度地震资料联合反演与CFP 2AV P 分析理论与方法的教学与研究工作。张建磊 硕士,1978年生;2001年毕业于石油大学(华东) 地
球探测与信息技术学院,2004年获得石油大学(华东) 地球探测与信息技术专业硕士学位, 现在东方地球物理公司物探技术研究中心从事地震正演和偏移方法研究。刘 炯 1980年生, 清华大学在读博士, 主要从事复杂介质
中地震波传播机理的研究。张丽华 工程师,1974年生;1996年毕业于长春地质学院,
现为吉林大学地球探测科学与技术专业博士研究生, 研究方向为石油地球物理测井数据处理与解释研究。柳永坡 博士生, 副教授,1971年生;1995年毕业于大庆石
油学院, 主要研究领域为:软件工程、测试技术、性能故障分析。王 月 1985年生; 现在中国石油大学(北京) 攻读硕士学
位, 主要从事石油物探市场分析及国际石油经济等方面的研究, 已发表论文4篇。