CN111753422A - 一种川西地区河道弯度的计算方法 - Google Patents

一种川西地区河道弯度的计算方法 Download PDF

Info

Publication number
CN111753422A
CN111753422A CN202010583180.6A CN202010583180A CN111753422A CN 111753422 A CN111753422 A CN 111753422A CN 202010583180 A CN202010583180 A CN 202010583180A CN 111753422 A CN111753422 A CN 111753422A
Authority
CN
China
Prior art keywords
river
river channel
channel
camber
width
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
Application number
CN202010583180.6A
Other languages
English (en)
Other versions
CN111753422B (zh
Inventor
张小菊
邓虎成
伏美燕
阮韵淇
赵爽
胡笑非
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chengdu Univeristy of Technology
Original Assignee
Chengdu Univeristy of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chengdu Univeristy of Technology filed Critical Chengdu Univeristy of Technology
Priority to CN202010583180.6A priority Critical patent/CN111753422B/zh
Publication of CN111753422A publication Critical patent/CN111753422A/zh
Application granted granted Critical
Publication of CN111753422B publication Critical patent/CN111753422B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Biology (AREA)
  • Computer Hardware Design (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Operations Research (AREA)
  • Probability & Statistics with Applications (AREA)
  • Evolutionary Computation (AREA)
  • Algebra (AREA)
  • Geometry (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了河道参数计算技术领域的一种川西地区河道弯度的计算方法，具体包括以下步骤：基于川西地区12口井所在河道的河道弯度和河道坡度数据，编制河道弯度和河道坡度的交汇图版，得出河道弯度和河道坡度的一个线性定量关系式，然后再此基础上，再根据现有提出的一系列参数公式，一步步反推(步骤2‑5)，得到最终的河道弯度、河道宽度、河道深度三者之间的定量关系式，可以直接从单井测井资料提取相关数据，由相关公式反演计算精确得到该单井所在位置的河道弯度指数，本发明能够减少传统经验方法的任意性，可以基于河道弯度更加精准地预测优质储层。

Description

一种川西地区河道弯度的计算方法
技术领域
本发明涉及河道参数计算技术领域，具体为一种川西地区河道弯度的计算方法。
背景技术
河道的形态是多种自然因素和力量共同作用的结果，包括河道结构成份，水流侵蚀能力,河道植物的密度和生长率，随河道沉淀物的构成及移动速度，沉淀物聚集范围，以及区域河流携带物的聚集和分解率等。河床的笔直和弯曲度直接影响河流速度，河道的形态主要包括河床的深潜度，平滑度及驳岸的坡度。因此河道的形态学对河道的物理形态和河道生态结构具有重大意义。为了进一步分析河道的特点及变迁规律，前人主要是选择河道弯度指数来对河道弯曲度进行定量表征及分析。
河道弯度指数(Ω)主要是指沿河道中心线的总长度(S/km)与河道两端点之间的直线距离(D/km)之比。
因为每口井又不可能都在同一条河道上，而且在同一条河道上的多口井，也是分布不同位置上。在具体研究河道弯度相关科学问题时，一般都是基于现有的单井资料来开展的，前人主要是根据沿河道中心线的总长度(S/km)与河道两端点之间的直线距离(D/km)之比来定义河道弯度指数，通过上述这种方法计算的，基于河道砂体的地震属性图来估算S和D的值(见图1)。假设在同一段河道分布的若干口井，若采用该套方法来计算，则计算出的河道弯度指数与某些井实际上所在位置的弯曲度大小存在明显的差异，若某口井刚好位于河道弯度极大处，则利用该方法可能计算出来的具体河道弯度指数会明显比实际的偏小，从而影响后续的研究分析。
前期研究表明在川西地区河道弯度控制了研究区砂岩储层的物性及砂体叠置类型，另外，侧积叠置均一型砂体叠置类型主要发育在河道弯度较大处(Ω：1.15-1.25)，且该处的砂岩储层物性较好。因此，在与川西地区相似的气田，可以基于河道弯度来预测优质储层。故优化计算河道弯度指数的方法对于河道弯度大小精确定量表征显得格外重要。
基于此，本发明设计了一种川西地区河道弯度的计算方法，以解决上述提到的问题。
发明内容
本发明的目的在于提供一种川西地区河道弯度的计算方法，客观的定量方法，统计理论基础较强，可以准确捕捉河道弯度、河道宽度及河道深度的内在交互关系，能够减少传统经验方法的任意性，在与川西地区相似的气田，可以基于河道弯度更加精准地来预测优质储层。
为实现上述目的，本发明提供如下技术方案：一种川西地区河道弯度的计算方法，具体包括以下步骤：
S1：基于川西地区12口井所在河道的河道弯度和河道坡度数据，编制河道弯度和河道坡度的交汇图版，得到该区的河道弯度和河道坡度的关系式：
Ω＝-0.3006*S+0.5326(R2＝0.4122) (1)
其中，Ω为河道弯度指数，S为河道弯度，R2为回归平方和占总误差平方和的比例，反映回归直线的拟合程度；
S2：根据河道弯度与年平均流量及沉积负载指数之间的关系，即：
S＝60M-0.38*Qm -0.32 (2)，
由公式(1)和公式(2)推导可得河道弯度与沉积负载指数及年平均流量之间的关系，即：
Ω＝-0.3006*(60M-0.38*Qm -0.32)+0.5326＝-18.036M-0.38*Qm -0.32+0.5326(R2＝0.4122)(3)；
其中，Ω为河道弯度指数，S为河道弯度，Qm为年平均流量，M为沉积负载指数，R2为回归平方和占总误差平方和的比例，反映回归直线的拟合程度；
S3：根据年平均流量与河道宽度及沉积负载指数之间的关系，即：
Qm 0.38＝W*M0.39/43.7 (4)，
由公式(3)和公式(4)推导可得河道弯度与河道宽度及沉积负载指数之间的关系，即：
Ω＝-18.036M-0.38*(W*M0.39/43.7)-0.32/0.38+0.5326＝-18.036/43.7(-0.32/0.38)*M(-0.2692/0.38)*W(-0.32/0.38)+0.5326(R2＝0.4122) (5)，
其中，Ω为河道弯度指数，W为河道宽度，M为沉积负载指数，R2为回归平方和占总误差平方和的比例，反映回归直线的拟合程度；
S4：根据河道宽度、河道深度及沉积负载指数之间的关系，即：
宽深比F＝W/D＝255M-1.08 (6)，
由公式(6)推导可得M＝(W/255D)(-1/1.08) (7)，
其中，W为河道宽度，M为沉积负载指数，D为河道深度；
S5：由公式(5)和(7)进一步地推导可得到河道弯度、河道宽度及河道深度之间的关系，即：
Ω＝[(-18.036)*255(-0.2692/0.4104)/43.7(-0.32/0.38)]*W(-0.0764/0.4104)*D(-0.2692/0.4104)+0.5326(R2＝0.4122) (8)，
其中，Ω为河道弯度指数，W为河道宽度，D为河道深度；
S6：根据满槽河流宽度及深度之间的关系，即：
W＝6.8D1.54 (9)
其中，D为河道深度，通过现有测井资料获取的不同井在不同河道的砂体厚度值，即为该口井所处河道位置的深度，根据公式(9)可以进一步计算出河道宽度；
S7：根据已有的河道深度和河道宽度数据代入到公式(8)中，即可得到不同井所在位置对应的河道弯度指数。
优选的，根据所述步骤S7得到不同井所在位置对应的河道弯度指数可以得到：河道弯度指数介于1.15-1.25之间的河道砂岩储层物性好且发育均一型砂体叠置类型；而介于1.00-1.04之间的河道砂岩储层物性差且发育复合型砂体叠置类型。
与现有技术相比，本发明的有益效果是：本发明提出的基于河道深度及河道宽度的川西地区河道弯度定量计算方法，可以直接从单井测井资料提取相关数据，由相关公式反演计算精确得到该单井所在位置的河道弯度指数，本发明能够减少传统经验方法的任意性，可以基于河道弯度更加精准地预测优质储层。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对实施例描述所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明地震振幅切片上河道形态的几何学参数示意图；
图2为本发明流程图；
图3为本发明川西地区河道弯度-河道坡度交汇图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。
实施例1
请参阅图2-3，本发明提供一种技术方案：一种川西地区河道弯度的计算方法，具体包括以下步骤：
S1：基于川西地区12口井所在河道的河道弯度和河道坡度数据，编制河道弯度和河道坡度的交汇图版，得到该区的河道弯度和河道坡度的关系式：
Ω＝-0.3006*S+0.5326(R2＝0.4122) (1)
其中，Ω为河道弯度指数，S为河道弯度，R2为回归平方和占总误差平方和的比例，反映回归直线的拟合程度；
S2：根据Schumm在1968年提出了河道弯度与年平均流量及沉积负载指数之间的关系，即：
S＝60M-0.38*Qm -0.32 (2)，
由公式(1)和公式(2)推导可得河道弯度与沉积负载指数及年平均流量之间的关系，即：
Ω＝-0.3006*(60M-0.38*Qm -0.32)+0.5326＝-18.036M-0.38*Qm -0.32+0.5326(R2＝0.4122)(3)；
其中，Ω为河道弯度指数，S为河道弯度，Qm为年平均流量，M为沉积负载指数，R2为回归平方和占总误差平方和的比例，反映回归直线的拟合程度；
S3：根据Schumm在1968年提出了年平均流量与河道宽度及沉积负载指数之间的关系，即：
Qm 0.38＝W*M0.39/43.7 (4)，
由公式(3)和公式(4)推导可得河道弯度与河道宽度及沉积负载指数之间的关系，即：
Ω＝-18.036M-0.38*(W*M0.39/43.7)-0.32/0.38+0.5326＝-18.036/43.7(-0.32/0.38)*M(-0.2692/0.38)*W(-0.32/0.38)+0.5326(R2＝0.4122) (5)，
其中，Ω为河道弯度指数，W为河道宽度，M为沉积负载指数，R2为回归平方和占总误差平方和的比例，反映回归直线的拟合程度；
S4：根据Schumm在1960年提出了河道宽度、河道深度及沉积负载指数之间的关系，即：
宽深比F＝W/D＝255M-1.08 (6)，
由公式(6)推导可得M＝(W/255D)(-1/1.08) (7)，
其中，W为河道宽度，M为沉积负载指数，D为河道深度；
S5：由公式(5)和(7)进一步地推导可得到河道弯度、河道宽度及河道深度之间的关系，即：
Ω＝[(-18.036)*255(-0.2692/0.4104)/43.7(-0.32/0.38)]*W(-0.0764/0.4104)*D(-0.2692/0.4104)+0.5326(R2＝0.4122) (8)，
其中，Ω为河道弯度指数，W为河道宽度，D为河道深度；
S6：满槽河流深度数据可以直接在测井柱子上读取研究层位的砂体厚度值即可，Leeder在1973年提出了满槽河流宽度及深度之间的关系，即：
W＝6.8D1.54 (9)
其中，D为河道深度，通过现有测井资料获取的不同井在不同河道的砂体厚度值，即为该口井所处河道位置的深度，根据公式(9)可以进一步计算出河道宽度；
S7：根据已有的河道深度和河道宽度数据代入到公式(8)中，即可得到不同井所在位置对应的河道弯度指数。
其中，根据步骤S7得到不同井所在位置对应的河道弯度指数可以得到：河道弯度指数介于1.15-1.25之间的河道砂岩储层物性好且发育均一型砂体叠置类型；而介于1.00-1.04之间的河道砂岩储层物性差且发育复合型砂体叠置类型。
在本说明书的描述中，参考术语“一个实施例”、“示例”、“具体示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
以上公开的本发明优选实施例只是用于帮助阐述本发明。优选实施例并没有详尽叙述所有的细节，也不限制该发明仅为所述的具体实施方式。显然，根据本说明书的内容，可作很多的修改和变化。本说明书选取并具体描述这些实施例，是为了更好地解释本发明的原理和实际应用，从而使所属技术领域技术人员能很好地理解和利用本发明。本发明仅受权利要求书及其全部范围和等效物的限制。

Claims (2)

1.一种川西地区河道弯度的计算方法，其特征在于：具体包括以下步骤：
S1：基于川西地区12口井所在河道的河道弯度和河道坡度数据，编制河道弯度和河道坡度的交汇图版，得到该区的河道弯度和河道坡度的关系式：
Ω＝-0.3006*S+0.5326(R2＝0.4122) (1)
其中，Ω为河道弯度指数，S为河道弯度，R2为回归平方和占总误差平方和的比例，反映回归直线的拟合程度；
S2：根据河道弯度与年平均流量及沉积负载指数之间的关系，即：
S＝60M-0.38*Qm -0.32 (2)，
由公式(1)和公式(2)推导可得河道弯度与沉积负载指数及年平均流量之间的关系，即：
Ω＝-0.3006*(60M-0.38*Qm -0.32)+0.5326＝-18.036M-0.38*Qm -0.32+0.5326(R2＝0.4122)(3)；
其中，Ω为河道弯度指数，S为河道弯度，Qm为年平均流量，M为沉积负载指数，R2为回归平方和占总误差平方和的比例，反映回归直线的拟合程度；
S3：根据年平均流量与河道宽度及沉积负载指数之间的关系，即：
Qm 0.38＝W*M0.39/43.7 (4)，
由公式(3)和公式(4)推导可得河道弯度与河道宽度及沉积负载指数之间的关系，即：
Ω＝-18.036M-0.38*(W*M0.39/43.7)-0.32/0.38+0.5326＝-18.036/43.7(-0.32/0.38)*M(-0.2692/0.38)*W(-0.32/0.38)+0.5326(R2＝0.4122) (5)，
其中，Ω为河道弯度指数，W为河道宽度，M为沉积负载指数，R2为回归平方和占总误差平方和的比例，反映回归直线的拟合程度；
S4：根据河道宽度、河道深度及沉积负载指数之间的关系，即：
宽深比F＝W/D＝255M-1.08 (6)，
由公式(6)推导可得M＝(W/255D)(-1/1.08) (7)，
其中，W为河道宽度，M为沉积负载指数，D为河道深度；
S5：由公式(5)和(7)进一步地推导可得到河道弯度、河道宽度及河道深度之间的关系，即：
Ω＝[(-18.036)*255(-0.2692/0.4104)/43.7(-0.32/0.38)]*W(-0.0764/0.4104)*D(-0.2692/0.4104)+0.5326(R2＝0.4122) (8)，
其中，Ω为河道弯度指数，W为河道宽度，D为河道深度；
S6：根据满槽河流宽度及深度之间的关系，即：
W＝6.8D1.54 (9)
其中，D为河道深度，通过现有测井资料获取的不同井在不同河道的砂体厚度值，即为该口井所处河道位置的深度，根据公式(9)可以进一步计算出河道宽度；
S7：根据已有的河道深度和河道宽度数据代入到公式(8)中，即可得到不同井所在位置对应的河道弯度指数。
2.根据权利要求1所述的一种川西地区河道弯度的计算方法，其特征在于：根据所述步骤S7得到不同井所在位置对应的河道弯度指数可以得到：河道弯度指数介于1.15-1.25之间的河道砂岩储层物性好且发育均一型砂体叠置类型；而介于1.00-1.04之间的河道砂岩储层物性差且发育复合型砂体叠置类型。
CN202010583180.6A 2020-06-23 2020-06-23 一种川西地区河道弯度的计算方法 Active CN111753422B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010583180.6A CN111753422B (zh) 2020-06-23 2020-06-23 一种川西地区河道弯度的计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010583180.6A CN111753422B (zh) 2020-06-23 2020-06-23 一种川西地区河道弯度的计算方法

Publications (2)

Publication Number Publication Date
CN111753422A true CN111753422A (zh) 2020-10-09
CN111753422B CN111753422B (zh) 2022-03-22

Family

ID=72676929

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010583180.6A Active CN111753422B (zh) 2020-06-23 2020-06-23 一种川西地区河道弯度的计算方法

Country Status (1)

Country Link
CN (1) CN111753422B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105607146A (zh) * 2015-09-10 2016-05-25 中国海洋石油总公司 一种曲流河砂体规模的定量表征方法
CN106777440A (zh) * 2015-11-24 2017-05-31 中国石油化工股份有限公司 三角洲前缘水下分流河道宽度预测方法
CN110927817A (zh) * 2018-09-20 2020-03-27 中国石油化工股份有限公司 一种薄层叠置河道水平井轨迹优化设计及控制的方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105607146A (zh) * 2015-09-10 2016-05-25 中国海洋石油总公司 一种曲流河砂体规模的定量表征方法
CN106777440A (zh) * 2015-11-24 2017-05-31 中国石油化工股份有限公司 三角洲前缘水下分流河道宽度预测方法
CN110927817A (zh) * 2018-09-20 2020-03-27 中国石油化工股份有限公司 一种薄层叠置河道水平井轨迹优化设计及控制的方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
XIAOJU ZHANG 等: "River channel pattern controls on the quality of sandstone reservoirs: A case study from the Jurassic Shaximiao formation of western Sichuan Basin, China", 《JOURNAL OF PETROLEUM SCIENCE AND ENGINEERING》 *
XU YI 等: "Pn-wave velocity and anisotropy of the western Sichuan and Longmen Mountain region,china", 《SCIENCE CHINA EARTH SCIENCES 》 *
刘畅: "川中龙岗地区侏罗系沙溪庙组一段河道砂体特征及沉积过程重建", 《中国优秀硕士学位论文全文数据库 基础科学辑》 *
卜淘: "川西坳陷东坡侏罗系沙溪庙组三角洲河道砂体构型", 《断块油气田》 *
张小菊: "川西东坡沙溪庙组砂体构型与气水分布研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 *

Also Published As

Publication number Publication date
CN111753422B (zh) 2022-03-22

Similar Documents

Publication Publication Date Title
CN110334431A (zh) 一种低渗透致密气藏单井控制储量计算及剩余气分析方法
CN105651676B (zh) 一种水平井规则开发井网下的储层非均质性表征方法
Jiongxin Evolution of mid‐channel bars in a braided river and complex response to reservoir construction: an example from the middle Hanjiang River, China
CN104631392B (zh) 一种基于河相关系的航道整治设计方法
Elliott Depositional architecture of a sand-rich, channelized turbidite system: the Upper Carboniferous Ross Sandstone Formation, western Ireland
CN101725346A (zh) 油藏井间动态连通性反演方法
CN106372341B (zh) 一种水库断面法水位库容曲线的修正方法
CN103939066A (zh) 一种一注多采井组定注水量确定油井产液量的方法
CN112989565B (zh) 一种基于水流挟沙能力的造床流量计算方法
CN115114834B (zh) 一种复杂情况下的压裂井试井模拟方法
CN101864950B (zh) 连续型致密砂岩气分布的预测方法
Allen The plan shape of current ripples in relation to flow conditions
CN107808068A (zh) 一种用于致密砂岩气藏高产富集评价的方法
CN114966886A (zh) 一种基于构型的致密气藏砂体连通性等级判定方法
CN111104746B (zh) 一种基于小波分析的河漫滩高程确定方法
CN111753422B (zh) 一种川西地区河道弯度的计算方法
Wu et al. A long-term morphological modeling study on the evolution of the Pearl River Delta, network system, and estuarine bays since 6000 yr BP
CN101592027B (zh) 一种依据修正毛管数曲线的复合驱油方法
CN106320255A (zh) 一种漫滩水流滩槽过流量的计算方法
CN110428497A (zh) 辫状河训练图像生成方法
Liu et al. The Control Theory and Application for Well Pattern Optimization of Heterogeneous Sandstone Reservoirs
CN112329192A (zh) 厚层砂岩油藏注聚均衡驱替设计方法
Gould et al. An analysis of infill drilling
CN114066166A (zh) 一种沙砾岩低渗透油藏co2混相驱筛选方法
RU2782712C1 (ru) Сетка переплетенных скважин для вытеснения, выполненная с возможностью введения текучей среды для разрабатывания особенно толстых или многопластовых нефтегазовых залежей

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