地质储量计算区块面积确定方法研究.pdf

第27卷 第2期 物探化探计算技术 2005年5月 收稿日期2004- 09- 24 文章编号10011749200502014703 地质储量计算区块面积确定方法研究 刘吉余,杨玉华,于润涛 大庆石油学院,黑龙江 大庆 163318 摘 要根据地质储量计算区块边界特征来计算含油面积,其核心技术是充分利用现有的石油 勘探、 开发数据库资料,根据储量计算区块边界各顶点坐标,由计算机来计算储量计算区块的面 积。该方法解决了用求积仪等常规方法求取面积时所存在的精度偏低、 操作繁琐、 工作量大等缺 点,用这种方法在面积确定过程中自动化程度高,具有方便、 灵活、 实用的特点。目前此方法已在 生产中进行了推广应用。 关键词地质储量;容积法;求积仪;边界;面积 中图分类号TE155 文献标识码 A 0 引言 石油、 天然气地质储量是石油和天然气在地下 的蕴藏量,它是油、 气田勘探开发成果的综合反映。 无论是勘探阶段还是开发阶段,油气储量一直是石 油工作者追寻的主要目标,是油气田勘探、 开发过 程中的一项极为重要的工作任务。 油气储量是指导 油气田勘探、 开发,确定投资规模的重要依据。因 此,石油、 天然气地质储量计算的准确与否至关重 要。 在油田开发过程中,常常需要用容积法计算不 同区块、 不同井组、 不同沉积单元的地质储量,为采 取各种开发措施提供地质储量上的依据。 容积法地 质储量计算参数主要包括含油面积、 有效厚度、 有 效孔隙度、 含油饱和度、 地面原油密度、 原油体积系 数等[1, 2]。 用容积法计算储量的可靠程度,取决于上 述影响因素的数量和质量。 随着勘探开发工作的不 断深入,部分油田储量与产量之间的不匹配现象日 益突出。对1984年以后申报的56个油田94个断 块地质储量复算容积法结果统计显示[3],复算 前、 后地质储量变化很大,不同油田的误差范围分 布主要集中在6～72之间。经分析表明,含油 面积对地质储量变化的影响占35 ,产生误差的 原因包括① 基础资料不全、 不准;② 由于研究 工作不够;③ 确定参数方法不合理。 对于某一具体储量计算区块,含油面积是储量 计算的首要参数,对其测定精度的高低,直接影响 到储量计算结果的可靠性。但在实际测定过程中, 一些常见方法往往存在精度偏低、 操作繁琐等缺 点。为此,我们对储量计算区块的含油面积确定方 法进行了研究,开发了一套计算储量区块面积的方 法,该方法具有计算速度快、 操作简单方便、 精度高 等优点,值得进一步推广应用。 1 传统确定地质储量计算区块面积 方法简述 传统确定储量计算区块面积的方法有在储量 计算区块井位图上,用求积仪法、 几何图形法、 厘米 方格纸法、 曲线仪法及称重法来确定含油面积,其 中求积仪法是最常用的一种方法。 在实际应用过程 中发现,上述方法均在某种程度上存在精度偏低、 计算及操作繁琐等缺点,而且往往需要用另一种方 法进行检查验证。例如,生产上常用求积仪法来确 定地质储量计算区块的面积,使用时要求用求积仪 在两个不同位置上各测一次,然后取其算术平均 值,并要求每次测定结果与算术平均值之差不大于 其值的1 200 或两次测定结果误差小于 2 。另 外,此法要求的工作条件很高如桌面图纸绝对的 平整 , 操作特别仔细,否则每次测定结果与算术平 均值之差就大于允许误差。 用求积仪法确定面积后 还需要用几何图形法或别的方法进行检查验证,操 作过程就更加繁琐,工作量巨大。 2 储量计算区块边界特征 储量计算区块边界特征不外乎有两种,一是以 相带为边界;另一种是以油水井连线为边界,油 水井井距之半连线为边界,或任意点连线为边 界。以相带为边界的线是圆滑的弧线,而以油水 井连线为边界、 油水井井距之半连线为边界、 或 任意点连线为边界的线是折线。 储量计算区块面积 的确定实际上就是确定圆滑弧线构成的封闭面积 或折线构成的多边形面积。 3 储量计算区块面积的计算原理 确定圆滑弧线构成的封闭面积可采用文献[4] 所述的方法,而多边形的面积计算有多种算法,根 据储量计算区块面积特征,我们采用多边形各顶点 坐标的算法来计算其面积。 设有按一定顺序顺时针或逆时针排列的n 个井点x1,y1、x2,y2、⋯⋯、xn,yn围成一个 没有交叉的多边形见图1,则该多边形的面积为 图1 储量计算区块顶点坐标五边形 计算面积计算示意图 Fig.1 The sketch map of calculating coordinatesof the acmes pentagonof the block boundary according to reserves and area S 1 2 6 n i 1 xiyi 1-xi 1yi1 式中 xn 1x1;yn 1y1 多边形各顶点坐标实际上就是井位坐标油 水井连线为边界 , 两口井的井距之半处坐标油 水井井距之半连线为边界 , 或任意点处坐标任 意点连线为边界。 井位坐标是已知的,两口井井距之半处坐标和 任意点处坐标可通过已知的井位坐标计算,并通过 实际坐标系统与屏幕坐标系统进行相互转换。 主要 代码如下 int m_nV iewW idth; 视图区宽度 int m_nV iewHeight; 视图区高度 int m_nCurV iewO rgX; 当前可视范围 的原点X int m_nCurV iewO rgY; 当前可视范围 的原点Y int m_nCurV iewExtX; 当前可视范围 终止X int m_nCurV iewExtY; 当前可视范围 终止Y 初始化图纸 void CTri A reaDoc∷IniPaperCRect rect { m_nV iewW idth rect. W idth ; m_nV iewHeight rect. Height ; } 实际坐标转换到视图坐标 void CTri A reaDoc∷WorldToV iewCPoint point. y m_nV iewHeight3point. y - m_nCurV iewO rgYm_nCurV iewExtY -m_ nCurV iewO rgY ; point. y m_nV iewHeight- point. y; } 视图坐标转换到实际坐标 void CTri A reaDoc∷V iewToWorldCPoint point. y m_nV iewHeight2point. y; point. y m_nCurV iewO rgY m_ nCurV iewExtY- m_nCurV iewO rgY 3 point. ym_ nV iewHeight; } 841物探化探计算技术27卷 4 地质储量计算区块面积的计算 软件 根据地质储量计算区块边界特征及任意面积 的计算原理,我们研制了地质储量计算区块面积的 计算系统。该系统是在W indow s 2000环境下,采 用M icrosoft V isualC 6. 0开发,系统流程及界面 如图所示见图2、 图 3 。 5 应用结果及效果分析 储量计算区块面积的计算系统开发完成以后, 在大庆某油田北一二排西S˚10- S10组 “上返 方案” 储量的精细计算中进行了应用。 北一二排西S˚10- S10组在垂向上共分 20个沉积单元,平面上划分了86个不规则小区 块,分不同厚度等级、 不同沉积单元计算地质储量。 同时,在北一二排西S˚10- S10组 “上返方案” 设计87个井组储量计算中进行了应用。 结果表明, 该方法具有精度高、 操作简单、 方便等特点,提高工 作效率达30倍以上,完全可以取代传统的面积计 算方法求积仪法、 几何图形法、 厘米方格纸法、 曲 线仪法、 称重法等。 目前该方法已在生产中进行了 推广应用。 图2 地质储量计算区块面积确定流程图 Fig. 2 The flow chart of geologic reserve calculation and block area definition 图3 地质储量计算区块面积计算系统常用工具及某一区块计算结果 Fig. 3 The common tools of geologic reserves calculation and block area system and calculating result of a certain block 参考文献 [1 ] 黎文清,李世安.油气田开发地质基础第二版 [M ].北京石油工业出版社, 1993. [2] 吴元燕,陈碧珏.油矿地质学第二版 [M ].北京石 油工业出版社, 1996. [3] 萧德铭,毕海滨.利用油基泥浆取心资料确定砂岩油 藏原始含油饱和度若干问题探讨[J ].中国海上油气 地质, 2003, 174 252. [4] 刘吉余,许洪东,王长生,等.任意面积储量计算方法 研究[J ].物探化探计算技术, 2003, 251 75. 作者简介刘吉余1961- ,男,博士,副教授,主 要从事石油开发、 地质、 油气资源定量化方面的研 究。 9412期刘吉余等地质储量计算区块面积确定方法研究 EX PLORA T ION, 2005, 272 0141 By analyzing the accessory element such as M n, Sr, M g, etc, and isotope such as 13 C, 18 O, 87Sr86Sr etc, the authors consider that themost i m 2 portant mechanism of the carbonate reservoir of Lower Paleozoic are dolom itization and karstifica2 tion.The reservoirs of Yeli- L iang jianshan a2 tion are valuable in ZhuangHai zone and theirs for2 mation was mainly concerned w ith dolom itizatin during the buried processing, especially concerned w ith the volume contraction of grains,and they mainly distributed in diagenetic dolom ites .Thema2 rine carbonate near the uncon itywasmainly re2 constructed by strong atmogenic freshwater during the Indo -China and Yanshan M ovement, being better reservoir .The reconstructed records behave as follow ing the ratio for 87Sr86Sr in the different carbonatem ineralsor structural components is high2 er than carbonate matrix, while the∆ 13C and ∆ 18O value in mass of samples in the different carbonate m inerals or structural components are lower than carbonatematrix.The geochem ical section of single well show s that the depth influenced by atmogenic freshwater is about 150 meter which is under the Caledonian affiliation plane of the uncon ity and at the same ti me, the freshwater reconstructed the dolom ite reservoirwhich had existed before the for2 mer. Key words karstification; dolom itization; ing mechanism; carbonate; zhuanghai zone; jiyang sag RESEARCH OF THE CALCULATINGM ETHOD OFGEOLOGI CRESERVES ANDDEFINING BLOCK AREA L iu Ji2yu, YAN G Yu2hua, YU Run2taoDaqing PetroleumInstitute, Daqing Heilongjiang 163318, China1COM PU T IN GT ECH N IQU ESFOR GEOPH YS ICALAN DGEOCH EM ICALEX2 PLORA T ION, 2005, 272 0147 The paper mainly discusses a of calcu2 lating oil-bearing area based on block boundary characters calculated by geological reserves .The key technique of the is to make full use of the ination from the present oilprospecting and exploitation database and to calculate the coordinates of the acmes of the block boundary according to re2 serves and block area and reservesw ith computer . The resolved the defects caused by common s w ith plani meter on calculating area, which produce low precision and require com2 plicatedoperationsandmuchworkload.The w ith high automation degree in the course of defining area has been w idely used and spread in production at present for its convenience, flexibility and practice. Key words reserves; volumetric ; plani me2 ter; boundary; area THE ACHIEVE M ENT OF THE L′TH I M PEND OF FRACTALINTERPOLATED SURFACE OF CHE M ISTRY ELE M ENT IN ORE BELT USING MATLAB N IXiu2jing, L IJian, HAN Ze2huaSchool of in2 ation adm inistration,Chengdu U niversity of Technology,Chengdu 610059, China1COM2 PU T IN G T ECH N IQU ES FOR GEOPH YS ICAL AN D GEOCH EM ICALEX PLORA T ION, 2005, 272 0150 M any problem s of geosciences are very compli2 cated,non2linearity and irregular and they may have better solution if non2linearity sciences are used for them.This paper introduces the interpolation function for fractal interpolated surface and discuss2 es the realization about chem istry element′sL th i m2 pend of fractal interpolated surface and the algo2 rithm sw ithMA TLAB. Based on the algorithm and the MA TLAB procedure,the fractal interpolated surfaces are created according to chem istry element data in practical ore belt2ore area of Gegonglong of Gongjue county of XizangTibetA utonomousRe2 gion.From the procedure of achievement, we can fully realize the convenience and high2efficiency of MA TLAB in solving practicalproblem s . Key wordsfractal geometry; fractal interpolated surface; MA TLAB A M ULTI-MAP SYSTE M I C ERROR CORREC- TI ONM ETHOD BASED ON BORDERL I M ITING PRINCIPLE J I Hong2jin,SH I Yan2xiang,DA I Yong2gang, HAO L i2bo, LU Ji2longCollege of Geo2Explo2 ration Science and Technology,Jilin U niversity, Changchun 130026, China . COM PU T IN G T E2 CH N IQU ES FOR GEOPH YS ICAL AN D GEO2 CH EM ICALEX PLORA T ION, 2005, 27 2 0154 To solve the problem about multi- map sys2 tem ic error correction in regional geochem ical explo2 ration,a computing based on the least squaress is developed.The basic theorem is the border li m iting principle developed in partition standardization .It alsoshow s thatthe is a extension of partition standardization , and its correction result is the same as lat2 ter in the case of two maps . Key words geochem icalmap; system ic error; cor2 rection; border li m iting principle 4ABSTRACTS Vol . 27 No. 2