致密储层中油气运移的数值模拟方法及装置的制造方法
【技术领域】
[0001] 本发明设及石油技术领域，特别设及一种致密储层中油气运移的数值模拟方法及 装置。
【背景技术】
[0002] 作为一种有效的接替资源，致密油在全球得W成功勘探和开发。该是自上个世纪 80年代W来，致密砂岩气藏之后，石油领域的又一重要突破。
[0003] 油气的运移、聚集是油气勘探的重点内容，深刻影响着油气勘探的方法及思路。作 为发生在地质历史时期的动态过程，致密储层内油气运移不可能直接观察；受实验技术条 件的限制，大尺度、长时间、可W考虑多种因素影响的实验往往也无法实现。因此数学模拟 技术正在成为油气运移研究的重要手段。
[0004] 由于研究的尺度、动力学条件或者基本物理条件的侧重不同，相应的数学模型也 有很大差异。从应用的角度来看，经过近几十年的完善，人们的认识也在不断深化，针对不 同盆地类型、地质条件、孔隙介质、流体性质的数学模型不断被提出。
[0005] 孔隙介质内多相渗流的理论最早被应用于油气运移的分析中，油气的运移过程被 视为等在宏观上均一的多相流动过程。因而可W用多相渗流的达西定律对流动过程进行分 析，多孔介质的达西流和流迹是目前较为成熟的模拟方法。超压作为油气向低渗透致密储 层中充注运移的主要动力，是致密储层石油成藏的决定性因素。超压控制下的流体流动具 有不连续性及流体活动的多期次性和周期性的特点，并W幕式排放的样式进行充注。油气 运动样式属于非达西流，该样使得基于渗流理论的模型在建立可靠的连续性方程时会有很 大的难度，受渗透率非均质的控制，运移路径难W有效预测。
【发明内容】
[0006] 为了提高优势通道模拟的预测精度，本发明提供了一种致密储层中油气运移的数 值模拟方法，所述方法包括：
[0007] S1 ;获取所述致密储层的构造格架图像及姪源岩数据，对所述构造格架图像进行 网格化处理；
[000引S2;根据所述姪源岩数据在所述构造格架图像中确定姪源岩网格，根据所述姪源 岩网格的生姪强度计算各姪源岩网格的排姪概率；
[0009]S3;将预设数量的人工妈蚁依次放置于按照所述排姪概率随机选择的姪源岩网格 处；
[0010]S4 ;计算各人工妈蚁所在网格与其相邻网格之间的动力差；
[0011]S5;对各人工妈蚁进行遍历，判断遍历到的当前人工妈蚁所在网格与其相邻网格 之间的动力差是否均小于等于预设阔值，若所述动力差均小于等于预设阔值，则将所述当 前人工妈蚁所在网格记录至富集区集合，并将所述当前人工妈蚁的路径保存至路径集合 中，直接执行步骤S7,否则执行步骤S6 ;
[0012] S6;根据网格的信息素及所述动力差计算所述当前人工妈蚁从其所在网格运移至 各相邻网格的运移概率，将所述当前人工妈蚁从其所在网格运移至根据运移概率从所述相 邻网格中随机选择的一个相邻网格，并将所述当前人工妈蚁运移至该相邻网格记录到所述 当前人工妈蚁的路径中；
[0013] S7 ;判断是否所有人工妈蚁所在网格与其相邻网格之间的动力差均小于等于预设 阔值，若是，则执行步骤S8,否则返回步骤S5 ;
[0014] S8 ;根据所述路径集合对各网格的信息素进行更新；
[0015] S9 ;返回步骤S3,直至返回步骤S3的次数满足预设次数后，执行步骤S10 ;
[0016] S10 ;将所述路径集合及富集区集合作为模拟结果。
[0017] 其中，步骤S4中，通过下式计算各人工妈蚁所在网格与其相邻网格之间的动力 差，
[001 引
[0019] 其中，i为人工妈蚁所在网格；j为所述人工妈蚁所在网格的相邻网格；?u(t)为 人工妈蚁所在网格i与相邻网格j之间的动力差；AP为人工妈蚁所在网格i与该相邻网 格j之间的剩余压差；曰为油水两相界面张力；Rj为该相邻网格j内的岩石孔隙半径；Ri为 人工妈蚁所在网格i内的岩石孔隙半径。
[0020] 其中，步骤S6中，根据所述网格的信息素及所述动力差通过下式计算所述当前妈 蚁从其所在网格运移至各相邻网格的运移概率，
[0021]
[002引其中，Pu(t)为从当前人工妈蚁所在网格i运移至相邻网格j的运移概率；Tu(t) 为该相邻网格j当前的信息素；a为继承启发因子；0为动力启发因子；allowed为与所述 当前人工妈蚁所在网格相邻的所有网格。
[0023] 其中，步骤S8中，根据所述路径集合通过下式对各网格的信息素进行更新，
[0024] Tu(t) = (1-丫）Tu(t-l)+M?AT
[002引其中，Ty(t-l)为网格j在上一次更新时的信息素；丫为信息数挥发比例；M为从 上一次更新到当前时刻之间，所述路径集合增加的路径中人工妈蚁运移至网格j的次数； AT为一次转运过程中信息素的增加量。
[0026] 其中，步骤S10之后还包括：
[0027]S11 ;通过下式计算所述路径集合中各条路径的通道系数，
[002引
[0029] 其中，C为路径的通道系数，T为该路径中的网格数，E〇z(t)为该路径中各相邻 网格之间的动力差的总和。
[0030] 本发明还公开了一种致密储层中油气运移的数值模拟装置，所述装置包括：
[0031] 图像获取单元，用于获取所述致密储层的构造格架图像及姪源岩数据，对所述构 造格架图像进行网格化处理；
[0032] 概率计算单元，用于根据所述姪源岩数据在所述构造格架图像中确定姪源岩网 格，根据所述姪源岩网格的生姪强度计算各姪源岩网格的排姪概率；
[0033] 姪源岩选择单元，用于将预设数量的人工妈蚁依次放置于按照所述排姪概率随机 选择的姪源岩网格处；
[0034] 动力差计算单元，用于计算各人工妈蚁所在网格与其相邻网格之间的动力差；
[0035] 路径保存单元，用于对各人工妈蚁进行遍历，判断遍历到的当前人工妈蚁所在网 格与其相邻网格之间的动力差是否均小于等于预设阔值，若所述动力差均小于等于预设阔 值，则将所述当前人工妈蚁所在网格记录至富集区集合，并将所述当前人工妈蚁的路径保 存至路径集合中；
[0036] 运移记录单元，用于根据网格的信息素及所述动力差计算所述当前人工妈蚁从其 所在网格运移至各相邻网格的运移概率，将所述当前人工妈蚁从其所在网格运移至根据运 移概率从所述相邻网格中随机选择的一个相邻网格，并将所述当前人工妈蚁运移至该相邻 网格记录到所述当前人工妈蚁的路径中；
[0037] 动力差判断单元，用于判断是否所有人工妈蚁所在网格与其相邻网格之间的动力 差均小于等于预设阔值；
[003引信息素更新单元，用于根据所述路径集合对各网格的信息素进行更新；
[0039] 调用判断单元，用于调用所述姪源岩选择单元，直至调用所述姪源岩选择单元的 次数满足预设次数后，调用系数计算单元；
[0040] 结果确定单元，用于将所述路径集合及富集区集合作为模拟结果。
[0041] 其中，所述动力差计算单元通过下式计算各人工妈蚁所在网格与其相邻网格之间 的动力差，
[0042]
[0043] 其中，i为人工妈蚁所在网格；j为所述人工妈蚁所在网格的相邻网格；?u(t)为 人工妈蚁所在网格i与相邻网格j之间的动力差；AP为人工妈蚁所在网格i与该相邻网 格j之间的剩余压差；曰为油水两相界面张力；Rj为该相邻网格j内的岩石孔隙半径；Ri为 妈蚁所在网格i内的岩石孔隙半径。
[0044] 其中，所述运移记录单元根据所述路径集合及所述动力差通过下式计算所述当前 人工妈蚁从其所在网格运移至各相邻网格的运移概率，
[0045]
[0046] 其中，pu(t)为从当前人工妈蚁所在网格i运移至相邻网格j的运移概率；Tu(t) 为该相邻网格j当前的信息素；a为继承启发因子；0为动力启发因子；allowed为与所述 当前人工妈蚁所在网格相邻的所有网格。
[0047] 其中，所述信息素更新单元根据所述路径集合通过下式对各网格的信息素进行更 新，
[0048] Tu(t)=(1-丫）Tu(t-l)+M?AT
[0049]其中，Tu(t-l)为网格j在上一次更新时的信息素；丫为信息数挥发比例；M为从 上一次更新到当前时刻之间，所述路径集合增加的路径中人工妈蚁运移至网格j的次数； AT为一次转运过程中信息素的增加量。
[0化0] 其中，所述装置还包括：
[0化1] 系数计算单元，用于通过下式计算所述路径集合中各条路径的通道系数，
[0化2]
[005引其中，C为路径的通道系数，T为该路径中的网格数，E0z(t)为该路径中各相邻 网格之间的动力差的总和。
[0化4]本发明将姪源岩网格等效为蚁窝，根据生姪强度从多个姪源岩释放代表油气的人 工妈蚁，并对各人工妈蚁的路径及停止运转的网格（即富集区集合中的网格）进行记录， 获得将所有人工妈蚁的路径及富集区集合，方式本身相对而言比基于达西渗流方式模拟的 方法更为简便，效果更为接近现实条件，预测精度高，同时，由于方法本身较为简便，计算迅 速。
【附图说明】