模拟退火算法数字岩心建模方法改进.pdfVIP

————弋、查全兰塑竺苎篓查查堕圭璺堡垒堡竺苎墨i 模拟退火算法数字岩心建模的方法改进 赵秀才“ 姚军1 衣艳静2 (1．中国石油大学(华东)石油工程学院，山东东营，257061 2．中油国际海外研究中心，北京，100083) 摘要：系统地介绍了模拟遇火算法构建教宇岩心的方法，包括建模参考函数、建模过程及重要参数的设 置方法。研兜发现建模规模的增大、建模参考函数的欠缺及教字岩心各向同性的假设导致算法生成的擞 字岩心与真实岩一0差异鞍大。通过在退火过程中适时地清除孤立的岩石颗粒可l；【较好地降低系统能量、 有效地优化孔暾空闻、加快退火速度。以计算机现有的运算速度，若考虑多个褒杂建模函最及系统的各向 异性。应用模担退^算法进行截字岩心建模特需要极大的时间开支，因此，该方法有待于进一步改进。 关键词：模拟退大；教字岩心；孤立岩石颗粒；系统能量 0引言 从岩石的微观结构出发预测其宏观性质是近几年国内外学者较为重视的一个研究领 域，数字岩心作为联系岩石微观结构和宏观性质的重要桥梁发挥着至关重要的作用。迄今 为止，用于构建数字岩心的方法有多种“]。模拟退火算法(simulatedAnnealingA190— rithm，下文以sAA简称)作为以岩心薄片分析为基础的图像重建方法之一早有学者研 究-2“J。由于它不仅能够考虑多个建模参考函数、丰富建模输入资料，还可适用于不同岩石 类型，因而较高斯场[53和基于过程模拟160的建模方法具有一定优势。本文详细介绍了sAA 的数字岩心建模方法，指出了其中的不足并进行了有效的改进。 l 基于SAA的数字岩心建模方法 1．1建模参考函数 本文所要建立的数字岩心为各向同性且具有统计均质性的两相系统(仅考虑孔隙和岩 石骨架)．建模时各个参考函数均以孔隙空间为分析对象，此时该数字岩心系统可由相函数 z(r)描述。以Z(r)为基础，岩心孔隙度≠可由统计平均值给出： zcj，=E：：凳熹 ㈤ *作者简介：赵秀才。男，1982年生，山东沂水县人，现为中国石油太学(华东)油气田开发专业在读博士 maLl：‘haoxiucan234@126．c口Tn， 《堪血世孤葬碜数每辱心谨画面万话甄再Z j囫l ≠一百两 (2) 式中，r为系统中的任意～点，一表示统计平均(下文与此相同)。除孔隙度外，两点概率函数 要函数，它们的定义如下： s(r)=Z百r—而i干万 (3) l 屯∈孔隙 L(r)=P(r，r+r)，P(r，，+r)= (4) O 其中，‘是线段[r，r+r]上的任意一点；r为间隔距离．为计算方便取，的单位为体素数。体 素数与体素棱长的乘积即物理长度，下文所有长度单位与此相同。对于二维图像，r的单位 是像素数。 分析公式(3)可知·函数s具有以下性质：s(o)一≠；对于具有统计均质性的系统，s(r) 值将随r的增大而逐步趋于稳定，反映出统计均质性。一般情况下，函数s曲线的长度截取 至曲线达到某一稳定值或水平波动不大时即可。函数L是描述多孔介质内同一相连通性能 的计算长度选取至其值减小到。的值即可。 1．2建模过程 本文以CT图像为基车资料，应用≠，s及L函数建模。≠．s，L均可通过用于建模的岩心 图像计算得到。设s。，L，分别为sAA建立的数字岩心系统中与s，L相对应的函数，则sAA 建模实际上是对系统不断优化直至S，，L，与s，L足够逼近为止。 首先随机生成孔隙度为≠的三维o，1(o表示岩石骨架，l表示孔隙)数字岩心系统，启动 退火进程。设某时刻计算出第^个系统的s，，L。，通过下式计算出此时系统(称为原系统)的 能量ET：