CN106529022B - 一种储层分选性精确划分方法 - Google Patents
一种储层分选性精确划分方法 Download PDFInfo
- Publication number
- CN106529022B CN106529022B CN201610982828.0A CN201610982828A CN106529022B CN 106529022 B CN106529022 B CN 106529022B CN 201610982828 A CN201610982828 A CN 201610982828A CN 106529022 B CN106529022 B CN 106529022B
- Authority
- CN
- China
- Prior art keywords
- particle
- point
- pore throat
- sorting
- survey line
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
Landscapes
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明涉及油气田勘探开发、矿产评价预测领域,尤其是一种储层分选性精确划分方法。本发明专利利用数据化软件对颗粒或孔隙的边界数值化,得到的颗粒或孔喉边界数据化点,设置充填步长,通过测线与颗粒或孔喉边界充填数据点相交次数的统计以及测线数据点的充填,求取颗粒或孔喉的等效直径,实现颗粒或孔喉的分选性定量计算,进而提出了一套颗粒或孔隙分选性的编程化计算方法。本发明专利提供了一种储层分选性精确划分方法,它实现了粒度分选性和孔喉分选性精确划分,不仅对储层分选适用,同样适用于岩石颗粒磨圆度、球度的计算。
Description
技术领域
本发明涉及油气田勘探开发、矿产评价预测领域,尤其是一种储层分选性精确划分方法。
背景技术
分选性是指碎屑颗粒或者孔喉的粗细均匀程度,大小均匀者,分选性好;大小混杂者,分选性差。岩石颗粒、孔喉的描述包括颗粒形态、尺寸、磨圆性、分选性,是影响岩石的传输特性以及力学特性的重要因素,在这其中储层的分选性是影响储层物性的重要因素,能够反映储层均质程度,同时也是水动力强弱的标志。储层的分选性分为粒度分选和孔喉分选,一般都是通过实验室分析来得到,前者是通过粒度分析的累积曲线计算得到,后者是通过压汞法的毛管压力曲线计算得到,主要图件有孔喉分布直方图、频率曲线及毛管压力曲线。这2种方法都是通过累积曲线来计算分选系数。对分选系数的计算,二者在形态上极具相似性,通常也都是采用福克和沃德的标准偏差公式来计算。在实际储层中,天然的颗粒以及孔喉大多是非规则形状的,很难准确定量表述其分选特征。
发明内容
本发明旨在解决上述问题,提供了一种储层分选性精确划分方法,它实现了粒度分选性和孔喉分选性精确划分。
本发明的技术方案为:一种储层分选性精确划分方法,具体步骤如下:
第一步利用数据化软件对颗粒或孔隙的边界数值化
边界数据化方法主要有两种:①利用MATLAB软件将图像去噪声干扰,使图片能够映储层孔喉、颗粒形状、大小,然后使用锐化函数对低通滤波后的图像进行锐化处理,以突出颗粒或孔隙的边界;通过颗粒或孔隙边缘检测,突出完整的颗粒或孔隙边界形态,最后将边界数值化;②直接利用数据化软件对孔喉或者颗粒的边界数据化。
第二步得到的颗粒或孔喉边界数据化点,设置充填步长b,即充填后的边界数据点的间距为b,颗粒的最小内接圆半径为r,时能够满足分选性的计算要求,充填后得到颗粒或孔喉边界数据化点集合Q。
第三步测线与颗粒或孔喉边界充填数据点的相交次数统计
在颗粒或孔喉分布的矩形区域内,所述的测线平行于矩形边界,所述的测线间距为d,测线条数为t,时,满足分选性计算的要求。统计每条测线与边界数据化点集合Q相交的交点个数ni,其中i为第i条测线。
对第i条测线插值得到间距为d的数据集合Gi,其中1≤i≤t,对于集合Gi中的任意一点P(Pxi,Pyi),当同时有集合Q中的两个点E(Exj,Eyj)、F(Fxk,Fyk)满足公式(1)条件时,认为P点是第i条测线与点集合Q的一个交点:
由公式(1)得到每条测线中交点的个数ni,交点坐标位置集合记为Ui(Uxm,Uym),其中m≤ni。
第四步测线数据点的充填
对于第i条测线的交点集合Ui(Uxm,Uym),交点的个数为ni,对交点Uw——Uw+1之间充填间距为d的数据点,得到该测线的数据充填点集Si。w满足条件:
如图3所示,lc、le为测线,对于测线lc,1、2交点满足公式(2)要求,对交点1、2充填间距为d的数据点,对于测线le,1、2交点以及3、4交点满足公式(2)要求,对交点1、2以及交点3、4之间充填间距为d的数据点。通过循环迭代,得到全部测线充填点集S(图4),其个数为node:
t为测线条数,Si为第i条测线的充填点集。
第五步颗粒或孔喉的等效直径ca
由测线充填点集S得到每个颗粒或孔喉的等效直径ca:
第六步颗粒或孔喉的分选性定量计算
采用福克和沃德所提出的标准偏差公式:
公式(5)-公式(6)中,Φ5、Φ16、Φ84、Φ95分别为颗粒累计为5%、16%、84%、95%对应的直径,对于粒度分选性计算,ca为颗粒等效直径;对于孔喉分选性计算,ca为孔喉等效直径。偏差σ越小,分选性越好;具体定性、定量分选性划分标准如表1所示。
表1颗粒或孔喉的分选性-定性、定量分选性划分标准表
标准偏差 | 分选等级 |
<0.35 | 分选极好 |
0.35~0.50 | 分选好 |
0.50~0.71 | 分选较好 |
0.71~1.00 | 分选中等 |
1.00~2.00 | 分选差 |
2.00~4.00 | 分选较差 |
>4.00 | 分选极差 |
本发明的有益效果是:本发明专利利用数据化软件对颗粒或孔喉的边界数值化,得到的颗粒或孔喉边界数据化点,设置充填步长,通过测线与颗粒或孔喉边界充填数据点的相交次数的统计以及测线数据点的充填,求取颗粒或孔喉的等效直径,实现颗粒或孔喉的分选性定量计算,进而提出了一套颗粒或孔喉分选性的编程化计算方法。本发明对于颗粒或孔喉的分选具有较高的实用价值,并且预测成本低廉、可操作性强,能大量减少人力、财力的支出,在油田实用性强;本发明专利不仅对储层分选适用,同样适用于岩石颗粒磨圆度、球度的计算。
附图说明
图1为一种储层分选性精确划分方法的流程图。
图2为颗粒或孔喉边界数值化示意图。
图3为测线与颗粒或孔喉边界充填数据点的相交次数的统计示意图。
图4为颗粒或孔喉面积求取示意图。
图3中,lc、le为测线,1、2、3、4分别为测线与颗粒或孔喉边界充填数据点。
具体实施方式
下面结合附图说明本发明的具体实施方式:
1)利用数据化软件对颗粒或孔喉的边界数值化,利用数据化软件对孔喉或者颗粒的边界数据化。
2)将得到的颗粒或孔喉边界数据化点(图2),设置充填步长b=10m,充填后得到颗粒或孔喉边界数据化点集合Q。
3)在颗粒或孔喉分布的矩形区域内,所设置d=100,对第i条测线,插值得到间距为d的数据集合Gi,对于集合Gi中的任意一点P(Pxi,Pyi),当同时有集合Q中的两个点满足公式(1)时,认为该点为测线与颗粒或孔喉边界充填数据点的交点,得到每条测线中交点的个数ni,交点坐标位置集合记为Ui(Uxm,Uym),其中m≤ni。
4)对于第i条测线的交点集合Ui(Uxm,Uym),交点的个数为ni,对交点Uw——Uw+1之间充填间距为d的数据点,得到该测线的数据充填点集Si。w满足公式(2)时,通过循环迭代,得到全部测线充填点集S(图4)。
5)通过公式(4)得到颗粒或孔喉的等效直径ca。
6)利用公式(5)-(6)实现颗粒或孔喉的分选性定量计算。
上面以举例方式对本发明进行了说明,但本发明不限于上述具体实施例,凡基于本发明所做的任何改动或变型均属于本发明要求保护的范围。
Claims (2)
1.一种储层分选性精确划分方法,计算的步骤如下:
1)利用MATLAB软件将图像去噪声干扰,使图片能够映储层孔喉、颗粒形状、大小,然后使用锐化函数对低通滤波后的图像进行锐化处理,以突出孔喉、颗粒的边界,通过孔喉、颗粒边缘检测,突出完整的孔喉、颗粒边界形态,最后将边界数值化,或者直接利用数据化软件对孔喉或者颗粒的边界数据化;
2)利用得到的颗粒或孔喉边界数据化点,设置充填步长b,所述的b满足其中r为颗粒的最小内接圆半径,充填后得到颗粒或孔喉边界数据化点集合记为Q;
3)测线与颗粒或孔喉边界充填数据点的相交次数统计,所述的测线平行于矩形统计单元边界,所述的测线的间距为d,测线条数为t;在颗粒或孔喉分布的矩形区域内,d满足时,即符合分选性计算的要求;对第i条测线插值得到间距为d的数据集合Gi,其中,1≤i≤t,对于集合Gi中的任意一点P(Pxi,Pyi),当同时有集合Q中的两个点E(Exj,Eyj)、F(Fxk,Fyk)满足公式(1)条件时,认为P点是第i条测线与点集合Q的一个交点:
由公式(1)得到每条测线中交点的个数ni,交点坐标位置集合记为Ui(Uxm,Uym),其中m≤ni;
4)测线与颗粒或孔喉边界充填数据点的交点之间进行点充填,得到全部测线充填点集S,其个数为node;
5)颗粒或孔喉的等效直径ca计算;由测线充填点集S得到每个颗粒或孔喉的等效直径ca的计算公式:
6)颗粒或孔喉的分选性定量计算;采用福克和沃德所提出的标准偏差公式:
公式(5)-公式(6)中,Φ5、Φ16、Φ84、Φ95分别为颗粒累计为5%、16%、84%、95%对应的直径,对于粒度分选性计算,ca为颗粒等效直径;对于孔喉分选性计算,ca为孔喉等效直径;偏差σ越小,分选性越好。
2.根据权利要求1所述的一种储层分选性精确划分方法,所述的测线与颗粒或孔喉边界充填数据点的交点间的点充填方法为:对于第i条测线的交点集合Ui(Uxm,Uym),交点的个数为ni,对交点Uw——Uw+1之间充填间距为d的数据点,得到该测线的数据充填点集Si;w满足条件:
通过循环迭代,得到全部测线充填点集S,其个数为node:
t为测线条数,Si为第i条测线的充填点集。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610982828.0A CN106529022B (zh) | 2016-11-09 | 2016-11-09 | 一种储层分选性精确划分方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610982828.0A CN106529022B (zh) | 2016-11-09 | 2016-11-09 | 一种储层分选性精确划分方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106529022A CN106529022A (zh) | 2017-03-22 |
CN106529022B true CN106529022B (zh) | 2019-09-06 |
Family
ID=58350407
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610982828.0A Expired - Fee Related CN106529022B (zh) | 2016-11-09 | 2016-11-09 | 一种储层分选性精确划分方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106529022B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3693794A (en) * | 1968-04-17 | 1972-09-26 | Buehler Ag Geb | Method and apparatus for separating granular goods |
CN101738402A (zh) * | 2008-11-12 | 2010-06-16 | 中国石油天然气股份有限公司 | 一种岩样二维ct图像分析系统 |
CN103778328A (zh) * | 2014-01-10 | 2014-05-07 | 中国石油大学(华东) | 基于标准偏差分析的储层敏感孔喉提取方法 |
CN106021788A (zh) * | 2016-05-31 | 2016-10-12 | 中国石油大学(华东) | 一种基于孔喉结构特征的致密储层分级评价标准划分方法 |
-
2016
- 2016-11-09 CN CN201610982828.0A patent/CN106529022B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3693794A (en) * | 1968-04-17 | 1972-09-26 | Buehler Ag Geb | Method and apparatus for separating granular goods |
CN101738402A (zh) * | 2008-11-12 | 2010-06-16 | 中国石油天然气股份有限公司 | 一种岩样二维ct图像分析系统 |
CN103778328A (zh) * | 2014-01-10 | 2014-05-07 | 中国石油大学(华东) | 基于标准偏差分析的储层敏感孔喉提取方法 |
CN106021788A (zh) * | 2016-05-31 | 2016-10-12 | 中国石油大学(华东) | 一种基于孔喉结构特征的致密储层分级评价标准划分方法 |
Non-Patent Citations (2)
Title |
---|
Modeling of Diffusion with Partitioning in Stratum Corneum Using a Finite Element Model;ANA M.BARBERO等;《Annals of Biomedical Engineering》;20050930;第33卷(第9期);第718-722页 |
基于VBA的储层评价数据处理方法;王高成 等;《物探与化探》;20091231;第33卷(第6期);第1281-1292页 |
Also Published As
Publication number | Publication date |
---|---|
CN106529022A (zh) | 2017-03-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110458119B (zh) | 一种非接触式测量的混凝土骨料级配快速识别方法 | |
CN104850682B (zh) | 基于位置的多点地质统计学建模方法 | |
JP2020502673A5 (zh) | ||
Leonard et al. | The challenges of probabilistic seismic‐hazard assessment in stable continental interiors: An Australian example | |
CN104483705B (zh) | 一种三维剩余静校正方法 | |
EP4242870A3 (en) | Interactive geographical map | |
CN104573705A (zh) | 一种建筑物激光扫描点云数据的聚类方法 | |
EP2296117A3 (en) | Ray-aggregation for ray-tracing during rendering of imagery | |
Zhang et al. | Development of a supervised software tool for automated determination of optimal segmentation parameters for ecognition | |
CN109712131A (zh) | 肺结节特征的量化方法、装置、电子设备及存储介质 | |
CN105224937A (zh) | 基于人体部件位置约束的细粒度语义色彩行人重识别方法 | |
CN104573333B (zh) | 利用聚类分析进行模型优选的方法 | |
CN108182440A (zh) | 一种基于渣片图像识别获取围岩类别的方法 | |
CN102819688A (zh) | 基于半监督分类的二维地震数据全层位追踪方法 | |
CN102073867B (zh) | 一种遥感图像分类方法及装置 | |
CN105354600B (zh) | 一种砂岩显微薄片的自动分类方法 | |
CN102938069A (zh) | 一种基于信息熵的纯净和混合像元自动划分方法 | |
CN102867183A (zh) | 一种车辆遗撒物检测方法、装置及智能交通监控系统 | |
CN110672072A (zh) | 一种用于裂隙岩体结构的信息多元解译与表征方法 | |
Kosarevych et al. | Image segmentation based on the evaluation of the tendency of image elements to form clusters with the help of point field characteristics | |
Ni et al. | Lineament length and density analyses based on the segment tracing algorithm: a case study of the gaosong field in gejiu tin mine, China | |
CN110501751B (zh) | 一种基于多分量梯度数据联合和深度加权的相关成像方法 | |
Han et al. | A trace map comparison algorithm for the discrete fracture network models of rock masses | |
US9299191B2 (en) | Adaptive artifact removal | |
CN106529022B (zh) | 一种储层分选性精确划分方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
TA01 | Transfer of patent application right | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20190807 Address after: 266580 Qingdao Changjiang Road, Huangdao District, Shandong, No. 66 Applicant after: China Petroleum University (East China) Address before: 100083 Xueyuan Road, Beijing, China University of Geosciences, (Beijing), 29 Applicant before: China University of Geosciences (Beijing) |
|
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190906 Termination date: 20201109 |