CN112883574A - 利用水力压裂微地震b值来优化页岩气井重复压裂方法 - Google Patents
利用水力压裂微地震b值来优化页岩气井重复压裂方法 Download PDFInfo
- Publication number
- CN112883574A CN112883574A CN202110188357.7A CN202110188357A CN112883574A CN 112883574 A CN112883574 A CN 112883574A CN 202110188357 A CN202110188357 A CN 202110188357A CN 112883574 A CN112883574 A CN 112883574A
- Authority
- CN
- China
- Prior art keywords
- value
- fracturing
- well
- magnitude
- earthquake
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 11
- 238000009826 distribution Methods 0.000 claims abstract description 23
- 235000009508 confectionery Nutrition 0.000 claims description 5
- 230000009466 transformation Effects 0.000 claims description 3
- 235000021185 dessert Nutrition 0.000 abstract description 5
- 238000012544 monitoring process Methods 0.000 abstract description 3
- 208000010392 Bone Fractures Diseases 0.000 description 14
- 206010017076 Fracture Diseases 0.000 description 14
- 230000000694 effects Effects 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 3
- 230000018109 developmental process Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000001965 increasing effect Effects 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/25—Methods for stimulating production
- E21B43/26—Methods for stimulating production by forming crevices or fractures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V9/00—Prospecting or detecting by methods not provided for in groups G01V1/00 - G01V8/00
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Fluid Mechanics (AREA)
- Geophysics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种利用水力压裂微地震b值来优化页岩气井重复压裂的方法,基于水力压裂微地震监测事件的数量分布及震级大小等数据,采用地震学中描述频率‑震级分布的古登堡‑里克特关系来计算压裂阶段微地震b值的大小,并利用b值的分布云图来识别区域的应力差相对大小,为重复压裂寻找改造甜点区提供参考。本方法可在一定程度上促使重复压裂形成更加复杂的缝网,从而有效提升页岩气产能。
Description
技术领域
本发明涉及页岩气开采技术领域。
背景技术
位于四川盆地的长宁-威远页岩气区块是我国国家级页岩气示范开发区,在水力压裂技术的增产作用下,页岩气产量日渐显著。但由于四川盆地较复杂的地质构造,断裂发育及高应力差导致初次压裂往往不能形成较理想的复杂缝网,压裂提产效果达不到预期,同时页岩气井在开采一段时间后产能衰减过快,这成为制约我国页岩气有效开发的瓶颈问题。通过开展重复压裂能够破碎更多潜在页岩储层,在重复压裂过程中形成更加复杂的缝网,可应对初次压裂效果不理想及产能下降过快等问题。
重复压裂前需进行甜点井段优选(“甜点”代指油气产量较高的区域),其中一个重要因素是甜点井段的水平地应力差。应力差是形成复杂缝网过程中不可忽视的一个关键因素,在较小的应力差条件下压裂通常更容易形成复杂缝网,有利于压裂提产。
现有技术尚缺乏对初次压裂后井周应力差分布的有效识别,致使重复压裂设计的最优甜点区域难以明确。
发明内容
本发明提供了一种利用水力压裂微地震b值来优化页岩气井重复压裂的方法,以弥补现有技术思路无法准确识别初次压裂后井周应力差的分布情况。
本发明的目的通过以下技术方案来实现:一种利用水力压裂微地震b值来优化页岩气井重复压裂的方法,将多个压裂段视作不同区域分别来计算b值,并采用插值手段将其外推到井周区域,从而得到水平井段的井周b值分布云图,利用b值的区域分布特征来对重复压裂的改造甜点区进行优选。
具体的,它包括以下步骤:
S1、确定对象井的压裂分段情况;
S2、确定压裂阶段计算微地震b值所需要的参数,包括微地震的数量、微地震震级大小等;
S3、利用地震学中的古登堡-里克特关系来计算压裂阶段微地震b值的大小;
S4、基于计算的各阶段微地震b值,绘制b值分布云图;
S5、利用b值云图识别井周应力差较小的改造区域;
S6、优化重复压裂设计,提升页岩气产能。
本发明具有以下优点:
(1)b值是基于实测微地震数据进行反演得到的描述地震频率与震级之间关系的一种分形维数,能够表征诱发微地震事件区域的主要应力机制和地震事件的破坏模式,采用b值来进行描述,可以更为真实地反映应力差的分布情况。
(2)较小的b值被认为是与逆冲断层相关的,而逆冲断层往往承受着较高的应力水平,其应力差相对更大,因此b值可作为区域应力差水平的一种度量。相比利用数值手段来识别区域应力差大小,采用b值无疑极大降低了计算复杂度,从而更加适用于现场。
附图说明
图1为示例井第1压裂阶段微地震数量及震级分布示意图;
图2为示例井第1压裂阶段微地震b值计算示意图;
图3为示例井井周b值分布示意图;
图4为示例井利用b值识别最优改造区域的示意图;
图5为本发明的流程图;
具体实施方式
下面结合附图对本发明做进一步的描述,本发明的保护范围不局限于以下所述:
在地震领域中,b值的计算往往只针对某一片特定区域,而油气开采领域中页岩气井水平井段的压裂施工被划分为多个子井段(压裂段)分别进行,本发明将多个压裂段视作不同区域分别来计算b值,并采用插值手段将其外推到井周区域,从而得到水平井段的井周b值分布云图,利用b值的区域分布特征来对重复压裂的改造甜点区进行优选。
S1、确定对象井的压裂分段情况,选取长宁区块某井进行示例分析,该井采取了分段压裂方式进行改造,水平井段共被划分为18个压裂段。
S2、确定计算微地震b值的基本参数,包括每个压裂段的微地震信号数量N及震级大小M。本实施例井的第1压裂段的微地震数据如图1所示,矩形方块代表示例井第1段压裂产生的微地震信号,星形符号代表震级较大的微地震信号。
S3、在地震领域中,b值可以作为区域应力差水平的一种度量,通常较小的b值可以指示应力差相对较大的区域,反之较大的b值则可以指示应力差更小的区域。那么基于b值的此种特征,将其运用到工程领域就可以很好地帮助识别重复压裂的改造甜点区域。计算第1段压裂微地震b值,b值计算公式为地震学中描述频率-震级分布的古登堡-里克特关系:
log10 N(m≥M)=a-b*M---------------------------------------------(1)
式中,N(m≥M)表示震级m大于或等于M的地震数量,b值是频率-震级分布线性部分斜率的绝对值,a值表征该区域地震的活性,当M=0时,就表示了该区域所有的地震数目,即10a个。如图2所示的指频率-震级图,N对应步骤S2中确定的第1压裂段微地震数量,m对应步骤S2中确定的第1压裂段微地震震级,M指代频率-震级图中的x轴。利用最大曲率解可求得第1段压裂微地震b值约为0.72,如图2所示。
S4、基于步骤S1确定的压裂分段情况,按分段数重复步骤S2、S3,依次计算后续压裂段的b值,并采用插值手段将其外推到井周区域,绘制示例井所有压裂段的井周b值分布云图,如图3所示。从图3中可以看到该示例井微地震监测的大震级信号均处在b值较低的区域,大震级信号常对应断层激活的剪切破碎带,且断层在应力差较大的条件下易于发生剪切破坏,说明b值与应力差存在明显相关性。此外,观察均匀分布的小震级信号可以发现,其绝大部分处在b值较高的区域,在应力差较小的条件下,压裂更容易形成均匀分布的复杂缝网,因此b值很好地反映了井周应力差的分布情况。
S5、目前现有技术尚缺乏对初次压裂后井周应力差分布的有效识别,本发明创新性地提出了利用微地震数据计算的b值来识别井周应力差的分布,以弥补现有技术的不足,进一步提高重复压裂的产能。基于步骤S4绘制的井周b值分布云图,将图中b值较大的区域圈出,该区域即为重复压裂的最优改造区域。图4展示了如何利用b值分布云图来识别重复压裂最优改造区域。从微地震的监测情况可知,该示例井的初次压裂效果并不理想,尤其是中间部分几乎没有微地震信号出现,说明初次压裂未实现对中间部分储层的有效改造。从该示例井的b值分布可看出,中间部分靠近井根的部位有一片b值较大的区域(b>1),说明该区域的应力差相对较小,是重复压裂改造的甜点井段区域。
S6、针对上述确定的重复压裂甜点井段进行重复压裂设计和优化,提高重复压裂效果,提升页岩气产能。
Claims (2)
1.一种利用水力压裂微地震b值来优化页岩气井重复压裂的方法,其特征是,将多个压裂段视作不同区域分别来计算b值,并采用插值手段将其外推到井周区域,从而得到水平井段的井周b值分布云图,利用b值的区域分布特征来对重复压裂的改造甜点区进行优选。
2.如权利要求1所述的方法,其特征是,具体的,它包括以下步骤:
S1、确定对象井的压裂分段情况;
S2、确定压裂阶段计算微地震b值所需要的参数,包括微地震的数量、微地震震级大小;
S3、利用地震学中的古登堡-里克特关系来计算压裂阶段微地震b值的大小;
S4、基于计算的各阶段微地震b值,绘制b值分布云图;
S5、利用b值云图识别井周应力差较小的改造区域;
S6、优化重复压裂设计,提升页岩气产能。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110188357.7A CN112883574B (zh) | 2021-02-18 | 2021-02-18 | 利用水力压裂微地震b值来优化页岩气井重复压裂方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110188357.7A CN112883574B (zh) | 2021-02-18 | 2021-02-18 | 利用水力压裂微地震b值来优化页岩气井重复压裂方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112883574A true CN112883574A (zh) | 2021-06-01 |
CN112883574B CN112883574B (zh) | 2022-08-09 |
Family
ID=76056514
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110188357.7A Active CN112883574B (zh) | 2021-02-18 | 2021-02-18 | 利用水力压裂微地震b值来优化页岩气井重复压裂方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112883574B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106845043A (zh) * | 2017-04-07 | 2017-06-13 | 东方宝麟科技发展(北京)有限公司 | 一种页岩气水平井重复压裂工艺流程及设计方法 |
CN108952691A (zh) * | 2017-05-17 | 2018-12-07 | 中国石油化工股份有限公司 | 一种水平井分段压裂多段裂缝形态评价方法 |
CN109883992A (zh) * | 2019-01-02 | 2019-06-14 | 西安电子科技大学 | 一种仿真流场图像rgb与折射率场分布的转化方法 |
CN110888161A (zh) * | 2018-09-07 | 2020-03-17 | 中国石油化工股份有限公司 | 微地震震级-频度关系中的b值计算方法及*** |
-
2021
- 2021-02-18 CN CN202110188357.7A patent/CN112883574B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106845043A (zh) * | 2017-04-07 | 2017-06-13 | 东方宝麟科技发展(北京)有限公司 | 一种页岩气水平井重复压裂工艺流程及设计方法 |
CN108952691A (zh) * | 2017-05-17 | 2018-12-07 | 中国石油化工股份有限公司 | 一种水平井分段压裂多段裂缝形态评价方法 |
CN110888161A (zh) * | 2018-09-07 | 2020-03-17 | 中国石油化工股份有限公司 | 微地震震级-频度关系中的b值计算方法及*** |
CN109883992A (zh) * | 2019-01-02 | 2019-06-14 | 西安电子科技大学 | 一种仿真流场图像rgb与折射率场分布的转化方法 |
Non-Patent Citations (1)
Title |
---|
李政等: "微地震方法的裂缝监测与储层评价", 《地球物理学报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN112883574B (zh) | 2022-08-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110704888B (zh) | 一种非常规油气藏加密井体积压裂施工参数优化设计方法 | |
US11371339B2 (en) | Method for determining hydraulic fracture orientation and dimension | |
US11009620B2 (en) | Method for determining favorable time window of infill well in unconventional oil and gas reservoir | |
Lorenzoni et al. | Structural health monitoring of the Roman Arena of Verona, Italy | |
CN104331745B (zh) | 油气藏内天然裂缝的分期、分成因预测评价方法 | |
CN112364422B (zh) | 基于mic-lstm的盾构施工地表变形动态预测方法 | |
Urbancic et al. | Effective constraint of RTA models utilizing microseismicity derived flow attributes | |
CN111123374A (zh) | 一种基于匹配滤波的探地雷达全波形反演方法 | |
CN113530536A (zh) | 致密砂岩气藏水平井压裂缝储层动用效率评价方法及*** | |
Brinkley et al. | An Eagle Ford case study: Monitoring fracturing propagation through sealed wellbore pressure monitoring | |
CN112883574B (zh) | 利用水力压裂微地震b值来优化页岩气井重复压裂方法 | |
CN109298448B (zh) | 一种致密气压裂工程甜点的预测方法和装置 | |
CN113534247B (zh) | 基于叠后地震数据的裂缝定量化预测方法及装置 | |
CN105467469A (zh) | 一种预测致密低渗透非均质储层中构造裂缝优势方向及其密度的方法 | |
CN112836442A (zh) | 一种页岩油水平井网老井水力裂缝注液量的确定方法 | |
CN106844858B (zh) | 地层裂缝发育区带预测方法及装置 | |
CN111663935B (zh) | 一种用于评价井筒环空气密性的方法及*** | |
CN112100796A (zh) | 防止页岩气井中套管受损的钻井轨迹确定方法及装置 | |
Yong et al. | Development Strategy Optimization of Ning201 Block Longmaxi Shale Gas | |
CN115828636B (zh) | 一种页岩气井组压裂的防压窜施工参数优化方法 | |
Manuello Bertetto et al. | Acoustic Emission Monitoring and Thrust Network Analysis of the Central Nave Vaults of the Turin Cathedral | |
Yang et al. | Adit deformation prediction based on ICEEMDAN dispersion entropy and LSTM-BP | |
Morton et al. | Self-similarity in rock fracturing and the behaviour of large-scale faults in the mining environment | |
Potsch et al. | Determination of rock mass behaviour types, a case study | |
CN118296891A (zh) | 一种岩石裂缝起裂预测方法、***、设备及存储介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |