CN115930847B - 一种三维结构面粗糙度评价指标的定量确定方法 - Google Patents
一种三维结构面粗糙度评价指标的定量确定方法 Download PDFInfo
- Publication number
- CN115930847B CN115930847B CN202211215898.5A CN202211215898A CN115930847B CN 115930847 B CN115930847 B CN 115930847B CN 202211215898 A CN202211215898 A CN 202211215898A CN 115930847 B CN115930847 B CN 115930847B
- Authority
- CN
- China
- Prior art keywords
- structural surface
- evaluation index
- calculating
- roughness
- sampling interval
- 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
Landscapes
- Length Measuring Devices By Optical Means (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
一种三维结构面粗糙度评价指标的定量确定方法,包括如下步骤:将岩石结构面切割成试样;采用激光扫描仪对结构面进行扫描,获取其数字高程分布图;计算或直接从激光扫描仪上获取扫描的取样间距d;计算结构面数字高程分布的标准差σ;通过沿垂直x轴分别取N条长度为L毫米的切线(N≥3),对于每一条切线分别计算结构面数字高程分布的相关距离τ;通过得到的取样间距、标准差、相关距离定量计算三维结构面粗糙度的评价指标仅需进行代数运算即可评价和比较结构面的粗糙度,避免了繁琐的测量和计算,为三维结构面粗糙度评价指标的评价和比较提供了理论与技术支撑。
Description
技术领域
本发明属于地质、水利、矿山、交通等岩土工程技术领域,具体涉及一种三维结构面粗糙度评价指标的定量确定方法,尤其适用于结构面发育区域的结构面粗糙度评价指标的定量计算。
背景技术
结构面粗糙度对岩体物理力学性质以及水文力学特征具有显著影响。Barton提出的十条结构面粗糙度曲线作为结构面粗糙度评估的基础,被广泛应用于工程实践和科学研究中,但是该方法存在人为估计的主观缺陷,且该方法主要用于结构面二维粗糙度的确定,而在工程实践中,由于结构面存在各向异性,不同的取样曲线往往会得到不同的粗糙度指标,使得结果具有一定的局限性。因此,确定一种三维结构面粗糙度评价指标的定量确定方法是岩体工程领域亟待突破的工程和科研瓶颈。
现有技术中,例如RS表征法用结构面表面实际面积与投影面积的比值作为粗糙度评价指标,其值主要取决于结构面形貌特征,但是无法反映取样间距对该评价指标的影响;法用剪切方向有效接触面积比/>和剪切方向倾角θ*来表征结构面粗糙度指标,该方法具有严密的数学推理,物理意义明确,但是在概念和确定方法上较为繁琐,且无法考虑取样间距对该指标的影响。
发明内容
本发明要解决的是现有技术中的三维结构面粗糙度评价指标的定量确定方法存在的概念和确定方法上较为繁琐,且无法考虑取样间距对该指标的影响的问题。
针对上述问题,本发明采用的技术方案是:一种三维结构面粗糙度评价指标的定量确定方法,包括如下步骤:
A、将岩石结构面切割成尺寸为L*L毫米的正方形块状结构面试样;
B、采用激光扫描仪对结构面进行扫描,获取其数字高程分布图;
C、计算或直接从激光扫描仪上获取扫描的取样间距d;
D、计算结构面数字高程分布的标准差σ;
E、通过沿垂直x轴分别取N条L毫米的切线(N≥3),对于每一条切线分别计算结构面数字高程分布的相关距离τ;
F、通过得到的取样间距、标准差、相关距离定量计算三维结构面粗糙度的评价指标
进一步的,所述步骤A中试样为标准试样时,尺寸为100毫米*100毫米。
进一步的,所述步骤C中计算或直接从激光扫描仪上获取扫描的取样间距d的具体方法为:
d=xi+1-xi=yi+1-yi (11)
其中xi+1和xi分别表示第i+1和第i个点的横坐标,yi+1和yi分别表示第i+1和第i个点的纵坐标;或者直接从激光扫描仪的说明手册上获取该扫描仪的扫描取样间距。
进一步的,所述步骤D中计算结构面数字高程分布的标准差σ的具体方法为:
n=(L/d+1)2 (12)
其中n表示激光扫描获取的总点数,L表示试样的尺寸,对于标准试样为100毫米,表示结构面高程平均值,zi表示结构面各点的高程。
进一步的,所述步骤E中计算结构面数字高程分布的相关距离τ的具体方法为:
沿垂直x轴分别取N条长度为L的切线N≥3),对于每一条切线进行计算:
Rz(t)=E[zy*zy+t] (15)
其中zy和zy+t表示切线上第y和y+t个点的高程,t为相对距离,为取样间距d的整数倍,为切线上所有点高程的平均值,Rz(t)表示zy*zy+t的均值,ρz(t)表示当相对距离为t时的相关系数;通过改变t的值即可得到不同相对距离条件下的相关系数;
对获取的(t,ρz(t))散点,采用相关系数函数进行拟合;
拟合结果τ即为该切线下的相关距离,对N条切线下的相关距离求平均值即为该结构面的相关距离。
进一步的,所述步骤E中相关系数函数包括指数型相关系数函数和高斯型相关系数函数ρ(t),分别为:
ρ(t)=exp(-2*|t/τ|) (17)
ρ(t)=exp(-2*t/τ)2 (18)
*表示乘号。
进一步的,所述步骤F通过得到的取样间距、标准差、相关距离定量计算三维结构面粗糙度评价指标的方法为包括非标准试样计算方法和标准试样计算方法。
进一步的,所述非标准试样计算方法的计算公式为:
其中a,b,c,e为待拟合参数,*表示乘号。
进一步的,所述标准试样计算方法的计算公式为:
式中*表示乘号。
本发明的有益效果和特点是:
1、本发明的三维结构面粗糙度评价指标的定量确定方法,考虑了取样设备特征、结构面高程分布特征、结构面形貌特征三个方面因素的影响,给出的计算公式仅需进行代数运算即可评价和比较结构面的粗糙度,避免了繁琐的测量和计算,为三维结构面粗糙度评价指标的评价和比较提供了理论与技术支撑。
2、本发明的三维结构面粗糙度评价指标的定量确定方法,具有物理意义明确、操作简单、计算精度高、适用性强,可用于结构面粗糙度的评价和不同结构面粗糙度的比较。
附图说明
图1为本发明较佳实施例的粗糙度评价方法的流程图;
图2为本发明较佳实施例的100毫米*100毫米的岩石结构面标准试样;
图3为本发明较佳实施例结构面的数字高程分布图;
图4为本发明较佳实施例结构面高程统计直方图。
图5为本发明较佳实施例结构面相关距离的拟合结果。
具体实施方式
下面将结合附图对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明同时考虑结构面分布的高程特征和取样间距特征,提出了一种三维结构面粗糙度评价指标的定量确定方法;具体实施例方案如下:
请参考图1-2,本发明涉及一种三维结构面粗糙度的定量计算方法,通过激光扫描仪获取其数字高程分布图后,计算取样间距、标准差和相关距离,并通过代数运算即可获取结构面粗糙度的评价指标。
下面结合附图对本发明进行进一步说明:
一种三维结构面粗糙度评价指标的定量确定方法,包括如下步骤:
A、将岩石结构面切割成尺寸为L*L毫米的正方形块状结构面试样;当试样为标准试样时,尺寸为100毫米*100毫米。
B、采用激光扫描仪对结构面进行扫描,获取其数字高程分布图;
C、计算或直接从激光扫描仪上获取扫描的取样间距d;所述步骤C中计算扫描的取样间距d的具体方法为:
d=xi+1-xi=yi+1-yi (21)
其中xi+1和xi分别表示第i+1和第i个点的横坐标,yi+1和yi分别表示第i+1和第i个点的纵坐标;或者直接从激光扫描仪的说明手册上获取该扫描仪的扫描取样间距。
D、计算结构面数字高程分布的标准差σ具体方法为:
n=(L/d+1)2 (22)
其中n表示激光扫描获取的总点数,L表示试样的尺寸,对于标准试样为100毫米,表示结构面高程平均值,zi表示结构面各点的高程。
E、通过沿垂直x轴分别取N条长度为L毫米的切线,对于每一条切线分别计算结构面数字高程分布的相关距离τ;具体方法为:
沿垂直x轴分别取N条长度为L的切线(N≥3),对于每一条切线进行计算:
Rz(t)=E[zy*zy+t] (25)
其中zy和zy+t表示切线上第y和y+t个点的高程;t为相对距离,为取样间距d的整数倍;Rz(t)表示zy*zy+t的均值,ρz(t)表示当相对距离为t时的相关系数;通过改变t的值即可得到不同相对距离条件下的相关系数;
对获取的(t,ρz(t))散点,采用相关系数函数进行拟合;相关系数函数包括指数型相关系数函数和高斯型相关系数函数ρ(t),分别为:
ρ(t)=exp(-2*|t/τ|) (27)
ρ(t)=exp(-2*t/τ)2 (28)
式中*表示乘号。
拟合结果τ即为该切线下的相关距离,对N条切线下的相关距离求平均值即为该结构面的相关距离。
F、通过得到的取样间距、标准差、相关距离定量计算三维结构面粗糙度的评价指标具体包括非标准试样计算方法和标准试样计算方法。
其中,非标准试样计算方法的计算公式为:
其中a,b,c,e为待拟合参数,*表示乘号。
其中,标准试样计算方法的计算公式为:
式中*表示乘号。
下面提供了一个具体算例:其实施流程如图1所示,具体实施步骤是:
A、将岩石结构面切割成100毫米*100毫米的标准试样,如图2所示;
B、采用激光扫描仪对结构面进行扫描,获取其数字高程分布图,如图3所示;
C、计算或直接从激光扫描仪上获取扫描的取样间距d;
d=xi+1-xi=yi+1-yi (31)
其中xi+1和xi分别表示第i+1和第i个点的横坐标,yi+1和yi分别表示第i+1和第i个点的纵坐标;或者直接从激光扫描仪的说明手册上获取该扫描仪的扫描取样间距。
为验证本发明在不同取样间距下的有效性,取样间距分别设置为0.5毫米,1.0毫米,2.0毫米,4.0毫米;
D、计算结构面数字高程分布的标准差σ;
n=(L/d+1)2 (32)
其中n表示激光扫描获取的总点数,L表示试样的尺寸,对于标准试样为100毫米,表示结构面高程平均值,zi表示结构面各点的高程;其高程统计直方图如图4所示;
E、计算结构面数字高程分布的相关距离τ;
沿垂直x轴分别取N条100毫米的切线(N=3),对于每一条切线进行计算:
Rz(t)=E[zy*zy+t] (5)
其中zy和zy+t表示切线x上第y和y+t个点的高程,t为相对距离,为取样间距d的整数倍,为切线上所有点高程的平均值,Rz(t)表示zy*zy+t的均值,ρz(t)表示当相对距离为t时的相关系数;通过改变t的值即可得到不同相对距离条件下的相关系数。
对获取的(t,ρz(t))散点,采用相关系数函数进行拟合,常用的相关系数函数有指数型相关系数函数和高斯型相关系数函数,分别为:
指数型ρ(t)=exp(-2*|t/τ|) (7)
高斯型ρ(t)=exp(-2*|t/τ|)2 (8)
拟合结果τ即为该切线下的相关距离,对N条切线下的相关距离求平均值即为该结构面的相关距离;
本次采用指数型相关系数函数进行拟合,三条切线的拟合结果如图5所示,三次拟合所得的相关距离分别为10.5毫米、13.5毫米、9.8毫米,并对三次的结果取平均值得最终的相关距离为11.27毫米。
F、将取样间距、标准差、相关距离代入公式计算三维结构面粗糙度的评价指标
计算结果如表1所示,将本发明的粗糙度评价指标与传统方法的粗糙度评价指标Z2S进行比较,其中Z2S的计算公式为:
表1:采用本发明计算结果与采用传统方法计算结果比较
结果表明,采用本发明所提出方法计算所得的粗糙度评价指标与采用传统方法所得的结果一致,进而验证了本发明的正确性,而本发明的计算公式更为方便简洁,从而体现了本发明的优越性。
以上显示和描述了本发明的基本原理和主要特征及优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的结构关系及原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。
Claims (8)
1.一种三维结构面粗糙度评价指标的定量确定方法,其特征在于:包括如下步骤:
A、将岩石结构面切割成尺寸为L*L毫米的正方形块状结构面试样;
B、采用激光扫描仪对结构面进行扫描,获取其数字高程分布图;
C、计算或直接从激光扫描仪上获取扫描的取样间距d;
D、计算结构面数字高程分布的标准差σ;
E、通过沿垂直x轴分别取N条长度为L毫米的切线,对于每一条切线分别计算结构面数字高程分布的相关距离τ;计算结构面数字高程分布的相关距离τ的具体方法为:
沿垂直x轴分别取N条长度为L的切线,N≥3,对于每一条切线进行计算:
Rz(t)=E[zy*zy+t] (1)
其中zy和zy+t表示切线上第y和y+t个点的高程;t为相对距离,为取样间距d的整数倍;为切线上所有点高程的平均值,Rz(t)表示zy*zy+t的均值,E表示取均值,ρz(t)表示当相对距离为t时的相关系数;通过改变t的值即可得到不同相对距离条件下的相关系数;
对获取的(t,ρz(t))散点,采用相关系数函数进行拟合;
拟合结果τ即为该切线下的相关距离,对N条切线下的相关距离求平均值即为该结构面的相关距离;
F、通过得到的取样间距、标准差、相关距离定量计算三维结构面粗糙度的评价指标
2.根据权利要求1所述的三维结构面粗糙度评价指标的定量确定方法,其特征在于:所述步骤A中试样为标准试样时,尺寸为100毫米*100毫米。
3.根据权利要求1所述的三维结构面粗糙度评价指标的定量确定方法,其特征在于:所述步骤C中计算扫描的取样间距d的具体方法为:
d=xi+1-xi=yi+1-yi (3)
其中xi+1和xi分别表示第i+1和第i个点的横坐标,yi+1和yi分别表示第i+1和第i个点的纵坐标;或者直接从激光扫描仪的说明手册上获取该扫描仪的扫描取样间距。
4.根据权利要求1所述的三维结构面粗糙度评价指标的定量确定方法,其特征在于:所述步骤D中计算结构面数字高程分布的标准差σ的具体方法为:
n=(L/d+1)2 (4)
其中n表示激光扫描获取的总点数,L表示试样的尺寸,对于标准试样为100毫米,表示结构面高程平均值,zi表示结构面各点的高程。
5.根据权利要求4所述的三维结构面粗糙度评价指标的定量确定方法,其特征在于:所述相关系数函数包括指数型相关系数函数和高斯型相关系数函数ρ(t),分别为:
ρ(t)=exp(-2*|t/τ|) (7)
ρ(t)=exp(-2*|t/τ|)2 (8)
式中*表示乘号。
6.根据权利要求1所述的三维结构面粗糙度评价指标的定量确定方法,其特征在于:所述步骤F通过得到的取样间距、标准差、相关距离定量计算三维结构面粗糙度评价指标的方法为包括非标准试样计算方法和标准试样计算方法。
7.根据权利要求6所述的三维结构面粗糙度评价指标的定量确定方法,其特征在于:所述非标准试样计算方法的计算公式为:
其中a,b,c,e为待拟合参数,*表示乘号。
8.根据权利要求6所述的三维结构面粗糙度评价指标的定量确定方法,其特征在于:所述标准试样计算方法的计算公式为:
式中*表示乘号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211215898.5A CN115930847B (zh) | 2022-09-30 | 2022-09-30 | 一种三维结构面粗糙度评价指标的定量确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211215898.5A CN115930847B (zh) | 2022-09-30 | 2022-09-30 | 一种三维结构面粗糙度评价指标的定量确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115930847A CN115930847A (zh) | 2023-04-07 |
CN115930847B true CN115930847B (zh) | 2023-09-22 |
Family
ID=86649662
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211215898.5A Active CN115930847B (zh) | 2022-09-30 | 2022-09-30 | 一种三维结构面粗糙度评价指标的定量确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115930847B (zh) |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0183270A2 (en) * | 1984-11-30 | 1986-06-04 | Kawasaki Steel Corporation | Method of determining glossinesses of surface of body |
JP2011069680A (ja) * | 2009-09-25 | 2011-04-07 | Ngk Insulators Ltd | 表面粗さ測定装置及び表面粗さ測定方法 |
CN104677313A (zh) * | 2015-02-16 | 2015-06-03 | 太原理工大学 | 多角度智能化岩体结构面粗糙度系数测量仪及其测量方法 |
CN104834806A (zh) * | 2015-03-09 | 2015-08-12 | 绍兴文理学院 | 结构面粗糙度系数尺寸效应取样代表性评价方法 |
CN106482674A (zh) * | 2016-11-14 | 2017-03-08 | 绍兴文理学院 | 基于中智数函数的结构面粗糙度尺寸效应的近似表达方法 |
CN107036905A (zh) * | 2017-03-28 | 2017-08-11 | 中国地质大学(武汉) | 一种岩体结构面二维粗糙度评价方法及系统 |
CN107563087A (zh) * | 2017-09-13 | 2018-01-09 | 绍兴文理学院 | 最优采样间距条件下的结构面粗糙度系数统计方法 |
CN107656902A (zh) * | 2017-09-13 | 2018-02-02 | 绍兴文理学院 | 不同采样间距下的结构面粗糙度系数统计方法 |
CN107843215A (zh) * | 2017-09-13 | 2018-03-27 | 绍兴文理学院 | 基于最优采样间距条件下的粗糙度系数分形评价模型构建方法 |
CN109543236A (zh) * | 2018-10-29 | 2019-03-29 | 绍兴文理学院 | 基于变异系数级比分析的岩体结构面粗糙度统计样本数确定方法 |
CN110360963A (zh) * | 2019-07-22 | 2019-10-22 | 中国科学院武汉岩土力学研究所 | 一种钻孔摄像技术的岩体结构面三维粗糙度描述方法 |
CN114166656A (zh) * | 2021-05-25 | 2022-03-11 | 中国科学院武汉岩土力学研究所 | 基于三维形貌参数建立节理抗剪强度弱化本构模型的方法 |
CN114216795A (zh) * | 2021-11-10 | 2022-03-22 | 中国科学院武汉岩土力学研究所 | 一种岩体结构破坏模式的快速识别方法 |
CN114396895A (zh) * | 2021-12-20 | 2022-04-26 | 河海大学 | 一种隧洞衬砌混凝土管片表面粗糙度的测量方法 |
-
2022
- 2022-09-30 CN CN202211215898.5A patent/CN115930847B/zh active Active
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0183270A2 (en) * | 1984-11-30 | 1986-06-04 | Kawasaki Steel Corporation | Method of determining glossinesses of surface of body |
JP2011069680A (ja) * | 2009-09-25 | 2011-04-07 | Ngk Insulators Ltd | 表面粗さ測定装置及び表面粗さ測定方法 |
CN104677313A (zh) * | 2015-02-16 | 2015-06-03 | 太原理工大学 | 多角度智能化岩体结构面粗糙度系数测量仪及其测量方法 |
CN104834806A (zh) * | 2015-03-09 | 2015-08-12 | 绍兴文理学院 | 结构面粗糙度系数尺寸效应取样代表性评价方法 |
CN106482674A (zh) * | 2016-11-14 | 2017-03-08 | 绍兴文理学院 | 基于中智数函数的结构面粗糙度尺寸效应的近似表达方法 |
CN107036905A (zh) * | 2017-03-28 | 2017-08-11 | 中国地质大学(武汉) | 一种岩体结构面二维粗糙度评价方法及系统 |
CN107563087A (zh) * | 2017-09-13 | 2018-01-09 | 绍兴文理学院 | 最优采样间距条件下的结构面粗糙度系数统计方法 |
CN107656902A (zh) * | 2017-09-13 | 2018-02-02 | 绍兴文理学院 | 不同采样间距下的结构面粗糙度系数统计方法 |
CN107843215A (zh) * | 2017-09-13 | 2018-03-27 | 绍兴文理学院 | 基于最优采样间距条件下的粗糙度系数分形评价模型构建方法 |
CN109543236A (zh) * | 2018-10-29 | 2019-03-29 | 绍兴文理学院 | 基于变异系数级比分析的岩体结构面粗糙度统计样本数确定方法 |
CN110360963A (zh) * | 2019-07-22 | 2019-10-22 | 中国科学院武汉岩土力学研究所 | 一种钻孔摄像技术的岩体结构面三维粗糙度描述方法 |
CN114166656A (zh) * | 2021-05-25 | 2022-03-11 | 中国科学院武汉岩土力学研究所 | 基于三维形貌参数建立节理抗剪强度弱化本构模型的方法 |
CN114216795A (zh) * | 2021-11-10 | 2022-03-22 | 中国科学院武汉岩土力学研究所 | 一种岩体结构破坏模式的快速识别方法 |
CN114396895A (zh) * | 2021-12-20 | 2022-04-26 | 河海大学 | 一种隧洞衬砌混凝土管片表面粗糙度的测量方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115930847A (zh) | 2023-04-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110780307B (zh) | 基于电瓶车车载式激光点云移动测量系统获取道路横断面的方法 | |
Martin et al. | Data reduction for reverse engineering | |
CN108871241B (zh) | 一种定量评价粗糙裂缝酸刻蚀后粗糙程度变化的方法 | |
CN108986048A (zh) | 基于线激光扫描三维点云快速复合滤波处理方法 | |
CN108444413B (zh) | 陶瓷墙地砖平整度检测装置及方法 | |
CN105066912B (zh) | 酸刻蚀物理模拟实验中岩板表面扫描数据的步长标定方法 | |
CN105571511B (zh) | 一种船舶外板成型精度在线检测方法 | |
CN111751392A (zh) | 一种基于双极化探地雷达的钢筋锈蚀检测方法 | |
CN106841115A (zh) | 基于空间频率特征的海面溢油检测方法 | |
CN110379011A (zh) | 一种基于改进三次b样条曲线的水下地形点云孔洞修补方法 | |
CN107036561B (zh) | 基于中智数函数的结构面粗糙度各向异性的近似表达方法 | |
CN106482674A (zh) | 基于中智数函数的结构面粗糙度尺寸效应的近似表达方法 | |
CN103218837B (zh) | 一种基于经验分布函数的不等组距直方图的绘制方法 | |
CN115930847B (zh) | 一种三维结构面粗糙度评价指标的定量确定方法 | |
CN112857315B (zh) | 基于三维激光扫描的钢管柱垂直度测量方法 | |
CN109543236B (zh) | 岩体结构面粗糙度统计样本数确定方法 | |
CN107870894B (zh) | 一种计算乘潮水位及其延时的条件概率方法 | |
CN105092452A (zh) | 一种机制砂球体类似度的检测方法 | |
CN104297753A (zh) | 一种基于自适应缩减算子的导航雷达图像反演海面风向方法 | |
CN113759387B (zh) | 一种基于三维激光雷达的海岸防浪建筑物越浪量测量方法 | |
CN115876661A (zh) | 一种基于三维探地雷达图像的沥青路面空隙率评价方法 | |
CN115711587A (zh) | 结构面三维粗糙度评价的剖面线条带最小样本数确定方法 | |
Hu et al. | Adaptive sampling method for laser measuring free-form surface | |
CN104121864A (zh) | 一种岩石错动裂缝面的间隙评价方法 | |
CN112906244A (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 |