模拟退火算法在 AVO 反演中的应用

模拟退火算法在 AVO 反演中的应用中国论文网5V7_q%O1m/o*i
g_(P
王 瑜中国论文网@ID/e%@x,nfg;o4b
（哈尔滨石油学院）
C+wt#Qa%^6v`0[摘 要]震反射波振幅随炮检距变化（简称AVO技术）是最近几十年发展的新的一种寻找油漆的地震勘探技术。其基本原理是反射系数随入射角度的变化中国论文网2FR!YPx4D
而变化，即反射系数是入射角度的函数，理论的基础是Zoppritz方程及其近似式。模拟退火反演方法是寻求全局最优解的一种较好的方法之一，该方法综合利用了中国论文网 Aj9iX"LS#eK
石油地质，地球物理等进行约束反演，从效果上看该方法是可行的。
(W^kc4Z8nE+GV0[关键词]AVO正演模拟；AVO反演；非线性反演；模拟退火算法
l+F:bqad P0本文首先阐述了进行AVO反演的模拟退火算法的基本原理。然后介绍了中国论文网1\X(E5U)? R{ u9o
AVO反演的地球物理基础，包括Zoppritz方程及其简化形式，并对各种简化形
?tj8l-\0式的特点加以介绍。最后利用模拟退火算法，结合AVO正演模拟，进行了相关中国论文网0o
r`-@"jf{
的AVO参数反演。从效果上看，利用该算法进行反演的速度快效果好，说明本
iTKgW@Jj8Y0算法适用于该反演问题。
SQ$o[^_*{0第一章 模拟退火(SA)算法
;J)}q6el#mE%l0由于传统的单一算法无法适应解决大规模非线性规划问题的需要。而从用中国论文网r9I~`"[+Y@#o)}`
于统计力学的Monte Carlo方法上得到启发，模拟退火算法(Simulated an-中国论文网8mnTW U;X1vX)d4D+e
nealing algorithm，简称为SA)由Kirkpatrick等人于1983年提出，并成功的将
,{1wR#z fG0其应用于组合优化领域，从而为多维问题和非线性问题的解决提出了一种新的中国论文网8GM&b(ba"_
思路，特别是为NP完全组合优化问题提出了一种有效地近似算法。模拟退火算
Vd;JOO'A7LS0法源于统计热力物理学,模拟的是一个在熔融状态下物体逐渐冷却达到结晶状
6}1V;wHs M-KF}0态的物理过程；采用Metropolis准则进行判断来决定是否接受；并采用一组冷中国论文网9om/k$R@3f.t,B,C
却进度表的参数控制算法的过程，从而使模拟的热系统最终冷却结晶达到系统
'z5Y6{JE3G;CL0能量最小值，获得算法的近似最优解。
_+S}k7b4i%R01.1 模拟退火算法描述中国论文网$N'|)xU)R1l#r6g
1.1.1 经典的Metropolis退火算法(MSA算法)
QA/@ i[ Q7N{p0(1).随机产生一个初始解中国论文网k J0LpF)T
,令最优解中国论文网[@;j$U2m3bU'l~\
,并计算目标函数值

S(M P|T;q]0;中国论文网4f ]1al;j6]Q
(2).设置初始温度T(0)=To,迭代次数i=1；中国论文网Rqj'P4s}J8J
(3).do while T(i)>Tmin
Yl,w$XJ1z8t/wIX,n0(a).for j=1:k
y6lxb5F3zo3^m0(b).对当前最优解按某一邻域函数，产生一个新解中国论文网 oEo\-V+S3b.w*Z*D+E
.计算新的目标函数
%m0D c3Z ~S0值E(x2),并计算目标函数的增量。
2{"e?Ugy?0(c).若目标函数增量
{-`Bf#FL PE0<0,则
@-??_L-zy0。中国论文网W`0\c)u Q
UTe
(d).若目标函数增量中国论文网8f iq'D.a#L%h:h
>0,中国论文网&l*M.p$q:cu$zh_XR
;
MPpF$E%C*TZ0h;v0(通
W}4ZP"U7F0过概率来接受能够有效的跳出局部极小值的束缚）。中国论文网Q-Q` ?.z\R
(e).中国论文网$iY xiP-E
；否则有
"?yg._ yYM0。
^BY3m+Z"_HHLb0End for
o:_L S0V T5?L0(4).i=i+1;
-kE,D M
g%Q'S\4n0(5).end do中国论文网/\$J nWd7OP,an
(6).输出当前最优解，计算结束中国论文网8c:G{\f1}l
1.2 冷却进度表
'kL{X6b-pm0冷冷却进度表(cooling schedule)是模拟退火算法的核心,,它包括初始温
nfr{;`9?"j0度T及衰减函数、对应的Markov链长度以及停止准则等4个方面, 是模拟退火中国论文网g5\z4G,U5_we
算法应用于实际问题的关健所在。下面注意讨论冷却进度表的选取方法：
.c1Br\)zL01.初始温度中国论文网@7YzR8v}
的选取中国论文网 n#W&els3t"c
初初始温度需要选取的足够高,以防止模拟退火算法落入局部极优值的陷
^u3RZ5N[d+A0阱中，如果温度选得太小,则模拟退火算法一旦落入局部极优值陷阱中就很难中国论文网zf8WG$cb/EnL_'S
再跳出来,从而无法逃离。中国论文网] X_c*HAA
2．衰减函数的选取
]%Qc!Yh(?Y0衰减函数用于控制温度的退火速度。温度下降的太快会导致解的质量下中国论文网^2z FKbJSi
降，类似于固体降温太快会出现瑕疵，降温太慢又会导致太长的程序运行时间，
/H9L/EK6N2Cc4w0尤其在问题的规模相对较大时表现得尤为突出。常用的退火机制是：中国论文网%k2q3K#MJ
(1).指数下降型
+Gy*F]wzCo0（3-9）中国论文网 s#c.rq,m,C]
式中:k为迭代次数;c为衰减因子;N为状态的个数;T0为初始温度。中国论文网$Y%Utaa"]
(2).双曲线下降型
3D~/y6AYU
k!{s9b0(3-10)
aCGK5xdy:V#g/k0式中:k为迭代次数;中国论文网"e
BTn]W
为初始温度。初始温度中国论文网b[/Bua,L%Ib#I8z[
不能取得太高,否则增加
Y3vf:o2Fv}v0计算时间浪费机时;中国论文网%[0E%dOK V%z
也不能太低,否则模型选取不能遍及整个模型空间,只是中国论文网 j0@
c X h
dt%{q(v
在初始模型附近选取,不能进行全局寻优。所以中国论文网z(L Od6p'd4F \
的确定只有通过实验计算得中国论文网y5wihTTtX&zQ0@
到。
3I+WG7K}:?-TWZ4aX0(3).常系数型中国论文网
~J[&uu'u,N
（3-11）中国论文网 vK%nv0g!ca
式中a为衰减系数，是一个非常接近于1的常数，该方法形式简单，而且对只中国论文网4Y"\4C[+jC'FC
要求寻找近似最优解的问题也足够有效，因而成为目前常用的一种方法。
hz%z"v$D)uD Mn:y03.Markov链长度L的选取
ST-Q)u{Vu0原则是：在衰减参数T的衰减函数已经选定的前提下，L应该选在控制参数中国论文网8X2Me0x6\s&Xm+G
的每一取值上都能恢复准平衡。
5R2y!D#y|c4a04.终止条件
oGz]yS U7}$W0有很多不同的终止条件，而且不同的终止条件对解的质量和算法性能有很中国论文网1Hb!q#B#N.VP
大的影响。这里只介绍一种常用的方法：上一个解与最新的一个最优解之间的中国论文网7W!?9mPy~
p
差值小于某个容差，即可停止此次Markov链的迭代。中国论文网9b L8Z B|2P
一个有效的冷却进度表的判断依据：中国论文网!SfYc|5_v+^
(1).算法的收敛：主要取决于衰减函数和马尔科夫链的长度及停止准则。中国论文网D4by9F2xx/|5^
(2). 算法的实验性能：最终解的质量和CPU的时间。
q2t3S5OCZ0第二章 AVO三参数反演
[
U c)q%E9j+r8o02.1 AVO反演的基础
XQ9X+R F li0AV O技术的理论基础是描述平面波在水平分界面上反射和透射的

Cj"l@|0Zoeppritz方程。尽管该方程早在20世纪初就已经建立，但由于其数学上的复杂
x!`V9]N)YM;PFK0性和物理上的非直观性，一直没有得到直接的应用。尽管该方程早在20世纪初
|{;L:m
o0就已经建立，但由于其数学上的复杂性和物理上的非直观性，一直没有得到直中国论文网k ]Y/Qw5prh*I"A
接的应用。为了克服由Zoeppritz方程导出的反射系数形式复杂及不易进行数
)Hk:E}{w GiV0值计算的困难，许多学者对Zoeppritz方程进行了简化，并从不同的各个方面帮中国论文网w*ZM-B1G/F
助我们理解岩性参数，纵横波速度，密度和泊松比对于反射系数的影响。中国论文网Q*|3Q}6uLm
2.1.2 简化的Zoeppritz方程中国论文网 qN5lq:U*mV
AVO技术的核心思想是利用在不同介质中，反射系数随入射角的变化规
!g$y&r+R |^M0律来寻找油气层。因此，必须建立一个具有普遍意义的方程，将反射系数表示成
t/?lF q/P8y]*\(]8P%l0入射角和地层参数的函数。完全形式的Zoeppritz方程全面考虑了平面纵波和中国论文网6L-|*fK"|0Q&C
横波入射在平界面两侧产生的纵横波和透射能量之间的关系。当反射波地震勘

q,tNH2_%m#O H0探使用主要产生纵波的震源，接收的是反射纵波时，完全的Zoeppritz方程可以中国论文网b-e%e#cjR${*?IX
被大大的简化，即只考虑平面纵波入射产生的反射振幅随入射角的变化情况。
.d2|RWm0一方面可以减少计算量，另一方面更有利于AVO技术的研究和应用。在这方面，
6q1T-M3](k {-HN4b0前人做了大量的工作，从不同角度对Zoeppritz方程进行了简化，提出了不同纵
4S!A s$Z@0波反射振幅的近似表达式。中国论文网%G;sP:U%Z:])c@
第三章 结论与认识中国论文网 bT qD4_i
利用AVO反演来提取储层参数，是目前进行叠前反演的一种有效的手段中国论文网_2@J(uY9t%Q
之一。而通过AVO反演获取的储层参数，如泊松比等，可以直接用来预测油气。中国论文网FM)Vo\
本论文通过中国论文网-Q,iRf-e(`p_V
可以得到以下几点结论：中国论文网mO;L1I:o$h|DJ4Bn
1.本文利用了在AVO参数反演中常用的Zoeppritz简化方程Aki和Richard中国论文网'\;Q^:x4zS#I*E0h
方程来进行反演。通过模型数据的计算可以发现简化方程可以准确的反演各层中国论文网6K.ky,Cpq,K2r:K
的反射系数并可以达到很好的反演效果。中国论文网2wK1y N4QAv+` xn\
2.常规的最优化方法依赖于初始模型的选取，易陷于局部最优，无法满足中国论文网3I`xeggb7~N
地震反演问题的需要。而模拟退火算法避免常规线性反演方法的缺点，不依赖
5wH0N:GS0oz-C0于初始模型的选取，同时可以很好的跳出局部最优解，适合应用于解决相关的中国论文网/@'vg"P{ s$d8wJ'sB
地震反演问题。中国论文网/Tm {E+uje Gt
3.通过前人对于三类含油气岩石的AVO特性的总结，可以较方便从实际
r/`W4pFD#k[0的资料上选取AVO特性明显的目的层进行反演，可以减少计算量,而且还可以中国论文网-S)~i2b8up!D
利用含油气的AVO特性对目的层段的反演结果进行检验。中国论文网L\&iFK6f-Nv
4.在AVO参数的反演中，影响反演结果的因素很多。如Zoppritz方程的简中国论文网*ud6W.WZ[+}$}%k3["u)[
化方程，反演算法，AVO正演模拟等各个环节都会影响到结果的好坏。因此，可中国论文网vSbMb1l"h9\~
以通过对于这些环节的改进使反演结果更为精确。