CN112780252A - 一种适用于随钻测井正演模拟的地质导向模型及正演方法 - Google Patents
一种适用于随钻测井正演模拟的地质导向模型及正演方法 Download PDFInfo
- Publication number
- CN112780252A CN112780252A CN201911011932.5A CN201911011932A CN112780252A CN 112780252 A CN112780252 A CN 112780252A CN 201911011932 A CN201911011932 A CN 201911011932A CN 112780252 A CN112780252 A CN 112780252A
- Authority
- CN
- China
- Prior art keywords
- drilled
- model
- modeling
- well
- sequence
- 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 23
- 238000005553 drilling Methods 0.000 title claims abstract description 12
- 238000004364 calculation method Methods 0.000 claims description 10
- 239000011435 rock Substances 0.000 claims description 5
- 230000015572 biosynthetic process Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
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
-
- 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
-
- G—PHYSICS
- G09—EDUCATION; CRYPTOGRAPHY; DISPLAY; ADVERTISING; SEALS
- G09B—EDUCATIONAL OR DEMONSTRATION APPLIANCES; APPLIANCES FOR TEACHING, OR COMMUNICATING WITH, THE BLIND, DEAF OR MUTE; MODELS; PLANETARIA; GLOBES; MAPS; DIAGRAMS
- G09B25/00—Models for purposes not provided for in G09B23/00, e.g. full-sized devices for demonstration purposes
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Business, Economics & Management (AREA)
- Geology (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Marine Sciences & Fisheries (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Health & Medical Sciences (AREA)
- General Business, Economics & Management (AREA)
- General Health & Medical Sciences (AREA)
- Animal Husbandry (AREA)
- Geophysics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- Agronomy & Crop Science (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Educational Administration (AREA)
- Educational Technology (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提出一种适用于随钻测井正演模拟的地质导向模型及正演方法，根据建好的地质导向模型，沿井轨迹提取出对应单元的属性后，还考虑了与该单元上下相邻的单元属性，以及利用该单元的顶点坐标预测了该单元上下层界面的倾角，并将这些信息与仪器参数相关联进行正演模拟预测，本发明提高预测精度。
Description
技术领域
本发明涉及随钻地质导向领域，具体涉及一种适用于随钻测井正演模拟的地质导向模型及正演方法。
背景技术
随着越来越多的油气田被勘探开发，我国油气资源呈现出整体分散、局部集中的分布，开发目标开始转向那些规模更小、油层更薄、裂缝油藏和物性差的油藏，这对能否成功钻遇油气层提出了更高的技术要求。地质导向技术正是在这样的大背景下得到了快速发展，该技术通过井下传感器测量的地质及工程参数实时识别地下构造及属性信息，指导钻头顺利中靶及最大限度的在目的层中钻进。
地质导向作业通常需要地质导向软件的辅助，该软件可以直观的将地质模型、沿井轨迹预测或实测得到的测井数据等相关信息展示出来，并能实时更新井轨迹和地质模型，以方便实现这种互动的作业方式。目前的沿井轨迹预测测井数据的方法均只是简单的将模型中的属性提取出来，直接利用该属性作为预测数据，并没有将该属性与仪器的探测特性相关联，致使预测的精度很低。
发明内容
针对现有技术存在的技术问题，本发明提出一种适用于随钻测井正演模拟的地质导向模型及正演方法，模型的属性与仪器参数的关联，提高预测精度，进而提高地质导向效果。
本发明一种适用于随钻测井正演模拟的地质导向模型，包括以下步骤：
步骤一、选定建模区域
在待钻井周围选定建模区域，建模区域包含待钻井和该区域内的已钻邻井；
步骤二、对地层进行分层
根据要正演模拟的曲线类型，选择相对应的邻井测井数据对地层进行分层；
步骤三、对建模区域进行xy平面正交网格剖分
构建直角坐标系，令z方向为垂直水平面向下，则xy平面与水平面平行；在建模区域内对xy平面进行正交网格剖分；
步骤四、将待钻井轨迹垂直投影到步骤三得到的平面网格上，并拾取出待钻井轨迹与网格的所有交点坐标x、y；
步骤五、利用所有邻井的层界面数据，选取插值方法，在步骤四中得到的交点坐标x、y处进行插值，得到所有交点位置处的各地层界面的z坐标值，即x、y、z坐标值；
步骤六、构建地质导向剖面
利用步骤五得到的所有空间点，按对应的层界面序号，将每个层界面的所有层界面空间点连接起来；
步骤七、根据邻井的属性数据，采用插值方法，得到二维剖面的各地层间每个单元格的属性值；
步骤八、将步骤七得到的各单元格属性值、序号及四个顶点坐标集成到一起，作为一个存储单元，实现对整个地质导向模型的存储。
进一步，所述步骤一中邻井至少要三口以上。
进一步，所述步骤二中正演模拟的曲线类型有伽马曲线或电阻率曲线，各邻井的分层序号一致。
进一步，所述步骤五、七中插值方法有距离反比插值、格林样条插值。
利用本发明模型进行正演的方法，包括以下步骤：
步骤一、对待钻井轨迹进行离散，离散后得到一个待进行正演模拟的序列，序列中每个元素均包含该元素的序号、三维坐标值以及该位置处的井斜角；
步骤二、将步骤一中的离散序列后每个元素的位置对应到模型中建好的地质导向剖面中，拾取出每个元素位置处小单元格的属性和四个顶点坐标，并以此单元格为基础，向上下扩展出相邻单元格的属性，根据该单元四个顶点坐标，计算出上面两点与水平面的角度，作为上界面倾角，同样，计算出下面两点与水平面的角度，作为下界面倾角；
步骤三、通过步骤一、二得到了井轨迹序列中每个元素对应的地层模型参数，则将每个元素信息和对应的地层模型参数信息组合成一个独立的模型单元，组合成一个正演模型序列；
步骤四、在实现正演计算时采用的是逐点运算，全部计算完成后，得到一个正演响应的序列，即为对应待钻井轨迹的正演响应结果。
进一步，所述步骤三中每个序列单元包含如下信息：井轨迹坐标、井斜角、井轨迹所在层及上下围岩属性信息、上下层界面位置及倾角。
进一步，所述步骤四中的逐点运算为将每个单元作为快速正演输入，通过调用正演计算模块，输出该正演模型单元下的仪器响应。
进一步，所述步骤四中待钻井轨迹的正演响应结果包含的信息有该元素的待钻井轨迹位置信息、正演响应值。
有益效果：
1、本发明根据建好的地质导向模型，沿待钻井轨迹提取出对应单元的属性后，还考虑了与该单元上下相邻的单元属性，以及利用该单元的顶点坐标预测了该单元上下层界面的倾角，并将这些信息与仪器参数相关联进行正演模拟预测，这样使得到的预测曲线与只利用原始属性值相比更接近仪器的真实响应，进而提高地质导向效果。
2、本发明给出的建模方法实现简单，且方便从模型中拾取出快速正演模型所需的参数，方便实现快速正演模拟，能有效增强地质导向的应用效果。
附图说明
图1在建模区域对xy平面进行剖分得到的平面网格
图2井轨迹垂直投影与水平网格交点示意图
图3沿井轨迹的二维垂直剖面以及层界面位置示意图
图4快速正演模型
具体实施方式
下面结合附图对本发明做进一步的描述。
一、地质导向模型构建实现过程
步骤一、选定建模区域
在待钻井周围选定建模区域，建模区域包含待钻井和该区域内的已钻邻井。一般情况下，至少要保证有三口以上的邻井资料。
步骤二、对地层进行分层
根据要正演模拟的曲线类型，选择相对应的邻井测井数据对地层进行分层。比如要正演模拟伽马曲线，则需根据邻井的伽马测井数据进行地层分层，要正演模拟电阻率曲线，则需根据邻井的电阻率测井数据进行地层分层，且要求各邻井的分层序号一致。
步骤三、对建模区域进行xy平面正交网格剖分
构建直角坐标系，令z方向为垂直水平面向下，则xy平面与水平面平行。选择建模区域一般为立方体形状，在建模区域内对xy平面进行正交网格剖分，该网格疏密程度可控制，得到如下图1所示的平面网格。邻井数量越多建模精度越高，但计算量越大。
步骤四、将待钻井轨迹垂直(沿z方向)投影到步骤三得到的平面网格上，并拾取出待钻井轨迹与网格的所有交点坐标(x，y)。见如下示意图2。
步骤五、利用所有邻井的层界面数据，选取插值方法，如距离反比插值、格林样条插值等，在步骤四中得到的交点坐标(x1，y1)，(x2，y2)，…，(xn，yn)处进行插值，得到所有交点位置处的各地层界面的z坐标值。
比如在(x1，y1)交点处插值可得到该位置处各地层界面点的z坐标值，假设有n个地层界面，则这些地层界面点可以用一些列空间点的坐标(x1，y1，z1)，(x1，y1，z2)，…，(x1，y1，zn)表示，很明显，这些空间点和待钻井轨迹在一个垂直剖面上。
步骤六、构建地质导向剖面，利用步骤五得到的所有空间点，按对应的层界面序号，将每个层界面的所有层界面空间点连接起来(即层界面曲线)，并显示在沿井轨迹的二维垂直切片上，如图3所示，其中水平坐标为井轨迹沿水平方向的位移，垂直坐标为垂直于水平面向下的z方向(相当于实际钻井作业中的垂深)。
步骤七、在二维剖面的各地层内，纵向上只有一个单元格，如单元格M，横向上可以根据需要划分多个单元格，根据邻井的属性数据，采用插值方法，如距离反比加权、克里金等，得到各层内每个小单元格的属性值，如伽马、电阻率等。
步骤八、将步骤七得到的各单元格属性值、序号及四个顶点(a、b、c、d)坐标集成到一起，作为一个存储单元，用来实现对整个地质导向模型的存储。序号为单元格的位置，模型数据集成位置、属性、坐标。
二、地质导向模型与正演模拟相结合的实现过程
步骤一、对待钻井轨迹进行离散，离散点的疏密程度可控，离散后得到一个待进行正演模拟的序列，序列中每个元素均包含该元素的序号、三维坐标值以及该位置处的井斜角。
步骤二、得到步骤一中的离散序列后，将序列中每个元素的位置对应到建好的地质导向剖面中，拾取出每个元素位置处小单元格的属性和四个顶点坐标，并以此单元格为基础，向上下扩展出相邻单元格的属性，则该单元格和上下相邻单元格的属性即可作为对应元素位置的目的层和上下围岩属性。根据该单元四个顶点坐标，计算出上面两点与水平面的角度，作为上界面倾角，同样，计算出下面两点与水平面的角度，作为下界面倾角。
为实现快速正演模拟，需要对地层模型进行简化，通常可将多层地层模型简化成如下三层模型，即只考虑仪器所在的目的层和上下围岩三层地层，并且可以考虑仪器与地层的相对倾角。
步骤三、通过步骤一、二得到了井轨迹序列中每个元素对应的地层模型参数，则将每个元素信息和对应的地层模型参数信息组合成一个独立的模型单元，这样可以组合成一个正演模型序列，每个序列单元包含如下信息：井轨迹坐标、井斜角、井轨迹所在层及上下围岩属性信息、上下层界面位置及倾角。
步骤四、由于正演模型序列的每个单元是独立的，在实现正演计算时采用的是逐点运算，即将每个单元作为快速正演输入，通过调用正演计算模块，输出该正演模型单元下的仪器响应，全部计算完成后，则可以得到一个正演响应的序列，该序列就是对应井轨迹的正演响应结果。该序列每个元素包含的信息如下：该元素的井轨迹位置信息、正演响应值。
上述具体实施方式仅用于解释和说明本发明的技术方案，但并不能构成对权利要求保护范围的限定。本领域技术人员应当清楚，在本发明技术方案的基础上做任何简单的变形和替换得到的新的技术方案，均落入本发明的保护范围。
Claims (8)
1.一种适用于随钻测井正演模拟的地质导向模型，其特征在于包括以下步骤：
步骤一、选定建模区域
在待钻井周围选定建模区域，建模区域包含待钻井和该区域内的已钻邻井；
步骤二、对地层进行分层
根据要正演模拟的曲线类型，选择相对应的邻井测井数据对地层进行分层；
步骤三、对建模区域进行xy平面正交网格剖分
构建直角坐标系，令z方向为垂直水平面向下，则xy平面与水平面平行；在建模区域内对xy平面进行正交网格剖分；
步骤四、将待钻井轨迹垂直投影到步骤三得到的平面网格上，并拾取出待钻井轨迹与网格的所有交点坐标x、y；
步骤五、利用所有邻井的层界面数据，选取插值方法，在步骤四中得到的交点坐标x、y处进行插值，得到所有交点位置处的各地层界面的z坐标值，即x、y、z坐标值；
步骤六、构建地质导向剖面
利用步骤五得到的所有空间点，按对应的层界面序号，将每个层界面的所有层界面空间点连接起来；
步骤七、根据邻井的属性数据，采用插值方法，得到二维剖面的各地层间每个单元格的属性值；
步骤八、将步骤七得到的各单元格属性值、序号及四个顶点坐标集成到一起，作为一个存储单元，实现对整个地质导向模型的存储。
2.如权利要求1所述的一种适用于随钻测井正演模拟的地质导向模型，其特征在于，所述步骤一中邻井至少要三口以上。
3.如权利要求1所述的一种适用于随钻测井正演模拟的地质导向模型，其特征在于，所述步骤二中正演模拟的曲线类型有伽马曲线或电阻率曲线，各邻井的分层序号一致。
4.如权利要求1所述的一种适用于随钻测井正演模拟的地质导向模型，其特征在于，所述步骤五、七中插值方法有距离反比插值、格林样条插值。
5.利用权利要求1所述的模型进行正演的方法，其特征在于包括以下步骤：
步骤一、对待钻井轨迹进行离散，离散后得到一个待进行正演模拟的序列，序列中每个元素均包含该元素的序号、三维坐标值以及该位置处的井斜角；
步骤二、将步骤一中的离散序列后每个元素的位置对应到模型中建好的地质导向剖面中，拾取出每个元素位置处小单元格的属性和四个顶点坐标，并以此单元格为基础，向上下扩展出相邻单元格的属性，根据该单元四个顶点坐标，计算出上面两点与水平面的角度，作为上界面倾角，同样，计算出下面两点与水平面的角度，作为下界面倾角；
步骤三、通过步骤一、二得到了井轨迹序列中每个元素对应的地层模型参数，则将每个元素信息和对应的地层模型参数信息组合成一个独立的模型单元，组合成一个正演模型序列；
步骤四、在实现正演计算时采用的是逐点运算，全部计算完成后，得到一个正演响应的序列，即为对应待钻井轨迹的正演响应结果。
6.如权利要求5所述的正演方法，其特征在于，所述步骤三中每个序列单元包含如下信息：井轨迹坐标、井斜角、井轨迹所在层及上下围岩属性信息、上下层界面位置及倾角。
7.如权利要求5所述的正演方法，其特征在于，所述步骤四中的逐点运算为将每个单元作为快速正演输入，通过调用正演计算模块，输出该正演模型单元下的仪器响应。
8.如权利要求5所述的正演方法，其特征在于，所述步骤四中待钻井轨迹的正演响应结果包含的信息有该元素的待钻井轨迹位置信息、正演响应值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911011932.5A CN112780252B (zh) | 2019-10-23 | 2019-10-23 | 一种适用于随钻测井正演模拟的地质导向模型及正演方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911011932.5A CN112780252B (zh) | 2019-10-23 | 2019-10-23 | 一种适用于随钻测井正演模拟的地质导向模型及正演方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112780252A true CN112780252A (zh) | 2021-05-11 |
| CN112780252B CN112780252B (zh) | 2022-08-05 |
Family
ID=75747209
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911011932.5A Active CN112780252B (zh) | 2019-10-23 | 2019-10-23 | 一种适用于随钻测井正演模拟的地质导向模型及正演方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112780252B (zh) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120232871A1 (en) * | 2011-03-10 | 2012-09-13 | Ivan Priezzhev | Method for 3-d gravity forward modeling and inversion in the wavenumber domain |
| CN103774988A (zh) * | 2013-06-05 | 2014-05-07 | 中国石油大学(华东) | 一种实时随钻伽马正演地质导向钻井方法 |
| CN104847264A (zh) * | 2015-05-05 | 2015-08-19 | 中国海洋石油总公司 | 一种实现地质导向的方法和装置 |
| CN105317431A (zh) * | 2014-07-30 | 2016-02-10 | 中国石油化工股份有限公司 | 一种用于解释和评价水平井测井参数的方法 |
| AU2016391066A1 (en) * | 2016-02-05 | 2018-06-28 | Halliburton Energy Services, Inc. | Optimized geosteering using real-time geological models |
| US20180252101A1 (en) * | 2017-03-06 | 2018-09-06 | Baker Hughes, A Ge Company, Llc | Creation of structural earth formation models |
| CN109657346A (zh) * | 2018-12-18 | 2019-04-19 | 中国石油大学(华东) | 基于积分计算的随钻方位伽马正演方法、装置及设备 |
-
2019
- 2019-10-23 CN CN201911011932.5A patent/CN112780252B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120232871A1 (en) * | 2011-03-10 | 2012-09-13 | Ivan Priezzhev | Method for 3-d gravity forward modeling and inversion in the wavenumber domain |
| CN103774988A (zh) * | 2013-06-05 | 2014-05-07 | 中国石油大学(华东) | 一种实时随钻伽马正演地质导向钻井方法 |
| CN105317431A (zh) * | 2014-07-30 | 2016-02-10 | 中国石油化工股份有限公司 | 一种用于解释和评价水平井测井参数的方法 |
| CN104847264A (zh) * | 2015-05-05 | 2015-08-19 | 中国海洋石油总公司 | 一种实现地质导向的方法和装置 |
| AU2016391066A1 (en) * | 2016-02-05 | 2018-06-28 | Halliburton Energy Services, Inc. | Optimized geosteering using real-time geological models |
| US20180252101A1 (en) * | 2017-03-06 | 2018-09-06 | Baker Hughes, A Ge Company, Llc | Creation of structural earth formation models |
| CN109657346A (zh) * | 2018-12-18 | 2019-04-19 | 中国石油大学(华东) | 基于积分计算的随钻方位伽马正演方法、装置及设备 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112780252B (zh) | 2022-08-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103279986B (zh) | 三维水平地质剖面图制作方法及其用途 | |
| CA2713019C (en) | Gridless geological modeling | |
| CA2985743C (en) | Geologic stratigraphy via implicit and jump functions | |
| EP3293552B1 (en) | System and method for editing geological models by switching between volume-based models and surface-based structural models augmented with stratigraphic fiber bundles | |
| US10853533B2 (en) | Three-dimensional fracture abundance evaluation of subsurface formation based on geomechanical simulation of mechanical properties thereof | |
| US20180348401A1 (en) | Constructing Structural Models of the Subsurface | |
| EP2300945A1 (en) | Virtual petroleum system | |
| CA2724115A1 (en) | 3d visualization of 2d geophysical data | |
| GB2474740A (en) | Gridless geological modeling of a structural framework | |
| EP3371629A1 (en) | Representing structural uncertainty in a mesh representing a geological environment | |
| CN110244021B (zh) | 一种基于各向异性插值的地层分层方法 | |
| US20220163692A1 (en) | Modeling and simulating faults in subterranean formations | |
| CN105719346A (zh) | 山地三维地质体建模方法及综合地学信息演示系统 | |
| EP3455459A2 (en) | Three-dimensional fracture abundance evaluation of subsurface formations | |
| CN107564101B (zh) | 一种基于AutoCAD建立可视化三维地质模型的方法 | |
| CN114542056A (zh) | 水平井地层解释模型快速构建方法和装置 | |
| US10650107B2 (en) | Three-dimensional subsurface formation evaluation using projection-based area operations | |
| CN112780252B (zh) | 一种适用于随钻测井正演模拟的地质导向模型及正演方法 | |
| CN115035258A (zh) | 一种高效的基于cad钻孔柱状图的城市三维地质建模方法 | |
| US11175423B2 (en) | Real time deformation of seismic slices using programmable shaders | |
| CN112292714B (zh) | 基于断层辐射的网格隔分 | |
| US11353622B2 (en) | Systems and methods for hydrocarbon reservoir three dimensional unstructured grid generation and development | |
| CN110322563B (zh) | 构建地质模型的补充钻孔选址方法及系统 | |
| US20220171091A1 (en) | Hydrocarbon flow simulation | |
| CN117272481A (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 |