CN116578895A - 一种串状溶洞地质钻孔探测与弹性波ct扫描的数据聚类分析方法 - Google Patents

一种串状溶洞地质钻孔探测与弹性波ct扫描的数据聚类分析方法 Download PDF

Info

Publication number
CN116578895A
CN116578895A CN202310407986.3A CN202310407986A CN116578895A CN 116578895 A CN116578895 A CN 116578895A CN 202310407986 A CN202310407986 A CN 202310407986A CN 116578895 A CN116578895 A CN 116578895A
Authority
CN
China
Prior art keywords
elastic wave
data
drilling
detection
wave
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.)
Pending
Application number
CN202310407986.3A
Other languages
English (en)
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.)
CCCC Fourth Harbor Engineering Co Ltd
First Engineering Co of CCCC Fourth Harbor Engineering Co Ltd
Original Assignee
CCCC Fourth Harbor Engineering Co Ltd
First Engineering Co of CCCC Fourth Harbor Engineering Co Ltd
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 CCCC Fourth Harbor Engineering Co Ltd, First Engineering Co of CCCC Fourth Harbor Engineering Co Ltd filed Critical CCCC Fourth Harbor Engineering Co Ltd
Priority to CN202310407986.3A priority Critical patent/CN116578895A/zh
Publication of CN116578895A publication Critical patent/CN116578895A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/23Clustering techniques
    • G06F18/232Non-hierarchical techniques
    • G06F18/2321Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
    • G06F18/23213Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
    • EFIXED CONSTRUCTIONS
    • E21EARTH DRILLING; MINING
    • E21BEARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00Testing 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
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment of water resources

Abstract

本发明公开了一种串状溶洞地质钻孔探测与弹性波CT检测的数据聚类分析方法,该方法包括:首先选择串状溶洞工程地点适当的位置进行钻孔,并将钻孔数据转化为不同地质对应的弹性波标准波速;随后,进行弹性波CT检测,获取弹性波CT扫描图像;提取图像信息中各点的CT值;将钻孔数据转化的弹性波标准波速混入提取的CT值中,获取混合数据;对混合数据中钻孔数据转化的弹性波标准波速进行标记处理,并对标记处理后的混合数据进行聚类分析;通过聚类分析的结果,判定钻孔结果和弹性波CT检测的结果的一致性;本发明将标记以后的钻孔探测数据和弹性波CT数据同时进行聚类分析;通过聚类分析结果,能够直观的分辨钻孔探测和弹性波CT检测的结果的一致性。

Description

一种串状溶洞地质钻孔探测与弹性波CT扫描的数据聚类分析 方法
技术领域
本发明涉及工程探测技术,具体而言,涉及一种串状溶洞地质钻孔探测与弹性波CT扫描的数据聚类分析方法。
背景技术
串状溶洞是一种地下岩溶地貌,通常由多个相互连接的小型洞穴组成,形成了复杂的三维空间结构;为了更好地了解该地貌的内部结构和特征,需要进行地质钻孔探测。
钻孔探测是一种常见的地质勘探方法,通过钻取地下岩石或土壤样本,获取其物理性质和化学成分等信息,以了解地下构造、矿产资源、工程地质情况等;弹性波CT检测是一种基于声波传播原理的无损检测技术,可以对地下介质的物理性质进行非破坏性测试,并获得高分辨率的图像数据;与传统的X射线CT检测相比,弹性波CT检测具有更好的局部探测能力,可以在不破坏地下介质的情况下获取更为详细的信息。
将钻孔探测和弹性波CT检测相结合,对工程地点的地质进行探测;并通过分析,将钻孔探测和弹性波CT的分析结果结合起来,验证弹性波CT检测是否和钻孔探测一致;分析方法通常是通过单独的对弹性波CT检测中收集到的数据进行聚类分析,然后再对比钻孔探测数据,是近年来对钻孔探测和弹性波CT进行结合分析的主要手段;但该方法需要拿钻孔探测和弹性波CT的数据专门对比,操作起来比较麻烦,且容易出错;本发明通过将钻孔探测数据转化为弹性波CT数据,并进行标记,将标记以后的钻孔探测数据和弹性波CT数据同时进行聚类分析;通过聚类分析结果,能够直观的给出钻孔探测和弹性波CT检测的结果的一致性。
发明内容
本发明的目的在于克服现有技术存在的以上问题,以及在原有技术基础上大大提高其技术效果;为此,本发明提供了一种串状溶洞地质钻孔探测与弹性波CT扫描的数据聚类分析方法,该方法包括:
选择串状溶洞工程地点适当的位置进行钻孔,并对钻孔数据进行处理,使得处理以后,将钻孔数据转化为不同地质对应的弹性波标准波速;
进行弹性波CT检测,获取弹性波CT扫描图像;提取弹性波CT图像信息,提取各点的CT值;所述CT值对应CT图像中各地质中弹性波波速,CT图像所对应的数据以矩阵的形式储存在计算机中;
将钻孔数据转化的弹性波标准波速混入提取的CT值中,获取混合数据;
对混合数据中钻孔数据转化的弹性波标准波速进行标记处理;并对标记处理后的混合数据进行聚类分析;通过聚类分析的结果,判定钻孔结果和弹性波CT检测的结果是否一致。
所述对钻孔数据进行处理包括:所述钻孔数据指钻孔的结果,钻孔的结果指工程地点不同位置对应的土质材料;根据土质情况,查找钻探数据表,获取不同土质材料对应的弹性波波速;所述对应的弹性波波速是一个范围值,取所述范围值最中间的值作为不同土质材料对应的弹性波标准波速。
所述进行弹性波CT检测包括:使用弹性波检测仪器对工程地点的地质情况进行检测,生成二维弹性波CT图像,提取CT图像中各点的CT值,所述CT值对应CT图像中各地质中弹性波波速。
所述将钻孔数据转化的弹性波标准波速混入提取的CT值中包括:将钻孔数据转化的弹性波标准波速混入到CT图像对应数据的矩阵中。
所述并对标记处理后的混合数据进行聚类分析包括:通过K均值聚类算法判别混合数据中每一点的数据所属的范围,所述每一点的数据所属的范围为波速范围;所属K均值聚类算法是根据离散分布的样本集D={x1,x2,…xn},x1、x2、直到xn分别代表单个样本,通过不断更新中心对象的方法,获得聚类所得簇划分C={C1,C2,…Ck},C1、C2、直到Ck分别代表单个簇,其中每个簇中包括多个处于同一簇的单个样本;聚类标准是使得平方误差最小化,所述平方误差公式为:其中,ui是簇Ci的均值向量,/>可以看出,Ed的值越小,从整体上讲,簇内样本x的相似度越高。
所述并对标记处理后的混合数据进行聚类分析还包括:通过k均值聚类算法对混合数据进行聚类分析,将混合数据分为k簇,即k类;每类中包括被标记的钻孔数据转化的弹性波标准波速和弹性波CT图像提取的各点CT值;分析每个分类中的被标记的钻孔数据转化的弹性波标准波速对应的土质材料与k均值聚类算法所得所属分类的土质材料是否相同,判定钻孔探测与弹性波CT检测结果是否一致。
本发明的有益效果是:
本发明提出了一种串状溶洞地质钻孔探测与弹性波CT扫描的数据聚类分析方法;根据钻孔探测得到的土层结构和弹性波在各土层结构的传播速度,将钻孔探测数据转化为弹性波CT数据,并进行标记;将标记以后的钻孔探测数据和弹性波CT数据同时进行聚类分析;通过聚类分析结果,直观的给出钻孔探测和弹性波CT检测的结果的一致性。
附图说明
图1:本发明的一种串状溶洞地质钻孔探测与弹性波CT扫描的数据聚类分析方法的流程图。
具体实施方式
以下结合附图对本发明的具体实施例进行详细说明;应当理解的是此处所给出的具体实施例仅用于说明和解释本发明,并不能用来限制本发明。
如图1所示,根据本发明的实施例的一种串状溶洞地质钻孔探测与弹性波CT扫描的数据聚类分析方法,该方法包括:步骤S200,选择串状溶洞工程地点适当的位置进行钻孔,并对钻孔数据进行处理,使得处理以后,将钻孔数据转化为不同地质对应的弹性波标准波速;步骤S201,进行弹性波CT检测,获取弹性波CT扫描图像;提取弹性波CT图像信息,提取各点的CT值;所述CT值对应CT图像中各地质中弹性波波速,CT图像所对应的数据以矩阵的形式储存在计算机中;步骤S202,将钻孔数据转化的弹性波标准波速混入提取的CT值中,获取混合数据;步骤S203,对混合数据中钻孔数据转化的弹性波标准波速进行标记处理;并对标记处理后的混合数据进行聚类分析;通过聚类分析的结果,判定钻孔结果和弹性波CT检测的结果是否一致。
具体地,本发明首先选择串状溶洞工程地点适当的位置进行钻孔,并将钻孔数据转化为不同地质对应的弹性波标准波速;随后,进行弹性波CT检测,获取弹性波CT扫描图像;提取图像信息中各点的CT值;将钻孔数据转化的弹性波标准波速混入提取的CT值中,获取混合数据;对混合数据中钻孔数据转化的弹性波标准波速进行标记处理,并对标记处理后的混合数据进行聚类分析;通过聚类分析的结果,判定钻孔结果和弹性波CT检测的结果的一致性。
步骤S200,选择串状溶洞工程地点适当的位置进行钻孔,并对钻孔数据进行处理,使得处理以后,将钻孔数据转化为不同地质对应的弹性波标准波速;具体地,根据钻孔探测的原则,在串状溶洞工程地点选择适当的位置进行钻孔探测;根据钻孔探测结果,将钻孔数据转化为不同地质对应的弹性波标准波速。
在上述实施例中,具体地,所述钻孔数据指钻孔的结果,钻孔的结果指工程地点不同位置对应的土质材料;根据土质情况,查找钻探数据表,获取不同土质材料对应的弹性波波速;所述对应的弹性波波速是一个范围值,取所述范围值最中间的值作为钻孔探测结果中不同土质材料对应的弹性波标准波速。
步骤S201,进行弹性波CT检测,获取弹性波CT扫描图像;提取弹性波CT图像信息,提取各点的CT值;所述CT值对应CT图像中各地质中弹性波波速,CT图像所对应的数据以矩阵的形式储存在计算机中;具体地,采用地震勘探仪器、电阻率仪器或雷达探测仪器等对工程地点适当位置进行弹性波探测;通过接收的结果,生成弹性波CT图像,提取CT图像中各点对应的CT值。
在上述实施例中,具体地,根据弹性波在不同地质结构的传播强度不同,进行弹性波CT图像中的CT值的提取;所述弹性波在不同地质结构的传播强度不同包括:波速、振幅和波长等;本发明通过提取CT图像中各点的CT值,提取的CT值对应CT图像中各地质中弹性波波速。
步骤S202,将钻孔数据转化的弹性波标准波速混入提取的CT值中,获取混合数据;具体地,根据弹性波在不同土质结构的传播速度不同的这一特性,将钻孔数据转化为对应频率的弹性波标准波速,并将转化以后的弹性波标准波速混入通过CT图像得到的CT值中。
在上述实施例中,具体地,将钻孔数据转化的弹性波标准波速混入提取的CT值中的方法是:将钻孔数据转化的弹性波标准波速混入到CT图像对应数据的矩阵中。
步骤S203,对混合数据中钻孔数据转化的弹性波标准波速进行标记处理;并对标记处理后的混合数据进行聚类分析;通过聚类分析的结果,判定钻孔结果和弹性波CT检测的结果是否一致;具体地,将钻孔数据转化的弹性波标准波速混入提取的CT值后,对混合数据中的钻孔数据转化的弹性波标准波速进行荧光标记;对标记处理以后的混合数据进行聚类分析,根据聚类分析结果,判断钻孔探测结果和弹性波CT检测结果是否一致。
在上述实施例中,具体地,聚类分析的方法为:通过K均值聚类算法判别混合数据中每一点的数据所属的范围,所述每一点的数据所属的范围为波速范围;所述K均值聚类算法是根据离散分布的样本集D={x1,x2,…xn},x1、x2、直到xn分别代表单个样本,通过不断更新中心对象的方法,获得聚类所得簇划分C={C1,C2,…Ck},C1到Ck分别代表单个簇,其中每个簇中包括多个处于同一簇的单个样本;聚类标准是使得平方误差最小化,所述平方误差公式为:其中,ui是簇Ci的均值向量,/>可以看出,Ed的值越小,从整体上讲,簇内样本x的相似度越高。
在上述实施例中,具体地,通过k均值聚类算法对混合数据进行聚类分析,将混合数据分为k簇,即k类;每类中包括被标记的钻孔数据转化的弹性波标准波速和弹性波CT图像提取的各点CT值;分析每个分类中的被标记的钻孔数据转化的弹性波标准波速对应的土质材料与k均值聚类算法所得所属分类的土质材料是否相同,若相同,则证明钻孔探测结果与弹性波CT检测结果一致。
需要理解的是,上述实施例为本发明的一个或多个实施例,基于本发明还有很多其他实施例及其变形;本行业的普通技术人员在没有作出开拓性的创新的时候,通过本发明进行的变形和修改,均属于本发明的保护范围。

Claims (6)

1.一种串状溶洞地质钻孔探测与弹性波CT检测的数据聚类分析方法,其特征在于,所述方法包括如下步骤:
选择串状溶洞工程地点适当的位置进行钻孔,并对钻孔数据进行处理,使得处理以后,将钻孔数据转化为不同地质对应的弹性波标准波速;
进行弹性波CT检测,获取弹性波CT扫描图像;提取弹性波CT图像信息,提取各点的CT值;所述CT值对应CT图像中各地质中弹性波波速,CT图像所对应的数据以矩阵的形式储存在计算机中;
将钻孔数据转化的弹性波标准波速混入提取的CT值中,获取混合数据;
对混合数据中钻孔数据转化的弹性波标准波速进行标记处理;并对标记处理后的混合数据进行聚类分析;通过聚类分析的结果,判定钻孔结果和弹性波CT检测的结果是否一致。
2.根据权利要求1所述的一种串状溶洞地质钻孔探测与弹性波CT扫描的数据聚类分析方法,其特征在于,所述对钻孔数据进行处理包括:所述钻孔数据指钻孔的结果,钻孔的结果指工程地点不同位置对应的土质材料;根据土质情况,查找钻探数据表,获取不同土质材料对应的弹性波波速;所述对应的弹性波波速是一个范围值,取所述范围值最中间的值作为不同土质材料对应的弹性波标准波速。
3.根据权利要求1所述的一种串状溶洞地质钻孔探测与弹性波CT扫描的数据聚类分析方法,其特征在于,所述进行弹性波CT检测包括:使用弹性波检测仪器对工程地点的地质情况进行检测,生成二维弹性波CT图像,提取CT图像中各点的CT值,所述CT值对应CT图像中各地质中弹性波波速。
4.根据权利要求1所述的一种串状溶洞地质钻孔探测与弹性波CT扫描的数据聚类分析方法,其特征在于,所述将钻孔数据转化的弹性波标准波速混入提取的CT值中包括:将钻孔数据转化的弹性波标准波速混入到CT图像对应数据的矩阵中。
5.根据权利要求1所述的一种串状溶洞地质钻孔探测与弹性波CT扫描的数据聚类分析方法,其特征在于,所述并对标记处理后的混合数据进行聚类分析包括:通过K均值聚类算法判别混合数据中每一点的数据所属的范围,所述每一点的数据所属的范围为波速范围;所属K均值聚类算法是根据离散分布的样本集D={x1,x2,…xn},x1、x2、直到xn分别代表单个样本,通过不断更新中心对象的方法,获得聚类所得簇划分C={C1,C2,…Ck},C1、C2、直到Ck分别代表单个簇,其中每个簇中包括多个处于同一簇的单个样本;聚类标准是使得平方误差最小化,所述平方误差公式为:其中,ui是簇Ci的均值向量,可以看出,Ed的值越小,从整体上讲,簇内样本x的相似度越高。
6.根据权利要求5所述的一种串状溶洞地质钻孔探测与弹性波CT扫描的数据聚类分析方法,其特征在于,所述并对标记处理后的混合数据进行聚类分析还包括:通过k均值聚类算法对混合数据进行聚类分析,将混合数据分为k簇,即k类;每类中包括被标记的钻孔数据转化的弹性波标准波速和弹性波CT图像提取的各点CT值;分析每个分类中的被标记的钻孔数据转化的弹性波标准波速对应的土质材料与k均值聚类算法所得所属分类的土质材料是否相同,判定钻孔探测与弹性波CT检测结果是否一致。
CN202310407986.3A 2023-04-17 2023-04-17 一种串状溶洞地质钻孔探测与弹性波ct扫描的数据聚类分析方法 Pending CN116578895A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310407986.3A CN116578895A (zh) 2023-04-17 2023-04-17 一种串状溶洞地质钻孔探测与弹性波ct扫描的数据聚类分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310407986.3A CN116578895A (zh) 2023-04-17 2023-04-17 一种串状溶洞地质钻孔探测与弹性波ct扫描的数据聚类分析方法

Publications (1)

Publication Number Publication Date
CN116578895A true CN116578895A (zh) 2023-08-11

Family

ID=87543648

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310407986.3A Pending CN116578895A (zh) 2023-04-17 2023-04-17 一种串状溶洞地质钻孔探测与弹性波ct扫描的数据聚类分析方法

Country Status (1)

Country Link
CN (1) CN116578895A (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101893800B1 (ko) * 2017-08-09 2018-09-04 제주대학교 산학협력단 암석 물리학상 구축을 통한 퇴적환경 해석 방법
CN108646298A (zh) * 2018-04-09 2018-10-12 陕西铁道工程勘察有限公司 多井多参数约束下隧道波速反演成像及围岩分级的方法
CN112257789A (zh) * 2020-10-26 2021-01-22 西南交通大学 对围岩等级进行识别的方法
CN113960695A (zh) * 2021-09-18 2022-01-21 山东大学 一种复杂城市环境富水岩溶精细探查方法
CN115494560A (zh) * 2022-10-10 2022-12-20 广东省工程勘察院 一种基于ct和钻探数据相结合的建模新方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101893800B1 (ko) * 2017-08-09 2018-09-04 제주대학교 산학협력단 암석 물리학상 구축을 통한 퇴적환경 해석 방법
CN108646298A (zh) * 2018-04-09 2018-10-12 陕西铁道工程勘察有限公司 多井多参数约束下隧道波速反演成像及围岩分级的方法
CN112257789A (zh) * 2020-10-26 2021-01-22 西南交通大学 对围岩等级进行识别的方法
CN113960695A (zh) * 2021-09-18 2022-01-21 山东大学 一种复杂城市环境富水岩溶精细探查方法
CN115494560A (zh) * 2022-10-10 2022-12-20 广东省工程勘察院 一种基于ct和钻探数据相结合的建模新方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
袁利永 等: "一种改进的半监督K-Means聚类算法", 《计算机工程与科学》, vol. 33, no. 6, 30 June 2011 (2011-06-30), pages 1 - 6 *

Similar Documents

Publication Publication Date Title
CN103026202A (zh) 获取多孔介质的相容和综合物理性质的方法
Wodecki et al. Process Monitoring in Heavy Duty Drilling Rigs—Data Acquisition System and Cycle Identification Algorithms
Koike et al. Spatial correlation structures of fracture systems for deriving a scaling law and modeling fracture distributions
Moustafa et al. A quantitative site-specific classification approach based on affinity propagation clustering
Mohammadpour et al. Automatic lineament extraction method in mineral exploration using CANNY algorithm and hough transform
Xu et al. Integration of NMR and conventional logs for vuggy facies classification in the Arbuckle Formation: A machine learning approach
CN109298464A (zh) 致密砂岩储层成岩相测井识别方法及装置
TWI626622B (zh) System and method for stereoscopic imaging of underground rock formation characteristics
Dong et al. Fracture identification in tight reservoirs by multiple kernel Fisher discriminant analysis using conventional logs
Chandan et al. An intelligent model for indian soil classification using various machine learning techniques
Timmer et al. Applying deep learning for identifying bioturbation from core photographs
Syambas An approach for predicting the shape and size of a buried basic object on surface ground penetrating radar system
Hu et al. A sample enhancement method based on simple linear iterative clustering superpixel segmentation applied to multibeam seabed classification
CN116578895A (zh) 一种串状溶洞地质钻孔探测与弹性波ct扫描的数据聚类分析方法
CN108732625B (zh) 一种井周岩石非均匀各向异性的识别方法及系统
CN110806444B (zh) 基于浅地层剖面仪与svm的海底底质识别与分类方法
Song et al. Pre-stack-texture-based reservoir characteristics and seismic facies analysis
Wang et al. Characteristic parameters extraction method of hidden Karst Cave from Borehole radar signal
Tang et al. Automatic extraction of rock discontinuities from the point cloud using dynamic DBSCAN algorithm
Min et al. Strategy of oversampling geotechnical parameters through geostatistical, SMOTE, and CTGAN methods for assessing susceptibility of landslide
Zaman et al. Developing a low-cost device for preliminary geotechnical site characterization using multi-channel analysis of surface wave
Green et al. Rise of the Machines: Improving the identification of possible graves in GPR data with interactive survey guidance and machine learning
US11668847B2 (en) Generating synthetic geological formation images based on rock fragment images
Dai et al. Study of an Automatic Picking Method for Multimode Dispersion Curves of Surface Waves Based on an Improved U-Net
CN114383959B (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