CN110751405B - 一种水驱特征曲线快速拟合方法及系统 - Google Patents
一种水驱特征曲线快速拟合方法及系统 Download PDFInfo
- Publication number
- CN110751405B CN110751405B CN201911019328.7A CN201911019328A CN110751405B CN 110751405 B CN110751405 B CN 110751405B CN 201911019328 A CN201911019328 A CN 201911019328A CN 110751405 B CN110751405 B CN 110751405B
- Authority
- CN
- China
- Prior art keywords
- coordinate
- water
- characteristic curve
- line segment
- curve
- 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 166
- 238000000034 method Methods 0.000 title claims abstract description 34
- 230000008859 change Effects 0.000 claims abstract description 47
- 230000014509 gene expression Effects 0.000 claims abstract description 33
- 238000004519 manufacturing process Methods 0.000 claims abstract description 15
- 238000000605 extraction Methods 0.000 claims abstract description 12
- 238000012417 linear regression Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000005065 mining Methods 0.000 claims description 3
- 238000012937 correction Methods 0.000 description 8
- 230000001186 cumulative effect Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000007796 conventional method Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- -1 hydrogen Chemical class 0.000 description 2
- 239000001257 hydrogen Substances 0.000 description 2
- 108010059993 Vancomycin Proteins 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 125000001421 myristyl group Chemical group [H]C([*])([H])C([H])([H])C([H])([H])C([H])([H])C([H])([H])C([H])([H])C([H])([H])C([H])([H])C([H])([H])C([H])([H])C([H])([H])C([H])([H])C([H])([H])C([H])([H])[H] 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- MYPYJXKWCTUITO-LYRMYLQWSA-N vancomycin Chemical compound O([C@@H]1[C@@H](O)[C@H](O)[C@@H](CO)O[C@H]1OC1=C2C=C3C=C1OC1=CC=C(C=C1Cl)[C@@H](O)[C@H](C(N[C@@H](CC(N)=O)C(=O)N[C@H]3C(=O)N[C@H]1C(=O)N[C@H](C(N[C@@H](C3=CC(O)=CC(O)=C3C=3C(O)=CC=C1C=3)C(O)=O)=O)[C@H](O)C1=CC=C(C(=C1)Cl)O2)=O)NC(=O)[C@@H](CC(C)C)NC)[C@H]1C[C@](C)(N)[C@H](O)[C@H](C)O1 MYPYJXKWCTUITO-LYRMYLQWSA-N 0.000 description 1
- 229960003165 vancomycin Drugs 0.000 description 1
- MYPYJXKWCTUITO-UHFFFAOYSA-N vancomycin Natural products O1C(C(=C2)Cl)=CC=C2C(O)C(C(NC(C2=CC(O)=CC(O)=C2C=2C(O)=CC=C3C=2)C(O)=O)=O)NC(=O)C3NC(=O)C2NC(=O)C(CC(N)=O)NC(=O)C(NC(=O)C(CC(C)C)NC)C(O)C(C=C3Cl)=CC=C3OC3=CC2=CC1=C3OC1OC(CO)C(O)C(O)C1OC1CC(C)(N)C(O)C(C)O1 MYPYJXKWCTUITO-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06315—Needs-based resource requirements planning or analysis
-
- 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
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- Data Mining & Analysis (AREA)
- Strategic Management (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Operations Research (AREA)
- Tourism & Hospitality (AREA)
- Pure & Applied Mathematics (AREA)
- Marketing (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Life Sciences & Earth Sciences (AREA)
- Agronomy & Crop Science (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- General Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Software Systems (AREA)
- General Health & Medical Sciences (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Primary Health Care (AREA)
- Probability & Statistics with Applications (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Game Theory and Decision Science (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供一种水驱特征曲线快速拟合方法及系统，该方法包括：根据水驱特征曲线类型，分别确定各水驱特征曲线待定参数；将待定参数代入含水率变化关系式，计算采出程度和含水率变化曲线；根据实际生产数据，选取适合描述油田含水变化规律的对应类型的水驱特征曲线；设定第一含水率值对应的坐标为拟合线段的初始坐标，设定第二含水率值对应的坐标为末端坐标，分别确定对应类型的水驱特征曲线的初始坐标和末端坐标的表达式；将待定参数估值代入初始坐标和末端坐标的表达式中，确定拟合线段的初始点和末端点并绘制拟合线段。通过该方案解决了工程软件制图效率低的问题，实现对水驱特征曲线的快速拟合，方便水驱特征曲线的绘制。
Description
技术领域
本发明涉及油藏工程领域，尤其涉及一种水驱特征曲线快速拟合方法及系统。
背景技术
水驱特征曲线是注水开发油田效果分析、水驱可采储量标定中最常用的方法之一。随着水驱油理论的逐步创新与完善，已形成近80种水驱特征曲线，其中，在矿场上应用较多的有马克西莫夫-童宪章曲线(简称甲型)、沙卓洛夫曲线(简称乙型)、西帕切夫曲线(简称丙型)、纳扎洛夫曲线(简称丁型)等12种，这些方法常见于油藏工程软件中，用于拟合水驱特征曲线。
一般在水驱特征曲线拟合过程中，利用算法编程可以快速对待定参数进行辨识，但具体确定拟合线段时，通常采用给定曲线某一纵坐标值，再通过水驱特征曲线预估模型计算横坐标的方式，对于单个曲线的绘制，常规方法较为适用，而具体应用到工程软件中，常规方法制图效率较低，使用较为不便。
发明内容
有鉴于此，本发明实施例提供了一种水驱特征曲线快速拟合方法及系统，以解决工程软件上水驱特征曲线拟合效率低的问题。
在本发明实施例的第一方面，提供了一种水驱特征曲线快速拟合方法，包括：
根据水驱特征曲线类型，分别确定各水驱特征曲线待定参数；
将所述待定参数代入含水率变化关系式，计算采出程度和含水率变化曲线；
根据实际生产数据，选取适合描述油田含水变化规律的对应类型的水驱特征曲线；
设定第一含水率值对应的坐标为拟合线段的初始坐标，设定第二含水率值对应的坐标为末端坐标，分别确定对应类型的水驱特征曲线的初始坐标和末端坐标的表达式；
将待定参数估值代入初始坐标和末端坐标的表达式中，确定拟合线段的初始点和末端点并绘制拟合线段。
在本发明实施例的第一方面，提供了一种水驱特征曲线快速拟合系统，包括：
确定模块：用于根据水驱特征曲线类型，分别确定各水驱特征曲线待定参数；
计算模块：用于将所述待定参数代入含水率变化关系式，计算采出程度和含水率变化曲线；
选取模块：用于根据实际生产数据，选取适合描述油田含水变化规律的对应类型的水驱特征曲线；
设定模块：用于设定第一含水率值对应的坐标为拟合线段的初始坐标，设定第二含水率值对应的坐标为末端坐标，分别确定对应类型的水驱特征曲线的初始坐标和末端坐标的表达式；
拟合模块：用于将待定参数估值代入初始坐标和末端坐标的表达式中，确定拟合线段的初始点和末端点并绘制拟合线段。
本发明实施例中，确定不同类型水驱特征曲线的待定参数，将待定参数代入含水率变化关系式，计算采出程度和含水率变化曲线，再根据实际生产数据，选取适合描述油田含水变化规律的对应类型的水驱特征曲线，分别确定对应类型的水驱特征曲线的初始坐标和末端坐标的表达式；将待定参数估值代入初始坐标和末端坐标的表达式中，确定拟合线段的初始点和末端点并绘制拟合线段。解决了传统水驱特征曲线拟合效率低的问题，可以借助工程软件，及水驱特征曲线中首末坐标表达式，实现对含水率区间内不同类型特征曲线的快速拟合，方便水驱特征曲线在矿场上的应用。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单介绍，显而易见地，下面描述的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获取其他附图。
图1为本发明的一个实施例提供的水驱特征曲线快速拟合方法的流程示意图；
图2为本发明的一个实施例提供的水驱特征曲线快速拟合系统的结构示意图。
具体实施方式
为使得本发明的发明目的、特征、优点能够更加的明显和易懂，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，下面所描述的实施例仅仅是本发明一部分实施例，而非全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围以下结合附图对本发明的原理和特征进行描述，所举实例只用于解释本发明，并非用于限定本发明的范围。
本发明的说明书或权利要求书及上述附图中的术语“包括”以及其他相近意思表述，意指覆盖不排他的包含，如包含一系列步骤或单元的过程、方法或系统、设备没有限定于已列出的步骤或单元。已列出的步骤或单元。此外，“第一”“第二”用于区分不同对象，并非用于描述特定顺序。
请参阅图1，本发明实施例提供的一种水驱特征曲线快速拟合方法的流程示意图，包括：
S101、根据水驱特征曲线类型，分别确定各水驱特征曲线待定参数；
对于众多不同类型的水驱特征曲线，可以直接转化为多元线性回归求解，如不带修正参数的甲型、乙型、丙型、丁型等水驱特征曲线，直接利用多元线性求得待定参数。
示例性的，如马克西莫夫-童宪章曲线(甲型)：
水驱特征曲线：Np＝a+blnWp；
对应含水率变化曲线：
令y＝Np，x＝lnWp，则未校正甲型水驱特征曲线转化为y＝a+bx。
如沙卓洛夫曲线(乙型)：
水驱特征曲线：Np＝a+blnLp；
对应含水率变化曲线：
令y＝Np，x＝lnLp，则未校正乙型水驱特征曲线转化为y＝a+bx。
如西帕切夫曲线(丙型)：
水驱特征曲线：
对应含水率变化曲线：
令y＝Lp/Np，x＝Lp，则未校正丙型水驱特征曲线转化为y＝a+bx。
如纳扎洛夫曲线(丁型)：
水驱特征曲线：
对应含水率变化曲线：
令y＝Lp/Np，x＝Wp，则未校正丁型水驱特征曲线转化为y＝a+bx。
如广义西帕切夫曲线(校正丙型)：
水驱特征曲线：
对应含水率变化曲线：
令y＝Lp/Np，x1＝Lp，x2＝1/Np，a0＝a，a1＝b，a2＝-C，则校正丙型水驱特征曲线转化为y＝a0+a1x1+a2x2。
如广义纳扎洛夫曲线(校正丁型)：
水驱特征曲线：
对应含水率变化曲线：
令y＝Lp/Np，x1＝Wp，x2＝1/Np，a0＝a，a1＝b，a2＝-C，则校正丁型水驱特征曲线转化为y＝a0+a1x1+a2x2。
对于一些其他类型的水驱特征曲线，还可以直接转化为含单变量的最小二乘法求解，如带修正参数的甲型、乙型、丙型、丁型等水驱特征曲线，借助试凑法或单变量求解法确定待定参数。
示例性的，如马克西莫夫-童宪章曲线(校正甲型)：
水驱特征曲线：Np＝a+bln(Wp+C)；
对应含水率变化曲线：
令y＝Np，x＝ln(Wp+C)，则校正甲型水驱特征曲线转化为含单变量C的y＝a+bx。
如沙卓洛夫曲线(校正乙型)：
水驱特征曲线：Np＝a+bln(Lp+C)；
对应含水率变化曲线：
令y＝Np，x＝ln(Lp+C)，则校正乙型水驱特征曲线转化为含单变量C的y＝a+bx。
如万吉业S～凸型曲线：
水驱特征曲线：ln(1-Np/N)＝a-bln(Lp+C)；
对应含水率变化曲线：
令y＝ln(1-Np/N)，x＝ln(Lp+C)，则万吉业S～凸型水驱特征曲线转化为含单变量C的y＝a+bx。
如布雷吉曲线(双对数曲线)：
水驱特征曲线：lnNp＝a+bln(Wp+C)；
对应含水率变化曲线：
令y＝lnNp，x＝ln(Wp+C)，则双对数型水驱特征曲线转化为含单变量C的y＝a+bx。
如万吉业超凸型：
水驱特征曲线：lnNp＝a+bln(Lp+C)；
对应含水率变化曲线：
令y＝lnNp，x＝ln(Lp+C)，则万吉业超凸型水驱特征曲线转化为含单变量C的y＝a+bx。
如陈元千曲线：
水驱特征曲线：
对应含水率变化曲线：
令y＝Np，
则陈元千型水驱特征曲线转化为含单变量C的y＝a+bx。
利用多元线性回归法或含单变量的最小二乘法可以依次确定待定参数为
或
S102、将所述待定参数代入含水率变化关系式，计算采出程度和含水率变化曲线；
将待定参数
或
代入含水变化关系式，可以得到不同类型的含水率变化曲线预估模型。
示例性的，甲型或校正甲型水驱特征曲线预估模型和对应含水率变化曲线预估模型分别为：
利用R＝Np/N，代入(2)式，可得
将(2)式代入(1)式，得
可以理解的是，R和fw分别对应采出程度和含水率变化曲线。
S103、根据实际生产数据，选取适合描述油田含水变化规律的对应类型的水驱特征曲线；
根据历史油藏开采数据，得到常用水驱特征曲线的拟合结果，其中包括待定参数估值、相关系数和适应性等信息。
S104、设定第一含水率值对应的坐标为拟合线段的初始坐标，设定第二含水率值对应的坐标为末端坐标，分别确定对应类型的水驱特征曲线的初始坐标和末端坐标的表达式；
由水驱特征曲线特点可知，横轴为采出程度，纵轴为含水率，则拟合的曲线线段坐标通式可由采出程度函数和含水率函数表示。
示例性的，甲型或校正甲型的横轴为Np，纵轴为lnWp(或ln(Wp+C))，则拟合线段通式为：
在实际中，一般油藏含水小于2％时为无水期开发，含水率为98％时，对应累积产油量为油藏水驱可采储量。由此，可定义含水率为2％时，对应坐标为拟合线段初端坐标，含水率为98％时，对应坐标为拟合线段的末端坐标。
则甲型或校正甲型曲线拟合线段初端坐标为
甲型或校正甲型曲线拟合线段末端坐标为：
进一步的，沙卓洛夫曲线(乙型或校正乙型)拟合线段各点坐标通式为：
则乙型或校正乙型曲线拟合线段初端坐标为：
乙型或校正乙型曲线拟合线段末端坐标为：
西帕切夫曲线(丙型)拟合线段各点坐标通式为：
则丙型曲线拟合线段初端坐标为：
丙型曲线拟合线段末端坐标为：
纳扎洛夫曲线(丁型或校正丁型)拟合线段各点坐标通式为：
则纳扎洛夫曲线曲线拟合线段初端坐标为：
纳扎洛夫曲线曲线拟合线段末端坐标为：
校正丙型拟合线段各点坐标通式为：
则校正丙型拟合线段初端坐标为：
校正丙型拟合线段末端坐标为：
万吉业S～凸型曲线拟合线段各点坐标通式为：
则万吉业S～凸型曲线拟合线段初端坐标为：
万吉业S～凸型曲线拟合线段末端坐标为：
布雷吉曲线(双对数曲线)拟合线段各点坐标通式为：
则布雷吉曲线(双对数曲线)拟合线段初端坐标为：
布雷吉曲线(双对数曲线)拟合线段末端坐标为：
万吉业超凸型拟合线段各点坐标通式为：
则万吉业超凸型拟合线段初端坐标为：
万吉业超凸型拟合线段末端坐标为：
陈元千水驱特征曲线对应的含水率变化曲线比较复杂，属超越方程，无法转化为用含水率函数直接表示累计产油量，因此，该水驱特征曲线拟合线段各点坐标通式时，应按下列步骤确定：
1)先利用
计算(试凑法)出含水率为fw时对应的累计产油量Np *；
2)其次，将Np *代入陈元千水驱特征曲线预估模型
中，计算出纵坐标
3)最后，确定出陈元千水驱特征曲线拟合线段各点坐标通式为：
按上述通式，可得到陈元千水驱特征曲线拟合线段初端坐标为：
式中，Np1 *是利用
计算(试凑法)出含水率为2％时对应的无水期累计产油量。
同样，可得到陈元千水驱特征曲线拟合线段末端坐标为：
式中，Np2 *是利用
计算(试凑法)出含水率为98％时对应的最终水驱可采储量。
具体的，根据采出程度和含水变化曲线，确定拟合线段的坐标通式；将第一含水率值和第二含水率值代入所述坐标通式中，得到初始坐标和末端坐标的表达式。
S105、将待定参数估值代入初始坐标和末端坐标的表达式中，确定拟合线段的初始点和末端点并绘制拟合线段。
根据拟合线段的坐标通式，线段首末端点可以得到一段拟合的水驱特征曲线，工程软件基于以上方法可以快速完成水驱特征曲线绘制，方便为油田勘探开发提供参考。
可以理解的是，以上各式中fw表示含水率，R表示采出程度，N表示地质储量，Np表示累计产油量，Lp表示累计产液量，Wp表示累计产水量。
通过本实施例提供的方法，可以在工程软件上快速完成水驱特征曲线的拟合绘制，基于已给出的12种水驱特征曲线拟合直线段初端坐标和末端坐标计算公式，方便了水驱特征曲线在矿场上应用。
应理解，上述实施例中各步骤的序号大小并不意味着执行顺序的先后，各过程的执行顺序应以其功能和内在逻辑确定，而不应对本发明实施例的实施过程构成任何限定。
图2为本发明实施例提供的一种水驱特征曲线快速拟合系统的结构示意图，该系统包括：
确定模块210：用于根据水驱特征曲线类型，分别确定各水驱特征曲线待定参数；
可选的，所述根据水驱特征曲线类型，分别确定各水驱特征曲线待定参数估值还包括：
通过多元线性回归或含单变量的最小二乘法求解所述待定参数。
计算模块220：用于将所述待定参数代入含水率变化关系式，计算采出程度和含水率变化曲线；
选取模块230：用于根据实际生产数据，选取适合描述油田含水变化规律的对应类型的水驱特征曲线；
设定模块240：用于设定第一含水率值对应的坐标为拟合线段的初始坐标，设定第二含水率值对应的坐标为末端坐标，分别确定对应类型的水驱特征曲线的初始坐标和末端坐标的表达式；
可选的，所述分别确定初始坐标和末端坐标的表达式包括：
根据采出程度和含水变化曲线，确定拟合线段的坐标通式；
将第一含水率值和第二含水率值代入所述坐标通式中，得到初始坐标和末端坐标的表达式。
拟合模块250：用于将待定参数估值代入初始坐标和末端坐标的表达式中，确定拟合线段的初始点和末端点并绘制拟合线段。
在上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有详述或记载的部分，可以参见其它实施例的相关描述。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成，所述的程序可以存储于一计算机可读取存储介质中，该程序在执行时，包括步骤S101至S105，所述的存储介质包括如：ROM/RAM、磁碟、光盘等。
以上所述，以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (3)
1.一种水驱特征曲线快速拟合方法，其特征在于，包括：
根据水驱特征曲线类型，分别确定各水驱特征曲线待定参数；
其中，通过多元线性回归或含单变量的最小二乘法求解所述待定参数；
将所述待定参数代入含水率变化关系式，计算采出程度和含水率变化曲线；
根据实际生产数据，选取适合描述油田含水变化规律的对应类型的水驱特征曲线；
设定第一含水率值对应的坐标为拟合线段的初始坐标，设定第二含水率值对应的坐标为末端坐标，分别确定对应类型的水驱特征曲线的初始坐标和末端坐标的表达式；
将待定参数估值代入初始坐标和末端坐标的表达式中，确定拟合线段的初始点和末端点并绘制拟合线段。
2.根据权利要求1所述的方法，其特征在于，所述分别确定初始坐标和末端坐标的表达式包括：
根据采出程度和含水变化曲线，确定拟合线段的坐标通式；
将第一含水率值和第二含水率值代入所述坐标通式中，得到初始坐标和末端坐标的表达式。
3.一种水驱特征曲线快速拟合系统，其特征在于，包括：
确定模块：用于根据水驱特征曲线类型，分别确定各水驱特征曲线待定参数；
其中，通过多元线性回归或含单变量的最小二乘法求解所述待定参数；
计算模块：用于将所述待定参数代入含水率变化关系式，计算采出程度和含水率变化曲线；
选取模块：用于根据实际生产数据，选取适合描述油田含水变化规律的对应类型的水驱特征曲线；
设定模块：用于设定第一含水率值对应的坐标为拟合线段的初始坐标，设定第二含水率值对应的坐标为末端坐标，分别确定对应类型的水驱特征曲线的初始坐标和末端坐标的表达式；
拟合模块：用于将待定参数估值代入初始坐标和末端坐标的表达式中，确定拟合线段的初始点和末端点并绘制拟合线段。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911019328.7A CN110751405B (zh) | 2019-10-24 | 2019-10-24 | 一种水驱特征曲线快速拟合方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911019328.7A CN110751405B (zh) | 2019-10-24 | 2019-10-24 | 一种水驱特征曲线快速拟合方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110751405A CN110751405A (zh) | 2020-02-04 |
| CN110751405B true CN110751405B (zh) | 2022-07-15 |
Family
ID=69279716
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911019328.7A Active CN110751405B (zh) | 2019-10-24 | 2019-10-24 | 一种水驱特征曲线快速拟合方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110751405B (zh) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104680244A (zh) * | 2014-12-24 | 2015-06-03 | 中国石油天然气股份有限公司 | 一种油气田产能预测方法及装置 |
| CN104912528A (zh) * | 2014-03-10 | 2015-09-16 | 中国石油化工股份有限公司 | 一种特高含水期水驱特征曲线拐点判识方法 |
| CN104915512A (zh) * | 2015-06-25 | 2015-09-16 | 陕西延长石油(集团)有限责任公司研究院 | 一种预测油田采出程度与含水率的方法 |
| CN107944126A (zh) * | 2017-11-21 | 2018-04-20 | 中国石油天然气股份有限公司 | 一种确定水驱油藏含水率的方法及装置 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2011133930A (ja) * | 2009-12-22 | 2011-07-07 | Fujitsu Ltd | 形状最適化プログラム、方法及び装置 |
-
2019
- 2019-10-24 CN CN201911019328.7A patent/CN110751405B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104912528A (zh) * | 2014-03-10 | 2015-09-16 | 中国石油化工股份有限公司 | 一种特高含水期水驱特征曲线拐点判识方法 |
| CN104680244A (zh) * | 2014-12-24 | 2015-06-03 | 中国石油天然气股份有限公司 | 一种油气田产能预测方法及装置 |
| CN104915512A (zh) * | 2015-06-25 | 2015-09-16 | 陕西延长石油(集团)有限责任公司研究院 | 一种预测油田采出程度与含水率的方法 |
| CN107944126A (zh) * | 2017-11-21 | 2018-04-20 | 中国石油天然气股份有限公司 | 一种确定水驱油藏含水率的方法及装置 |
Non-Patent Citations (1)
| Title |
|---|
| 水驱特征曲线法对油田进行动态预测;李忠臣;《内江科技》;20180625(第06期);79-81页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110751405A (zh) | 2020-02-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107944599B (zh) | 油气水平井产量的预测方法 | |
| CN111144663B (zh) | 计及出力波动过程的海上风电场超短期风功率预测方法 | |
| CN110685651B (zh) | 一种多层合采气井产量劈分方法及系统 | |
| CN111852463B (zh) | 气井产能评价方法及设备 | |
| CN108444884B (zh) | 基于Harris模型的油相相对渗透率曲线表征方法 | |
| Zhang et al. | When doesn't cokriging outperform kriging? | |
| CN110084431B (zh) | 一种页岩气井产量分析预测方法及系统 | |
| CN113486538B (zh) | 一种非常规油气井产量预测及压裂效果评价方法 | |
| CN112392478B (zh) | 一种快速预测低渗致密油藏经济可采储量的方法 | |
| CN108959739B (zh) | 一种抽水蓄能电站过渡过程压力脉动的分析方法及装置 | |
| CN104653170A (zh) | 特高含水阶段可采储量标定方法 | |
| CN110751405B (zh) | 一种水驱特征曲线快速拟合方法及系统 | |
| CN103809510B (zh) | 一种面向高精加工的自由曲面往复式刀具轨迹规划方法 | |
| CN109818377B (zh) | 一种基于振幅变化的自动发电控制性能评估方法及系统 | |
| CN110610288A (zh) | 一种油气井生产数据智能系统分析方法 | |
| CN113445988A (zh) | 一种低渗碳酸盐岩气藏气井产能评价方法 | |
| CN116629686A (zh) | 一种评估企业能耗数据的方法及装置 | |
| CN108051860B (zh) | 一种利用趋势面形态生成误差校正面的方法 | |
| CN115456260A (zh) | 客服话务量预测方法 | |
| CN112489148B (zh) | 一种砂岩型铀矿平米铀量等值线图成图方法 | |
| Delou et al. | A Robust Adaptive MPC coupled with Kalman Filter for an Electric Submersible Pump System: A multi-model approach | |
| CN113392372A (zh) | 气井全生命周期无阻流量确定方法及系统 | |
| CN111985146A (zh) | 一种措施潜力井的确定方法和确定装置 | |
| CN118049239B (zh) | 一种暗挖工程数据处理和隧道施工质量控制系统 | |
| CN115906591B (zh) | 一种基于XGBoost网络的油井动液面计算方法 |
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 |