CN112084181A - 一种压力恢复测试早期数据校正及复原方法 - Google Patents
一种压力恢复测试早期数据校正及复原方法 Download PDFInfo
- Publication number
- CN112084181A CN112084181A CN202010942476.2A CN202010942476A CN112084181A CN 112084181 A CN112084181 A CN 112084181A CN 202010942476 A CN202010942476 A CN 202010942476A CN 112084181 A CN112084181 A CN 112084181A
- Authority
- CN
- China
- Prior art keywords
- data
- early
- pressure
- slope
- pressure recovery
- 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.)
- Withdrawn
Links
- 238000011084 recovery Methods 0.000 title claims abstract description 60
- 238000012360 testing method Methods 0.000 title claims abstract description 49
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000004140 cleaning Methods 0.000 claims abstract description 12
- 238000009499 grossing Methods 0.000 claims abstract description 12
- 230000002159 abnormal effect Effects 0.000 claims abstract description 10
- 238000001914 filtration Methods 0.000 claims description 4
- 230000000630 rising effect Effects 0.000 claims description 4
- 238000010586 diagram Methods 0.000 description 7
- 239000003129 oil well Substances 0.000 description 6
- 238000004519 manufacturing process Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
Images
Classifications
-
- 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/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- 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
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/21—Design, administration or maintenance of databases
- G06F16/215—Improving data quality; Data cleansing, e.g. de-duplication, removing invalid entries or correcting typographical errors
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/24—Querying
- G06F16/245—Query processing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Databases & Information Systems (AREA)
- Data Mining & Analysis (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mathematical Physics (AREA)
- General Engineering & Computer Science (AREA)
- Geophysics (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Geochemistry & Mineralogy (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Computational Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Mathematical Optimization (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Computational Linguistics (AREA)
- Algebra (AREA)
- Software Systems (AREA)
- Complex Calculations (AREA)
Abstract
本发明属于数据处理领域,具体涉及一种压力恢复测试早期数据校正及复原方法。包括以下步骤:①数据清洗,清洗早期异常形态数据,包括上翘、下凹等形态;②数据光滑,对于整个压力恢复曲线,由于地层、机器及人为等各种因素导致的噪音数据进行光滑,以获得光滑的数据;③数据复原,采用牛顿插值算法,复原压力恢复曲线中缺失的早期数据。本方法可以应对复杂的压力恢复测试数据,对早期缺失的数据进行复原,而且设计的数据校正及复原方法简单、有效。
Description
技术领域
本发明属于数据处理领域,具体涉及一种压力恢复测试早期数据校正及复原方法。
背景技术
试井方法是通过油气藏压力动态反应来研究油气藏特性的方法。试井方法是通过给油气藏一系列“信号”,这个信号通常是改变井的工作制度而产生的,测试在给定信号作用下,油气藏的压力动态反应。试井方法主要可分为两大类,一类是不稳定试井方法,连续测试由于井的工作制度发生变化而引起的井底压力随时间的变化过程,通过这一压力变化过程的特征来研究井和储层的特性参数的试井方法,主要包括压力恢复试井、压力降落试井、变流量试井、干扰试井和脉冲试井等;另一类是稳定试井,也叫产能试井,通过测试井在几个不同稳定工作制度下,井底压力与产量之间的关系来研究井的生产(注入)能力的试井方法,主要包括回压试井、等时试井和一点法试井等。不稳定试井方法中,压力恢复试井是最常用、最有效的方法。
压力恢复测试指油井生产稳定后下存储式压力计到井内预定位置,关井录取井底压力随时间的变化,反映了油井在改变工作制度后油层压力重新分布的过程。
压力恢复测试的方法是先将油井或气井以不变的产量生产一定时间以后,关井使其井底压力不断恢复,并且用井底压力计连续地记录关井时间与不断上升的井底压力。部分低压低渗油井压力恢复测试由于受矿场条件、测试成本等限制,经常采用井口关井的方式,而在关井之后到将压力计放入井底预定位置期间的压力恢复数据是无法获得的,因此造成早期的压力恢复数据数据缺失,压力恢复曲线形态异常,无法准确解释相关地层参数,即无法通过恢复曲线来研究油气藏特性,影响矿场的生产及措施调整。
发明内容
本发明要解决的技术问题是提供一种压力恢复测试早期数据校正及复原方法,其一方面为了解决现有压力恢复测试数据噪点多,数据曲线异常问题;另一方面为了解决压力恢复测试早期数据缺失导致的试井解释难以进行的问题。
一种压力恢复测试早期数据校正及复原方法,其特征在于:包括以下步骤:
①数据清洗,清洗早期异常形态数据,包括上翘、下凹等形态;
②数据光滑,对于整个压力恢复曲线,由于地层、机器及人为等各种因素导致的噪音数据进行光滑,以获得光滑的数据;
③数据复原,采用牛顿插值算法,复原压力恢复曲线中缺失的早期数据。
本发明具体操作步骤为:
①首先,通过利用早期压力数据点间的斜率进行数据清洗,其斜率存在极大值,且该极大值远超出平均斜率,判断其为不合理的早期压力数据;获得该极大值点的位置,并剔除前面的压力数据;
其次,还存在斜率较大的数据,通过分析测试数据的压力数据点间的斜率,其早期数据斜率存在极小值,获得该极小值点的位置,并剔除前面的压力数据;
最后,若存在压力数据存在先下降后上升的异常情况,其斜率为负,因此剔除早期斜率为负的数据;并且结合早期数据存在非第一个位置的极小值,对其进行处理,若不存在上述情况,则不对其处理;
②采用均值滤波算法进行平滑处理,即在压力数据中,对待处理的数据点给定一个大小的窗口,以窗口内所有数据的均值来代替原来的数据的方法;
其中,n为窗口大小;
③设压力恢复函数p(t)在[ta,tb]上有定义,对于n+1个节点(tj,pj),j=0,1,…,n,其中tj∈[ta,tb]互不相同,满足p(tj)=pj,j=0,1,…,n;对任意t∈[ta,tb],由一阶差商的定义得
p(t)=p(t0)+p[t,t0](t-t0)
类似地,由二阶差商至n阶差商的定义得到下列方程
p(t,t0)=p(t0,t1)+p[t,t0,t1](t-t1)
p(t,t0,t1)=p(t0,t1,t2)+p[t,t0,t1,t2](t-t2)
………
p(t,t0,t1,L,tn-1)=p(t0,t1,t2,L,tn)+p[t,t0,t1,t2,L,tn](t-tn)
整理得:
p(t)=p(t0)+p[t0,t1](t-t0)+p[t0,t1,t2](t-t0)(t-t1)+L+
p[t0,t1,t2,L,tn](t-t0)(t-t1)L(t-tn-1)+
p[t,t0,t1,t2,L,tn](t-t0)(t-x1)L(t-tn-1)(t-tn)
即为带余项的牛顿插值公式;
αn(t)=p(t0)+p[t0,t1](t-t0)+p[t0,t1,t2](t-t0)(t-t1)+L+
令p[t0,t1,t2,L,tn](t-t0)(t-t1)L(t-tn-1)
βn(t)=p[t,t0,t1,t2,L,tn](t-t0)(t-t1)L(t-tn-1)(t-tn)
则可以表示为
p(t)=αn(t)+βn(t)
其中,αn(t)为次数不高于n次的t的多项式,可以验证αn(ti)=p(ti)(i=0,1,Ln),αn(t)称为过n+1个节点的αn(t)的n阶牛顿插值多项式,βn(t)称为牛顿插值多项式的余项。
本发明的有益效果:
结合数据清洗、数据平滑和数据插值等步骤,解决了压力恢复测试数据噪点多,数据曲线异常问题,并将压力恢复测试的早期缺失数据进行了插值构建和复原,为提高压力恢复曲线解释结果的准确性,有助于指导矿场的生产及措施调整。
附图说明
图1为本发明的一种压力恢复测试早期数据校正及复原方法的一优选实施例的步骤①中的第一步处理前的压力恢复数据及斜率示意图。
图2为本发明的一种压力恢复测试早期数据校正及复原方法的一优选实施例的步骤①中的第一步处理后的压力恢复数据及斜率示意图。
图3为本发明的一种压力恢复测试早期数据校正及复原方法的一优选实施例的步骤①中的第二步处理前的压力恢复数据及斜率示意图。
图4为本发明的一种压力恢复测试早期数据校正及复原方法的一优选实施例的步骤①中的第二步处理后的压力恢复数据及斜率示意图。
图5为本发明的一种压力恢复测试早期数据校正及复原方法的一优选实施例的步骤①中的第三步处理后的压力恢复数据及斜率示意图。
图6为本发明的一种压力恢复测试早期数据校正及复原方法的一优选实施例的均值滤波示意图。
图7为本发明的一种压力恢复测试早期数据校正及复原方法的一优选实施例的压力恢复测试原始数据及早期数据复原结果对比示意图。
具体实施方式
下面结合附图对本发明作进一步的说明。
一种压力恢复测试早期数据校正及复原方法,包括以下步骤:
①数据清洗,清洗早期异常形态数据,包括上翘、下凹等形态;
首先,通过利用早期压力数据点间的斜率进行数据清洗,其斜率存在极大值,且该极大值远超出平均斜率,判断其为不合理的早期压力数据,如图1所示;获得该极大值点的位置,并剔除前面的压力数据,结果如图2所示;
其次,还存在斜率较大的数据,通过分析测试数据的压力数据点间的斜率,其早期数据斜率存在极小值,如图3所示,获得该极小值点的位置,并剔除前面的压力数据,结果如图4所示;
最后,若存在压力数据存在先下降后上升的异常情况,其斜率为负,因此剔除早期斜率为负的数据;并且结合早期数据存在非第一个位置的极小值,对其进行处理。若不存在上述情况,则不对其处理,如图5所示。
②数据光滑,对于整个压力恢复曲线,由于地层、机器及人为等各种因素导致的噪音数据进行光滑,以获得光滑的数据;
采用均值滤波算法进行平滑处理,即在压力数据中,对待处理的数据点给定一个大小的窗口,以窗口内所有数据的均值来代替原来的数据的方法,如图6所示。
其中,n为窗口大小。
③数据复原,采用牛顿插值算法,复原压力恢复曲线中缺失的早期数据,结果如图7所示。
设压力恢复函数p(t)在[ta,tb]上有定义,对于n+1个节点(tj,pj),j=0,1,…,n,其中tj∈[ta,tb]互不相同,满足p(tj)=pj,j=0,1,…,n。对任意t∈[ta,tb],由一阶差商的定义得
p(t)=p(t0)+p[t,t0](t-t0)
类似地,由二阶差商至n阶差商的定义得到下列方程
p(t,t0)=p(t0,t1)+p[t,t0,t1](t-t1)
p(t,t0,t1)=p(t0,t1,t2)+p[t,t0,t1,t2](t-t2)
………
p(t,t0,t1,L,tn-1)=p(t0,t1,t2,L,tn)+p[t,t0,t1,t2,L,tn](t-tn)
整理得:
p(t)=p(t0)+p[t0,t1](t-t0)+p[t0,t1,t2](t-t0)(t-t1)+L+
p[t0,t1,t2,L,tn](t-t0)(t-t1)L(t-tn-1)+
p[t,t0,t1,t2,L,tn](t-t0)(t-x1)L(t-tn-1)(t-tn)
即为带余项的牛顿插值公式。
αn(t)=p(t0)+p[t0,t1](t-t0)+p[t0,t1,t2](t-t0)(t-t1)+L+
令p[t0,t1,t2,L,tn](t-t0)(t-t1)L(t-tn-1)
βn(t)=p[t,t0,t1,t2,L,tn](t-t0)(t-t1)L(t-tn-1)(t-tn)
则可以表示为
p(t)=αn(t)+βn(t)
其中,αn(t)为次数不高于n次的t的多项式,可以验证αn(ti)=p(ti)(i=0,1,Ln)。αn(t)称为过n+1个节点的αn(t)的n阶牛顿插值多项式,βn(t)称为牛顿插值多项式的余项。
以上对本发明及其实施方式进行了描述,这种描述没有限制性,附图中所示的也只是本发明的实施方式之一,实际的结构并不局限于此。总而言之如果本领域的普通技术人员受其启示,在不脱离本发明创造宗旨的情况下,不经创造性的设计出与该技术方案相似的结构方式及实施例,均应属于本发明的保护范围。
Claims (2)
1.一种压力恢复测试早期数据校正及复原方法,其特征在于:包括以下步骤:
①数据清洗,清洗早期异常形态数据,包括上翘、下凹等形态;
②数据光滑,对于整个压力恢复曲线,由于地层、机器及人为等各种因素导致的噪音数据进行光滑,以获得光滑的数据;
③数据复原,采用牛顿插值算法,复原压力恢复曲线中缺失的早期数据。
2.根据权利要求1所述的一种压力恢复测试早期数据校正及复原方法,其特征在于:具体操作步骤为:
①首先,通过利用早期压力数据点间的斜率进行数据清洗,其斜率存在极大值,且该极大值远超出平均斜率,判断其为不合理的早期压力数据;获得该极大值点的位置,并剔除前面的压力数据;
其次,还存在斜率较大的数据,通过分析测试数据的压力数据点间的斜率,其早期数据斜率存在极小值,获得该极小值点的位置,并剔除前面的压力数据;
最后,若存在压力数据存在先下降后上升的异常情况,其斜率为负,因此剔除早期斜率为负的数据;并且结合早期数据存在非第一个位置的极小值,对其进行处理,若不存在上述情况,则不对其处理;
②采用均值滤波算法进行平滑处理,即在压力数据中,对待处理的数据点给定一个大小的窗口,以窗口内所有数据的均值来代替原来的数据的方法;
其中,n为窗口大小;
③设压力恢复函数p(t)在[ta,tb]上有定义,对于n+1个节点(tj,pj),j=0,1,…,n,其中tj∈[ta,tb]互不相同,满足p(tj)=pj,j=0,1,…,n;对任意t∈[ta,tb],由一阶差商的定义得
p(t)=p(t0)+p[t,t0](t-t0)
类似地,由二阶差商至n阶差商的定义得到下列方程
p(t,t0)=p(t0,t1)+p[t,t0,t1](t-t1)
p(t,t0,t1)=p(t0,t1,t2)+p[t,t0,t1,t2](t-t2)
………
p(t,t0,t1,L,tn-1)=p(t0,t1,t2,L,tn)+p[t,t0,t1,t2,L,tn](t-tn)
整理得:
p(t)=p(t0)+p[t0,t1](t-t0)+p[t0,t1,t2](t-t0)(t-t1)+L+p[t0,t1,t2,L,tn](t-t0)(t-t1)L(t-tn-1)+p[t,t0,t1,t2,L,tn](t-t0)(t-x1)L(t-tn-1)(t-tn)
即为带余项的牛顿插值公式;
αn(t)=p(t0)+p[t0,t1](t-t0)+p[t0,t1,t2](t-t0)(t-t1)+L+
令p[t0,t1,t2,L,tn](t-t0)(t-t1)L(t-tn-1)
βn(t)=p[t,t0,t1,t2,L,tn](t-t0)(t-t1)L(t-tn-1)(t-tn)
则可以表示为
p(t)=αn(t)+βn(t)
其中,αn(t)为次数不高于n次的t的多项式,可以验证αn(ti)=p(ti)(i=0,1,Ln),αn(t)称为过n+1个节点的αn(t)的n阶牛顿插值多项式,βn(t)称为牛顿插值多项式的余项。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010942476.2A CN112084181A (zh) | 2020-09-09 | 2020-09-09 | 一种压力恢复测试早期数据校正及复原方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010942476.2A CN112084181A (zh) | 2020-09-09 | 2020-09-09 | 一种压力恢复测试早期数据校正及复原方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112084181A true CN112084181A (zh) | 2020-12-15 |
Family
ID=73732972
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010942476.2A Withdrawn CN112084181A (zh) | 2020-09-09 | 2020-09-09 | 一种压力恢复测试早期数据校正及复原方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112084181A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113027532A (zh) * | 2021-03-31 | 2021-06-25 | 山东思科赛德矿业安全工程有限公司 | 一种基于综采支架工作阻力数据挖掘来压信息的方法 |
CN114704247A (zh) * | 2022-06-07 | 2022-07-05 | 西安钰雪石油科技有限公司 | 油气井液面智能监测方法及装置 |
-
2020
- 2020-09-09 CN CN202010942476.2A patent/CN112084181A/zh not_active Withdrawn
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113027532A (zh) * | 2021-03-31 | 2021-06-25 | 山东思科赛德矿业安全工程有限公司 | 一种基于综采支架工作阻力数据挖掘来压信息的方法 |
CN113027532B (zh) * | 2021-03-31 | 2023-01-24 | 山东思科赛德矿业安全工程有限公司 | 一种基于综采支架工作阻力数据挖掘来压信息的方法 |
CN114704247A (zh) * | 2022-06-07 | 2022-07-05 | 西安钰雪石油科技有限公司 | 油气井液面智能监测方法及装置 |
CN114704247B (zh) * | 2022-06-07 | 2022-08-16 | 西安钰雪石油科技有限公司 | 油气井液面智能监测方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112084181A (zh) | 一种压力恢复测试早期数据校正及复原方法 | |
CN110321595B (zh) | 一种测井提取静态品质系数的断层封闭性评价方法 | |
CN109057786B (zh) | 用于非常规油气藏可采储量估算的自适应模型的确定方法 | |
CN115345378B (zh) | 一种基于机器学习的页岩气井产量评价方法 | |
CN108614902B (zh) | 一种生产制度变化的页岩气井生产数据分析方法 | |
Morris et al. | Exponentiation Bias in Production and Reserves Estimates Made by Use of Exponential-Decline Models | |
CN109187744B (zh) | 一种基于静力触探试验的剪切波速估算方法 | |
CN116644284A (zh) | 一种地层分类特征因子确定方法、系统、电子设备及介质 | |
CN115030708A (zh) | 一种基于不关井干扰试井的干扰井识别系统及方法 | |
CN114021620B (zh) | 基于bp神经网络特征提取的电潜泵故障诊断方法 | |
CN109184661B (zh) | 用于底水油藏水平井高产液位置识别的监测方法及其系统 | |
CN113762614A (zh) | 一种基于静态油气田数据的页岩气井产能预测的方法 | |
CN112324419B (zh) | 一种基于压降曲线分析的裂缝参数反演及评价方法 | |
CN112613171B (zh) | 一种煤储层渗透率确定方法及装置 | |
CN111338310B (zh) | 一种工业过程稳态工况识别与分类方法 | |
WO2017024533A1 (zh) | 一种评价致密气藏储层应力敏感性方法 | |
CN110991363B (zh) | 不同采煤工序下煤矿安全监控系统co涌出特征提取方法 | |
CN112145166B (zh) | 一种页岩气水平井压裂过程中井下情况识别预判方法 | |
CN116028865A (zh) | 基于特征扩增的航空发动机故障诊断方法 | |
CN109948291A (zh) | 一种砂体不连续界线方向自适应识别方法 | |
CN114060018B (zh) | 一种储层动态储量确定方法、系统、设备及可读存储介质 | |
Mohamed et al. | State of the Art in Characterization of Frac Stage Geometry, Skin and Conductivity Using Pressure Leakoff | |
CN111985146B (zh) | 一种措施潜力井的确定方法和确定装置 | |
CN112001055B (zh) | 一种基于微构造的低幅稀油油藏含水率预测方法 | |
CN112360449B (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 | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20201215 |
|
WW01 | Invention patent application withdrawn after publication |