CN112302639B - 一种适于裂缝复杂渗流空间体的曲率表征方法 - Google Patents

一种适于裂缝复杂渗流空间体的曲率表征方法 Download PDF

Info

Publication number
CN112302639B
CN112302639B CN202011068616.4A CN202011068616A CN112302639B CN 112302639 B CN112302639 B CN 112302639B CN 202011068616 A CN202011068616 A CN 202011068616A CN 112302639 B CN112302639 B CN 112302639B
Authority
CN
China
Prior art keywords
point
curvature
complex
curved surface
average
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
Application number
CN202011068616.4A
Other languages
English (en)
Other versions
CN112302639A (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 Institute Of Urban Safety And Emergency Management
Original Assignee
Chengdu Institute Of Urban Safety And Emergency Management
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 Institute Of Urban Safety And Emergency Management filed Critical Chengdu Institute Of Urban Safety And Emergency Management
Priority to CN202011068616.4A priority Critical patent/CN112302639B/zh
Publication of CN112302639A publication Critical patent/CN112302639A/zh
Application granted granted Critical
Publication of CN112302639B publication Critical patent/CN112302639B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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
    • 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
    • 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)
  • Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Mining & Mineral Resources (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Fluid Mechanics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Probability & Statistics with Applications (AREA)
  • Operations Research (AREA)
  • Evolutionary Biology (AREA)
  • Algebra (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明公开了一种适于裂缝复杂渗流空间体的曲率表征方法，涉及油气开发技术领域，包括以下步骤：获取岩心裂缝的表面A和表面B的三维数据；将表面A和表面B构成的复杂空间体等效为一个光滑平面A’和一个复杂曲面B’所构成的空间体；取所述复杂曲面B’上的任一点G(i，j)，建立7×7标准网格，获得点G(i，j)四个方向上的平均曲率k(i,j)，进而获得所述复杂曲面B’的曲率，从而获得复杂渗流空间体的曲率。本发明采用7×7标准网格来计算复杂曲面B’的三维数据点四个方向的曲率，并对各方向上每一三维数据点不同距离上的曲率赋予不同的权重来计算各三维数据点的平均曲率，提高了复杂曲面B’中三维数据点曲率表征的精度。

Description

一种适于裂缝复杂渗流空间体的曲率表征方法
技术领域
本发明涉及油气田开发技术领域，特别涉及一种适于裂缝复杂渗流空间体的曲率表征方法。
背景技术
碳酸盐岩储层油气储量规模大，约蕴藏着全球48％的原油储量及28％的天然气储量，单井产量高，在世界油气工业中受到的关注度日益增加。但是碳酸盐岩储层非均质性极强，油气储集不连续，通过井眼沟通储集体的概率较小。因此利用酸压工艺造缝，有效沟通天然裂缝及远井地带储集体已经成为碳酸盐岩油气藏开采、建产、稳产和增产的关键技术之一，对于发现油气田、核实储量、稳定和提高油气田产量以及提高采收率和经济高效地发开油田发挥着举足轻重的作用。
酸蚀裂缝导流能力是衡量酸压效果的关键参数之一，通过工艺手段实现较高的、稳定的酸蚀裂缝导流能力，对于目前深层、高温高压的碳酸盐岩储层改造来说至关重要。而储层的岩石力学性质、泵注工艺和闭合应力等因素直接决定导流能力的形成及变化规律，那么正确评价酸蚀裂缝的导流能力的变化规律，也是制定合理酸压工艺方案的基础。
无论是人工裂缝还是天然裂缝，均会在原地应力作用下逐渐趋于闭合。研究表明，闭合应力对导流能力影响的关键因素在于裂缝趋于闭合时的裂缝表面形态。酸液对裂缝壁面的刻蚀形态及闭合应力作用下裂缝壁面的接触状态是影响着导流能力大小的重要因素，一方面由于裂缝面上微凸体尺寸、分布范围和方向都不尽相同，使得研究裂缝面接触问题变得十分复杂。另一方面，微观角度上，酸压裂缝主要依靠裂缝表面形成非均匀刻蚀或者裂缝滑移错动等方式，形成能使裂缝保持张开的非均质性支撑裂缝，进而表现出在宏观上酸液的溶蚀沟槽等刻蚀形态，使得裂缝在闭合时形成有效的渗流空间。因此，酸蚀裂缝渗流空间的形态与大小直接决定导流能力的高低，准确地表征酸蚀裂缝渗流空间对评价酸蚀裂缝的导流能力具有重要的意义。但目前来看，酸蚀裂缝的空间形貌表征参数缺乏，裂缝导流能力研究依旧停留在裂缝单面的形态之上等。
发明内容
本申请的目的在于克服现有技术中酸蚀裂缝这一复杂空间体曲率难以精确表征的问题，提供一种适于裂缝复杂渗流空间体的曲率表征方法，实现复杂空间体的曲率精确表征。
为了实现上述发明目的，本申请提供了以下技术方案：一种适于裂缝复杂渗流空间体的曲率表征方法，包括以下步骤：
获取岩心裂缝的表面A和表面B的三维数据；
将所述表面A和表面B分别等效为一个光滑平面A’和一个复杂曲面B’，使所述表面A和表面B构成的复杂空间体被等效为一个光滑平面A’和一个复杂曲面B’所构成的相对规则的空间体，得到复杂曲面B’的三维数据；
取所述复杂曲面B’上的任一点G(i，j)，并以所述点G(i，j)为中心建立7×7标准网格，并根据曲率表达式，获得点G(i，j)的曲率k(i,j)，进而获得所述复杂曲面B’的三维数据每一点的曲率；
根据所述复杂曲面B’中每一三维数据点的曲率得到所述复杂曲面B’的平均曲率，即为所述表面A和表面B构成的复杂空间体的平均曲率。
在上述技术方案中，将裂缝上下两个复杂曲面等效为一个复杂曲面和一个平面，以形成相对规则的空间体，简化酸蚀裂缝等这一类复杂空间体的复杂程度，进而提高了复杂空间体表征的可能性，降低了表征复杂空间体的难度和计算量。此外，本技术方案采用7×7标准网格来计算复杂曲面B’的三维数据点的曲率，提高了复杂曲面B’中三维数据点曲率表征的精度，有利于复杂曲面B’的起伏特性的准确表征。
进一步地，所述表面A和表面B构成的复杂空间体等效为一个光滑平面A’和一个复杂曲面B’所构成的相对规则的空间体，具体包括以下步骤：
以所述表面A和所述表面B的最低点所在的面为基准面，获得获取所述表面A的三维数据点(xa、ya、za)和所述表面B(xa、ya、zb)的三维数据点；
以所述表面A和所述表面B上对应点的za和zb的差值等效为复杂曲面B’对应的三维数据点的z值，即可获得复杂曲面B’对应的三维数据点，进而获得复杂曲面B’的三维数据。
进一步地，所述点G(i，j)的曲率k(i,j)为该点在水平方向上的平均曲率kx(i,j)、竖直方向上的平均曲率ky(i,j)、斜向上45°方向上的平均曲率kd(i,j)、斜向下45°方向上的平均曲率ku(i,j)的平均值。
进一步地，所述点G(i，j)在水平方向/竖直方向/斜向上45°方向/斜向下45°方向中任一方向上的平均曲率由三部分曲率的加权平均值构成。
即：所述点G(i，j)在水平方向的平均曲率为点G(i，j)在水平方向上的曲率、点G(i-1，j)在水平方向上的曲率，及点G(i+1，j)在水平方向上的曲率的加权平均值；所述点G(i，j)在竖直方向的平均曲率为点G(i，j)在竖直方向上的曲率、点G(i，j-1)在竖直方向上的曲率，及点G(i，j+1)在竖直方向上的曲率的加权平均值；所述点G(i，j)在斜向上45°方向的平均曲率为点G(i，j)在斜向上45°方向上的曲率、点G(i-1，j-1)在斜向上45°方向上的曲率，及点G(i+1，j+1)在斜向上45°方向上曲率的加权平均值；所述点G(i，j)在斜向下45°方向的平均曲率为点G(i，j)在斜向下45°方向上的曲率、点G(i+1，j-1)在斜向下45°方向上的曲率，及点G(i-1，j+1)在斜向下45°方向的曲率的加权平均值。
进一步地，所述点G(i，j)分别在水平方向/竖直方向/斜向下45°方向/斜向上45°方向的曲率的权重值均为2/4；所述点G(i-1，j)、点G(i，j-1)、点G(i+1，j-1)、点G(i-1，j-1)分别在水平方向/竖直方向/斜向下45°方向/斜向上45°方向的曲率的权重值均为1/4；所述点G(i+1，j)、点G(i，j+1)、点G(i-1，j+1)、点G(i+1，j+1)分别在水平方向/竖直方向/斜向下45°方向/斜向上45°方向的曲率的权重值均为1/4。
进一步地，所述点G(i，j)在水平方向/竖直方向/斜向下45°方向/斜向上45°方向的曲率根据所述曲率表达式计算；所述点G(i-1，j)、点G(i，j-1)、点G(i+1，j-1)、点G(i-1，j-1)分别在水平方向/竖直方向/斜向下45°方向/斜向上45°方向的曲率根据所述曲率表达式计算；所述点G(i+1，j)、点G(i，j+1)、点G(i-1，j+1)、点G(i+1，j+1)分别在水平方向/竖直方向/斜向下45°方向/斜向上45°方向的曲率根据所述曲率表达式计算。
进一步地，根据所述曲率表达式计算所述点G(i，j)在水平方向/竖直方向/斜向下45°方向/斜向上45°方向的曲率，根据所述曲率表达式计算点G(i-1，j)、点G(i，j-1)、点G(i+1，j-1)、点G(i-1，j-1)分别在水平方向/竖直方向/斜向下45°方向/斜向上45°方向的曲率，根据所述曲率表达式计算点G(i+1，j)、点G(i，j+1)、点G(i-1，j+1)、点G(i+1，j+1)分别在水平方向/竖直方向/斜向下45°方向/斜向上45°方向的曲率时，在7×7标准网格中对一阶求导和二阶求导中的差分计算分别赋权重值；
其中，所述点G(i，j)在水平方向的曲率在进行一阶求导和二阶求导时采用三阶差分法，每阶求导分为三个部分：第一部分为点G(i，j)左侧第一点、点G(i，j)右侧第一点，该部分权重为3/6；第二部分为点G(i，j)左侧第二点、点G(i，j)右侧第二点，该部分权重分别为2/6；第三部分为点G(i，j)左侧第三点、点G(i，j)右侧第三点，该部分权重分别为1/6；
其中，所述点G(i，j)在竖直方向的曲率在进行一阶求导和二阶求导时采用三阶差分法，每阶求导分为三个部分：第一部分为点G(i，j)上方第一点、点G(i，j)下方第一点，该部分权重为3/6；第二部分为点G(i，j)上方第二点、点G(i，j)下方第二点，该部分权重分别为2/6；第三部分为点G(i，j)上方第三点、点G(i，j)下方第三点，该部分权重分别为1/6；
其中，所述点G(i，j)在斜向下45°方向上的曲率在进行一阶求导和二阶求导时采用二阶差分法，每阶求导分为两个部分：第一部分为点G(i，j)左上第一点、点G(i，j)右下第一点，该部分权重为2/3；第二部分为G(i，j)左上第二点、点G(i，j)右下第二点，该部分权重为1/3。
其中，所述点G(i，j)在斜向上45°方向上的曲率在进行一阶求导和二阶求导时采用二阶差分法，每阶求导分为两个部分：第一部分为点G(i，j)右上第一点、点G(i，j)左下第一点，该部分权重为2/3；第二部分为G(i，j)右上第二点、点G(i，j)左下第二点，该部分权重为1/3。
进一步地，所述曲率表达式具体为：
Figure GDA0003055066580000041
进一步地，所述复杂曲面B’的平均曲率为所述复杂曲面B’中每一三维数据点的曲率之和与所述复杂曲面B’的面积的比值。
进一步地，所述复杂曲面B’的平均曲率为：
Figure GDA0003055066580000042
其中，m-表示复杂曲面B’的长度；n-表示复杂曲面B’的宽度；(i，j)表示曲面B’中的任一点，k(i，j)-表示复杂曲面B’中任一点的三维曲率。
由于所述表面A和表面B构成的复杂空间体等效为一个光滑平面A’和一个复杂曲面B’所构成的复杂空间体；因此所述表面A和表面B构成的复杂空间体的平均曲率即等效为所述复杂曲面B’的平均曲率，故所述表面A和表面B构成的复杂空间体的平均曲率等于所述复杂曲面B’的平均曲率。
与现有技术相比，本发明的具有以下有益效果：
本申请公开了一种适于裂缝复杂渗流空间体的曲率表征方法，将裂缝上下两个复杂曲面等效为一个复杂曲面和一个平面，以形成相对规则的空间体，简化酸蚀裂缝等这一类复杂空间体的复杂程度，进而提高了复杂空间体表征的可能性，降低了表征复杂空间体的难度和计算量。此外，本申请采用7×7标准网格来计算复杂曲面B’的三维数据点的曲率，提高了复杂曲面B’中三维数据点曲率表征的精度，有利于复杂曲面B’的起伏特性的准确表征。此外，本申请在每个三维数据点的计算时均计算了四个方向上的曲率值，并对于每一方向的不同点赋予不同的权重值，提高了每一三维数据点的曲率计算精度，使其计算得到的平均曲率能够准确反映该三维空间点的空间特征，为复杂空间体的表征提供更为可靠的方法与参数。
附图说明
图1是本发明一些实施例中公开的复杂曲面A和复杂曲面B等效为B’的示意图；
图2是本发明一些实施例中公开的7×7标准网格示意图；
图3是本发明一些实施例中公开的裂缝面a和b经酸蚀后的三维形貌图；
图4是本发明一些实施例中公开的从酸蚀后裂缝面a和b闭合后抽提复杂渗流空间体的等效三维图。
具体实施方式
下面结合试验例及具体实施方式对本发明作进一步的详细描述。但不应将此理解为本发明上述主题的范围仅限于以下的实施例，凡基于本发明内容所实现的技术均属于本发明的范围。
需要说明的是，虽然本申请的背景技术中记载了酸蚀裂缝这一复杂空间体表征的相关现有技术，但本领域技术人员应该知道，本申请公开的适于裂缝复杂渗流空间体的曲率表征方法也可以应用在其他复杂空间体的表征，如：封闭空间的定量表征、宏观与微观裂缝空间体的定量表征、三维地形图的定量表征、复杂曲面的定量表征等。
现有技术中，裂缝的两个面大多也均为粗糙面，即每一裂缝具有两个复杂曲面，两个复杂曲面构成的空间即为复杂空间体。但由于两个面都为复杂曲面，因此导致其形成复杂空间体难以表征。而常规的三维曲面曲率表征方法大多只能针对一个复杂曲面进行表征，其表征精度较低，因此在应用于裂缝形成的复杂空间体就会造成较大的误差，导致曲率的表征结果参考性较差。
为了解决上述技术问题，发明人在本申请中提出了一种适于裂缝复杂渗流空间体的曲率表征方法，参阅图1、图2，该方法具体种适于裂缝复杂渗流空间体的曲率表征方法，包括以下步骤：
获取岩心裂缝的表面A和表面B的三维数据；
将所述表面A和表面B构成的复杂空间体等效为一个光滑平面A’和一个复杂曲面B’所构成的相对规则的空间体，得到复杂曲面B’的三维数据；
以所述复杂曲面B’的三维数据的任一三维数据点k(i，j)为中心建立7×7标准网格，获得点k(i，j)的曲率，进而获得所述复杂曲面B’的三维数据每一点的曲率；
根据所述复杂曲面B’中每一三维数据点的曲率得到所述复杂曲面B’的平均曲率。
需要说明的是，所述岩心裂缝的表面A和表面B的三维数据采用TN 3DOMS.S三维形貌扫描仪获得，其精度优选为0.05mm。
具体实施过程中，可采用以下步骤进行扫描：
扫描裂缝面需将裂缝面烘干；
将裂缝面放置于加持台并保持水平，然后将加持台置于扫描托盘中央。
打开三维形貌扫描仪电源及拍照镜头盖，检查电脑、图像采集卡及相应传感器线对应连接是否正常；
打开扫描软件，调节曝光、增益参数至能清晰观察到裂缝面；
点击扫描按钮进行扫描，扫描结束后保存相应的点云数据；
打开Geomagic软件，导入保存的裂缝面点云数据，通过人工剪裁去除点云噪点和无效区域，保存剪裁后的点云数据；
使用surfer软件导入保存的剪裁后点云数据，生成裂缝面的三维形貌图，完成裂缝面的数字化提取。
需要说明的是，所述表面A和表面B构成的复杂空间体等效为一个光滑平面A’和一个复杂曲面B’所构成的相对规则的空间体，具体包括以下步骤：
以所述表面A和所述表面B的最低点所在的面为基准面，获得获取所述表面A的三维数据点(xa、ya、za)和所述表面B的三维数据点(xa、ya、zb)；
以所述表面A和所述表面B上对应点的za和zb的差值等效为复杂曲面B’对应的三维数据点的z值，即可获得复杂曲面B’对应的三维数据点，进而获得复杂曲面B’的三维数据。
即：取所述表面A的任一三维数据点C(xa、ya、za)和所述表面B的对应的三维数据点D(xa、ya、zb)，获得点C和点D在高度差d＝za-zb，等效时将所述表面B的最低点所在水平面作为基准面，将高度差d作为裂缝空间的等效高度坐标，即可获得复杂曲面B’的这一对应点的三维数据(xa、ya、d)，即为所述表面A和表面B上这一点的等效点；计算所有所述表面A和表面B的对应点，即可获得所述复杂曲面B’。
根据曲率定义可获得曲率表达式：
Figure GDA0003055066580000061
其中，θ—表示二维曲线的某点在曲线上的倾角；
s—表示这一点的弧长；
定义复杂曲面B’上的任一点G(i，j)的曲率为k(i,j)，并以所述点G(i，j)为中心建立7×7标准网格，其7×7标准网格如图2所示。
所述点G(i，j)的曲率k(i,j)均为该点在水平方向上的平均曲率kx(i,j)、竖直方向上的平均曲率ky(i,j)、斜向上45°方向上的平均曲率kd(i,j)、斜向下45°方向上的平均曲率ku(i,j)的平均值。
即：复杂曲面B’上的任一点G(i，j)的曲率为k(i,j)计算如下：
Figure GDA0003055066580000071
其中，所述点G(i，j)在每个方向上的平均曲率由三部分组成，在水平方向的平均曲率为点G(i，j)在水平方向上的曲率且权值2/4、点G(i-1，j)在水平方向上的曲率且权值1/4，及点G(i+1，j)在水平方向上的曲率且权值1/4的加权平均值；所述点G(i，j)在竖直方向的平均曲率为点G(i，j)在竖直方向上的曲率且权值2/4、点G(i，j-1)在竖直方向上的曲率且权值1/4，及点G(i，j+1)在竖直方向上的曲率且权值1/4的加权平均值；所述点G(i，j)在斜向上45°方向的平均曲率为点G(i，j)在斜向上45°方向上的曲率且权值2/4、点G(i-1，j-1)在斜向上45°方向上的曲率且权值1/4，及点G(i+1，j+1)在斜向上45°方向上曲率且权值1/4的加权平均值；所述点G(i，j)在斜向下45°方向的平均曲率为点G(i，j)在斜向下45°方向上的曲率且权值2/4、点G(i+1，j-1)在斜向下45°方向上的曲率且权值1/4，及点G(i-1，j+1)在斜向下45°方向的曲率且权值1/4的加权平均值。
根据式(1)计算点G(i，j)在不同方向曲率过程中的求导步骤时，在7×7标准网格中对一阶求导和二阶求导中的差分计算分别赋权重值，其中，点G(i，j)水平方向的曲率在进行一阶求导和二阶求导时采用三阶差分法，每阶求导分为三个部分：第一部分为点G(i，j)左侧第一点、点G(i，j)右侧第一点，该部分权重为3/6；第二部分为点G(i，j)左侧第二点、点G(i，j)右侧第二点，该部分权重分别为2/6；第三部分为点G(i，j)左侧第三点、点G(i，j)右侧第三点，该部分权重分别为1/6。点G(i，j)在竖直方向的曲率在进行一阶求导和二阶求导时采用三阶差分法，每阶求导分为三个部分：第一部分为点G(i，j)上方第一点、点G(i，j)下方第一点，该部分权重为3/6；第二部分为点G(i，j)上方第二点、点G(i，j)下方第二点，该部分权重分别为2/6；第三部分为点G(i，j)上方第三点、点G(i，j)下方第三点，该部分权重分别为1/6。点G(i，j)在斜向下45°方向上的曲率在进行一阶求导和二阶求导时采用二阶差分法，每阶求导分为两个部分：第一部分为点G(i，j)左上第一点、点G(i，j)右下第一点，该部分权重为2/3；第二部分为G(i，j)左上第二点、点G(i，j)右下第二点，该部分权重为1/3。点G(i，j)在斜向上45°方向上的曲率在进行一阶求导和二阶求导时采用二阶差分法，每阶求导分为两个部分：第一部分为点G(i，j)右上第一点、点G(i，j)左下第一点，该部分权重为2/3；第二部分为G(i，j)右上第二点、点G(i，j)左下第二点，该部分权重为1/3。
再次参阅图2，并根据式(2)可以得出，点G(4，4)的曲率k(4，4)为：
Figure GDA0003055066580000081
对于点G(4，4)的水平方向上，点G(4，4)水平方向上的平均曲率为：
Figure GDA0003055066580000082
并根据式(1)，获得：
Figure GDA0003055066580000083
根据所述7×7标准网格，分别对点G(4，4)在水平方向的三个部分进行一阶求导和二阶求导，并分别代入三部分的权重值，即可获得点G(4，4)在水平方向上的曲率，进而获得点G(4，4)的水平方向上的平均曲率，其一阶求导和二阶求导过程如下：
Figure GDA0003055066580000084
Figure GDA0003055066580000085
Figure GDA0003055066580000086
同理的，点G(4，4)竖直方向上的平均曲率为：
Figure GDA0003055066580000087
其中：
Figure GDA0003055066580000088
根据所述7×7标准网格，分别对点G(4，4)在竖直方向的三个部分进行一阶求导和二阶求导，并分别代入三部分的权重值，即可获得点G(4，4)在竖直方向上的曲率，进而获得点G(4，4)的竖直方向上的平均曲率，其一阶求导和二阶求导过程如下：
Figure GDA0003055066580000091
Figure GDA0003055066580000092
Figure GDA0003055066580000093
同理的，点G(4，4)斜向上45°方向上的平均曲率：
Figure GDA0003055066580000094
其中：
Figure GDA0003055066580000095
根据所述7×7标准网格，分别对点G(4，4)在斜向上45°方向上的两个部分进行一阶求导和二阶求导，并分别代入两部分的权重值，即可获得点G(4，4)在斜向上45°方向上的曲率，进而获得点G(4，4)的斜向上45°方向上的平均曲率，其一阶求导和二阶求导过程如下：
Figure GDA0003055066580000096
Figure GDA0003055066580000097
Figure GDA0003055066580000098
同理的，点G(4，4)斜向下45°方向上的平均曲率：
Figure GDA0003055066580000101
其中：
Figure GDA0003055066580000102
根据所述7×7标准网格，分别对点G(4，4)在斜向下45°方向上的两个部分进行一阶求导和二阶求导，并分别代入两部分的权重值，即可获得点G(4，4)在斜向下45°方向上的曲率，进而获得点G(4，4)的斜向下45°方向上的平均曲率，其一阶求导和二阶求导过程如下：
Figure GDA0003055066580000103
Figure GDA0003055066580000104
Figure GDA0003055066580000105
根据上述计算得到的点G(4，4)在水平方向上的平均曲率kx(4，4)、竖直方向上的平均曲率ky(4，4)、斜向下45°方向上的平均曲率kd(4，4)、斜向上45°方向上的平均曲率ku(4，4)，进而即可获得点G(4，4)的平均曲率k(4，4)。
需要说明的是，Zij为7×7标准网格中每一点的坐标；△y为在竖直方向相邻两点的距离；△x为在水平方向相邻两点的距离；△u在斜向上45°方向上相邻两点的距离；△d为在斜向下45°方向上相邻两点的距离。
计算所述复杂曲面B’中每一三维数据点的曲率，并根据式(23)计算所述复杂曲面B’的平均曲率Ka：
Figure GDA0003055066580000106
其中，m-表示复杂曲面B’的长度；n-表示复杂曲面B’的宽度。
以下以某气藏的灰岩为例，从该气藏的地层中采取长宽高尺寸分别为140mm×40mm×40mm的岩板作为标准岩板，并沿高度中心线劈裂得到自然断裂裂缝面a和b；采用20％稠化酸体系(20％HCl+0.8％胶凝剂+1％高温缓蚀剂+1％铁离子稳定剂+0.1％表面活性剂)对裂缝面a和b闭合后形成的裂缝空间进行酸蚀，酸蚀参数：酸液流速100mL/min，实验酸量200mL，实验温度110℃。酸蚀后裂缝面a和b即为表面A和表面B，酸蚀后裂缝面a和b经闭合后的复杂空间体即为表面A和表面B组成的复杂空间体。
采用三维光学扫描仪(TNS)对酸蚀后裂缝面a和b进行三维扫描，获得裂缝面a和b的三维点云数据信息，并通过人工剪裁去除点云噪点和无效区域，通过surfer软件生成裂缝面的三维形貌图，完成酸蚀裂缝面的数字化提取，如图3所示；其中图3(a)为裂缝面a经酸蚀后的三维形貌图，其中图3(b)为裂缝面b经酸蚀后的三维形貌图。
按照本申请所公开的：表面A和表面B构成的复杂空间体等效为一个光滑平面A’和一个复杂曲面B’所构成的相对规则的空间体，这一复杂空间体等效方式，将酸蚀后的裂缝面a和裂缝面b构成的复杂空间体等效成如图4所示的复杂空间体，然后根据式(23)计算该等效复杂空间体的曲率为0.0856mm-1；因此酸蚀后的裂缝面a和裂缝面b形成的复杂渗流空间体的曲率为0.0856mm-1。
以上详细描述了本发明的优选实施方式，但是，本发明并不限于上述实施方式中的具体细节，在本发明的技术构思范围内，可以对本发明的技术方案进行多种简单变型，这些简单变型均属于本发明的保护范围。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (9)

1.一种适于裂缝复杂渗流空间体的曲率表征方法，其特征在于，包括以下步骤：
获取岩心裂缝的表面A和表面B的三维数据；
将所述表面A和表面B构成的复杂空间体等效为一个光滑平面A’和一个复杂曲面B’所构成的复杂空间体，得到复杂曲面B’的三维数据；
取所述复杂曲面B’上的任一点G(i，j)，并以所述点G(i，j)为中心建立7×7标准网格，并根据曲率表达式，获得点G(i，j)的曲率k(i,j)，进而获得所述复杂曲面B’的三维数据每一点的曲率；其中，i代表网格图的行数，j代表网格图的列数；
根据所述复杂曲面B’中每一三维数据点的曲率得到所述复杂曲面B’的平均曲率，即为所述表面A和表面B构成的复杂空间体的平均曲率。
2.根据权利要求1所述的适于裂缝复杂渗流空间体的曲率表征方法，其特征在于，所述表面A和表面B构成的复杂空间体等效为一个光滑平面A’和一个复杂曲面B’所构成的复杂空间体具体包括以下步骤：
以所述表面A和所述表面B的最低点所在的面为基准面，获得获取所述表面A的三维数据点(xa、ya、za)和所述表面B(xa、ya、zb)的三维数据点；
以所述表面A和所述表面B上对应点的za和zb的差值等效为复杂曲面B’对应的三维数据点的z值，x＝xa，y＝ya，即可获得复杂曲面B’对应的三维数据点，进而获得复杂曲面B’的三维数据。
3.根据权利要求1所述的适于裂缝复杂渗流空间体的曲率表征方法，其特征在于，所述点G(i，j)的曲率k(i,j)均为该点在水平方向上的平均曲率kx(i,j)、竖直方向上的平均曲率ky(i,j)、斜向上45°方向上的平均曲率kd(i,j)、斜向下45°方向上的平均曲率ku(i,j)的平均值。
4.根据权利要求3所述的适于裂缝复杂渗流空间体的曲率表征方法，其特征在于，所述点G(i，j)在水平方向的平均曲率为点G(i，j)在水平方向上的曲率、点G(i-1，j)在水平方向上的曲率，及点G(i+1，j)在水平方向上的曲率的加权平均值；所述点G(i，j)在竖直方向的平均曲率为点G(i，j)在竖直方向上的曲率、点G(i，j-1)在竖直方向上的曲率，及点G(i，j+1)在竖直方向上的曲率的加权平均值；所述点G(i，j)在斜向上45°方向的平均曲率为点G(i，j)在斜向上45°方向上的曲率、点G(i-1，j-1)在斜向上45°方向上的曲率，及点G(i+1，j+1)在斜向上45°方向上曲率的加权平均值；所述点G(i，j)在斜向下45°方向的平均曲率为点G(i，j)在斜向下45°方向上的曲率、点G(i+1，j-1)在斜向下45°方向上的曲率，及点G(i-1，j+1)在斜向下45°方向的曲率的加权平均值。
5.根据权利要求4所述的适于裂缝复杂渗流空间体的曲率表征方法，其特征在于，所述点G(i，j)在水平方向/竖直方向/斜向下45°方向/斜向上45°方向的曲率的权重值均为2/4；所述点G(i-1，j)、点G(i，j-1)、点G(i+1，j-1)、点G(i-1，j-1)分别在水平方向/竖直方向/斜向下45°方向/斜向上45°方向的曲率的权重值均为1/4；所述点G(i+1，j)、点G(i，j+1)、点G(i-1，j+1)、点G(i+1，j+1)分别在水平方向/竖直方向/斜向下45°方向/斜向上45°方向的曲率的权重值均为1/4。
6.根据权利要求4所述的适于裂缝复杂渗流空间体的曲率表征方法，其特征在于，所述点G(i，j)在水平方向/竖直方向/斜向下45°方向/斜向上45°方向的曲率根据所述曲率表达式计算；所述点G(i-1，j)、点G(i，j-1)、点G(i+1，j-1)、点G(i-1，j-1)分别在水平方向/竖直方向/斜向下45°方向/斜向上45°方向的曲率根据所述曲率表达式计算；所述点G(i+1，j)、点G(i，j+1)、点G(i-1，j+1)、点G(i+1，j+1)分别在水平方向/竖直方向/斜向下45°方向/斜向上45°方向的曲率根据所述曲率表达式计算。
7.根据权利要求6所述的适于裂缝复杂渗流空间体的曲率表征方法，其特征在于，根据所述曲率表达式计算所述点G(i，j)在水平方向/竖直方向/斜向下45°方向/斜向上45°方向的曲率，根据所述曲率表达式计算点G(i-1，j)、点G(i，j-1)、点G(i+1，j-1)、点G(i-1，j-1)分别在水平方向/竖直方向/斜向下45°方向/斜向上45°方向的曲率，根据所述曲率表达式计算点G(i+1，j)、点G(i，j+1)、点G(i-1，j+1)、点G(i+1，j+1)分别在水平方向/竖直方向/斜向下45°方向/斜向上45°方向的曲率时，在7×7标准网格中对一阶求导和二阶求导中的差分计算分别赋权重值；
其中，所述点G(i，j)在水平方向的曲率在进行一阶求导和二阶求导时采用三阶差分法，每阶求导分为三个部分：第一部分为点G(i，j)左侧第一点、点G(i，j)右侧第一点，该部分权重为3/6；第二部分为点G(i，j)左侧第二点、点G(i，j)右侧第二点，该部分权重分别为2/6；第三部分为点G(i，j)左侧第三点、点G(i，j)右侧第三点，该部分权重分别为1/6；
其中，所述点G(i，j)在竖直方向的曲率在进行一阶求导和二阶求导时采用三阶差分法，每阶求导分为三个部分：第一部分为点G(i，j)上方第一点、点G(i，j)下方第一点，该部分权重为3/6；第二部分为点G(i，j)上方第二点、点G(i，j)下方第二点，该部分权重分别为2/6；第三部分为点G(i，j)上方第三点、点G(i，j)下方第三点，该部分权重分别为1/6；
其中，所述点G(i，j)在斜向下45°方向上的曲率在进行一阶求导和二阶求导时采用二阶差分法，每阶求导分为两个部分：第一部分为点G(i，j)左上第一点、点G(i，j)右下第一点，该部分权重为2/3；第二部分为G(i，j)左上第二点、点G(i，j)右下第二点，该部分权重为1/3；
其中，所述点G(i，j)在斜向上45°方向上的曲率在进行一阶求导和二阶求导时采用二阶差分法，每阶求导分为两个部分：第一部分为点G(i，j)右上第一点、点G(i，j)左下第一点，该部分权重为2/3；第二部分为G(i，j)右上第二点、点G(i，j)左下第二点，该部分权重为1/3。
8.根据权利要求1所述的适于裂缝复杂渗流空间体的曲率表征方法，其特征在于，所述曲率表达式为：
Figure FDA0003055066570000031
其中，x、y为空间任一点的坐标；
θ—表示二维曲线的某点在曲线上的倾角；
s—表示这一点的弧长。
9.根据权利要求1所述的适于裂缝复杂渗流空间体的曲率表征方法，其特征在于：所述表面A和表面B构成的复杂空间体的平均曲率为所述复杂曲面B’中每一三维数据点的平均曲率之和与所述复杂曲面B’的面积的比值；其具体表达式为：
Figure FDA0003055066570000041
其中，m-表示复杂曲面B’的长度；n-表示复杂曲面B’的宽度；(i，j)表示曲面B’中的任一点，k(i，j)-表示曲面B’中任一点的三维曲率。
CN202011068616.4A 2020-10-09 2020-10-09 一种适于裂缝复杂渗流空间体的曲率表征方法 Active CN112302639B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011068616.4A CN112302639B (zh) 2020-10-09 2020-10-09 一种适于裂缝复杂渗流空间体的曲率表征方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011068616.4A CN112302639B (zh) 2020-10-09 2020-10-09 一种适于裂缝复杂渗流空间体的曲率表征方法

Publications (2)

Publication Number Publication Date
CN112302639A CN112302639A (zh) 2021-02-02
CN112302639B true CN112302639B (zh) 2021-08-03

Family

ID=74487962

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011068616.4A Active CN112302639B (zh) 2020-10-09 2020-10-09 一种适于裂缝复杂渗流空间体的曲率表征方法

Country Status (1)

Country Link
CN (1) CN112302639B (zh)

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6662111B2 (en) * 2001-06-29 2003-12-09 Exxonmobil Upstream Research Company Method for analyzing reflection curvature in seismic data volumes
US9816365B2 (en) * 2013-08-23 2017-11-14 Halliburton Energy Services, Inc. Fracturing treatments in subterranean formations using reducible materials
CN104422959A (zh) * 2013-09-05 2015-03-18 中国石油化工股份有限公司 一种检测储层边界的曲率属性方法
CN106154326A (zh) * 2016-08-18 2016-11-23 中国石油天然气集团公司 一种纵弯褶皱裂缝密度评价的方法及装置
CN107843923A (zh) * 2016-09-19 2018-03-27 中国石油化工股份有限公司 基于地震属性体的曲率分析方法
CN110954958A (zh) * 2018-09-26 2020-04-03 中国石油化工股份有限公司 裂缝与断层的预测方法及系统

Also Published As

Publication number Publication date
CN112302639A (zh) 2021-02-02

Similar Documents

Publication Publication Date Title
CN107146280B (zh) 一种基于切分的点云建筑物重建方法
CN103955939B (zh) 三维扫描系统中点云拼接用边界特征点配准方法
CN105136054B (zh) 基于地面三维激光扫描的构筑物精细变形监测方法及系统
CN112257722B (zh) 基于抗差非线性高斯-赫尔默特模型的点云拟合方法
CN112669359B (zh) 一种三维点云配准方法、装置、设备及存储介质
CN112017220A (zh) 一种基于抗差约束最小二乘算法的点云精确配准方法
CN105066912B (zh) 酸刻蚀物理模拟实验中岩板表面扫描数据的步长标定方法
CN110060342A (zh) 一种三维曲面拟合方法
CN109871818B (zh) 基于法向量分布直方图和协方差描述子的人脸识别方法
CN114037706A (zh) 一种基于三维点云模型的预制梁段表面平整度检测方法
CN111476900B (zh) 一种基于Voronoi图和高斯分布的离散裂缝网络模型构建方法
CN112302639B (zh) 一种适于裂缝复杂渗流空间体的曲率表征方法
CN108334692A (zh) 预测增材制造零件变形的方法
Li et al. Differential evolution algorithm-based range image registration for free-form surface parts quality inspection
CN106495089B (zh) 湿法刻蚀石英晶体少量晶面获取全晶面刻蚀速率的方法
CN107945272A (zh) 一种基于高清钻孔彩电的岩体结构面的搜索方法
CN108733621A (zh) 基于bfgs算法的概率积分模型中参数的反演方法
CN115564820B (zh) 基于贪婪投影三角化的体积确定方法、系统、设备及介质
CN114581556B (zh) 一种区域地质调查中的数字填图方法
CN116051771A (zh) 一种基于无人机倾斜摄影模型的光伏bim屋顶自动建模方法
CN116309026A (zh) 一种基于统计局部特征描述与匹配的点云配准方法及系统
CN113963114A (zh) 一种基于多边形生长的离散边界点追踪方法
CN108875163B (zh) 一种评价三维裂缝网络连通性的方法和系统
Li et al. Building boundaries extaction from point clouds using dual-threshold Alpha Shapes
CN113933909B (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