CN104992067B - 一种基于裂缝产状的裂缝分析方法 - Google Patents
一种基于裂缝产状的裂缝分析方法 Download PDFInfo
- Publication number
- CN104992067B CN104992067B CN201510413045.6A CN201510413045A CN104992067B CN 104992067 B CN104992067 B CN 104992067B CN 201510413045 A CN201510413045 A CN 201510413045A CN 104992067 B CN104992067 B CN 104992067B
- Authority
- CN
- China
- Prior art keywords
- crack
- msub
- mrow
- fracture
- inclination angle
- 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.)
- Active
Links
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明涉及一种基于裂缝产状的裂缝分析方法,其包括以下步骤:加载裂缝的深度、倾角和方位角信息,根据每条裂缝的产状,将井中的裂缝显示在SMIT图中;根据裂缝的倾角和方位角的范围,在SMIT图中对裂缝进行分组,分成低角度裂缝和不同走向的高角度裂缝;分组统计不同组裂缝的倾角和方位角信息;利用选择的裂缝计算不同组裂缝的线密度、面密度和体密度。本发明的有益效果为:能够在分组的基础上对裂缝进行分组分析,计算不同组不同产状裂缝的裂缝密度,从而提高裂缝分析的精度,使裂缝分析更好地服务于生产,提高裂缝油藏的研究精度,提高裂缝油藏的勘探开发精度,帮助提高裂缝油气的开发效率。
Description
技术领域
本发明属于石油勘探开发领域,具体涉及一种基于裂缝产状的裂缝分析方法。
背景技术
裂缝对油气藏勘探开发至关重要,裂缝分析是裂缝研究的重要技术环节。不同产状的裂缝其开度和传递率不同,其对油气藏渗流的贡献也不同,因此裂缝研究和裂缝建模需要按照产状对裂缝进行分组,在分组的基础上计算不同组裂缝的裂缝密度。然而现有技术只能对裂缝进行整体分析,计算总的裂缝密度,不能进行分组分析。
发明内容
为了解决现有技术存在的上述问题,本发明提供了一种基于裂缝产状的裂缝分析方法,该方法能够对裂缝进行分组分析。
本发明所采用的技术方案为:一种基于裂缝产状的裂缝分析方法,其包括以下步骤:1)加载裂缝的深度、倾角和方位角信息,根据每条裂缝的产状,将井中的裂缝显示在SMIT图中;2)根据裂缝的倾角和方位角的范围,在SMIT图中对裂缝进行分组,分成低角度裂缝和不同走向的高角度裂缝;3)分组统计不同组裂缝的倾角和方位角信息;4)利用选择的裂缝计算不同组裂缝的线密度、面密度和体密度。
所述步骤2)中在SMIT图中对裂缝进行分组,其具体包括:①确定裂缝的倾角和方位角的范围;手动输入裂缝的倾角和方位角的范围或者通过鼠标事件产生的扇形区来确定裂缝的倾角和方位角的范围;②利用倾角和方位角的范围对裂缝进行分组。
所述步骤4)中,不同组裂缝的线密度、面密度和体密度分别为:线密度:
式中,FLine表示线密度;m1表示裂缝步长,其取值范围为1~10米;n1表示m1步长内的裂缝条数;面密度:
式中,m2表示井段的厚度,其取值范围为1-10米,n2表示厚度为m2的井段的裂缝条数,Si表示每条裂缝的面积,S=π×d×arcsin(α),d表示井筒的半径,α表示裂缝倾角;体密度:
式中,vi表示每条裂缝的体积,v=π×d2×arcsin(α)×A,A表示裂缝开度
由于采用以上技术方案,本发明的有益效果为:1、本发明能够在裂缝SMIT图上按照产状对裂缝进行分组,在分组的基础上对裂缝进行分组分析,计算不同组不同产状裂缝的裂缝密度,从而提高裂缝分析的精度,使裂缝分析更好地服务于生产,提高裂缝油藏的研究精度,提高裂缝油藏的勘探开发精度。2、裂缝分析研究能提高裂缝油气的开发效率,帮助裂缝油气藏实现稳油控水,能够为石油公司和服务公司带来巨大的经济效益和社会效益。
附图说明
图1是本发明基于裂缝产状的裂缝分析方法的流程图;
图2是井中裂缝在SMIT图中的显示效果图;
图3是裂缝在SMIT图中的分组效果图。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
如图1所示,本发明提供了一种基于裂缝产状的裂缝分析方法,其具体包括以下步骤:
5)加载裂缝的深度、倾角和方位角等数据信息,根据每条裂缝的产状,将井中的裂缝显示在如图2所示的SMIT图中。
SMIT图中正北N处为0°,顺时针一周表示0°~360°。SMIT图的中心点倾角为0°,其最外沿倾角为90°,将每条裂缝以小三角形的形式显示在图2中。
2)根据裂缝的倾角和方位角的范围,在SMIT图中对裂缝进行分组,分成低角度裂缝和不同走向的高角度裂缝,其具体过程为:
①确定裂缝的倾角和方位角的范围;
裂缝的倾角和方位角的范围可以通过手动输入,也可以通过鼠标拖拽产生的扇形区来确定。
②利用倾角和方位角的范围对裂缝进行分组。
如图3所示的低角度裂缝,其倾角范围为0°~30°,其方位角范围为30°~330°。
3)分组统计不同组裂缝的倾角和方位角信息。
4)利用选择的裂缝计算不同组裂缝的线密度、面密度和体密度,其具体过程为:
线密度:
式中,FLine表示线密度。m1表示裂缝步长,其取值范围为1~10米;如果裂缝密度大,m1值取1~5米,如果裂缝密度小,m1值取5~10米。n1表示m1步长内的裂缝条数。
面密度:
首先,计算得到单条裂缝的面积:
S=π×d×arcsin(α) (2)
式中,d表示井筒的半径,利用井径曲线可以得到d值。α表示裂缝倾角。
其次,计算得到面密度:
式中,m2表示井段的厚度,其取值范围为1-10米,n2表示厚度为m2的井段的裂缝条数,Si表示每条裂缝的面积。
体密度:
首先,计算单条裂缝的体积:
v=π×d2×arcsin(α)×A (4)
式中,A表示裂缝开度。
其次,计算得到体密度:
式中,vi表示每条裂缝的体积。
本发明不局限于上述最佳实施方式,任何人在本发明的启示下都可得出其他各种形式的产品,但不论在其形状或结构上作任何变化,凡是具有与本申请相同或相近似的技术方案,均落在本发明的保护范围之内。
Claims (1)
1.一种基于裂缝产状的裂缝分析方法,其包括以下步骤:
1)加载裂缝的深度、倾角和方位角信息,根据每条裂缝的产状,将井中的裂缝显示在SMIT图中;
2)根据裂缝的倾角和方位角的范围,在SMIT图中对裂缝进行分组,分成低角度裂缝和不同走向的高角度裂缝;
在SMIT图中对裂缝进行分组,其具体包括:
①确定裂缝的倾角和方位角的范围;
手动输入裂缝的倾角和方位角的范围或者通过鼠标事件产生的扇形区来确定裂缝的倾角和方位角的范围;
②利用倾角和方位角的范围对裂缝进行分组;
3)分组统计不同组裂缝的倾角和方位角信息;
4)利用选择的裂缝计算不同组裂缝的线密度、面密度和体密度;
不同组裂缝的线密度、面密度和体密度分别为:
线密度:
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>L</mi>
<mi>i</mi>
<mi>n</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>n</mi>
<mn>1</mn>
</msub>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>,</mo>
</mrow>
式中,FLine表示线密度;m1表示裂缝步长,其取值范围为1~10米;n1表示m1步长内的裂缝条数;
面密度:
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>a</mi>
<mi>r</mi>
<mi>e</mi>
<mi>a</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mn>1</mn>
<msub>
<mi>n</mi>
<mn>2</mn>
</msub>
</munderover>
<msub>
<mi>S</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<mi>&pi;</mi>
<mo>&times;</mo>
<msub>
<mi>m</mi>
<mn>2</mn>
</msub>
<mo>&times;</mo>
<msup>
<mi>d</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
式中,m2表示井段的厚度,其取值范围为1-10米,n2表示厚度为m2的井段的裂缝条数,Si表示每条裂缝的面积,S=π×d×arcsin(α),d表示井筒的半径,α表示裂缝倾角;
体密度:
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>v</mi>
<mi>o</mi>
<mi>l</mi>
<mi>u</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mn>1</mn>
<msub>
<mi>n</mi>
<mn>2</mn>
</msub>
</munderover>
<msub>
<mi>v</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<mi>&pi;</mi>
<mo>&times;</mo>
<msub>
<mi>m</mi>
<mn>2</mn>
</msub>
<mo>&times;</mo>
<msup>
<mi>d</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
式中,vi表示每条裂缝的体积,v=π×d2×arcsin(α)×A,A表示裂缝开度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510413045.6A CN104992067B (zh) | 2015-07-14 | 2015-07-14 | 一种基于裂缝产状的裂缝分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510413045.6A CN104992067B (zh) | 2015-07-14 | 2015-07-14 | 一种基于裂缝产状的裂缝分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104992067A CN104992067A (zh) | 2015-10-21 |
CN104992067B true CN104992067B (zh) | 2018-02-09 |
Family
ID=54303880
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510413045.6A Active CN104992067B (zh) | 2015-07-14 | 2015-07-14 | 一种基于裂缝产状的裂缝分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104992067B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106291704B (zh) * | 2016-09-20 | 2018-02-27 | 中国地质大学(北京) | 一种不同尺度裂缝面密度预测方法 |
CN110927819B (zh) * | 2018-09-20 | 2021-08-31 | 中国石油化工股份有限公司 | 一种裂缝发育程度表征方法 |
CN114384584B (zh) * | 2020-10-20 | 2024-05-28 | 中国石油天然气股份有限公司 | 裂缝建模方法及装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103278864A (zh) * | 2013-05-10 | 2013-09-04 | 中国石油天然气股份有限公司 | 洞缝型储层的地质特征参数及分布的测定方法及装置 |
CN103278436A (zh) * | 2013-02-01 | 2013-09-04 | 西安石油大学 | 特低渗透双重介质砂岩油藏微观孔隙结构的定量表征方法 |
-
2015
- 2015-07-14 CN CN201510413045.6A patent/CN104992067B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103278436A (zh) * | 2013-02-01 | 2013-09-04 | 西安石油大学 | 特低渗透双重介质砂岩油藏微观孔隙结构的定量表征方法 |
CN103278864A (zh) * | 2013-05-10 | 2013-09-04 | 中国石油天然气股份有限公司 | 洞缝型储层的地质特征参数及分布的测定方法及装置 |
Non-Patent Citations (4)
Title |
---|
储层岩心裂缝与试件裂缝定量描述方法研究;黄辅琼等;《测井技术》;19971231;第21卷(第5期);356-360 * |
储集层构造裂缝描述与定量预测;汪必峰;《万方学位论文库》;20080923;第10,65-84以及图5-1,5-2,5-19,表5-7 * |
储集层裂缝预测分析的多参数判据法;闫相祯等;《石油勘探与开发》;20091231;第36卷(第6期);749-755 * |
构造裂缝多参数定量计算模型;季宗镇等;《中国石油大学学报(自然科学版)》;20100228;第34卷(第1期);24-27 * |
Also Published As
Publication number | Publication date |
---|---|
CN104992067A (zh) | 2015-10-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107301306B (zh) | 用于致密砂岩气藏压裂水平井的动态无阻流量预测方法 | |
CN102339339B (zh) | 一种分析缝洞型油藏剩余油分布的方法 | |
CN105239990B (zh) | 一种超低渗致密油藏拟本井侧向驱替水平井布井方法 | |
CN104992067B (zh) | 一种基于裂缝产状的裂缝分析方法 | |
CN108868748A (zh) | 一种页岩气水平井重复压裂裂缝开启压力的计算方法 | |
CN108691526B (zh) | 实验确定碳酸盐岩油气藏粗糙裂缝酸刻蚀导流能力的方法 | |
CN104018828B (zh) | 一种基于演化进程的曲流河砂体储层建筑结构分析方法 | |
CN107577831A (zh) | 一种缝洞型碳酸盐岩油气藏溶洞规模大小计算方法 | |
CN107220493A (zh) | 基于微地震事件的页岩气水平井网络裂缝建模方法 | |
CN102707333A (zh) | 页岩气资源量/储量的测量方法 | |
CN104747180A (zh) | 一种用于注水开发的缝洞型油藏分析方法及其应用 | |
CN103400020A (zh) | 一种测算多条相交离散裂缝流动状况的油藏数值模拟方法 | |
CN113027409B (zh) | 模拟水平井压裂缝网内支撑剂运移的实验装置 | |
CN109236273B (zh) | 油田开发生产动态数据处理方法 | |
CN107387051A (zh) | 低渗透非均质油藏多段压裂水平井重复压裂选井的方法 | |
CN107288626A (zh) | 一种页岩气地质甜度与工程甜度计算方法 | |
CN109209307A (zh) | 一种定量分析低渗透油藏注水开发效果的方法 | |
CN105350961A (zh) | 低渗透非均质应力敏感储层体积压裂水平井产量预测方法 | |
CN108661619A (zh) | 一种确定碳酸盐岩油气藏酸蚀填砂裂缝铺砂浓度的方法 | |
CN112069654B (zh) | 一种碳酸盐岩酸化数值模拟方法 | |
CN114186440B (zh) | 一种地质-工程“双轨制”页岩可压性综合评价方法 | |
CN106897531A (zh) | 一种低渗透石灰岩储层渗透率的定量评价方法 | |
CN114154430A (zh) | 一种压裂油藏co2驱油流动模拟方法 | |
CN104727789B (zh) | 中高渗砂岩油藏水驱波及系数及过水倍数动态描述方法 | |
Hu et al. | Simulation of immiscible water-alternating-CO2 flooding in the Liuhua oilfield offshore Guangdong, China |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |