利用多点地质统计学方法模拟岩相分布
PPT
汇报提纲
研究现状 目标与内容 过程论述 结论
研究现状
多点地质统计学是目前储层地质随机建模的研 究热点方向。相对于传统的两点地质统计学方 法，多点地质统计学在进行储层地质建模时能 够利用更多空间点数据（n>2）之间的相关性 进行统计计算，因此建立的模型更加满足实际 的情况。 多点地质统计学算法snesim是其中一种基于 概率统计的方法，由于snesim算法的参数设 置较为复杂，因此该算法的各种参数对模型的 控制需要进行研究和分析，从而指导如何使用 该算法建立高质量的储层模型。
硬数据条件化
向网格中加载硬数据：
硬数据条件化模拟结果：
非平稳模拟
旋转角控制河道局部的方向 仿射性控制河道宽度
手绘旋转与仿射性图像：
1.3
0.8 0.5 15°
75°
45°
格式转换：
参数文件复制：
参数设置：
运行得到实现：
涠洲11-2油田4井区岩相模拟
研究区属于辫状河三角洲前缘亚相,微相种类较多,面积较大，但 是钻井较少且分布集中,为了减少没有井控制的地方模拟的不确 定性,本次采用录井解释的砂泥岩相进行岩相模拟，相的种类减 少，有助于减少模拟的不确定性.
利用不同油组的平面沉积微相图作为训练图像：
序贯指示法模拟结果：
snesim算法模拟结果：
用多点法模拟的岩相展布图较为连续且具有一定的真实性
结论
设置不同参数对程序的影响： 目标比例能有效控制模拟结果相比例 目标相比例的设置应与训练图像边际 概率相近，否则要考虑更换训练图像 旋转角控制河道局部的方向，仿射性 控制河道宽度 用多点法模拟的岩相展布图较为连续 且具有一定的真实性
多点地质统计学原理、方法及应用__概述及解释说明
多点地质统计学原理、方法及应用概述及解释说明1. 引言1.1 概述本文旨在探讨多点地质统计学的原理、方法及应用，为读者提供一个全面了解该领域的概述。
多点地质统计学是一门研究如何有效地利用多变量数值以及空间数据进行地质分析和预测的学科。
它通过综合多种数据，包括物理测量数据、遥感图像数据和野外调查数据等，来实现对不同地质现象和过程的建模与研究。
1.2 文章结构本文按照以下结构组织内容：首先介绍多点地质统计学的基本原理，包括其定义与概念、基本假设以及原理解释。
随后，针对多点地质统计学的方法进行详细阐述，探讨数据收集与预处理、变量选择和缺失值处理以及统计模型拟合与优化算法应用等关键步骤。
接下来，我们将通过具体案例研究来展示多点地质统计学在矿产资源评估与勘探、地下水资源管理与保护以及石油勘探与开发中的应用实践。
最后，在结论部分对全文进行概括总结，并展望未来多点地质统计学研究的发展方向。
1.3 目的本文旨在全面介绍多点地质统计学的原理、方法及应用，以帮助读者对该领域有一个清晰的认识。
通过阐述基本原理和方法，读者可以了解多点地质统计学在地质分析和预测中的重要性。
此外，通过具体案例的引入，读者将能够更好地理解多点地质统计学在实际问题中的应用价值和潜力。
最后，通过对未来研究方向的展望，读者可以获得一些启示，并为自己在该领域开展研究提供参考。
2. 多点地质统计学原理2.1 定义与概念多点地质统计学是一种广泛应用于地质科学领域的统计学方法。
它通过对多个地点上的地质数据进行收集、分析和解释，旨在揭示地下资源的分布规律和空间变异性。
多点地质统计学基于一系列假设和方法，能够提供可靠的预测结果和决策依据。
2.2 基本假设在多点地质统计学中，存在几个基本假设：- 空间自相关假设：相邻位置上的地质现象存在关联性，即一个位置的观测值可能受到相邻位置观测值的影响。
- 空间平稳假设：在整个研究区域内，不同位置上的地质变量具有类似的变异性。
由两点到多点的地质统计学储层建模 (2)
断块油气田2012年9月断块油气田FAULT -BLOCK OIL ＆GAS FIELD 由两点到多点的地质统计学储层建模陈培元1，姜楠1，杨辉廷1，刘学利2（1.西南石油大学资源与环境学院，四川成都610500；2.中国石化西北油田分公司，新疆乌鲁木齐830011）基金项目：“十二五”国家科技重大专项“大型油气田及煤层气开发”子课题“塔里木盆地大型碳酸盐岩油气田勘探开发示范工程”（2011ZX05049-04）摘要传统的两点地质统计学建模方法，以象元为空间赋值单元、变差函数为工具建立确定性的模型，或者应用各种随机模拟方法建立可选的模型，在精确表征复杂的空间结构及目标体几何形态方面有一定的局限性。
有别于两点地质统计学的多点地质统计学，可有效地解决更广泛的地质模拟问题。
然而，在实际应用过程中，受岩-相模型及与之相对应的训练图像可靠性的影响，结果变得比较复杂。
因此，选择合适的训练图像及恰当的算法可有助于提高储层建模的精度和效率。
以××油田曲流河沉积为例，采用两点和多点统计学方法构建模型。
对比发现，基于多点地质统计学的地质建模方法真实可再现河流相的沉积形态，还降低随机建模的不确定性。
尽管模拟结果与井点真实数据之间存在误差，但通过调整随搜索半径、训练图像大小及概率计算中临近点个数限制，可显著提高模型精度。
关键词两点地质统计学；多点地质统计学；储层；随机模拟中图分类号：TE319文献标志码：A收稿日期：2012-04-01；改回日期：2012-07-10。
作者简介：陈培元，男，1984年生，博士，主要从事油藏描述和油藏地质建模研究。
E -mail ：swpua409@ 。
引用格式：陈培元，姜楠，杨辉廷，等.由两点到多点的地质统计学储层建模［J ］.断块油气田，2012，19（5）：596-599.Chen Peiyuan ，Jiang Nan ，Yang Huiting ，et al.Reservoir stochastic modeling using geostatistics from two -point to multiple -point ［J ］.Fault -Block Oil &Gas Field ，2012，19（5）：596-599.Reservoir stochastic modeling using geostatistics from two -point to multiple -pointChen Peiyuan 1,Jiang Nan 1,Yang Huiting 1,Liu Xueli 2(1.School of Resources and Environment,Southwest Petroleum University,Chengdu 610500,China;2.Northwest Oilfield Company,SINOPEC,Urumqi 830011,China)Abstract:Conventional two -point geostatistics modeling mainly uses the pixel -based method and variogram to establish the deterministic model,or uses all kinds of stochastic simulation methods to establish optional model.But it can not fully reflect the variability of the space structure and geometric shape of object.Unlike the two -point geostatistics,the multiple -point geostatistics can solve the problems of geologic simulation widely.In actual application process,due to the lithofacies and the reliability of corresponding training image,it is necessary choosing the suitable training image and appropriate algorithm to improve the accuracy and efficiency of simulation.Taking the meandering river sedimentation of some oilfield as an example,the reservoir model is built by two -point and multiple -point geostatistics parison results of two models show that the method based on the multiple -point not only represents the real sedimentary form of fluvial facies,but also reduces the uncertainty of stochastic modeling effectively and improves the modeling accuracy at maximum.Although the simulation results do not agree with the real well date completely,the prediction accuracy can be improved through adjusting the search radius,the size of training image and the number limit of nearest -neighbor points used in probability calculation.Key words:two -point geostatistics;multiple -point geostatistics;reservoir;stochastic simulation近年来，随着储层建模技术的不断发展，人们对储层地质模型的要求越来越高。
储层多点地质统计学随机建模方法
储层多点地质统计学随机建模方法摘要：多点地质统计学使用训练图像代替变差函数，将更多的地质资料整合到储层建模过程中，使得最终模型更加符合地质认识。
随着研究的不断深入，越来越多的地质工作人员开始熟悉这一方法，凭借自身的独特优势，多点地质统计学将在储层建模领域占得重要的一席。
关键词：多点地质统计学训练图像储层建模一、多点地质统计学与训练图像基于变差函数的传统地质统计学随机模拟是目前储层非均质性模拟的常用方法。
然而，变差函数只能建立空间两点之间的相关性，难于描述具有复杂空间结构和几何形态的地质体的连续性和变异性。
针对这一问题，多点地质统计学方法应运而生。
该方法着重表达空间中多点之间的相关性，能够有效克服传统地质统计学在描述空间形态较复杂的地质体方面的不足。
多点地质统计学的基本工具是训练图像，其地位相当于传统地质统计学中的变差函数。
对于沉积相建模而言，训练图像相当于定量的相模式，实质上就是一个包含有相接触关系的数字化先验地质模型，其中包含的相接触关系是建模者认为一定存在于实际储层中的。
二、地质概念模型转换成图像训练地质工作人员擅于根据自己的先验认识、专业知识或现有的类比数据库来建立储层的概念模型。
当地质工作人员认为某些特定的概念模型可以反映实际储层的沉积微相接触关系时，这些概念模型就可以转换或直接作为训练图像来使用。
利用训练图像整合先验地质认识，并在储层建模过程中引导井间相的预测，是多点地质统计学模拟的一个突破性贡献。
可以将训练图像看作是一个显示空间中相分布模式的定量且直观的先验模型。
地质解释成果图、遥感数据或手绘草图都可以作为训练图像或建立训练图像的要素来使用。
理想状态下，应当建立一个训练图像库，这样一来建模人员就可以直接选取和使用那些包含目标储层典型沉积模式的训练图像，而不需要每次都重新制作训练图像。
三、多点模拟原理进行多点模拟，需要使用地质统计学中的序贯模拟。
但是，多点模拟与传统的基于变差函数的两点模拟是不同的。
利用多点地质统计学方法进行垦西71断块沉积微相建模
[ 6] 尹艳树 , 吴胜和 , 张昌民, 等 164
Hale Waihona Puke 基于储层骨架的多点地质统计学方法 [ J]
[ 7] 张伟 , 林承焰 , 董春梅 多点地质统计学在秘鲁 D 油田地质建模中的应用 [ J] 24~ 28 [ 8] 郭秀蓉 , 程守田 , 刘星 55~ 58
油藏描述中的小层划分与对比 ! ! ! 以垦西油田 K 71 断块东营组为例 [ J ]
2
2 1
垦西 71 断块沉积微相建模
垦西 71 断块东二段地质概况 垦西 71 断块东营组砂体为一套三角洲前缘亚相沉积。根据各沉积时间单元的岩性组合特征、砂体
形态、厚度及电性特征等 , 东二段 ( Ed2 ) 主要发育了水下分流河道沉积、分流河道边缘沉积和分流间 湾沉积等 3 种微相[ 8] 。
[ 收稿日期 ] 2010 03 01 [ 基金项目 ] 国家 863 计划项目 ( 2007A A 060500) 。 [ 作者简介 ] 杨宏伟 ( 1967 ) , 男 , 1989 年大学毕业 , 高级工程师 , 博士生 , 现主要从事综合解释及油藏建模方面的研究工作。
第 32 卷第 6 期
杨宏伟 : 利用多点地质统计学方法进行垦西 71 断块沉积微相建模
∃ 225 ∃
2 2
训练图像的建立
1) 水下分流河道宏观统计特征 沉积微相研究表明垦西 71 断块东 二段沉积时 期, 河 道砂 体呈弯 曲条 带状 , 剖面 呈顶 平底 凸 透镜 状, 河 道宽一般 100 ~ 300m, 砂体厚 度一 般 4~ 5m, 河 道弯曲 度小。河 道边 缘砂呈透镜 状, 单 一规 模的河 道边 缘砂规模一般不超过 200m 。物源方 向为北东向。 2) 水下分流河道及河道边缘微 相的地质统计 对 59 口井的数据分
不同类型河道砂体构型要素划分方法及三维表征技术研究
累积概率（％）
0.35 0.3
0.25 0.2
0.15 0.1
0.05 0
1～2
2～3
分支河道砂体厚度概率图
3～4
4～5
5～6
砂体厚度（m）
6～7
1.2 1 0.8 0.6 0.4 0.2 0 7～8
35 30 25 20 15 10
5 0
50～100
主河道河道砂体宽度概率图
101～150
151～200 201～250 河道砂体宽度（m）
建立知识库思路
搭资以建 料 知架的识构收库，集架确和构定学为研习基究阶础内段，容，通。 结调过合研数沉国据积内库学外软知件识完，库成总研知结究识野成入外果库， 露提实头供现知搭网识建络库知功，识能并库，对架指比构导油的地田思质实路基， 际总础区结研块野究解外工剖露作，头。充和实现知代识沉库积内资 容料。，为野外露头知识库的建 立作准备。
三维表征是目前进行构型要素表征的主要手段 之一，从地质统计学角度来看，可分为三个方面：
(1) 基于象元的两点统计学算法 (2)基于目标的示性点过程算法 (3)多点地质统计学算法
该方法相对简单、成熟，一般是采 用变差函数的方式模拟，变差函数分 析结对果于的基优于劣目将标会的直示接性影点响过到程最算终法模而 型言的，好需坏要。填在写进地行质变体差的函几数何分特析征时参，数， 往不影这能质于面多往可响些否条河比点都避，参建件道传地采免表数立的砂统质用地征需起地体的统人会结要一质构随计机对果从套知型机学交变存地完识表建具互差在质备库征模有模函不知的，、方独拟数确识、将条法特方产定库符会件有的式生性中合直化较优，人。得实接数大势因为到际影据提，此的。地响方高对，， 到但最算终法模仍型在的完好善坏之。中，特别是计算机
复杂断块油藏的三维地质建模方法--以国外某油田为例
复杂断块油藏的三维地质建模方法--以国外某油田为例金春玉;宋扬【摘要】随着油气勘探开发的不断深入，一些具有复杂断块构造特征的油藏已逐渐成为开发的主体目标。
建立精准的复杂断块油藏三维地质模型对指导油田开发有着重要的意义。
三维地质建模技术能够更细致、准确地研究地下的油藏，为油藏的后续开发提供可靠的地质依据。
以国外某油田断块构造发育区块为研究对象，针对研究区地质条件复杂、断层发育的特点，以地质、地震、测井资料为基础，搭建构造框架模型，应用地质统计学理论建立储层岩相及物性参数模型，揭示构造和储层空间分布特征，最终建立一个三维定量的油藏地质模型，在储层计算中各断块误差均小于5%，符合精度要求。
该模型将为该区块数值模拟和井位设计等后续工作提供可靠的地质依据，同时也为同类复杂断块油藏的三维地质建模提供借鉴。
%With the development of oil and gas exploration,some complex fault-block reservoirs have be-come the main objective.3D geologic modeling of complex fault-block reservoirs is an important significance to guide oilfield development.Study on underground reservoirs,3D geologic modeling technology is more particular and accurate,which supplies reliable data for future development.Study on an overseas oilfield which has the complex fault-block,the geologic condition of the complex fault-block is complex.Based on geological,seismic and logging data,the structural modeling had been built.The lithofacies modeling and physical property modeling had been built by geostatistics theory which reveal the structural attitude and spatial distribution.3D geologic modeling of the block had been built at last.The results showed that the reserves error ofthese fault-blocks were within 5%,which meet the requirement.The 3D geologic modeling of the block will provide reliable geological basis forthe numerical simulation and the location and so on, and provide some useful reference for 3D geologic modeling technology of similar complex fault-block res-ervoirs.【期刊名称】《河北联合大学学报（自然科学版）》【年(卷),期】2015(000)001【总页数】6页(P15-20)【关键词】地质建模;复杂断块;构造模型【作者】金春玉;宋扬【作者单位】河北联合大学，河北唐山 063009;大庆油田有限责任公司勘探开发研究院，黑龙江大庆 163712【正文语种】中文【中图分类】TE12随着油气勘探开发的不断深入，一些具有复杂断块构造特征的油藏已逐渐成为开发的主体目标。
苏里格气田致密砂岩气藏剩余气分布特征及其挖潜
第44卷 第5期 新 疆 石 油 地 质Vol. 44，No.52023年10月 XINJIANG PETROLEUM GEOLOGY Oct. 2023文章编号：1001-3873（2023）05-0554-08 DOI ：10.7657/XJPG20230506苏里格气田致密砂岩气藏剩余气分布特征及其挖潜石耀东1，王丽琼1，臧苡澄2，张吉1，3，李鹏2，李旭1（1.中国石油 长庆油田分公司 第四采气厂，内蒙古 鄂尔多斯 017300；2.中国石油 长庆油田分公司 勘探开发研究院，西安 710018；3.低渗透油气田勘探开发国家工程实验室，西安 710018）摘 要：苏里格气田中区苏36-11区块已开发17年，开发程度和储量动用程度均高，储集层非均质性强，储量动用不均衡，剩余气分布复杂，剩余气分布的确定及挖潜是气田稳产的关键。
通过储集层构型精细表征，明确剩余气分布的主要影响因素，确定不同类型剩余气分布规律，提出对应的挖潜对策。
研究结果表明：研究区含气砂体主要分布在4级构型单元心滩坝与点坝中，整体规模小，宽度为150~500 m ，长度为300~800 m ，连通性差，受各级次渗流屏障影响大，区块北东—南西向主砂带开发程度最高，地层压力低，剩余气主要分布在区块西北部盒8段下亚段；剩余气分布主要受储集层非均质与开采非均匀影响，可分为井网未控制型、复合砂体阻流带型、水平井未动用次产层型、直定向井未射开气层型和投产未采出型5类；提出井间加密、老井侧钻、查层补孔和老井挖潜4种动用措施，调整方案后，预测可稳产7年，采收率可达45%。
关键词：苏里格气田；致密砂岩；储集层构型；剩余气储量评价；剩余气分布；挖潜对策；开发中—后期；开发调整方案中图分类号：TE122 文献标识码：A©2018 Xinjiang Petroleum Geology. Creative Commons Attribution-NonCommercial 4.0 International License 收稿日期：2022-11-12 修订日期：2023-04-13基金项目：国家科技重大专项（2016ZX05050）；中国石油科技重大专项（2016E-0509）第一作者：石耀东（1973-），男，陕西靖边人，高级工程师，气田开发与生产管理，（Tel ）************（E-mail ）syd_cq@通讯作者：王丽琼（1989-），女，甘肃华池人，高级工程师，硕士，油气田开发，（Tel ）************（E-mail ）wangliqiong12_cq@petrochina..Distribution and Potential Tapping Strategies of Remaining Gasin Tight Sandstone Gas ReservoirsSHI Yaodong 1，WANG Liqiong 1，ZANG Yicheng 2，ZHANG Ji 1，3，LI Peng 2，LI Xu 1（1.No.4 Gas Production Plant, Changqing Oilfield Company, PetroChina, Ordos, Inner Mongolia 017300, China;2.Research Institute of Exploration and Development, Changqing Oilfield Company, PetroChina, Xi ’an, Shaanxi 710018, China ；3.National Engineering Laboratory for Exploration and Development of Low Permeability Oil and Gas Fields, Xi ’an, Shaanxi 710018, China ）Abstract ：The Su 36⁃11 block in the central area of Sulige gas field has been developed for 17 years, with high degrees of development and reserves producing. The strong reservoir heterogeneity in this block leads to uneven producing of reserves and complex distribution of re⁃maining gas. Distribution determination and potential tapping of the remaining gas are crucial for maintaining stable production in the gas field. By accurately characterizing the reservoir architecture, the main factors influencing remaining gas distribution were identified, the distribution patterns of different types of remaining gas were determined, and corresponding strategies for recovering the remaining gas were proposed. The research results show that the gas⁃bearing sand bodies in the study area are mainly distributed in the 4th⁃order architec⁃ture units, such as channel bar and point bar, these sand bodies are significantly affected by various levels of flow barriers, with small over⁃all scale, poor connectivity, width of 150-500 m and length of 300-800 m. The main NE⁃SW sand belt in the block has been developed the most, with low formation pressure, and the remaining gas is mainly distributed in the lower He 8 member in the northwestern part of the block. Remaining gas, whose distribution is mainly influenced by reservoir heterogeneity and uneven development, can be divided into five types: gas uncontrolled by well pattern, gas in composite sand body flow barrier, gas in secondary pay zone unexploited by horizontal well, gas in unperforated gas⁃bearing layer in vertical well, and gas unproduced. Four potential tapping measures were proposed, including well infilling, reperforation, sidetracking and potential tapping in exsisting wells. According to the adjusted development plan, it is predicted that stable production can be maintained for 7 years with the recovery efficiency reaching 45%.Keywords ：Sulige gas field; tight sandstone; reservoir architecture; remaining gas reserves evaluation; remaining gas distribution; potential tapping; middle-late development stage; adjusted development plan中国致密气资源总量及开发潜力巨大，约占全球资源量的十分之一，主要分布在鄂尔多斯盆地、四川盆地、塔里木盆地等区域。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
收稿日期:2005-01-18基金项目:国家教育部高等学校优秀青年教师教学科研奖励计划(T RA POYT )资助作者简介:冯国庆(1974-),男,山东荷泽人,博士,主要从事油藏描述和油藏数值模拟研究.文章编号:1673-064X(2005)05-0009-03利用多点地质统计学方法模拟岩相分布Stochastic simulation of lithofacies distribution using multi -point g eostatistics冯国庆1,陈浩2,张烈辉1,李允1(1.西南石油学院油气藏地质及开发工程国家重点实验室,四川成都610500;2.新疆油田分公司陆梁油田作业区,新疆克拉玛依834000)摘要:文中介绍了多点地质统计学的基本原理,并利用该算法模拟了我国东部某砂岩油藏的岩相分布.通过对储层非均质性模拟研究,揭示储层在岩性、物性和含油气性的各向异性或非均质性,揭示砂体展布、连通程度以及在横向、纵向上的变化规律,有利于开发方案的制定及注采井网的布署.模拟结果表明,多点地质统计学模拟方法可以很好地再现砂体局部的变异性和非均质性.由此可得出以下结论:snesim 模拟算法可以快速灵活地进行多点地质统计模拟,多点地质统计学比两点地质统计学更具有优势.关键词:多点地质统计学;随机模拟;训练图像中图分类号:TE19;T E121.1+5 文献标识码:A 传统的地质统计学是利用变差函数描述地质变量的相关性和变异性,通过建立在某个方向上两点之间的地质变量的变化关系来描述空间的变化特性.但是,建立在两点统计关系上的变差函数本身在描述储层非均质性上有很大的不足,它反映的仅仅是空间两点之间的相关性,不能充分描述复杂几何形状砂体如河道砂体和冲积扇砂体空间的连续性和变异性.当井资料较少时,用于计算实验变差函数的点对很少,它也就不能正确反应空间两点之间的相关性.建立在两点基础上的变差函数在储层地质建模中存在一定的不足,而多点地质统计学[1-4]是建立在多个点的相关关系上,所以它在解决描述空间变量的连续性和变异性方面得到越来越广泛的应用.斯坦福大学的Journel 教授曾指出多点地质统计学是今后地质统计学发展的方向,它的优势已越来越显著.1 多点地质统计学的原理在阐述多点地质统计学之前,首先回顾一下变差函数的地质统计学方法是如何模拟储层岩相分布的,以序贯指示模拟算法为例进行说明.该方法的基本原理[5]简述如下:假设在模拟区域有k 种岩相s 1,s 2,,s k ,对于模拟目标区域内的每一相,定义指示变量:I (u )=1 Z(u)I s k0 Z(u)|s k,对于任一待模拟点,其出现第k 种相的概率为:P (I k =1|Z(u A )=s A ,P A ),A 为待估点所包括的条件区域,利用两点地质统计学方法计算该概率是采用克里格方法:P (I k =1|Z(u A )=s A ,P A )=E (I k )+E nA =1KA [1-E (I k )].2005年9月第20卷第5期西安石油大学学报(自然科学版)Journal of Xi c an Shiy ou U niversity(N atural Science Edition)Sep.2005V ol.20No.5其中,K A 为克里格方法确定的权系数,它通过求解由变差函数或协方差函数建立的克里格方程组来确定.多点地质统计学与两点地质统计学的主要区别在于上面的概率的确定方法不同,它首先引入一训练图像,通过在训练图像中寻找与待估点内条件数据分布完全相同的事件的个数来确定概率分布,因此它可以反映出多个位置的联合变异性.例如,计算图1(a)中u 点的概率时,相应的条件数据场为d n ={z (u 1),z (u 2),z (u 3),z (u 4)},其基本方法是首先要在训练图像(b)中寻找与图(a)中数据分布完全相同的事件的个数,即要在训练图像中找出与图(a)几何完全相同的区域,同时在该区域中相同的位置处z (u 1),z (u 2),z (u 3),z (u 4)的值完全相同.在训练图像中一共找到4个既能满足条件数据u 1,u 2,u 3,u 4数值,同时又能满足它们分空间几何形状的事件,在这4个事件中,3个事件的u 点的值为0,只有1个事件中u 点值为1,因此u 点岩相为1的条件概率为P {u =1|d n }=1/4,而P {u =0|d n }=3/4,这样便可求出了u 点的条件概率.图1 多点统计方法估计条件概率的原理上述方法不仅考虑了区域内条件数据的值,而且也考虑了条件数据的几何形状.而两点地质统计学只是依靠z (u 1),z (u 2),z (u 3),z (u 4)的值及各点与u 点距离通过求解克里格方程组来确定u 点的概率,并没有考虑d n 的几何形状和各条件数据的配位关系.基于上述原理,Sebastien Strebelle 提出了snes -im 模拟算法[2],利用该算法可以快速、灵活地模拟岩相分布.该方法的具体步骤为:(1)利用非条件模拟建立三维训练图像;(2)定义通过所有待估结点的随机路径;(3)对随机路径中的任意待估点l(=1,2,,,l):¹定义查找范围内的条件数据;º保留邻区的数据点;»在训练图像中寻找与该区域内条件数据完全相同的事件,计算该点岩相的分布概率.¼由Monto Carlo 法得到位置处的一个模拟值;½将模拟结果归入条件指示数据集中.(4)重复上一步模拟,直到所有的点全被模拟.训练图像既可以通过非条件模拟求出,也可以通过该地区的地质露头资料分析得出.对训练图像的条件非模拟可以选择非条件的布尔模拟方法,其方法和原理参见文献[6].2 实例分析对于开发中后期的砂岩油藏储层参数模拟采用两阶段模拟方法可以较为准确地反映储层的非均质性,而/两步建模0的第一步就是要建立储层结构或流动单元模型,模拟沉积体在空间排列的复杂性;利用多点统计学模拟方法可以较好地完成砂体骨架模拟.模拟区域选择我国东部某砂岩油藏第15小层,在该层一共有64口井,测井资料解释结果表明有26口井钻遇砂体,另外38口井钻遇泥岩,砂体比例为40%.对岩相进行编码,砂岩为1,泥岩为0,图2为该层井位分布图.图2 第15小层井位和井点岩相图采用上述方法模拟砂体的分布.首先建立训练图像,运用布尔模拟方法,把砂体比例40%输入,为保证训练图像数据充足,网格划分为250@250@1,一共由62500个模拟数据组成,布尔模拟结果见图3.把条件数据和布尔模拟生成的训练图像,输入到snesim 模拟算法中进行模拟.根据该区域的特点,椭圆最大搜索半径选为300m,搜索半径内最多)10)西安石油大学学报(自然科学版)的条件数据设为30,搜索主方向选择物源方向5b ,得到该层的砂泥岩分布(图4).从模拟结果看出,它很好地满足了条件数据,即在各井点处的模拟结果与数据相一致,这表明该方法为条件模拟.同时,模拟的砂体展布方向和趋势与依靠地质经验手工绘制的砂体展布图(图5)比较吻合,在模拟的左下角与左上角砂体的展布与手工勾绘的几乎完全一致,但该方法在局部区域表现出砂体展布的非均质性和不确定性,与手工勾画砂体展布的平滑而唯一的表现是具有一定差别的,它充分体现了砂体局部的变异性和非均质性.图3布尔模拟生成的训练图像图4 多点地质统计方法模拟的砂泥岩展布图图5 手工绘制的该小层的砂体展布图3 结 论(1)多点地质统计学是今后地质统计学发展的主要方向,它可以联合反映空间多个位置点的几何形状和相互配位关系;在模拟具有复杂形状地质体分布时,它比两点地质统计学方法具有更大的优势.(2)利用snesim 模拟算法可以快速灵活地进行多点地质统计模拟,模拟的岩相展布图具有一定的真实性,它为储层参数的两阶段模拟奠定了基础.参考文献:[1] Strebelle S.Condit ional simulation of complex geologicalstructures using multiplepoint statist ics[J].M athematical G eolog y,2002,34(1):1-22.[2] Strebelle S,Journel A G.Reser voir modeling using mult-i ple -point stat i stics[C].SPE71324,2001.[3] Jef Caers.Stochast ic Reser voir Characterization U singM ultiple -po int statistics [A ].In:Proceedings o f the I AM G 99,F ifth Annual Conference of the I nternational Associat ion for M at hemat ical Geology [C].L ippars S G,et al ,1999.467-472.[4] Jef Caers.G eostatistical history matching under tr aining -image based geo logical model constraints [A ].SP E77429,2002.[5] 冯国庆,李允.应用序贯指示模拟方法模拟沉积微相[J].西南石油学院学报,2001,23(2):1-4.[6] Deutsch ,Journel.GSL IB:G eostatistical Softw ar e Libraryand U ser .s Guide,2nd ed[M ].New Yor k:Ox ford U n-iversity Press,1998.104-109.编辑:国伍玲)11)冯国庆等:利用多点地质统计学方法模拟岩相分布ABSTRACTS OF THE PRESENT ISSUE(JXSYU ISSN1673-064X)Determination of the geologic age of the bottom interface of Quaternary of Caidamu Basin Abstract:Up to now there are different viewpoints on the geolog ic age the bottom interface of Quaternary. T he determination of the geolog ic age of the bottom interface of Quaternary of Caidamu Basin is not only favor-able to the theoretical study of Quaternary but also of importance to the observation of biog as reservoirs.The ge-ologic age is determined as2.5M a B.P.based on the variation of sporopollen assemblages,the evolution of palaeoclimate,the unconformity formed by intense tectonic movement,the contact state of sediment and the v ariation of archaeom agnetism.The bottom interface is identical w ith that of Quaternary loess in the northwest of China determined by Liu Tun-sheng of our country and that of Quaternary g lacial age beginning determ ined by Shakleto N.T.based on the profile of DSDP site552A hole in North Atlantic Ocean.It is a unconformity surface in w hich the basin g raduates into conformity contact.The palaeoclimate becomes cold suddenly across it.A lot of strong first appears here.Key words:Caidam u Basin;fluvia-l lacustrine sedimentation;bottom interface of QuaternaryZ H U X iao-m in1,K AN G A n2,K AN G Qiang3(1.Key Laboratory of Education M inistry for Hydrocarbon Accumulation M echanism,China University of Petroleum,Beijing,102249,China; 2.China Offshore Petroleum Exploration Corporation,Beijing,100021,China;3.China Oil&Gas Exploration Development Cor-poration,Beijing,100011,China)JXSYU2005V.20N.5p.1-4Effects of volcanic activities on the distribution of accommodation space during the development of volcanic se-quence in Xinyang DepressionAbstract:The effects of the volcanic activities on the formation and the distribution of accommodation space among regressive system tracts of volcanic sequence are discussed w ith Eogene volcanic sequence in Xinyang De-pression being taken as an ex ample.The features of volcanic sequences and the patterns of the accommodation space in regressive system tracts show that,the local accommodation space of the fault-depression lake basin w ith volcanic cones is nonhomogeneous,w hich is different from common fault-depression lake basins.By studying se-quence stratigraphy in Xinyang Depression,it is realized that the volcanic cones in reg ressive system tracts lead to the nonhomogeneous distribution of the accom modation space in this depression,w hich enriches the content of terrigenous sequence stratigraphy.Key words:Xinyang Depression,volcanic activity;sequence;accommodation space;nonhomog eneity L I Ji-shan1,JIANG Zai-x ing2(1.Institute of Porous Flow&Fluid M echanics,Chinese Academy of Sc-i ences,Langfang065007,H ebei,China;2.Faculty of Earth Resource and Information,China U niversity of Petroleum,Dongying257061,Shandong,China)JXSYU2005V.20N.5p.5-8Stochastic simulation of lithofacies distribution using mult-i point geostatisticsAbstract:The research of reservoir heterogeneity is one of the main contents of reservoir description,the heterogeneity of lithology,the heterogeneity of the physical parameters of reservoirs,and the distribution,con-nectivity and vertical and horizontal variations of sandbodies can be revealed by the simulation of reservoir hetero-geneity,which is useful for deploy ing reasonable development plan and arranging w ell netw ork.But traditional g eostatistics describes the correlation and variability of geologic variables by m eans of variate difference function, and it can not fully reflect the correlation and variability of the geologic bodies with com plicated geometric shape since it is only based on two-point geostatistics.M ult-i point geostatistics is based on multiple-point correlation relationship.It has great advantages in reflecting the joint variability of the geolog ical variables.The theory of multiple-point geostatistics is introduced in this paper and a case is also g iven.T he lithofacies distribution of the sandstone reservoir in a oilfield in eastern China is simulated,and the result shows that using multiple-point geo-statistics can reveal the local heterogeneity and variability of sandbodies.Key words:multiple-point geostatistics;stochastic simulation;training im ageÔFEN G G uo-qing1,C H EN H ao2,Z H ANG L ie-hui1,LI Yun1(1.State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation,Southw est Petroleum Institute,Chengdu610500,Sichuan,China; 2.Lu-liang Oilfeild Operational Zone,Xinjiang Oilfield Company,Kelamayi834000,Xijiang,China)JXSYU2005 V.20N.5p.9-11Inorganic geochemical explanation of abnormal vitrinite reflectance:taking Niu-38well in Dongying Sag, Jiyang Depression as an exampleAbstract:The v itrinite reflectance of N iu-38well is restrained,it decreases from0.6%of the upper to 0.45%of the low er of cored interval.T he causes leading to the abnormality of the vitrinite reflectance are ana-lyzed from organic geochemistry and inorg anic geochemistry separately.The explanation from organic g eochem-istry is that,vitrinite reflectance decreases w ith the increase of the contents of sapropelic components,and hy-drogen-rich org anic matters restrain the increase of the vitrinite reflectance.From inorganic geochem istry,the study results in this paper show that,the information w ater in the studied interval becomes fresh and its acidity becomes a little strong from the low er to the upper.The acidic condition is favorable to the m aturation of organic matters.In the studied interval,the contents of K,Na,Mg and Li vary a little,therefore,there are no correla-tions between them and vitrinite reflectance;there are neg ative correlations betw een the contents of Ca,Mn,Sr and B and vitrinite reflectance,they maybe restrain the increase of vitrinite reflectance;there are positive corre-lations betw een the contents of Fe,Ba,Co,Zn,Ni and Rb and vitrinite reflectance,they maybe promote the increase of vitrinite reflectance.Key words:Dongy ing Sag;vitrinite reflectance abnormality;inorganic geochemistry;therm al evolution of organic matterXIE Ming-j u1,2,QI U N an-sheng1,2(1.Key Laboratory of Education Ministry for H ydrocarbon Accumula-tion M echanism,China U niversity of Petroleum,Beijing102249,China; 2.Basin&Reservoir Research Cen-ter,China University of Petroleum,Beijing102249,China)JXSYU2005V.20N.5p.12-16Application of nonlinearly constrained reservoir inversion technique in W eicheng Oilfield Abstract:Linear seismic inversion can only predict the distribution feature of thick sandbodies,but it can. t obtain the know ledge of the thin layer such as substratum to meet the demands of progressive exploration. Based on processing the seismic and log data of target,the distinct know ledg e of the structure of the reservoir can be gained by means of the nonlinear seismic inversion constrained by structure and log conditions,w hich can g uide the adjustment of oilfield deployment.T his study method is used in the progressive ex ploration of Weibei area in Weicheng Oilfield.The target processing of the seism ic and log data improves the resolution of the old data,and then the thin sandbodies are clearly recognized by nonlinearly constrained inversion,w hich are verified by drilling data.The fact demonstrates that nonlinearly constrained inversion has important function to the know ledge of reservoirs and the interpretation and verification of structures,and it is also an important tool to g uide progressive exploration.Key words:nonlinear constraint;reservoir inversion;w avelet frequency-division;logg ing curve reconstruc-tion;prog ressive exploration and developmentL I Zong-tian1,LI U Wei2(1.Department of Petroleum Engineering,China University of Geosciences, Beijing100081,China; 2.No.3Production Plant,Zhong yuan Oilfield Company,Puyang457001,Henan, China)JXSYU2005V.20N.5p.17-21Method of computing the mechanical parameters of formation and its application in Sulige area of Changqing OilfieldAbstract:A set of models for calculating three formation pressure profiles(form ation pore pressure,caving pressure and fracture pressure)and other formation mechanical parameters are established based on the study of borehole wall destabilization,logging data and ordinary empirical formulas.T he calculated m echanical parame-ters can help us design the safe density lim it of drilling fluid to make sure that the borehole w all is stable duringÕ。