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

模拟退火算法在 AVO 反演中的应用
,jd5[V3Ap;H0王 瑜
|[3ww7R0（哈尔滨石油学院）中国论文网qe7?cIg
[摘 要]震反射波振幅随炮检距变化（简称AVO技术）是最近几十年发展的新的一种寻找油漆的地震勘探技术。其基本原理是反射系数随入射角度的变化中国论文网v!T+Ir+iq
而变化，即反射系数是入射角度的函数，理论的基础是Zoppritz方程及其近似式。模拟退火反演方法是寻求全局最优解的一种较好的方法之一，该方法综合利用了中国论文网9G.`x P_y5T:bng
石油地质，地球物理等进行约束反演，从效果上看该方法是可行的。
&] {8hj cU0I#[ zE0[关键词]AVO正演模拟；AVO反演；非线性反演；模拟退火算法中国论文网$SX4yz|Kb$@
o-lC
本文首先阐述了进行AVO反演的模拟退火算法的基本原理。然后介绍了中国论文网4@.m!@1aq
AVO反演的地球物理基础，包括Zoppritz方程及其简化形式，并对各种简化形
Ih-I T1y/D0式的特点加以介绍。最后利用模拟退火算法，结合AVO正演模拟，进行了相关中国论文网O]@O Dmv
的AVO参数反演。从效果上看，利用该算法进行反演的速度快效果好，说明本中国论文网$mkBwQ| w#xg
算法适用于该反演问题。
5^!\$mK In
H1G0第一章 模拟退火(SA)算法

Xxf"Fr{$NU+}0由于传统的单一算法无法适应解决大规模非线性规划问题的需要。而从用
8liH;U}R0于统计力学的Monte Carlo方法上得到启发，模拟退火算法(Simulated an-中国论文网,pG?\+fN4U$Y-@
nealing algorithm，简称为SA)由Kirkpatrick等人于1983年提出，并成功的将
.vfcGoU'~&f0其应用于组合优化领域，从而为多维问题和非线性问题的解决提出了一种新的
$b2nV];~M vG"?0S0思路，特别是为NP完全组合优化问题提出了一种有效地近似算法。模拟退火算中国论文网_7Zf(yPA8G|t
法源于统计热力物理学,模拟的是一个在熔融状态下物体逐渐冷却达到结晶状
L&n3x[6R1d~"G0态的物理过程；采用Metropolis准则进行判断来决定是否接受；并采用一组冷
O#eFY)_;c)d0却进度表的参数控制算法的过程，从而使模拟的热系统最终冷却结晶达到系统
4E8eA)b?
^4OR$U @0能量最小值，获得算法的近似最优解。中国论文网`w1L oS k!f-w+~W:~
1.1 模拟退火算法描述中国论文网7l-z8n(X:g\
1.1.1 经典的Metropolis退火算法(MSA算法)
vX(y%P(D-h:~0(1).随机产生一个初始解
0kRnT2X Z1y0f0,令最优解中国论文网"Q PK8N])Ra
,并计算目标函数值
'd
{$X{!u/}FD;L0;中国论文网"?N,f7lk&gh&r
(2).设置初始温度T(0)=To,迭代次数i=1；中国论文网0` d*F J8rQb9R7qQ
(3).do while T(i)>Tmin
'{.W9w.WkAz'@0(a).for j=1:k
,F!Tl+tI|8?]:y0(b).对当前最优解按某一邻域函数，产生一个新解
_ r2D2S,k;B7T(dA0.计算新的目标函数中国论文网&eubl4X8}
值E(x2),并计算目标函数的增量。

Z3|)y
myv m2M'Qg0(c).若目标函数增量中国论文网;C,Ir5wR
<0,则中国论文网]-Z\_b1fB
。中国论文网1ID
]1w})XK
(d).若目标函数增量
N1s
aA.r_*g0>0,中国论文网4ucqH(s"V` X
;
Bt$G8Hx i-j0(通中国论文网Y"fW*O9rW3_ah
过概率来接受能够有效的跳出局部极小值的束缚）。中国论文网I`A+K9G%I
(e).中国论文网iD+~)`!Eq-Fpd
；否则有
?DX0kZ0。
5r-V_ Jy;iC(N0End for
L Cw:u8h-J$QH-[#G0(4).i=i+1;中国论文网q|7nd'm+{)I8v8v
(5).end do中国论文网^+IO1twe"W
(6).输出当前最优解，计算结束
Mh^"g2DP01.2 冷却进度表
XD5UC"OXX$w0冷冷却进度表(cooling schedule)是模拟退火算法的核心,,它包括初始温
m.iW2}H9I G0度T及衰减函数、对应的Markov链长度以及停止准则等4个方面, 是模拟退火
gek/Zr
|;rF0算法应用于实际问题的关健所在。下面注意讨论冷却进度表的选取方法：中国论文网z ^n6GGK
1.初始温度
"u0Cp#SXR]p^0的选取中国论文网2B|` {K0Jm1_A
初初始温度需要选取的足够高,以防止模拟退火算法落入局部极优值的陷
.G.Nzm`6}([4~(V0阱中，如果温度选得太小,则模拟退火算法一旦落入局部极优值陷阱中就很难中国论文网 bJC"@7?f]6k]
再跳出来,从而无法逃离。中国论文网3R,[Pn \^s
2．衰减函数的选取
e&}vy\t(C Na0衰减函数用于控制温度的退火速度。温度下降的太快会导致解的质量下
K:I'grPt0降，类似于固体降温太快会出现瑕疵，降温太慢又会导致太长的程序运行时间，
$y*\'aEIX;p]!L0尤其在问题的规模相对较大时表现得尤为突出。常用的退火机制是：中国论文网ZU1F
v(vW8Z
(1).指数下降型
"B9o Q u.PA*?!}xS0（3-9）
C&K D.@NDN0式中:k为迭代次数;c为衰减因子;N为状态的个数;T0为初始温度。中国论文网K]lD7g
(2).双曲线下降型中国论文网o0I x Z}#`
(3-10)
5@1Ee.hkX4D0式中:k为迭代次数;
Kxg,w'F*K!m b0为初始温度。初始温度
${qD(_vT/PN.n3J0不能取得太高,否则增加
jWD$d*B2z0计算时间浪费机时;
{8lB\&vZ^0也不能太低,否则模型选取不能遍及整个模型空间,只是
4h-RjFX Fh0在初始模型附近选取,不能进行全局寻优。所以中国论文网7oh~$X-Q8y2`
的确定只有通过实验计算得中国论文网#x.}:R)? Ukef
到。中国论文网2YYS@ g6G3u/U
(3).常系数型
#crM!v%~f][.S0（3-11）
-h3W`GH9tR0式中a为衰减系数，是一个非常接近于1的常数，该方法形式简单，而且对只

k{T2`qakf#x t0要求寻找近似最优解的问题也足够有效，因而成为目前常用的一种方法。中国论文网.I&S:fFd^6^"s"n
3.Markov链长度L的选取中国论文网V k_Tm`wEpz
原则是：在衰减参数T的衰减函数已经选定的前提下，L应该选在控制参数中国论文网Q4ew'c8Z9B9`5fD
的每一取值上都能恢复准平衡。中国论文网"sar `gd]g
4.终止条件
-ux9D,Z T|0有很多不同的终止条件，而且不同的终止条件对解的质量和算法性能有很中国论文网V&BK7GN!k%T{{
大的影响。这里只介绍一种常用的方法：上一个解与最新的一个最优解之间的
aeg6F-rq|h0差值小于某个容差，即可停止此次Markov链的迭代。中国论文网5[u:C^$c ff@5?
一个有效的冷却进度表的判断依据：中国论文网/l z^!E$OD6T Vc/?
(1).算法的收敛：主要取决于衰减函数和马尔科夫链的长度及停止准则。中国论文网6YRl1U2L4H(j
(2). 算法的实验性能：最终解的质量和CPU的时间。
3C[!w"DVK0第二章 AVO三参数反演
q d1F-w"ay X02.1 AVO反演的基础
x+gCU(J5}0AV O技术的理论基础是描述平面波在水平分界面上反射和透射的中国论文网:| ?,u-F/xV
Zoeppritz方程。尽管该方程早在20世纪初就已经建立，但由于其数学上的复杂
"Ja%e:P^}0性和物理上的非直观性，一直没有得到直接的应用。尽管该方程早在20世纪初
xA6~Q*d/nL0就已经建立，但由于其数学上的复杂性和物理上的非直观性，一直没有得到直中国论文网 ^]B(h}'?|
接的应用。为了克服由Zoeppritz方程导出的反射系数形式复杂及不易进行数中国论文网O8z-K
t H
值计算的困难，许多学者对Zoeppritz方程进行了简化，并从不同的各个方面帮中国论文网(rnz3Hz
s
助我们理解岩性参数，纵横波速度，密度和泊松比对于反射系数的影响。中国论文网D)fgS9u0m5u$f[
2.1.2 简化的Zoeppritz方程中国论文网gF3W'l$nT
AVO技术的核心思想是利用在不同介质中，反射系数随入射角的变化规
-TX#`&eqaY0律来寻找油气层。因此，必须建立一个具有普遍意义的方程，将反射系数表示成中国论文网t.?%pX^ T?b
入射角和地层参数的函数。完全形式的Zoeppritz方程全面考虑了平面纵波和中国论文网q7G-gn"O"v}I
横波入射在平界面两侧产生的纵横波和透射能量之间的关系。当反射波地震勘
}l%g-ICZG#Wpi0探使用主要产生纵波的震源，接收的是反射纵波时，完全的Zoeppritz方程可以
Z~ aAMrH-A{0被大大的简化，即只考虑平面纵波入射产生的反射振幅随入射角的变化情况。
)l5EQLJuyO7kn0一方面可以减少计算量，另一方面更有利于AVO技术的研究和应用。在这方面，
"C9mz$o
}fLH0前人做了大量的工作，从不同角度对Zoeppritz方程进行了简化，提出了不同纵中国论文网/`Au/LL4V,F;fD;|
波反射振幅的近似表达式。
.d)S"L&M%r0第三章 结论与认识中国论文网 UzB}
r#w @b
利用AVO反演来提取储层参数，是目前进行叠前反演的一种有效的手段中国论文网+d/y4d0EMFbi4R
之一。而通过AVO反演获取的储层参数，如泊松比等，可以直接用来预测油气。
l#R N&uW5jQjb0本论文通过
#Vl,D Pt vk#@_0可以得到以下几点结论：
9?OuSb3y
|01.本文利用了在AVO参数反演中常用的Zoeppritz简化方程Aki和Richard
C+B?hy9w,l3t x0方程来进行反演。通过模型数据的计算可以发现简化方程可以准确的反演各层
L2?g0icT `0的反射系数并可以达到很好的反演效果。
b$_q&W,m)zU{02.常规的最优化方法依赖于初始模型的选取，易陷于局部最优，无法满足中国论文网w&?1V%k;E4M
地震反演问题的需要。而模拟退火算法避免常规线性反演方法的缺点，不依赖
7^F `#i4N0于初始模型的选取，同时可以很好的跳出局部最优解，适合应用于解决相关的
mE3m&`*d w(L*K.~%\zhJ0地震反演问题。
+S&c0F La'^3v}03.通过前人对于三类含油气岩石的AVO特性的总结，可以较方便从实际
9dq`4Z*U.U0的资料上选取AVO特性明显的目的层进行反演，可以减少计算量,而且还可以中国论文网S7Q$z#Zi;]\~
利用含油气的AVO特性对目的层段的反演结果进行检验。
ge"r9hm u"u.Q9I04.在AVO参数的反演中，影响反演结果的因素很多。如Zoppritz方程的简中国论文网XG~
G.sU2aI8}
化方程，反演算法，AVO正演模拟等各个环节都会影响到结果的好坏。因此，可中国论文网&Lq0k.w8CGU]
以通过对于这些环节的改进使反演结果更为精确。