CN110056336A - 一种页岩气缝网压裂施工压力曲线自动诊断方法 - Google Patents
一种页岩气缝网压裂施工压力曲线自动诊断方法 Download PDFInfo
- Publication number
- CN110056336A CN110056336A CN201910468722.2A CN201910468722A CN110056336A CN 110056336 A CN110056336 A CN 110056336A CN 201910468722 A CN201910468722 A CN 201910468722A CN 110056336 A CN110056336 A CN 110056336A
- Authority
- CN
- China
- Prior art keywords
- pressure
- net
- net pressure
- follows
- parameter
- 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 39
- 238000003745 diagnosis Methods 0.000 title claims abstract description 22
- 238000004364 calculation method Methods 0.000 claims abstract description 29
- 239000012530 fluid Substances 0.000 claims abstract description 29
- 238000003909 pattern recognition Methods 0.000 claims abstract description 15
- 238000006243 chemical reaction Methods 0.000 claims abstract description 11
- 230000008569 process Effects 0.000 claims abstract description 11
- 238000004458 analytical method Methods 0.000 claims abstract description 8
- 238000013517 stratification Methods 0.000 claims description 18
- 238000004088 simulation Methods 0.000 claims description 12
- 239000007788 liquid Substances 0.000 claims description 11
- 230000003068 static effect Effects 0.000 claims description 11
- 238000010276 construction Methods 0.000 claims description 9
- 238000000418 atomic force spectrum Methods 0.000 claims description 8
- 230000008859 change Effects 0.000 claims description 8
- 230000004888 barrier function Effects 0.000 claims description 4
- 230000005484 gravity Effects 0.000 claims description 2
- 238000009877 rendering Methods 0.000 claims description 2
- 238000013459 approach Methods 0.000 claims 1
- 208000010392 Bone Fractures Diseases 0.000 description 36
- 206010017076 Fracture Diseases 0.000 description 36
- 239000004576 sand Substances 0.000 description 6
- 238000011161 development Methods 0.000 description 4
- 241000208340 Araliaceae Species 0.000 description 3
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 3
- 235000003140 Panax quinquefolius Nutrition 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 3
- 239000003795 chemical substances by application Substances 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 235000008434 ginseng Nutrition 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 239000011435 rock Substances 0.000 description 3
- 239000000243 solution Substances 0.000 description 3
- 238000002405 diagnostic procedure Methods 0.000 description 2
- 238000005065 mining Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013211 curve analysis Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
Classifications
-
- 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
- E21B47/00—Survey of boreholes or wells
- E21B47/06—Measuring temperature or pressure
-
- 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
-
- 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
- E21B43/267—Methods for stimulating production by forming crevices or fractures reinforcing fractures by propping
-
- 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
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Algebra (AREA)
- Geophysics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Databases & Information Systems (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Excavating Of Shafts Or Tunnels (AREA)
- Force Measurement Appropriate To Specific Purposes (AREA)
Abstract
本发明公开了一种页岩气缝网压裂施工压力曲线自动诊断方法,包括以下步骤:S1、运用流体动力学理论建立页岩气缝网压裂井底净压力折算模型;S2、运用数值分析理论建立井底净压力动态拟合模型;S3、运用裂缝延伸理论建立裂缝延伸模式识别模型;S4、结合步骤S1、S2、S3,通过数值计算流程实现对页岩气缝网压裂施工压力曲线进行自动诊断。本发明实现了对页岩气缝网压裂施工压力曲线的自动、快速、准确诊断分析,解决了页岩气缝网压裂施工压力曲线特征复杂多变,缺乏针对性的压力曲线自动诊断方法的问题。
Description
技术领域
本发明涉及页岩气开发技术领域,具体的涉及一种页岩气缝网压裂施工压力曲线自动诊断方法。
背景技术
水力压裂施工压力曲线诊断是压裂工艺技术理论的一个重要环节。施工压力是裂缝动态扩展力学行为的外在表现,而裂缝动态扩展与储层岩石力学性质、压裂流体性质、地应力分布、天然裂缝发育状况等地层特性紧密相关。因此,通过施工压力曲线分析,可以识别水力裂缝延伸行为,判断地层特性,为压裂优化设计和现场施工参数调整提供依据和指导。
目前,国内外学者研究主要围绕常规压裂开展了施工压力曲线特征研究。1981年,Nolte和Smith将压力双对数曲线归纳为四种典型情况:Ⅰ、斜率约为1/4的阶段:表示裂缝在长度方向延伸,这是正常的施工曲线;Ⅱ、斜率约为0的阶段:表示地层内天然微裂隙张开,使得滤失量与注入量持平,另外还可能是压力超过上覆地层应力,形成T型缝;Ⅲ、斜率约为1的阶段:表示裂缝延伸受阻,缝内压力急剧上升;如果斜率大于1则表示裂缝内发生了堵塞,这种情况下应合理控制施工砂比和施工排量,以保证施工顺利进行,而对于端部脱砂压裂施工,则希望支撑剂在一定缝长时形成砂堵,然后通过控制排量和砂比,使裂缝完全填满;对于常规的加砂压裂,出现斜率为1或大于1时,应立即采取措施,以免井筒内发生砂卡;Ⅳ、斜率为负的阶段:表示裂缝穿过低应力层,裂缝发生不稳定增长,也可能是沟通了天然裂缝,使滤失量大大增加,此结果会导致裂缝内砂堵,压力又将很快上升。随后,杨悦(2012)、曲冠政(2013)等国内学者将压力双对数曲线进一步细分为六大类型,并对施工早期压力下降进行了解释。Martinez(1993)、Fan(1995)、曹学军(2002)、刘先灵(2003)、王正茂(2004)等学者基于常规压裂施工压力曲线特征研究成果,通过建立裂缝延伸反演模型,提出了多种常规压裂施工压力曲线分析方法。此外,王腾飞(2009)、张喆(2014)等学者分别借助分形几何、神经网络等理论,对常规压裂施工压力曲线分析方法进行了改进。
近年来,国内外页岩气开发规模迅速扩大。然而,对于具有低孔、超低渗特征的页岩气藏,常规压裂已无法满足其商业化开发的需求,需采用水平井分段分簇缝网压裂方式进行增产改造。由于页岩储层非均质性较强、天然裂缝、层理缝等力学结构弱面较为发育、压裂液滤失量较大等特点,水力裂缝通常以非均匀、非连续性的形式进行延伸扩展形成裂缝网络,导致页岩压裂过程中的压力曲线与常规压裂的单一特征曲线不同,而往往呈现出复杂的多阶段、多变化特征,因此常规的压力曲线诊断方法并不适用于页岩气缝网压裂。但是,目前国内外尚未针对页岩气缝网压裂提出系统的施工压力曲线诊断方法。
因此,亟需建立一种适用于页岩气缝网压裂施工压力曲线自动诊断方法,将有助于提高页岩气缝网压裂设计的科学性和现场施工参数调整的针对性,进一步挖掘页岩气藏增产潜力。
发明内容
针对上述问题,本发明提供一种页岩气缝网压裂施工压力曲线自动诊断方法,目的在于,解决页岩气缝网压裂施工压力曲线特征复杂多变,缺乏针对性的压力曲线自动诊断方法的问题。
本发明采用下述的技术方案:
一种页岩气缝网压裂施工压力曲线自动诊断方法,包括以下步骤:
S1、运用流体动力学理论建立页岩气缝网压裂井底净压力折算模型;
S2、运用数值分析理论建立井底净压力动态拟合模型;
S3、运用裂缝延伸理论建立裂缝延伸模式识别模型;
S4、结合步骤S1、S2、S3,通过数值计算流程实现对页岩气缝网压裂施工压力曲线进行自动诊断。
优选的,所述步骤S1中建立页岩气缝网压裂井底净压力折算模型包括以下步骤:
A、收集以下数据,包括:井身结构、井筒参数、射孔参数、压裂施工参数、压裂液参数、支撑剂参数、地应力条件数据、天然裂缝参数、层理缝参数;
B、利用步骤A中的数据,计算井筒内液体流动压降、射孔孔眼压降、携砂液柱静压力、井底净压力,并绘制井底净压力曲线。
优选的,所述井身结构包括井筒长度、井筒垂直深度,井筒参数包括井筒直径、井筒管壁绝对粗糙度,射孔参数包括孔眼数量、孔眼直径、孔眼流量系数,压裂施工参数包括井口压力、泵注排量、支撑剂比例,压裂液参数包括粘度、密度,支撑剂参数为支撑剂密度,地应力条件包括最小水平主应力、最大水平主应力、垂向应力、储层与上下隔层应力差,天然裂缝参数包括逼近角、倾角、内聚力、摩擦系数、抗张强度,层理缝参数为抗张强度。
优选的,所述井筒内液体流动压降计算公式如下:
式中:Δpwf—井筒内液体流动压降,Pa;λ—水力摩阻系数,无量纲;L—井筒长度,m;D—井筒直径,m;v—压裂液流速,m/s;—压裂液密度,kg/m3;
所述压裂液流速v为:
式中:q—压裂液流量,即泵注排量,m3/s;
所述水力摩阻系数(λ)的取值与管流的流动状态有关;
①当Re<2000时,流动状态为层流,水力摩阻系数为:
λ=64/Re (3)
②当2000<Re<59.7/ε8/7时,流动状态为水力光滑紊流,水力摩阻系数为:
λ=0.3164/Re0.25 (4)
③当59.7/ε8/7<Re<(665-765lgε)/ε时,流动状态为混合摩擦紊流,水力摩阻系数为:
④当Re>(665-765lgε)/ε时,流动状态为粗糙紊流,水力摩阻系数为:
λ=1/[2lg(3.7D/R)]2 (6)
式中:Re—雷诺数,无量纲;ε—井筒管壁相对粗糙度,m;R—井筒管壁绝对粗糙度,m;
所述管壁相对粗糙度(ε)和雷诺数(Re)分别为:
ε=2R/D (7)
式中:μ—压裂液粘度,Pa·s;
所述射孔孔眼压降,计算公式如下:
式中:npf—射孔孔眼数量,个;dpf—射孔孔眼直径,m;α—孔眼流量系数,一般取0.8~0.85,无量纲;
井筒所述内携砂液柱静压力计算公式如下:
pwh=[(1-Vp)ρl+Vpρp]gZ (10)
式中:pwh—井筒内携砂液柱静压力,Pa;Vp—支撑剂比例,%;ρp—支撑剂密度,kg/m3;g—重力加速度,m/s2;Z—井筒垂直深度,m;
由于井筒直径和携砂液密度随井筒深度变化而变化,因此井筒内携砂液柱静压力(pwh)和井筒内液体流动压降(Δpwf)需要分别取井筒垂直深度和长度微元计算并积分求得;所述井底净压力计算公式如下:
式中:p—井底净压力,Pa;ph—井口压力,Pa;σhmin—地层最小水平主应力,Pa。
优选的,所述步骤S2中建立井底净压力动态拟合模型包括以下步骤:
a、建立井底净压力数据序列p1,p2,p3,…,…,pN以及相对应的时间数据序列t1,t2,t3,…,tN,并依次提取数据序列进行计算;
b、将净压力数据序列和时间数据序列中的当前数据分别设为参考压力数据和参考时间数据,即把参考序号r赋值为当前序号;
c、取下一个净压力和时间数据为当前数据,设当前序号为i;
d、计算出当前净压力指数拟合值和斜率拟合值,当前净压力指数平均值和斜率平均值,当前净压力拟合相对误差值;
e、若当前净压力拟合相对误差值>10%,则返回步骤b,将当前净压力和时间设为参考净压力和时间,即把参考序号r重新赋值为当前序号i;若当前净压力拟合相对误差值≤10%,则返回步骤c,继续进行下一个数据的计算;
f、当所有数据计算完毕时,即i=N时,绘制净压力指数平均值曲线。
优选的,所述净压力指数拟合值的计算公式如下:
式中:n—净压力指数拟合值,无量纲;p—井底净压力,Pa;t—时间,s;下标:i—当前净压力和当前时间数据序号,无量纲;r—参考净压力和参考时间数据序号,无量纲;
所述净压力斜率拟合值的计算公式如下:
式中:k—净压力斜率拟合值,Pa/(s)e;
所述净压力指数平均值的计算公式如下:
式中:—净压力指数平均值,无量纲;
所述净压力斜率平均值的计算公式如下:
式中:—净压力斜率平均值,Pa/(s)n;
所述净压力拟合相对误差值的计算公式如下:
式中:χ—净压力拟合相对误差值,无量纲。
优选的,所述步骤S3中,页岩气缝网压裂过程中水力裂缝延伸模式包括:缝网延伸模式、延伸受阻模式、正常延伸模式、层理缝延伸模式、缝高方向延伸模式、快速滤失模式,各类模式对应着不同的压力曲线特征识别方程,建立裂缝延伸模式识别模型包括以下步骤:
Ⅰ、利用井底净压力与净压力指数平均值数据序列,依次判断是否为缝网延伸模式、延伸受阻模式、正常延伸模式、层理缝延伸模式、缝高延伸模式、快速滤失模式,若数据满足其中一个模式,则继续判断下一个数据点;
Ⅱ、当所有数据点处理完毕后,停止计算,输出裂缝延伸模式识别曲线。
优选的,所述缝网延伸模式的识别方程如下:
其中:
式中:σv—地层垂向应力,Pa;σHmax—地层最大水平主应力,Pa;θnf—地层天然裂缝逼近角,°;—地层天然裂缝倾角,°;Stnf—地层天然裂缝抗张强度,Pa;τo—地层天然裂缝内聚力,Pa;knf—地层天然裂缝摩擦系数,无量纲;σnf—地层天然裂缝壁面所受法向应力值,Pa;στ—地层天然裂缝壁面所受切向应力值,Pa;ΔS—储层与上下隔层应力差,Pa;
所述延伸受阻模式的识别方程如下:
所述正常延伸模式的识别方程如下:
所述层理缝延伸模式的识别方程如下:
式中:Stnf—层理缝抗张强度,Pa;
所述缝高方向延伸模式的识别方程如下:
所述快速滤失模式的识别方程如下:
优选的,所述步骤S4的数值计算流程包括以下步骤:
S901、参数输入:井身结构参数、井筒参数、射孔参数、压裂施工参数、压裂液参数、支撑剂参数、地应力条件、天然裂缝参数、层理缝参数;
S902、计算井筒内液体流动压降、射孔孔眼压降、携砂液柱静压力;将压裂施工井口压力折算为井底净压力;
S903、建立井底净压力数据与时间数据序列;对井底净压力进行动态拟合,并计算净压力指数平均值;
S904、基于井底净压力和净压力指数平均值,自动识别各段对应的裂缝延伸模式;
S905、数据输出:井底净压力值、净压力指数平均值、裂缝延伸模式;
S906、图像绘制:井底净压力曲线、净压力指数平均值曲线、裂缝延伸模式识别曲线,以及施工压力曲线诊断图。
本发明的有益效果是:
与现有技术相比,该方法专门针对页岩气藏储层特征,结合页岩气缝网压裂工艺特点,建立了页岩气缝网压裂井底净压力折算模型、井底净压力动态拟合模型、裂缝延伸模式识别模型,从而提出了一种适用于页岩气缝网压裂的施工压力曲线自动诊断方法。
由于该方法综合系统地考虑了页岩气缝网压裂过程中常见的六种裂缝延伸模式对施工压力的影响,为此,可对页岩气缝网压裂施工压力曲线进行自动、快速、准确的诊断分析,解决了页岩气缝网压裂施工压力曲线特征复杂多变,缺乏针对性的压力曲线自动诊断方法的问题。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例的附图作简单地介绍,显而易见地,下面描述中的附图仅仅涉及本发明的一些实施例,而非对本发明的限制。
图1为本发明页岩气缝网压裂现场施工曲线图;
图2为本发明方法的计算流程框图;
图3为本发明井底净压力曲线;
图4为本发明井底净压力指数平均值曲线;
图5为本发明裂缝延伸模式识别曲线;
图6为本发明施工压力曲线诊断图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例的附图,对本发明实施例的技术方案进行清楚、完整地描述。显然,所描述的实施例是本发明的一部分实施例,而不是全部的实施例。基于所描述的本发明的实施例,本领域普通技术人员在无需创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
除非另外定义,本公开使用的技术术语或者科学术语应当为本公开所属领域内具有一般技能的人士所理解的通常意义。本公开中使用的“包括”或者“包含”等类似的词语意指出现该词前面的元件或者物件涵盖出现在该词后面列举的元件或者物件及其等同,而不排除其他元件或者物件。“上”、“下”、“左”、“右”等仅用于表示相对位置关系,当被描述对象的绝对位置改变后,则该相对位置关系也可能相应地改变。
下面结合附图和实施例对本发明进一步说明。
如图1至图6所示,一种页岩气缝网压裂施工压力曲线自动诊断方法,包括以下步骤:
已知某一页岩气藏地质工程矿场实际参数如表1所示:
表1页岩气藏地质工程矿场实际参数
该页岩气藏内某水平井缝网压裂现场施工曲线如图1所示,井身结构参数如表2所示。
表2井身结构参数
根据本发明方法的计算流程(如图2所示)开展实例计算:
S1、运用流体动力学理论建立页岩气缝网压裂井底净压力折算模型,用于将压裂施工井口压力曲线折算为井底净压力曲线。
具体步骤为:
①利用计算公式(1)~(8)计算井筒内液体流动压降;
②利用公式(9)计算射孔孔眼压降;
③利用公式(10)计算携砂液柱静压力;
④利用公式(11)折算井底净压力,并绘制井底净压力曲线。
基于页岩气缝网压裂井底净压力折算模型,折算得到井底净压力曲线如图3所示。
S2、运用数值分析理论建立井底净压力动态拟合模型,用于计算井底净压力指数平均值,为下一步裂缝延伸模式识别提供依据。
具体步骤为:
a、建立井底净压力数据序列p1,p2,p3,…,…,pN以及相对应的时间数据序列t1,t2,t3,…,tN,并依次提取数据序列进行计算;
b、将净压力数据序列和时间数据序列中的当前数据分别设为参考压力数据和参考时间数据,即把参考序号r赋值为当前序号;
c、取下一个净压力和时间数据为当前数据,设当前序号为i;
d、计算出当前净压力指数拟合值和斜率拟合值,当前净压力指数平均值和斜率平均值,当前净压力拟合相对误差值;
e、若当前净压力拟合相对误差值>10%,则返回步骤b,将当前净压力和时间设为参考净压力和时间,即把参考序号r重新赋值为当前序号i;若当前净压力拟合相对误差值≤10%,则返回步骤c,继续进行下一个数据的计算;
f、当所有数据计算完毕时,即i=N时,绘制净压力指数平均值曲线。基于井底净压力动态拟合模型,计算得到井底净压力指数平均值曲线如图4所示。
S3、运用裂缝延伸理论建立裂缝延伸模式识别模型,用于根据井底净压力与净压力指数平均值识别裂缝延伸模式。
具体步骤为:
①依次处理井底净压力与净压力指数平均值数据序列;
②判断是否满足方程式(17)~(18),若满足,则识别为缝网延伸模式,并继续判断下一个数据点;
③判断是否满足方程式(19),若满足,则识别为延伸受阻模式,并继续判断下一个数据点;
④判断是否满足方程式(20),若满足,则识别为正常延伸模式,并继续判断下一个数据点;
⑤判断是否满足方程式(21),若满足,则识别为层理缝延伸模式,并继续判断下一个数据点;
⑥判断是否满足方程式(22),若满足,则识别为缝高延伸模式,并继续判断下一个数据点;
⑦判断是否满足方程式(23),若满足,则识别为快速滤失模式,并继续判断下一个数据点;
⑧当所有数据点处理完毕后,停止计算,输出裂缝延伸模式识别曲线。
基于裂缝延伸模式识别模型,判断页岩缝网压裂过程中不同阶段的裂缝延伸模式见表3,裂缝延伸模式识别曲线如图5所示。
表3裂缝延伸模式表
S4、结合步骤S1、S2、S3,通过数值计算流程实现对页岩气缝网压裂施工压力曲线进行自动诊断。
基于上述计算流程和结果,绘制施工压力曲线诊断图,如图6所示。可以看出,该页岩气藏水平井压裂施工过程中,裂缝延伸模式依次为缝高方向延伸→快速滤失→正常延伸→延伸受阻→缝网延伸→层理缝延伸→缝高方向延伸→快速滤失。其中,有利于压裂缝网形成的阶段(缝网延伸与层理缝延伸)总时长为43min,占总压裂施工时长的34.13%。
以上所述,仅是本发明的较佳实施例而已,并非对本发明作任何形式上的限制,虽然本发明已以较佳实施例揭露如上,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,当可利用上述揭示的技术内容作出些许更动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。
Claims (9)
1.一种页岩气缝网压裂施工压力曲线自动诊断方法,其特征在于,包括以下步骤:
S1、运用流体动力学理论建立页岩气缝网压裂井底净压力折算模型;
S2、运用数值分析理论建立井底净压力动态拟合模型;
S3、运用裂缝延伸理论建立裂缝延伸模式识别模型;
S4、结合步骤S1、S2、S3,通过数值计算流程实现对页岩气缝网压裂施工压力曲线进行自动诊断。
2.根据权利要求1所述的一种页岩气缝网压裂施工压力曲线自动诊断方法,其特征在于,所述步骤S1中建立页岩气缝网压裂井底净压力折算模型包括以下步骤:
A、收集以下数据,包括:井身结构、井筒参数、射孔参数、压裂施工参数、压裂液参数、支撑剂参数、地应力条件数据、天然裂缝参数、层理缝参数;
B、利用步骤A中的数据,计算井筒内液体流动压降、射孔孔眼压降、携砂液柱静压力、井底净压力,并绘制井底净压力曲线。
3.根据权利要求2所述的一种页岩气缝网压裂施工压力曲线自动诊断方法,其特征在于,所述井身结构包括井筒长度、井筒垂直深度,井筒参数包括井筒直径、井筒管壁绝对粗糙度,射孔参数包括孔眼数量、孔眼直径、孔眼流量系数,压裂施工参数包括井口压力、泵注排量、支撑剂比例,压裂液参数包括粘度、密度,支撑剂参数为支撑剂密度,地应力条件包括最小水平主应力、最大水平主应力、垂向应力、储层与上下隔层应力差,天然裂缝参数包括逼近角、倾角、内聚力、摩擦系数、抗张强度,层理缝参数为抗张强度。
4.根据权利要求2所述的一种页岩气缝网压裂施工压力曲线自动诊断方法,其特征在于,所述井筒内液体流动压降计算公式如下:
式中:Δpwf—井筒内液体流动压降,Pa;λ—水力摩阻系数,无量纲;L—井筒长度,m;D—井筒直径,m;v—压裂液流速,m/s;ρ—压裂液密度,kg/m3;
所述压裂液流速v为:
式中:q—压裂液流量,即泵注排量,m3/s;
所述水力摩阻系数(λ)的取值与管流的流动状态有关;
①当Re<2000时,流动状态为层流,水力摩阻系数为:
λ=64/Re (3)
②当2000<Re<59.7/ε8/7时,流动状态为水力光滑紊流,水力摩阻系数为:
λ=0.3164/Re0.25 (4)
③当59.7/ε8/7<Re<(665-765lgε)/ε时,流动状态为混合摩擦紊流,水力摩阻系数为:
④当Re>(665-765lgε)/ε时,流动状态为粗糙紊流,水力摩阻系数为:
λ=1/[2lg(3.7D/R)]2 (6)
式中:Re—雷诺数,无量纲;ε—井筒管壁相对粗糙度,m;R—井筒管壁绝对粗糙度,m;
所述管壁相对粗糙度(ε)和雷诺数(Re)分别为:
ε=2R/D (7)
式中:μ—压裂液粘度,Pa·s;
所述射孔孔眼压降,计算公式如下:
式中:npf—射孔孔眼数量,个;dpf—射孔孔眼直径,m;α—孔眼流量系数,一般取0.8~0.85,无量纲;
井筒所述内携砂液柱静压力计算公式如下:
pwh=[(1-Vp)ρl+Vpρp]gZ (10)
式中:pwh—井筒内携砂液柱静压力,Pa;Vp—支撑剂比例,%;ρp—支撑剂密度,kg/m3;g—重力加速度,m/s2;Z—井筒垂直深度,m;
由于井筒直径和携砂液密度随井筒深度变化而变化,因此井筒内携砂液柱静压力(pwh)和井筒内液体流动压降(Δpwf)需要分别取井筒垂直深度和长度微元计算并积分求得;所述井底净压力计算公式如下:
式中:p—井底净压力,Pa;ph—井口压力,Pa;σhmin—地层最小水平主应力,Pa。
5.根据权利要求1所述的一种页岩气缝网压裂施工压力曲线自动诊断方法,其特征在于,所述步骤S2中建立井底净压力动态拟合模型包括以下步骤:
a、建立井底净压力数据序列p1,p2,p3,…,…,pN以及相对应的时间数据序列t1,t2,t3,…,tN,并依次提取数据序列进行计算;
b、将净压力数据序列和时间数据序列中的当前数据分别设为参考压力数据和参考时间数据,即把参考序号r赋值为当前序号;
c、取下一个净压力和时间数据为当前数据,设当前序号为i;
d、计算出当前净压力指数拟合值和斜率拟合值,当前净压力指数平均值和斜率平均值,当前净压力拟合相对误差值;
e、若当前净压力拟合相对误差值>10%,则返回步骤b,将当前净压力和时间设为参考净压力和时间,即把参考序号r重新赋值为当前序号i;若当前净压力拟合相对误差值≤10%,则返回步骤c,继续进行下一个数据的计算;
f、当所有数据计算完毕时,即i=N时,绘制净压力指数平均值曲线。
6.根据权利要求5所述的一种页岩气缝网压裂施工压力曲线自动诊断方法,其特征在于,所述净压力指数拟合值的计算公式如下:
式中:n—净压力指数拟合值,无量纲;p—井底净压力,Pa;t—时间,s;下标:i—当前净压力和当前时间数据序号,无量纲;r—参考净压力和参考时间数据序号,无量纲;
所述净压力斜率拟合值的计算公式如下:
式中:k—净压力斜率拟合值,Pa/(s)e;
所述净压力指数平均值的计算公式如下:
式中:—净压力指数平均值,无量纲;
所述净压力斜率平均值的计算公式如下:
式中:—净压力斜率平均值,Pa/(s)n;
所述净压力拟合相对误差值的计算公式如下:
式中:χ—净压力拟合相对误差值,无量纲。
7.根据权利要求1所述的一种页岩气缝网压裂施工压力曲线自动诊断方法,其特征在于,所述步骤S3中,页岩气缝网压裂过程中水力裂缝延伸模式包括:缝网延伸模式、延伸受阻模式、正常延伸模式、层理缝延伸模式、缝高方向延伸模式、快速滤失模式,各类模式对应着不同的压力曲线特征识别方程,建立裂缝延伸模式识别模型包括以下步骤:
Ⅰ、利用井底净压力与净压力指数平均值数据序列,依次判断是否为缝网延伸模式、延伸受阻模式、正常延伸模式、层理缝延伸模式、缝高延伸模式、快速滤失模式,若数据满足其中一个模式,则继续判断下一个数据点;
Ⅱ、当所有数据点处理完毕后,停止计算,输出裂缝延伸模式识别曲线。
8.根据权利要求7所述的一种页岩气缝网压裂施工压力曲线自动诊断方法,其特征在于,所述缝网延伸模式的识别方程如下:
其中:
式中:σv—地层垂向应力,Pa;σHmax—地层最大水平主应力,Pa;θnf—地层天然裂缝逼近角,°;—地层天然裂缝倾角,°;Stnf—地层天然裂缝抗张强度,Pa;τo—地层天然裂缝内聚力,Pa;knf—地层天然裂缝摩擦系数,无量纲;σnf—地层天然裂缝壁面所受法向应力值,Pa;στ—地层天然裂缝壁面所受切向应力值,Pa;ΔS—储层与上下隔层应力差,Pa;
所述延伸受阻模式的识别方程如下:
所述正常延伸模式的识别方程如下:
所述层理缝延伸模式的识别方程如下:
式中:Stnf—层理缝抗张强度,Pa;
所述缝高方向延伸模式的识别方程如下:
所述快速滤失模式的识别方程如下:
9.根据权利要求1所述的一种页岩气缝网压裂施工压力曲线自动诊断方法,其特征在于,所述步骤S4的数值计算流程包括以下步骤:
S901、参数输入:井身结构参数、井筒参数、射孔参数、压裂施工参数、压裂液参数、支撑剂参数、地应力条件、天然裂缝参数、层理缝参数;
S902、计算井筒内液体流动压降、射孔孔眼压降、携砂液柱静压力;将压裂施工井口压力折算为井底净压力;
S903、建立井底净压力数据与时间数据序列;对井底净压力进行动态拟合,并计算净压力指数平均值;
S904、基于井底净压力和净压力指数平均值,自动识别各段对应的裂缝延伸模式;
S905、数据输出:井底净压力值、净压力指数平均值、裂缝延伸模式;
S906、图像绘制:井底净压力曲线、净压力指数平均值曲线、裂缝延伸模式识别曲线,以及施工压力曲线诊断图。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910468722.2A CN110056336B (zh) | 2019-05-31 | 2019-05-31 | 一种页岩气缝网压裂施工压力曲线自动诊断方法 |
US16/752,651 US10689972B1 (en) | 2019-05-31 | 2020-01-26 | Automatic diagnosis method for wellhead pressure curve of hydraulic fracturing in shale gas horizontal well |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910468722.2A CN110056336B (zh) | 2019-05-31 | 2019-05-31 | 一种页岩气缝网压裂施工压力曲线自动诊断方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110056336A true CN110056336A (zh) | 2019-07-26 |
CN110056336B CN110056336B (zh) | 2020-01-07 |
Family
ID=67325190
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910468722.2A Active CN110056336B (zh) | 2019-05-31 | 2019-05-31 | 一种页岩气缝网压裂施工压力曲线自动诊断方法 |
Country Status (2)
Country | Link |
---|---|
US (1) | US10689972B1 (zh) |
CN (1) | CN110056336B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110905472A (zh) * | 2019-10-29 | 2020-03-24 | 中国石油集团川庆钻探工程有限公司 | 确定基于复合暂堵体系的实时转向压裂参数的方法 |
CN111636855A (zh) * | 2020-06-10 | 2020-09-08 | 华美孚泰油气增产技术服务有限责任公司 | 一种t型缝形成风险现场判别的方法及相应压裂工艺 |
CN112832723A (zh) * | 2019-11-22 | 2021-05-25 | 中国石油化工股份有限公司 | 一种气井用自适应控水装置及其设计方法 |
CN112983379A (zh) * | 2021-04-06 | 2021-06-18 | 华美孚泰油气增产技术服务有限责任公司 | 一种控制页岩气井均衡压裂改造的方法 |
CN113107447A (zh) * | 2021-04-14 | 2021-07-13 | 中煤科工集团重庆研究院有限公司 | 一种井下多孔段并行动态压裂系统及其施工方法 |
CN113435059A (zh) * | 2021-07-13 | 2021-09-24 | 中国石油大学(北京) | 模型建立方法、裂缝起裂事件诊断方法和装置 |
CN114444414A (zh) * | 2022-01-26 | 2022-05-06 | 北京科技大学 | 一种确定致密储层中多段压裂水平井最大裂缝间距的方法 |
Families Citing this family (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020097060A2 (en) | 2018-11-05 | 2020-05-14 | Schlumberger Technology Corporation | Fracturing operations pump fleet balance controller |
CN111764892B (zh) * | 2020-07-03 | 2023-05-12 | 中国石油天然气集团有限公司 | 页岩气井变压变产阶段生产制度适应性判别方法 |
CN111915171B (zh) * | 2020-07-21 | 2024-02-13 | 中国石油大学(华东) | 一种深度堵调效果综合评价方法 |
CN111927421B (zh) * | 2020-08-18 | 2022-07-15 | 中国石油化工股份有限公司 | 一种基于页岩气不稳定线性流理论的重复压裂选井方法 |
CN111852431B (zh) * | 2020-08-19 | 2022-04-12 | 西南石油大学 | 一种缝网压裂中缝网结构参数的优化方法和装置 |
CN114086946B (zh) * | 2020-08-24 | 2023-08-22 | 中国石油天然气股份有限公司 | 压裂缝的缝高确定方法 |
CN111963149B (zh) * | 2020-08-26 | 2024-03-01 | 中国石油天然气集团有限公司 | 一种考虑滞地液量增压的压裂后地层压力求取方法 |
CN114198095B (zh) * | 2020-08-31 | 2023-08-22 | 中国石油天然气股份有限公司 | 井筒失效数据采集分析方法及其装置、存储介质 |
US20220082727A1 (en) * | 2020-09-11 | 2022-03-17 | Saudi Arabian Oil Company | Method and system for determining energy-based brittleness |
CN112031759B (zh) * | 2020-09-17 | 2023-12-01 | 中国石油天然气股份有限公司 | 一种直井多分支缝压裂现场施工控制方法 |
CN112145166B (zh) * | 2020-10-27 | 2024-04-30 | 中国石油天然气集团有限公司 | 一种页岩气水平井压裂过程中井下情况识别预判方法 |
CN112377178B (zh) * | 2020-11-24 | 2024-04-30 | 中国石油天然气股份有限公司 | 一种边水水侵前缘的定量诊断方法 |
CN114592840B (zh) * | 2020-12-04 | 2023-10-27 | 中国石油天然气股份有限公司 | 暂堵压裂方法及其应用 |
CN114718513B (zh) * | 2021-01-05 | 2024-05-28 | 中国石油天然气股份有限公司 | 煤层气的产气量预估方法及装置 |
CN113605879B (zh) * | 2021-08-02 | 2023-03-21 | 中国石油大学(北京) | 煤储层原始地层压力的计算方法及装置 |
CN115707856A (zh) * | 2021-08-18 | 2023-02-21 | 中国石油天然气股份有限公司 | 效果参数确定方法、装置、设备以及介质 |
CN113657056A (zh) * | 2021-08-30 | 2021-11-16 | 中国石油大学(华东) | 降低井底沉砂与提高射孔孔眼进砂效率的直井压裂加砂优化设计方法 |
CN113969773A (zh) * | 2021-10-21 | 2022-01-25 | 中国石油化工股份有限公司 | 一种粗糙天然裂缝压裂液滤失测试方法 |
CN114580100B (zh) | 2022-02-22 | 2022-09-16 | 西南石油大学 | 压裂水平井全井筒压力计算方法、设备和计算机可读储存介质 |
CN114925626B (zh) * | 2022-04-25 | 2024-03-19 | 中国石油大学(北京) | 考虑孔眼动态磨蚀的多裂缝扩展预测方法及设备 |
CN115659736B (zh) * | 2022-10-19 | 2023-11-03 | 西南石油大学 | 一种深层页岩气水平井缝内转向压裂缝网扩展计算方法 |
CN116070455A (zh) * | 2023-02-22 | 2023-05-05 | 西南石油大学 | 一种控制水力压裂支撑剂回流的实时油嘴尺寸计算方法 |
CN116029167B (zh) * | 2023-03-31 | 2023-07-04 | 中国石油大学(华东) | 一种水平井多簇压裂非均匀射孔的定量设计方法 |
CN116629154B (zh) * | 2023-05-24 | 2024-01-09 | 西南石油大学 | 分形复合气藏压裂井瞬态压力计算方法、系统和设备 |
CN116455946B (zh) * | 2023-06-19 | 2023-09-19 | 安徽井上天华科技有限公司 | 一种基于云端高频井口压力生产数据分析方法 |
CN117633409B (zh) * | 2024-01-25 | 2024-04-09 | 中国科学院地质与地球物理研究所 | 计算页岩油气储层压裂缝网渗流参数的方法、系统及设备 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107524432A (zh) * | 2016-06-21 | 2017-12-29 | 中国石油化工股份有限公司 | 一种压裂施工中的实时调参方法 |
CN108019203A (zh) * | 2016-10-31 | 2018-05-11 | 中国石油化工股份有限公司 | 一种分层压裂裂缝沟通识别方法 |
WO2018200735A1 (en) * | 2017-04-25 | 2018-11-01 | Borehole Seismic, Llc. | Non-fracturing restimulation of unconventional hydrocarbon containing formations to enhance production |
CN108829945A (zh) * | 2018-05-29 | 2018-11-16 | 西南石油大学 | 一种基于压裂施工压力实时评估增产改造体积方法 |
CN108868748A (zh) * | 2018-04-28 | 2018-11-23 | 中国石油化工股份有限公司江汉油田分公司石油工程技术研究院 | 一种页岩气水平井重复压裂裂缝开启压力的计算方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2179637C1 (ru) * | 2001-05-08 | 2002-02-20 | Чикин Андрей Егорович | Способ определения характеристик скважины, призабойной зоны и пласта и устройство для его осуществления |
US9879507B2 (en) * | 2015-10-22 | 2018-01-30 | Dennis W. Gilstad | Adaptive stimulation system |
-
2019
- 2019-05-31 CN CN201910468722.2A patent/CN110056336B/zh active Active
-
2020
- 2020-01-26 US US16/752,651 patent/US10689972B1/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107524432A (zh) * | 2016-06-21 | 2017-12-29 | 中国石油化工股份有限公司 | 一种压裂施工中的实时调参方法 |
CN108019203A (zh) * | 2016-10-31 | 2018-05-11 | 中国石油化工股份有限公司 | 一种分层压裂裂缝沟通识别方法 |
WO2018200735A1 (en) * | 2017-04-25 | 2018-11-01 | Borehole Seismic, Llc. | Non-fracturing restimulation of unconventional hydrocarbon containing formations to enhance production |
CN108868748A (zh) * | 2018-04-28 | 2018-11-23 | 中国石油化工股份有限公司江汉油田分公司石油工程技术研究院 | 一种页岩气水平井重复压裂裂缝开启压力的计算方法 |
CN108829945A (zh) * | 2018-05-29 | 2018-11-16 | 西南石油大学 | 一种基于压裂施工压力实时评估增产改造体积方法 |
Non-Patent Citations (2)
Title |
---|
张平等: "《水力压裂净压力拟合分析解释技术研究与应用》", 《油气井测试》 * |
逄仁德: "《鄂尔多斯盆地陆相页岩气缝网压裂技术应用分析》", 《中国石油勘探》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110905472A (zh) * | 2019-10-29 | 2020-03-24 | 中国石油集团川庆钻探工程有限公司 | 确定基于复合暂堵体系的实时转向压裂参数的方法 |
CN110905472B (zh) * | 2019-10-29 | 2021-10-22 | 中国石油集团川庆钻探工程有限公司 | 确定基于复合暂堵体系的实时转向压裂参数的方法 |
CN112832723A (zh) * | 2019-11-22 | 2021-05-25 | 中国石油化工股份有限公司 | 一种气井用自适应控水装置及其设计方法 |
CN112832723B (zh) * | 2019-11-22 | 2022-12-09 | 中国石油化工股份有限公司 | 一种气井用自适应控水装置及其设计方法 |
CN111636855A (zh) * | 2020-06-10 | 2020-09-08 | 华美孚泰油气增产技术服务有限责任公司 | 一种t型缝形成风险现场判别的方法及相应压裂工艺 |
CN112983379A (zh) * | 2021-04-06 | 2021-06-18 | 华美孚泰油气增产技术服务有限责任公司 | 一种控制页岩气井均衡压裂改造的方法 |
CN113107447A (zh) * | 2021-04-14 | 2021-07-13 | 中煤科工集团重庆研究院有限公司 | 一种井下多孔段并行动态压裂系统及其施工方法 |
CN113435059A (zh) * | 2021-07-13 | 2021-09-24 | 中国石油大学(北京) | 模型建立方法、裂缝起裂事件诊断方法和装置 |
CN113435059B (zh) * | 2021-07-13 | 2022-09-27 | 中国石油大学(北京) | 模型建立方法、裂缝起裂事件诊断方法和装置 |
CN114444414A (zh) * | 2022-01-26 | 2022-05-06 | 北京科技大学 | 一种确定致密储层中多段压裂水平井最大裂缝间距的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110056336B (zh) | 2020-01-07 |
US10689972B1 (en) | 2020-06-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110056336A (zh) | 一种页岩气缝网压裂施工压力曲线自动诊断方法 | |
CN109209337B (zh) | 一种考虑岩屑床的水平井钻井润滑性实验装置及实验方法 | |
CN103510893B (zh) | 一种利用流量监控实现井底压力控制的钻井装备与方法 | |
CN104790916B (zh) | 一种利用油套压力平衡法排除气井积液方法 | |
CN105672972B (zh) | 一种多级水力喷射压裂管柱性能评价方法 | |
CN202755927U (zh) | 一种利用流量监控实现井底压力控制的钻井装备 | |
CN108612508A (zh) | 一种利用高强度可降解暂堵球的页岩气水平井压裂新方法 | |
CN108468538A (zh) | 一种页岩水压裂缝扩展预测方法 | |
CN103899581B (zh) | 一种反循环生产、反循环起水力喷射泵泵芯方法及装置 | |
CN103527148B (zh) | 一种双套管单泵分层控压合层排采技术及装置 | |
CN105626036A (zh) | 一种确定油藏合理产液量油藏工程计算方法 | |
CN107843513A (zh) | 水合物开采井机械筛管冲蚀评价仿真系统及其仿真测试方法 | |
CN108894777A (zh) | 一种分层压裂多层合采油气藏储层物性及裂缝特性参数的确定方法 | |
CN110118069A (zh) | 一种超深井钻井压力控制设备及操作方法 | |
CN110761743A (zh) | 一种稠油掺稀气举工艺管柱及实施方法 | |
CN105715220B (zh) | 一种钻井控压系统 | |
CN111639434B (zh) | 一种页岩气水平井投球暂堵转向压裂的优化方法 | |
CN110259421A (zh) | 一种裂缝性的致密油藏注水补充能量方法 | |
CN204532261U (zh) | 一种可控型可溶球座多级压裂滑套 | |
CN109973059B (zh) | 一种缝洞型油藏注气体机理的快速分析方法 | |
CN113743037B (zh) | 一种低渗透油藏注水诱导动态裂缝变导流能力计算方法 | |
CN114135265B (zh) | 一种海上油田低渗储层低成本高效改造工艺方法 | |
CN109522619B (zh) | 一种低渗油田周期采油方法 | |
CN109931038B (zh) | 一种缝洞型油藏注氮气设计方法 | |
CN107339096A (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 | ||
OL01 | Intention to license declared | ||
OL01 | Intention to license declared |