CN109977589B - 不位于表层的等厚等强度叠层矿的储量确定方法 - Google Patents

不位于表层的等厚等强度叠层矿的储量确定方法 Download PDF

Info

Publication number
CN109977589B
CN109977589B CN201910277231.XA CN201910277231A CN109977589B CN 109977589 B CN109977589 B CN 109977589B CN 201910277231 A CN201910277231 A CN 201910277231A CN 109977589 B CN109977589 B CN 109977589B
Authority
CN
China
Prior art keywords
ore
laminated
medium
thickness
layer
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
CN201910277231.XA
Other languages
English (en)
Other versions
CN109977589A (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.)
Hunan Powerise Information Technology Co ltd
Central South University
Original Assignee
Hunan Powerise Information Technology Co ltd
Central South University
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 Hunan Powerise Information Technology Co ltd, Central South University filed Critical Hunan Powerise Information Technology Co ltd
Priority to CN201910277231.XA priority Critical patent/CN109977589B/zh
Publication of CN109977589A publication Critical patent/CN109977589A/zh
Application granted granted Critical
Publication of CN109977589B publication Critical patent/CN109977589B/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

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Length Measuring Devices Characterised By Use Of Acoustic Means (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)

Abstract

本发明公开了一种不位于表层的等厚等强度叠层矿的储量确定方法，其包括：获取等厚等强度叠层矿的矿层状态，以及叠层矿所处位置的上、下层介质的厚度和材料属性，其中矿层状态包括材料属性、层数和一层矿的厚度；根据叠层矿的矿层状态、上下层介质厚度以及上下层介质的材料属性，计算叠层矿的弧长；获取需测试的地表范围，并根据所述地表范围计算弧长的倍数；获取叠层矿的厚度，并计算叠层矿的储量，所述叠层矿的储量为弧长、倍数以及叠层矿的厚度三者的乘积。本发明可以简单便捷地计算出不位于表层的等厚等强度叠层矿的矿产储量。

Description

不位于表层的等厚等强度叠层矿的储量确定方法
技术领域
本发明涉及矿产勘查技术领域，尤其涉及一种不位于表层的等厚等强度叠层矿的储量确定方法。
背景技术
在地质研究中，矿产储量的估算十分重要，但是由于叠层矿褶皱后的矿床结构复杂，目前没有简单的针对该类矿床的储量估算方法。叠层矿存在等厚等强度叠层矿、不等厚不等强度叠层矿等多种情况，其中，等厚等强度叠层矿表示叠层矿中各矿层的厚度相等，弹性模量或粘度相同，且各矿层之间无摩擦。本发明仅针对不位于表层的等厚等强度叠层矿进行研究。
发明内容
本发明的主要目的在于提供一种不位于表层的等厚等强度叠层矿的储量确定方法，旨在解决目前没有针对褶皱后叠层矿的储量估算方法的问题。
为实现上述目的，本发明提供的一种不位于表层的等厚等强度叠层矿的储量确定方法，包括：
获取等厚等强度叠层矿的矿层状态，以及叠层矿所处位置的上、下层介质的厚度和材料属性，其中矿层状态包括材料属性、层数和一层矿的厚度；
根据叠层矿的矿层状态、上下层介质厚度以及上下层介质的材料属性，计算叠层矿的弧长；
获取需测试的地表范围，并根据所述地表范围计算弧长的倍数；
获取叠层矿的厚度，并计算叠层矿的储量，所述叠层矿的储量为弧长、倍数以及叠层矿的厚度三者的乘积。
优选地，所述根据叠层矿的矿层状态、上下层介质厚度以及上下层介质的材料属性，计算叠层矿的弧长的步骤包括：
当所述上层介质和所述下层介质的厚度有限相等，且二者均为偏弹性体，且弹性模量相等时，按照第一方程计算叠层矿的弧长，其中所述第一方程为：
Figure BDA0002020411680000021
S为弧长，h为上/下层介质的厚度，t为一层矿的厚度，n为叠层矿的层数，B为叠层矿的弹性模量，B0为上/下层介质的弹性模量。
优选地，所述根据叠层矿的矿层状态、上下层介质厚度以及上下层介质的材料属性，计算叠层矿的弧长的步骤包括：
当所述上层介质和所述下层介质的厚度有限相等，且二者均为偏弹性体，且弹性模量不相等时，按照第二方程计算叠层矿的弧长，其中所述第二方程为：
Figure BDA0002020411680000022
S为弧长，h为上/下层介质的厚度，t为一层矿的厚度，n为叠层矿的层数，B为叠层矿的弹性模量，B1为上层介质的弹性模量，B2为下层介质的弹性模量。
优选地，所述根据叠层矿的矿层状态、上下层介质厚度以及上下层介质的材料属性，计算叠层矿的弧长的步骤包括：
当所述上层介质和所述下层介质的厚度有限不相等，且二者均为偏弹性体，且弹性模量相等时，按照第三方程计算叠层矿的弧长，其中所述第三方程为：
Figure BDA0002020411680000023
S为弧长，h1为上层介质的厚度，h2为下层介质的厚度，t为一层矿的厚度，n为叠层矿的层数，B为叠层矿的弹性模量，B0为上/下层介质的弹性模量。
优选地，所述根据叠层矿的矿层状态、上下层介质厚度以及上下层介质的材料属性，计算叠层矿的弧长的步骤包括：
当所述上层介质和所述下层介质的厚度有限不相等，且二者均为偏弹性体，且弹性模量不相等时，按照第四方程计算叠层矿的弧长，其中所述第四方程为：
Figure BDA0002020411680000024
S为弧长，h1为上层介质的厚度，h2为下层介质的厚度，t为一层矿的厚度，n为叠层矿的层数，B为叠层矿的弹性模量，B1为上层介质的弹性模量，B2为下层介质的弹性模量。
优选地，所述根据叠层矿的矿层状态、上下层介质厚度以及上下层介质的材料属性，计算叠层矿的弧长的步骤包括：
当所述上层介质和所述下层介质的厚度有限相等，且二者均为偏粘滞性体，且粘度相等时，按照第五方程计算叠层矿的弧长，其中所述第五方程为：
Figure BDA0002020411680000031
S为弧长，h为上/下层介质的厚度，t为一层矿的厚度，n为叠层矿的层数，μ为叠层矿的粘度，μ0为上/下层介质的粘度。
优选地，所述根据叠层矿的矿层状态、上下层介质厚度以及上下层介质的材料属性，计算叠层矿的弧长的步骤包括：
当所述上层介质和所述下层介质的厚度有限相等，且二者均为偏粘滞性体，且粘度不相等时，按照第六方程计算叠层矿的弧长，其中所述第六方程为：
Figure BDA0002020411680000032
S为弧长，h为上/下层介质的厚度，t为一层矿的厚度，n为叠层矿的层数，μ为叠层矿的粘度，μ1为上层介质的粘度，μ2为下层介质的粘度。
优选地，所述根据叠层矿的矿层状态、上下层介质厚度以及上下层介质的材料属性，计算叠层矿的弧长的步骤包括：
当所述上层介质和所述下层介质的厚度有限不相等，且二者均为偏粘滞性体，且粘度相等时，按照第七方程计算叠层矿的弧长，其中所述第七方程为：
Figure BDA0002020411680000033
S为弧长，h1为上层介质的厚度，h2为下层介质的厚度，t为一层矿的厚度，n为叠层矿的层数，μ为叠层矿的粘度，μ0为上/下层介质的粘度。
优选地，所述根据叠层矿的矿层状态、上下层介质厚度以及上下层介质的材料属性，计算叠层矿的弧长的步骤包括：
当所述上层介质和所述下层介质的厚度有限不相等，且二者均为偏粘滞性体，且粘度不相等时，按照第八方程计算叠层矿的弧长，其中所述第八方程为：
Figure BDA0002020411680000034
S为弧长，h1为上层介质的厚度，h2为下层介质的厚度，t为一层矿的厚度，n为叠层矿的层数，μ为叠层矿的粘度，μ1为上层介质的粘度，μ2为下层介质的粘度。
本发明的原理和构思如下：
褶皱是野外广泛发育的构造现象之一，它是成层或面理化岩石由地壳运动而发生弯曲形成的，尽管它们在野外姿态百态，可以分为多种模型，但都是波状起伏的。如果我们假设岩层在受力变形褶皱到形成明显挠曲的褶皱过程中岩层本身缩短量不大，只是褶皱而已，则可以把任何时候的相当弧长近似当成初始波长，再应用初始波长公式进行计算是可以的，也就是说如果知道岩层及周围介质的弹性模量或粘度，就可以计算出最终弧长来。根据弧长和矿层厚度，那么在一定地表范围内就可以精确得出矿产储量。
本发明提出的技术方案中，先根据不同的边界条件，得出不同的弧长计算公式，再根据弧长计算公式计算获得的弧长和矿产厚度，地表范围，即可得到矿产储量，简单便捷。
附图说明
图1为本发明不位于表层的等厚等强度叠层矿的储量确定方法一实施例的流程示意图；
图2为不位于表层的等厚等强度叠层矿埋于一个厚度为h的介质中的褶皱示意图。
本发明目的的实现、功能特点及优点将结合实施例，参照附图做进一步说明。
具体实施方式
应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
请参照图1，图1为本发明一实施例提出的不位于表层的等厚等强度叠层矿的储量确定方法。在该实施例中，所述不位于表层的等厚等强度叠层矿的储量确定方法，包括：
步骤S100，获取等厚等强度叠层矿的矿层状态，以及叠层矿所处位置的上、下层介质的厚度和材料属性，其中矿层状态包括材料属性、层数和一层矿的厚度；
步骤S200，根据叠层矿的矿层状态、上下层介质厚度以及上下层介质的材料属性，计算叠层矿的弧长；
步骤S300，获取需测试的地表范围，并根据所述地表范围计算弧长的倍数；
步骤S400，获取叠层矿的厚度，并计算叠层矿的储量，所述叠层矿的储量为弧长、倍数以及叠层矿的厚度三者的乘积。
具体地，可以通过矿层对比图来获知叠层矿的层数、一层矿的厚度以及上下层介质的厚度，需要说明的是，矿层对比图是根据某一方向上的探矿工程资料(一般即为钻孔柱状图)，以缩小了比例尺的柱状图形式，将其排列在同一平面上，以说明沉积矿层的层位、层数、厚度、上下岩层、接触关系及它们的变化情况等。
材料属性主要为弹性模块或粘度，该属性可以通过数字岩心分析技术来获得。该方法首先利用扫描电子显微镜对矿层进行数字岩心成像，再使用二值图像分割法将数字岩心图像转换成表征空隙空间和颗粒结构的二值图，在二值图的基础上模拟岩石力学实验，计算弹性模量。
其中，所述根据叠层矿的矿层状态、上下层介质厚度以及上下层介质的材料属性，计算叠层矿的弧长的步骤包括：
当所述上层介质和所述下层介质的厚度有限相等，且二者均为偏弹性体，且弹性模量相等时，按照第一方程计算叠层矿的弧长，其中所述第一方程为：
Figure BDA0002020411680000057
S为弧长，h为上/下层介质的厚度，t为一层矿的厚度，n为叠层矿的层数，B为叠层矿的弹性模量，B0为上/下层介质的弹性模量，请参照图2。需要说明的是，本发明中同一符号表示同一含义。
具体地，第一方程的推导来源为：
总力矩∑BI＝nBI；
阻力矩
Figure BDA0002020411680000051
介质的弹性常数
Figure BDA0002020411680000052
临界载荷
Figure BDA0002020411680000053
令
Figure BDA0002020411680000054
解得：
Figure BDA0002020411680000055
以下通过实验验证了第一方程的正确性，实验材料：主层为橡胶，介质为海绵，盖层与基底为沙粘土，实验材料已知值见表1，测量数据及计算结果见表2。
表1实验材料已知值
Figure BDA0002020411680000056
Figure BDA0002020411680000061
表2测量数据及计算结果
Figure BDA0002020411680000062
表2验证了第一方程的正确性和实用性。
当所述上层介质和所述下层介质的厚度有限相等，且二者均为偏弹性体，且弹性模量不相等时，按照第二方程计算叠层矿的弧长，其中所述第二方程为：
Figure BDA0002020411680000063
S为弧长，h为上/下层介质的厚度，t为一层矿的厚度，n为叠层矿的层数，B为叠层矿的弹性模量，B1为上层介质的弹性模量，B2为下层介质的弹性模量。
当所述上层介质和所述下层介质的厚度有限不相等，且二者均为偏弹性体，且弹性模量相等时，按照第三方程计算叠层矿的弧长，其中所述第三方程为：
Figure BDA0002020411680000064
S为弧长，h1为上层介质的厚度，h2为下层介质的厚度，t为一层矿的厚度，n为叠层矿的层数，B为叠层矿的弹性模量，B0为上/下层介质的弹性模量。
具体地，第三方程的推导来源为：
介质的弹性常数
Figure BDA0002020411680000065
临界载荷
Figure BDA0002020411680000066
令
Figure BDA0002020411680000067
解得：
Figure BDA0002020411680000068
当所述上层介质和所述下层介质的厚度有限不相等，且二者均为偏弹性体，且弹性模量不相等时，按照第四方程计算叠层矿的弧长，其中所述第四方程为：
Figure BDA0002020411680000071
S为弧长，h1为上层介质的厚度，h2为下层介质的厚度，t为一层矿的厚度，n为叠层矿的层数，B为叠层矿的弹性模量，B1为上层介质的弹性模量，B2为下层介质的弹性模量。
当所述上层介质和所述下层介质的厚度有限相等，且二者均为偏粘滞性体，且粘度相等时，按照第五方程计算叠层矿的弧长，其中所述第五方程为：
Figure BDA0002020411680000072
S为弧长，h为上/下层介质的厚度，t为一层矿的厚度，n为叠层矿的层数，μ为叠层矿的粘度，μ0为上/下层介质的粘度。
当所述上层介质和所述下层介质的厚度有限相等，且二者均为偏粘滞性体，且粘度不相等时，按照第六方程计算叠层矿的弧长，其中所述第六方程为：
Figure BDA0002020411680000073
S为弧长，h为上/下层介质的厚度，t为一层矿的厚度，n为叠层矿的层数，μ为叠层矿的粘度，μ1为上层介质的粘度，μ2为下层介质的粘度。
当所述上层介质和所述下层介质的厚度有限不相等，且二者均为偏粘滞性体，且粘度相等时，按照第七方程计算叠层矿的弧长，其中所述第七方程为：
Figure BDA0002020411680000074
S为弧长，h1为上层介质的厚度，h2为下层介质的厚度，t为一层矿的厚度，n为叠层矿的层数，μ为叠层矿的粘度，μ0为上/下层介质的粘度。
以下通过实验验证了第七方程的正确性，实验材料为石膏、砂、粘土混合材料，实验材料已知值见表3，测量数据及计算结果见表4。
表3实验材料已知值
Figure BDA0002020411680000075
其中，主层厚度t，由于做模型过程中的误差，在0.45-0.7cm平均为0.65cm。(普遍较0.5cm大)。
表4测量数据及计算结果
Figure BDA0002020411680000081
从表4可以发现第七方程基本与实验结果相符合，误差在允许的范围内。
当所述上层介质和所述下层介质的厚度有限不相等，且二者均为偏粘滞性体，且粘度不相等时，按照第八方程计算叠层矿的弧长，其中所述第八方程为：
Figure BDA0002020411680000082
S为弧长，h1为上层介质的厚度，h2为下层介质的厚度，t为一层矿的厚度，n为叠层矿的层数，μ为叠层矿的粘度，μ1为上层介质的粘度，μ2为下层介质的粘度。
本发明提出的技术方案中，先根据不同的边界条件，得出不同的弧长计算公式，再根据弧长计算公式计算获得的弧长和矿产厚度，地表范围，即可得到矿产储量，简单便捷。
在本说明书的描述中，参考术语“一实施例”、“另一实施例”、“其他实施例”、或“第一实施例～第X实施例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料、方法步骤或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
需要说明的是，在本文中，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者装置不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者装置所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括该要素的过程、方法、物品或者装置中还存在另外的相同要素。
上述本发明实施例序号仅仅为了描述，不代表实施例的优劣。
通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到上述实施例方法可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件，但很多情况下前者是更佳的实施方式。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质(如ROM/RAM、磁碟、光盘)中，包括若干指令用以使得一台终端(可以是手机，计算机，服务器，空调器，或者网络设备等)执行本发明各个实施例所述的方法。
上面结合附图对本发明的实施例进行了描述，但是本发明并不局限于上述的具体实施方式，上述的具体实施方式仅仅是示意性的，而不是限制性的，本领域的普通技术人员在本发明的启示下，在不脱离本发明宗旨和权利要求所保护的范围情况下，还可做出很多形式，这些均属于本发明的保护之内。

Claims (9)

1.一种不位于表层的等厚等强度叠层矿的储量确定方法，其特征在于，包括：
获取等厚等强度叠层矿的矿层状态，以及叠层矿所处位置的上、下层介质的厚度和材料属性，其中矿层状态包括材料属性、层数和一层矿的厚度；
根据叠层矿的矿层状态、上下层介质厚度以及上下层介质的材料属性，计算叠层矿的弧长；
获取需测试的地表范围，并根据所述地表范围计算弧长的倍数；
获取叠层矿的厚度，并计算叠层矿的储量，所述叠层矿的储量为弧长、倍数以及叠层矿的厚度三者的乘积。
2.如权利要求1所述的不位于表层的等厚等强度叠层矿的储量确定方法，其特征在于，所述根据叠层矿的矿层状态、上下层介质厚度以及上下层介质的材料属性，计算叠层矿的弧长的步骤包括：
当所述上层介质和所述下层介质的厚度有限相等，且二者均为偏弹性体，且弹性模量相等时，按照第一方程计算叠层矿的弧长，其中所述第一方程为：
Figure FDA0002020411670000011
S为弧长，h为上/下层介质的厚度，t为一层矿的厚度，n为叠层矿的层数，B为叠层矿的弹性模量，B0为上/下层介质的弹性模量。
3.如权利要求1所述的不位于表层的等厚等强度叠层矿的储量确定方法，其特征在于，所述根据叠层矿的矿层状态、上下层介质厚度以及上下层介质的材料属性，计算叠层矿的弧长的步骤包括：
当所述上层介质和所述下层介质的厚度有限相等，且二者均为偏弹性体，且弹性模量不相等时，按照第二方程计算叠层矿的弧长，其中所述第二方程为：
Figure FDA0002020411670000012
S为弧长，h为上/下层介质的厚度，t为一层矿的厚度，n为叠层矿的层数，B为叠层矿的弹性模量，B1为上层介质的弹性模量，B2为下层介质的弹性模量。
4.如权利要求1所述的不位于表层的等厚等强度叠层矿的储量确定方法，其特征在于，所述根据叠层矿的矿层状态、上下层介质厚度以及上下层介质的材料属性，计算叠层矿的弧长的步骤包括：
当所述上层介质和所述下层介质的厚度有限不相等，且二者均为偏弹性体，且弹性模量相等时，按照第三方程计算叠层矿的弧长，其中所述第三方程为：
Figure FDA0002020411670000021
S为弧长，h1为上层介质的厚度，h2为下层介质的厚度，t为一层矿的厚度，n为叠层矿的层数，B为叠层矿的弹性模量，B0为上/下层介质的弹性模量。
5.如权利要求1所述的不位于表层的等厚等强度叠层矿的储量确定方法，其特征在于，所述根据叠层矿的矿层状态、上下层介质厚度以及上下层介质的材料属性，计算叠层矿的弧长的步骤包括：
当所述上层介质和所述下层介质的厚度有限不相等，且二者均为偏弹性体，且弹性模量不相等时，按照第四方程计算叠层矿的弧长，其中所述第四方程为：
Figure FDA0002020411670000022
S为弧长，h1为上层介质的厚度，h2为下层介质的厚度，t为一层矿的厚度，n为叠层矿的层数，B为叠层矿的弹性模量，B1为上层介质的弹性模量，B2为下层介质的弹性模量。
6.如权利要求1所述的不位于表层的等厚等强度叠层矿的储量确定方法，其特征在于，所述根据叠层矿的矿层状态、上下层介质厚度以及上下层介质的材料属性，计算叠层矿的弧长的步骤包括：
当所述上层介质和所述下层介质的厚度有限相等，且二者均为偏粘滞性体，且粘度相等时，按照第五方程计算叠层矿的弧长，其中所述第五方程为：
Figure FDA0002020411670000023
S为弧长，h为上/下层介质的厚度，t为一层矿的厚度，n为叠层矿的层数，μ为叠层矿的粘度，μ0为上/下层介质的粘度。
7.如权利要求1所述的不位于表层的等厚等强度叠层矿的储量确定方法，其特征在于，所述根据叠层矿的矿层状态、上下层介质厚度以及上下层介质的材料属性，计算叠层矿的弧长的步骤包括：
当所述上层介质和所述下层介质的厚度有限相等，且二者均为偏粘滞性体，且粘度不相等时，按照第六方程计算叠层矿的弧长，其中所述第六方程为：
Figure FDA0002020411670000031
S为弧长，h为上/下层介质的厚度，t为一层矿的厚度，n为叠层矿的层数，μ为叠层矿的粘度，μ1为上层介质的粘度，μ2为下层介质的粘度。
8.如权利要求1所述的不位于表层的等厚等强度叠层矿的储量确定方法，其特征在于，所述根据叠层矿的矿层状态、上下层介质厚度以及上下层介质的材料属性，计算叠层矿的弧长的步骤包括：
当所述上层介质和所述下层介质的厚度有限不相等，且二者均为偏粘滞性体，且粘度相等时，按照第七方程计算叠层矿的弧长，其中所述第七方程为：
Figure FDA0002020411670000032
S为弧长，h1为上层介质的厚度，h2为下层介质的厚度，t为一层矿的厚度，n为叠层矿的层数，μ为叠层矿的粘度，μ0为上/下层介质的粘度。
9.如权利要求1所述的不位于表层的等厚等强度叠层矿的储量确定方法，其特征在于，所述根据叠层矿的矿层状态、上下层介质厚度以及上下层介质的材料属性，计算叠层矿的弧长的步骤包括：
当所述上层介质和所述下层介质的厚度有限不相等，且二者均为偏粘滞性体，且粘度不相等时，按照第八方程计算叠层矿的弧长，其中所述第八方程为：
Figure FDA0002020411670000033
S为弧长，h1为上层介质的厚度，h2为下层介质的厚度，t为一层矿的厚度，n为叠层矿的层数，μ为叠层矿的粘度，μ1为上层介质的粘度，μ2为下层介质的粘度。
CN201910277231.XA 2019-04-08 2019-04-08 不位于表层的等厚等强度叠层矿的储量确定方法 Active CN109977589B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910277231.XA CN109977589B (zh) 2019-04-08 2019-04-08 不位于表层的等厚等强度叠层矿的储量确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910277231.XA CN109977589B (zh) 2019-04-08 2019-04-08 不位于表层的等厚等强度叠层矿的储量确定方法

Publications (2)

Publication Number Publication Date
CN109977589A CN109977589A (zh) 2019-07-05
CN109977589B true CN109977589B (zh) 2022-11-11

Family

ID=67083469

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910277231.XA Active CN109977589B (zh) 2019-04-08 2019-04-08 不位于表层的等厚等强度叠层矿的储量确定方法

Country Status (1)

Country Link
CN (1) CN109977589B (zh)

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11297561A (ja) * 1998-04-09 1999-10-29 Nec Corp 複合ペロブスカイト化合物を用いた積層セラミックコンデンサ
JP5862364B2 (ja) * 2012-02-24 2016-02-16 Tdk株式会社 誘電体積層薄膜
CN102967889A (zh) * 2012-11-09 2013-03-13 湖南创博龙智信息科技股份有限公司 一种矿床储量的测量方法
CN108957581B (zh) * 2017-05-22 2020-07-28 中国石油化工股份有限公司 一种分析页岩地层岩石物理各向异性的方法

Also Published As

Publication number Publication date
CN109977589A (zh) 2019-07-05

Similar Documents

Publication Publication Date Title
Leong et al. Universal soil shrinkage curve equation
CN104641259A (zh) 利用多相流表征单元体积确定的数字岩石分析系统与方法
Li et al. Fractal characteristics based on different statistical objects of process-based digital rock models
US20180120479A1 (en) Porosity prediction based on effective stress
Bradley Software for efficient static dislocation–traction calculations in fault simulators
CN110441209B (zh) 一种基于致密储层数字岩心计算岩石渗透率的方法
Liu et al. Reconstruction of 3-D digital cores using a hybrid method
Sain et al. How computational rock-physics tools can be used to simulate geologic processes, understand pore-scale heterogeneity, and refine theoretical models
Panteghini et al. A single numerically efficient equation for approximating the Mohr–Coulomb and the Matsuoka–Nakai failure criteria with rounded edges and apex
Nadimi et al. Numerical modelling of rough particle contacts subject to normal and tangential loading
CN109977589B (zh) 不位于表层的等厚等强度叠层矿的储量确定方法
Olafuyi et al. Spontaneous imbibition in small cores
Sweeney et al. Pore size distributions calculated from 3-D images of DEM-simulated powder compacts
CN110008596B (zh) 位于表层的等厚等强度叠层矿的储量确定方法
Xue et al. Effect of plate length on subduction kinematics and slab geometry: Insights from buoyancy‐driven analog subduction models
Fokker Subsidence prediction and inversion of subsidence data
US11879825B2 (en) Method for digitally characterizing the permeability of rock
Sun et al. Estimation of petrophysical parameters of heterogeneous carbonate rock sample with multi-scale CT images
WO2020131698A1 (en) Method for characterizing the porosity of rock
Liou et al. A new method developed for fractal dimension and topothesy varying with the mean separation of two contact surfaces
Al Balushi et al. Numerical generation of stress-dependent permeability curves
Isik et al. A shape parameter for soil particles using a computational method
Duhan et al. XFEM for multiphysics analysis of edge dislocations with nonuniform misfit strain: A novel enrichment implementation
CN111767650A (zh) 一种数字岩芯等效弹性参数估算方法及装置
Sajeva et al. Strain-accumulation mechanisms in sands under isotropic stress

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