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

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

db@f J0王 瑜
)q
QV$r*^b:v!VgF s0（哈尔滨石油学院）中国论文网1C1vN"zw9y+Z3_.w
[摘 要]震反射波振幅随炮检距变化（简称AVO技术）是最近几十年发展的新的一种寻找油漆的地震勘探技术。其基本原理是反射系数随入射角度的变化
*kq
A"Bm3X&{e)]H0而变化，即反射系数是入射角度的函数，理论的基础是Zoppritz方程及其近似式。模拟退火反演方法是寻求全局最优解的一种较好的方法之一，该方法综合利用了
?Gx5M!BE9s4D~0石油地质，地球物理等进行约束反演，从效果上看该方法是可行的。中国论文网u2@OeDq]5g"U,kZ
[关键词]AVO正演模拟；AVO反演；非线性反演；模拟退火算法中国论文网@p L'QC h)E
本文首先阐述了进行AVO反演的模拟退火算法的基本原理。然后介绍了
"\&i{8N%it3M$}0AVO反演的地球物理基础，包括Zoppritz方程及其简化形式，并对各种简化形
}O"y0\)QZ"V$A%?0式的特点加以介绍。最后利用模拟退火算法，结合AVO正演模拟，进行了相关中国论文网9rB7l|*@
的AVO参数反演。从效果上看，利用该算法进行反演的速度快效果好，说明本中国论文网 }%R"??7z*]2H)D*j
算法适用于该反演问题。
'U&}
^.?
N2`0第一章 模拟退火(SA)算法中国论文网_"EP9v^"|
由于传统的单一算法无法适应解决大规模非线性规划问题的需要。而从用

]&M eE8B.x0于统计力学的Monte Carlo方法上得到启发，模拟退火算法(Simulated an-

qM&S\`p1r0nealing algorithm，简称为SA)由Kirkpatrick等人于1983年提出，并成功的将
tfZ+JA$y!y*}%z/{v#z*C0其应用于组合优化领域，从而为多维问题和非线性问题的解决提出了一种新的

omMZ-fRI0思路，特别是为NP完全组合优化问题提出了一种有效地近似算法。模拟退火算

u Hqy*o0法源于统计热力物理学,模拟的是一个在熔融状态下物体逐渐冷却达到结晶状
*{.\$Cz#W0态的物理过程；采用Metropolis准则进行判断来决定是否接受；并采用一组冷中国论文网%s6jb{]G
却进度表的参数控制算法的过程，从而使模拟的热系统最终冷却结晶达到系统
Rm%z._ bJ!k0能量最小值，获得算法的近似最优解。中国论文网RiEmyWb9D!h-]9J
1.1 模拟退火算法描述
K!lI*g4f m|01.1.1 经典的Metropolis退火算法(MSA算法)
"q-d0u$v;G0(1).随机产生一个初始解
6R/V&Zq
v\
C i"t;fA{0,令最优解
0u7q(T0U1NhBj0A0,并计算目标函数值中国论文网XU,EZ
Eo
;中国论文网
EW:sH#Lk a Au
(2).设置初始温度T(0)=To,迭代次数i=1；中国论文网 ^4A2`A`
H3dh
(3).do while T(i)>Tmin中国论文网;O5~
The_
(a).for j=1:k中国论文网]_F_ N'V*C+g
(b).对当前最优解按某一邻域函数，产生一个新解中国论文网#je)R$od.u+OP%h
.计算新的目标函数
,m7m5F/gP*vJ0值E(x2),并计算目标函数的增量。中国论文网2g c%YRQ7\PW(r9OH?
(c).若目标函数增量
5p&b^%m3}i0<0,则中国论文网"Q4f/ATFsq,s
。
~"w+Q?*a6d:Jus0(d).若目标函数增量
bUYZJ0>0,
%{)N4n? ^GI3A)w0;
)?a8`{)A5_L|}0(通中国论文网;? VnB|#]
过概率来接受能够有效的跳出局部极小值的束缚）。中国论文网 oOQ+i:v?&[|\+Y,B r
(e).
#@;ZKZ"\Nz0；否则有
-_A7Ob:R L0。
|6Y%Bm-R,M0i0End for
Pb`|t%qH2v+Yd[#]3V0(4).i=i+1;中国论文网reM{ ~ F|
(5).end do中国论文网]cD6rK
(6).输出当前最优解，计算结束中国论文网J!A~(Nu
1.2 冷却进度表
yR3V7L(d3\ Da$b0冷冷却进度表(cooling schedule)是模拟退火算法的核心,,它包括初始温
)~0L5}d5^&Vf0度T及衰减函数、对应的Markov链长度以及停止准则等4个方面, 是模拟退火中国论文网8p#NN(gcQ
算法应用于实际问题的关健所在。下面注意讨论冷却进度表的选取方法：中国论文网 z`gj2^ J#oe
1.初始温度
3AA?];EIP6tM0的选取中国论文网 qQ4\h)k
初初始温度需要选取的足够高,以防止模拟退火算法落入局部极优值的陷
t2Lt'nO:wvM!R0阱中，如果温度选得太小,则模拟退火算法一旦落入局部极优值陷阱中就很难中国论文网
P0Q[,MH$Y%B6lXC
再跳出来,从而无法逃离。中国论文网'_#^@1B&p$j+_
2．衰减函数的选取
.X8Og:ll7J$Df`0衰减函数用于控制温度的退火速度。温度下降的太快会导致解的质量下
9x(_;[5P }w"n"I$x0降，类似于固体降温太快会出现瑕疵，降温太慢又会导致太长的程序运行时间，中国论文网W*N:ni,F3l*hMM9\'?
尤其在问题的规模相对较大时表现得尤为突出。常用的退火机制是：
5E GY"eB0(1).指数下降型
r4~"L*UF#z0（3-9）
Dg:d%W/^w)l2b+D*J,KT0式中:k为迭代次数;c为衰减因子;N为状态的个数;T0为初始温度。
T SRu2d
]}0(2).双曲线下降型中国论文网 Z.MPs+d_L
(3-10)中国论文网1m8v S9k1z#G,j
式中:k为迭代次数;中国论文网"Sh}&d"B
为初始温度。初始温度
oH1N.F f0不能取得太高,否则增加中国论文网^4NomB%N
计算时间浪费机时;中国论文网;ah ~:`}OsX
也不能太低,否则模型选取不能遍及整个模型空间,只是
e8xx-C'q V0在初始模型附近选取,不能进行全局寻优。所以中国论文网.J2`+G9w YT
的确定只有通过实验计算得中国论文网6hr;M,Cl`"Ux
到。
)LU&AJ#{eQn0(3).常系数型
!H`&aW[z.x2U'Xi ?0（3-11）中国论文网8l9T0Y#w-CC F$W*B
式中a为衰减系数，是一个非常接近于1的常数，该方法形式简单，而且对只
b_7n-H Hl#S@v%^'H0要求寻找近似最优解的问题也足够有效，因而成为目前常用的一种方法。中国论文网IH/Urq_8O
3.Markov链长度L的选取
;\M![,s%t%C-\0原则是：在衰减参数T的衰减函数已经选定的前提下，L应该选在控制参数
i'v(A? T.^O0的每一取值上都能恢复准平衡。
A|.?8i
X-@} B04.终止条件中国论文网:yBY%xmCn){r
有很多不同的终止条件，而且不同的终止条件对解的质量和算法性能有很
1B~\l.P"KY
B@0大的影响。这里只介绍一种常用的方法：上一个解与最新的一个最优解之间的中国论文网X-`3T{:C.lM
差值小于某个容差，即可停止此次Markov链的迭代。
#^1i.xc\0M#O0一个有效的冷却进度表的判断依据：
-f(pu.I!j8w`KoC0(1).算法的收敛：主要取决于衰减函数和马尔科夫链的长度及停止准则。中国论文网2z
Y@3Iu%]
(2). 算法的实验性能：最终解的质量和CPU的时间。
5Nf gzUT"_0第二章 AVO三参数反演
P sok{3?!\0W&`02.1 AVO反演的基础中国论文网$Qh5Qi_O4v [7p%d3G
AV O技术的理论基础是描述平面波在水平分界面上反射和透射的中国论文网(}%B,Ihx
Zoeppritz方程。尽管该方程早在20世纪初就已经建立，但由于其数学上的复杂中国论文网^t\)B"vWjj
性和物理上的非直观性，一直没有得到直接的应用。尽管该方程早在20世纪初中国论文网t,\j&ja#]
就已经建立，但由于其数学上的复杂性和物理上的非直观性，一直没有得到直
3}J;F#FK8Sv0接的应用。为了克服由Zoeppritz方程导出的反射系数形式复杂及不易进行数中国论文网 jHJGG/EINP
值计算的困难，许多学者对Zoeppritz方程进行了简化，并从不同的各个方面帮
$x-P"IP
D\d0助我们理解岩性参数，纵横波速度，密度和泊松比对于反射系数的影响。中国论文网1D7hg9@`
2.1.2 简化的Zoeppritz方程中国论文网 Th fv8\$P
AVO技术的核心思想是利用在不同介质中，反射系数随入射角的变化规
Qy D'r;v5K,K`0律来寻找油气层。因此，必须建立一个具有普遍意义的方程，将反射系数表示成中国论文网)y*O?Kr3^ }
入射角和地层参数的函数。完全形式的Zoeppritz方程全面考虑了平面纵波和中国论文网-@#Wv1F$Wc5_
横波入射在平界面两侧产生的纵横波和透射能量之间的关系。当反射波地震勘中国论文网 }&[0~j$al%V N K
探使用主要产生纵波的震源，接收的是反射纵波时，完全的Zoeppritz方程可以中国论文网Q$O%{[+X
被大大的简化，即只考虑平面纵波入射产生的反射振幅随入射角的变化情况。
#KI/Sye0一方面可以减少计算量，另一方面更有利于AVO技术的研究和应用。在这方面，中国论文网y;f'Fq$o1s
前人做了大量的工作，从不同角度对Zoeppritz方程进行了简化，提出了不同纵中国论文网oVn0zB!\'p,n
波反射振幅的近似表达式。
D9K${,@yw0第三章 结论与认识中国论文网(Z3OI W&H k;@L#i
利用AVO反演来提取储层参数，是目前进行叠前反演的一种有效的手段
x aC-n+G
Ycb)s0之一。而通过AVO反演获取的储层参数，如泊松比等，可以直接用来预测油气。中国论文网o/x u[,OPc
本论文通过
4t C;eD)e'DG W|0可以得到以下几点结论：中国论文网T%^!PF].P,]g |;[
1.本文利用了在AVO参数反演中常用的Zoeppritz简化方程Aki和Richard
1`L\`3m0方程来进行反演。通过模型数据的计算可以发现简化方程可以准确的反演各层中国论文网[ K!d"@9Zpa0fr
的反射系数并可以达到很好的反演效果。
6QT*g$dFRz02.常规的最优化方法依赖于初始模型的选取，易陷于局部最优，无法满足中国论文网y-H_k$Y+F
地震反演问题的需要。而模拟退火算法避免常规线性反演方法的缺点，不依赖
N R0zEz~ p0于初始模型的选取，同时可以很好的跳出局部最优解，适合应用于解决相关的中国论文网W$KZ;uo/H3lX
地震反演问题。中国论文网1l9A1n6S'|4sT.I:uR
3.通过前人对于三类含油气岩石的AVO特性的总结，可以较方便从实际
vL$ig A#wA^jb~0的资料上选取AVO特性明显的目的层进行反演，可以减少计算量,而且还可以中国论文网2TrXi,{
利用含油气的AVO特性对目的层段的反演结果进行检验。
g5j)Wl[KlkL9^04.在AVO参数的反演中，影响反演结果的因素很多。如Zoppritz方程的简中国论文网"`-E y-a}P2s1k
化方程，反演算法，AVO正演模拟等各个环节都会影响到结果的好坏。因此，可中国论文网;C|:{k!G'wmt
以通过对于这些环节的改进使反演结果更为精确。