CN118094359B - 基于磨料水射流截割的煤岩裂纹风险预测方法 - Google Patents
基于磨料水射流截割的煤岩裂纹风险预测方法 Download PDFInfo
- Publication number
- CN118094359B CN118094359B CN202410508380.3A CN202410508380A CN118094359B CN 118094359 B CN118094359 B CN 118094359B CN 202410508380 A CN202410508380 A CN 202410508380A CN 118094359 B CN118094359 B CN 118094359B
- Authority
- CN
- China
- Prior art keywords
- model
- coal
- cutting
- rock
- coal rock
- 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
- 239000011435 rock Substances 0.000 title claims abstract description 140
- 239000003245 coal Substances 0.000 title claims abstract description 133
- 238000005520 cutting process Methods 0.000 title claims abstract description 133
- 238000000034 method Methods 0.000 title claims abstract description 51
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 36
- 238000004458 analytical method Methods 0.000 claims abstract description 21
- 230000008859 change Effects 0.000 claims abstract description 18
- 230000006870 function Effects 0.000 claims description 44
- 238000012706 support-vector machine Methods 0.000 claims description 27
- 238000012360 testing method Methods 0.000 claims description 15
- 238000012549 training Methods 0.000 claims description 15
- 238000009826 distribution Methods 0.000 claims description 13
- 238000005457 optimization Methods 0.000 claims description 13
- 239000011148 porous material Substances 0.000 claims description 13
- 230000015572 biosynthetic process Effects 0.000 claims description 12
- 238000005516 engineering process Methods 0.000 claims description 10
- 238000012544 monitoring process Methods 0.000 claims description 9
- 239000013598 vector Substances 0.000 claims description 9
- 238000002591 computed tomography Methods 0.000 claims description 8
- 238000002790 cross-validation Methods 0.000 claims description 7
- 238000005259 measurement Methods 0.000 claims description 7
- 238000007476 Maximum Likelihood Methods 0.000 claims description 6
- 238000013507 mapping Methods 0.000 claims description 6
- 208000010392 Bone Fractures Diseases 0.000 claims description 5
- 206010017076 Fracture Diseases 0.000 claims description 5
- 230000000694 effects Effects 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000000704 physical effect Effects 0.000 claims description 4
- 230000004044 response Effects 0.000 claims description 4
- 238000007619 statistical method Methods 0.000 claims description 4
- 238000012795 verification Methods 0.000 claims description 4
- 238000012669 compression test Methods 0.000 claims description 3
- 230000009022 nonlinear effect Effects 0.000 claims description 3
- 238000010998 test method Methods 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 230000009466 transformation Effects 0.000 claims description 2
- 239000000463 material Substances 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 229910000831 Steel Inorganic materials 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000010219 correlation analysis Methods 0.000 description 2
- 238000005336 cracking Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 239000010959 steel Substances 0.000 description 2
- 238000010200 validation analysis Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000013499 data model Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 238000011478 gradient descent method Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000013179 statistical model Methods 0.000 description 1
- 230000003313 weakening effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2411—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on the proximity to a decision surface, e.g. support vector machines
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B24—GRINDING; POLISHING
- B24C—ABRASIVE OR RELATED BLASTING WITH PARTICULATE MATERIAL
- B24C1/00—Methods for use of abrasive blasting for producing particular effects; Use of auxiliary equipment in connection with such methods
- B24C1/04—Methods for use of abrasive blasting for producing particular effects; Use of auxiliary equipment in connection with such methods for treating only selected parts of a surface, e.g. for carving stone or glass
- B24C1/045—Methods for use of abrasive blasting for producing particular effects; Use of auxiliary equipment in connection with such methods for treating only selected parts of a surface, e.g. for carving stone or glass for cutting
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
- G01N23/046—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N9/00—Investigating density or specific gravity of materials; Analysing materials by determining density or specific gravity
- G01N9/24—Investigating density or specific gravity of materials; Analysing materials by determining density or specific gravity by observing the transmission of wave or particle radiation through the material
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/245—Classification techniques relating to the decision surface
- G06F18/2451—Classification techniques relating to the decision surface linear, e.g. hyperplane
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/27—Regression, e.g. linear or logistic regression
-
- 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
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0635—Risk analysis of enterprise or organisation activities
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/12—Classification; Matching
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Biology (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Analytical Chemistry (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Computation (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Operations Research (AREA)
- Mathematical Physics (AREA)
- Strategic Management (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Entrepreneurship & Innovation (AREA)
- Economics (AREA)
- Pure & Applied Mathematics (AREA)
- Mechanical Engineering (AREA)
- Dispersion Chemistry (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Pulmonology (AREA)
- Radiology & Medical Imaging (AREA)
- Development Economics (AREA)
- Probability & Statistics with Applications (AREA)
- Educational Administration (AREA)
Abstract
本发明涉及裂纹预测分析技术领域,具体涉及基于磨料水射流截割的煤岩裂纹风险预测方法,包括以下步骤:测量煤岩的初始物理特性,初始物理特性包括煤岩基本结构和强度特性,同时建立煤岩基本结构‑强度模型;获得切割路径横截面的煤岩基本结构,基于煤岩基本结构‑强度模型,分析切割路径上的抗压强度变化;在磨料水射流切割过程中,实时监测沿切割路径的声发射特征,以捕捉切割中局部切割位置的煤岩微裂纹发生情况,基于微裂纹发生情况验证基本结构‑强度模型;实时调整切割路径和/或磨料水射流的工作参数。本发明,不仅增强了煤岩切割过程中的风险管理能力,还提高了切割操作的安全性和可靠性。
Description
技术领域
本发明涉及裂纹预测分析技术领域,尤其涉及基于磨料水射流截割的煤岩裂纹风险预测方法。
背景技术
磨料水射流截割技术是一种广泛应用于硬质材料加工的高效方法,尤其在煤岩开采和加工领域中具有重要作用。这种技术利用高速水流和磨料颗粒的联合作用来切割煤岩,既可以实现精确切割,又能最小化热影响和机械应力,保持煤岩的结构完整性。然而,切割过程中煤岩的裂纹形成问题一直是技术和安全上的挑战。裂纹的产生不仅会降低煤岩的商业价值,也可能导致切割效率降低和潜在的安全风险。
传统上,对煤岩裂纹风险的预测主要依赖于经验判断和简单的理论模型,缺乏对煤岩内部复杂结构变化的深入理解。这种方法往往无法准确预测切割过程中的裂纹发展,限制了磨料水射流截割技术的优化和效率提升。此外,缺乏有效的实时监测手段来跟踪切割过程中的微裂纹形成,使得及时调整切割策略以避免裂纹扩展变得困难。
因此,急需一种能够精确预测煤岩在磨料水射流切割过程中裂纹风险的方法,并实时监测裂纹的发展,以指导切割参数的调整和优化。
发明内容
基于上述目的,本发明提供了基于磨料水射流截割的煤岩裂纹风险预测方法。
基于磨料水射流截割的煤岩裂纹风险预测方法,包括以下步骤:
S1:在磨料水射流切割开始前,测量煤岩的初始物理特性,初始物理特性包括煤岩基本结构和强度特性,其中,煤岩基本结构为密度、孔隙率,强度特性表现为抗压强度,同时建立煤岩基本结构-强度模型;
S2:基于S1的煤岩的初始物理特性测量,获得切割路径横截面的煤岩基本结构,基于煤岩基本结构-强度模型,分析切割路径上的抗压强度变化,从而预测裂纹风险;
S3:在磨料水射流切割过程中,实时监测沿切割路径的声发射特征,以捕捉切割中局部切割位置的煤岩微裂纹发生情况,基于微裂纹发生情况验证基本结构-强度模型;
S4:根据验证后基本结构-强度模型,以及S2中预测的裂纹风险,实时调整切割路径和/或磨料水射流的工作参数,包括改变切割速度、调整压力或流量,以优化切割过程并最小化裂纹风险。
进一步的,所述S1具体包括:
S11:采用X射线计算机断层扫描技术测量煤岩的整体密度和孔隙率,通过分析扫描图像,计算煤岩内部的孔隙结构和总体密度分布,评估煤岩的孔隙特性和整体密实度;
S12:利用标准化的岩石力学试验方法测量煤岩的抗压强度,包括单轴或多轴压缩测试,以获得煤岩在不同应力条件下的强度参数;
S13:结合S11和S13的测试结果,综合分析煤岩基本结构和强度特性,建立煤岩基本结构-强度模型,用于描述煤岩的密度、孔隙率与其抗压强度之间的关系。
进一步的,所述煤岩基本结构-强度模型基于非线性回归模型,具体包括:
收集大量煤岩样本数据,包括样本的密度、孔隙率和对应的抗压强度测量值;
对收集的数据进行统计分析,识别密度、孔隙率与抗压强度之间的基本关系和影响模式;
选择多项式回归模型来描述非线性关系,模型形式为:
,其中,是抗压强度,是煤岩的密度,是孔隙率,是模型参数,是描述非线性关系的指数,是误差项,、表示密度和孔隙率的非线性影响;
参数估计:使用最大似然估计方法来估计模型参数a,b,c,n,m;
模型验证:通过留一交叉验证技术验证模型的能力。
进一步的,所述最大似然估计方法具体包括:
基于多项式回归模型中的误差项遵循正态分布,构建似然函数,似然函数是给定参数时,观测到的数据出现概率的函数,似然函数表示为所有独立观测数据抗压强度的联合概率密度函数:
,其中分别是第个观测样本的抗压强度、密度和孔隙率,表示误差项的方差,反映模型中随机误差的大小;
最大化似然函数:参数的估计值是使似然函数最大化的值,通过对似然函数取对数(即对数似然函数)并应用数值优化方法来实现,因为对数转换后的乘积形式变为求和形式,表示为:
;
计算和优化:利用梯度下降法数值优化技术来求解对数似然函数的最大值,计算对数似然函数关于每个参数的导数,并迭代调整参数值以找到使对数似然函数最大化的参数组合。
进一步的,所述留一交叉验证技术表示为LOOCV,用于验证多项式回归模型能力,验证步骤如下:
准备数据集,数据集中包括若干数据点,每个数据点包括煤岩的密度、孔隙率和抗压强度,数据集中的每个样本都将轮流用作测试数据;
实施LOOCV:对于数据集中的每一个样本:
A.留出一个样本:从数据集中选取一个样本作为测试数据,其余的样本作为训练数据;
B.建立模型:使用训练数据建立多项式回归模型;
C.模型预测:使用留出的测试样本来评估模型的预测能力,计算预测值与实际值之间的差异;
计算平均误差:重复实施LOOCV,直到每个样本都被用作测试数据一次,记录每一次的预测误差,然后计算误差的平均值,平均值表示模型的总体预测能力;
分析结果:分析LOOCV的结果来评估多项式回归模型的稳定性和预测能力,低平均误差值表示模型具有良好预测性能。
进一步的,所述S2具体包括:
重建横截面图像:基于X射线计算机断层扫描技术扫描结果,通过图像重建技术生成煤岩的横截面图像,以显示预定切割路径区域的内部结构;
通过横截面图像,计算预定切割路径区域内煤岩密度,图像中,不同密度的区域会显示不同的灰度级别,可以利用这一特点对煤岩密度进行定量化分析,根据横截面图像中孔隙和裂缝的分布和大小,评估切割路径区域的孔隙率;
映射切割路径,将密度和孔隙率数据映射到预定的切割路径上,将具体的图像数据点与切割路径的几何位置对应起来,从而获得沿着切割路径的密度和孔隙率分布。
进一步的,所述S2中的预测裂纹风险包括对于切割路径上的每一点或某区域,将相应的密度和孔隙率数据代入基本结构-强度模型中计算预测的抗压强度,获得切割过程中各个点或区域的抗压强度估计;
分析沿切割路径的抗压强度预测值,关注低于临界抗压强度的区域,为高风险裂纹区域。
进一步的,所述S3具体包括:
安装声发射传感器在切割路径附近,用于实时捕捉煤岩在切割点的声发射信号,指示煤岩结构变化和裂纹形成活动;
收集传感器数据,并利用快速傅里叶变换FFT对振动数据进行频率分析,以及对声发射信号进行时频分析,确定切割中煤岩的动态响应特性;
结合实时监测数据和预定切割路径信息,使用支持向量机SVM模型识别煤岩微裂纹形成的早期迹象。
进一步的,所述快速傅里叶变换FFT将时域信号转换为频域信号,表示为:,其中,是频域中的第个频率分量,是时域中的第个样本,是样本总数;
所述时频分析用于理解声发射信号随时间的变化,通过短时傅里叶变换进行时频分析,以确定在切割过程中何时何地发生变化。
进一步的,所述支持向量机SVM模型具体包括:
基于频率分析以及时频分析结果中提取特征,包括峰值幅度、能量、频率内容和持续时间,形成特征向量用于SVM训练和预测;
SVM模型建立:通过一个超平面分隔不同类别的数据点,最大化不同类别之间的边界,对于微裂纹识别问题,将其视为二分类问题,即微裂纹存在,即正类,或不存在,即负类;
通过最小化目标函数来训练SVM模型:
,其中,是超平面的法向量,是正则化参数(控制误分类的惩罚),是松弛变量(用于处理不完全可分的情况),是训练样本数量,目标是找到超平面的法向量和偏差项使得目标函数最小化;
决策函数为:,当的值大于零时,被分类为正类,即微裂纹存在,否则为负类,即微裂纹不存在;
使用训练好的SVM模型对新收集的声发射数据进行分类,以判断煤岩是否在某个切割位置形成微裂纹。
本发明的有益效果:
本发明,通过结合煤岩基本结构-强度模型与实时监测技术,提供了高度精准的裂纹风险预测方法,煤岩的密度、孔隙率以及抗压强度等物理特性是裂纹形成的关键因素,利用X射线计算机断层扫描(CT)进行初始测量,可以获得煤岩的详细结构特性,为建立基本结构-强度模型提供准确数据,实时监测技术,特别是声发射监测和支持向量机(SVM)分析,能够在切割过程中捕捉微裂纹的形成和发展。这种实时性允许操作者在裂纹发展到危险程度之前及时调整切割参数,如改变切割速度或压力,从而有效控制裂纹风险。此外,实时监测的数据可以用来验证和优化建立的基本结构-强度模型,通过将模型预测的裂纹风险点与实时监测到的微裂纹发生情况进行对比,可以评估模型的准确性。这种精准的裂纹风险预测和实时监测相结合的方法,不仅增强了煤岩切割过程中的风险管理能力,还提高了切割操作的安全性和可靠性。通过有效预测和控制裂纹的形成,可以避免切割过程中的潜在危险,确保煤岩结构的完整性和切割作业的顺利进行。
本发明,利用实时数据分析和模型预测来指导磨料水射流切割参数的调整或者切割路径的优化,动态调整策略可以根据煤岩的实际响应优化切割过程,提高切割效率,并确保切割质量,通过减少不必要的停机时间和避免过度切割,显著提升整体的切割性能和经济效益。
附图说明
为了更清楚地说明本发明或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例的风险预测逻辑框图;
图2为本发明实施例的方法流程示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,对本发明进一步详细说明。
如图1-2所示,基于磨料水射流截割的煤岩裂纹风险预测方法,包括以下步骤:
S1:在磨料水射流切割开始前,测量煤岩的初始物理特性,初始物理特性包括煤岩基本结构和强度特性,其中,煤岩基本结构为密度、孔隙率,强度特性表现为抗压强度,同时建立煤岩基本结构-强度模型;
S2:基于S1的煤岩的初始物理特性测量,获得切割路径横截面的煤岩基本结构,基于煤岩基本结构-强度模型,分析切割路径上的抗压强度变化,从而预测裂纹风险;
S3:在磨料水射流切割过程中,实时监测沿切割路径的声发射特征,以捕捉切割中局部切割位置的煤岩微裂纹发生情况,基于微裂纹发生情况验证基本结构-强度模型;
S4:根据验证后基本结构-强度模型,以及S2中预测的裂纹风险,实时调整切割路径和/或磨料水射流的工作参数,包括改变切割速度、调整压力或流量,以优化切割过程并最小化裂纹风险。
S1具体包括:
S11:采用X射线计算机断层扫描技术测量煤岩的整体密度和孔隙率,通过分析扫描图像,计算煤岩内部的孔隙结构和总体密度分布,评估煤岩的孔隙特性和整体密实度;
S12:利用标准化的岩石力学试验方法测量煤岩的抗压强度,包括单轴或多轴压缩测试,以获得煤岩在不同应力条件下的强度参数;
S13:结合S11和S13的测试结果,综合分析煤岩基本结构和强度特性,建立煤岩基本结构-强度模型,用于描述煤岩的密度、孔隙率与其抗压强度之间的关系。
在水射流切割过程中,煤岩的密度和孔隙率通常不会发生本质性的变化,但切割作用会在局部区域内影响煤岩的结构特性,具体表现如下:
煤岩的总体密度在水射流切割过程中基本保持不变,然而,在切割接触区域,由于物料被去除,会造成局部密度的相对变化。
孔隙率会在切割区域发生轻微变化,水射流切割通过高速水流和磨料的冲击作用,可能会在煤岩表面或切割区域内部产生新的微小裂纹或扩大现有的孔隙,从而略微增加局部的孔隙率。这种变化通常局限于切割区域附近。
因此,虽然整体的密度和孔隙率变化不大,但在切割过程中,尤其是在水射流直接作用的区域,煤岩的微观结构会经历一些改变,得到煤岩基本结构-强度关系模型后,在后续预测裂纹时,只需在切割时实时监测分析密度、孔隙率的实时变化,即可预测抗压强度变化,抗压强度变化可用于预测是否会产生裂纹。
煤岩的密度、孔隙率与其抗压强度之间的关系通常是非线性的,这种非线性关系源于煤岩复杂的内部结构和不同因素之间相互作用的复杂性。通常认为,密度较高的煤岩因为结构更为紧密,具有更高的抗压强度。然而,这种关系并不总是线性的,因为煤岩的密度增加到一定程度后,其抗压强度的增加率可能会下降。孔隙率高的煤岩通常抗压强度较低,因为孔隙作为裂纹扩展的通道,减弱了煤岩的整体强度。这种关系往往是非线性的,特别是当孔隙率在某一范围内变化时,其对抗压强度的影响会更加显著。
因此,为了准确描述煤岩的物理特性与其抗压强度之间的关系,需要利用非线性模型来反映这种复杂的相互依赖性。这种非线性模型可以更准确地预测在不同条件下煤岩的行为,特别是在进行裂纹风险预测和磨料水射流切割参数优化时,因此,煤岩基本结构-强度模型基于非线性回归模型,具体包括:
收集大量煤岩样本数据,包括样本的密度、孔隙率和对应的抗压强度测量值;
对收集的数据进行统计分析,识别密度、孔隙率与抗压强度之间的基本关系和影响模式,将收集的煤岩样本数据整理成可分析的表格形式,每一行代表一个样本,每一列代表一个变量(如密度、孔隙率和抗压强度);
进行描述性统计分析,包括计算密度、孔隙率和抗压强度的均值、中位数、标准差、最大值和最小值,以获取数据的基本分布情况和特点。
进行相关性分析,计算密度、孔隙率与抗压强度之间的相关系数,通常使用皮尔逊相关系数进行线性相关性分析。相关系数的计算公式为:,其中,和分别是单个样本的密度(或孔隙率)和抗压强度的值,和分别是密度(或孔隙率)和抗压强度的平均值;
选择多项式回归模型来描述非线性关系,模型形式为:
,其中,是抗压强度,是煤岩的密度,是孔隙率,是模型参数,是描述非线性关系的指数,是误差项,、表示密度和孔隙率的非线性影响;
参数估计:使用最大似然估计方法来估计模型参数a,b,c,n,m;
模型验证:通过留一交叉验证技术验证模型的能力。
最大似然估计方法具体包括:
基于多项式回归模型中的误差项遵循正态分布,构建似然函数,似然函数是给定参数时,观测到的数据出现概率的函数,似然函数表示为所有独立观测数据抗压强度的联合概率密度函数:
,其中分别是第个观测样本的抗压强度、密度和孔隙率,表示误差项的方差,反映模型中随机误差的大小;
最大化似然函数:参数的估计值是使似然函数最大化的值,通过对似然函数取对数(即对数似然函数)并应用数值优化方法来实现,因为对数转换后的乘积形式变为求和形式,表示为:
;
计算和优化:利用梯度下降法数值优化技术来求解对数似然函数的最大值,计算对数似然函数关于每个参数的导数,并迭代调整参数值以找到使对数似然函数最大化的参数组合。
为了应用梯度下降法,需要计算对数似然函数相对于每个参数的梯度,对每个参数求偏导数:
,同样,对进行偏导以得到它们各自的梯度表达式。
选择初始参数:选择合适的初始值开始迭代,初值可基于经验、先前研究来设定。
迭代更新参数:在每一步迭代中,根据梯度方向和选择的学习率(步长)更新参数值。对于参数,更新公式为:
,其中是学习率,控制每次更新的步长;类似地,对其他参数进行更新。
迭代直到收敛:重复执行上述步骤,直到参数值收敛或达到预设的迭代次数。收敛的标准是梯度的变化小于某个阈值或对数似然函数的改变小于某个预定值。
留一交叉验证技术表示为LOOCV,用于验证多项式回归模型能力,验证步骤如下:
准备数据集,数据集中包括若干数据点,每个数据点包括煤岩的密度、孔隙率和抗压强度,数据集中的每个样本都将轮流用作测试数据;
实施LOOCV:对于数据集中的每一个样本:
A.留出一个样本:从数据集中选取一个样本作为测试数据,其余的样本作为训练数据;
B.建立模型:使用训练数据建立多项式回归模型;
C.模型预测:使用留出的测试样本来评估模型的预测能力,计算预测值与实际值之间的差异;
计算平均误差:重复实施LOOCV,直到每个样本都被用作测试数据一次,记录每一次的预测误差,然后计算误差的平均值,平均值表示模型的总体预测能力;
分析结果:分析LOOCV的结果来评估多项式回归模型的稳定性和预测能力,低平均误差值表示模型具有良好预测性能;
对于每一次留一交叉验证,计算预测值与实际值之间的差异用均方误差MSE来计算:
,其中,是第个留出样本的实际抗压强度,是模型预测的抗压强度,对所有求平均值得到总体预测误差:,其中,是数据集中样本的总数。
S2具体包括:
重建横截面图像:基于X射线计算机断层扫描技术扫描结果,通过图像重建技术生成煤岩的横截面图像,以显示预定切割路径区域的内部结构;
通过横截面图像,计算预定切割路径区域内煤岩密度,图像中,不同密度的区域会显示不同的灰度级别,可以利用这一特点对煤岩密度进行定量化分析,根据横截面图像中孔隙和裂缝的分布和大小,评估切割路径区域的孔隙率,使用图像分析软件来识别和计算孔隙空间,从而得出孔隙率;
映射切割路径,将密度和孔隙率数据映射到预定的切割路径上,将具体的图像数据点与切割路径的几何位置对应起来,从而获得沿着切割路径的密度和孔隙率分布。
S2中的预测裂纹风险包括对于切割路径上的每一点或某区域,将相应的密度和孔隙率数据代入基本结构-强度模型中计算预测的抗压强度,获得切割过程中各个点或区域的抗压强度估计;
分析沿切割路径的抗压强度预测值,关注低于临界抗压强度的区域,为高风险裂纹区域。
S3具体包括:
安装声发射传感器在切割路径附近,用于实时捕捉煤岩在切割点的声发射信号,指示煤岩结构变化和裂纹形成活动;
收集传感器数据,并利用快速傅里叶变换FFT对振动数据进行频率分析,以及对声发射信号进行时频分析,确定切割中煤岩的动态响应特性;
结合实时监测数据和预定切割路径信息,使用支持向量机SVM模型识别煤岩微裂纹形成的早期迹象。
快速傅里叶变换FFT将时域信号转换为频域信号,表示为:
,其中,是频域中的第个频率分量,是时域中的第个样本,是样本总数,这个转换揭示了信号中的主要频率成分,帮助识别由裂纹形成引起的声发射事件;
时频分析用于理解声发射信号随时间的变化,通过短时傅里叶变换进行时频分析,以确定在切割过程中何时何地发生变化。
支持向量机SVM模型具体包括:
基于频率分析以及时频分析结果中提取特征,包括峰值幅度、能量、频率内容和持续时间,形成特征向量用于SVM训练和预测;
SVM模型建立:通过一个超平面分隔不同类别的数据点,最大化不同类别之间的边界,对于微裂纹识别问题,将其视为二分类问题,即微裂纹存在,即正类,或不存在,即负类;
通过最小化目标函数来训练SVM模型:
,其中,是超平面的法向量,是正则化参数(控制误分类的惩罚),是松弛变量(用于处理不完全可分的情况),是训练样本数量,目标是找到超平面的法向量和偏差项使得目标函数最小化;
决策函数为:,当的值大于零时,被分类为正类,即微裂纹存在,否则为负类,即微裂纹不存在;
使用训练好的SVM模型对新收集的声发射数据进行分类,以判断煤岩是否在某个切割位置形成微裂纹。
基于获得的切割过程中各个点或区域的抗压强度估计,生成一个裂纹风险分布图,标出各个预测的高风险区域。对比基本结构-强度模型预测的裂纹风险区域与SVM模型实时监测分析到的实际微裂纹发生区域。关注预测为高风险并且也确实捕捉到微裂纹活动的区域。一致性高的区域表明模型预测的准确性。统计模型成功预测裂纹发生的次数相对于总的裂纹发生次数。计算模型的准确率、召回率指标,来量化模型的性能。
所属领域的普通技术人员应当理解:以上任何实施例的讨论仅为示例性的,并非旨在暗示本发明的范围被限于这些例子;在本发明的思路下,以上实施例或者不同实施例中的技术特征之间也可以进行组合,步骤可以以任意顺序实现,并存在如上所述的本发明的不同方面的许多其它变化,为了简明它们没有在细节中提供。
本发明旨在涵盖落入权利要求的宽泛范围之内的所有这样的替换、修改和变型。因此,凡在本发明的精神和原则之内,所做的任何省略、修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.基于磨料水射流截割的煤岩裂纹风险预测方法,其特征在于,包括以下步骤:
S1:在磨料水射流切割开始前,测量煤岩的初始物理特性,初始物理特性包括煤岩基本结构和强度特性,其中,煤岩基本结构为密度、孔隙率,强度特性表现为抗压强度,同时建立煤岩基本结构-强度模型,具体包括:
S11:采用X射线计算机断层扫描技术测量煤岩的整体密度和孔隙率,通过分析扫描图像,计算煤岩内部的孔隙结构和总体密度分布,评估煤岩的孔隙特性和整体密实度;
S12:利用标准化的岩石力学试验方法测量煤岩的抗压强度,包括单轴或多轴压缩测试,以获得煤岩在不同应力条件下的强度参数;
S13:结合S11和S13的测试结果,综合分析煤岩基本结构和强度特性,建立煤岩基本结构-强度模型,用于描述煤岩的密度、孔隙率与其抗压强度之间的关系;
所述煤岩基本结构-强度模型基于非线性回归模型,具体包括:
收集大量煤岩样本数据,包括样本的密度、孔隙率和对应的抗压强度测量值;
对收集的数据进行统计分析,识别密度、孔隙率与抗压强度之间的基本关系和影响模式;
选择多项式回归模型来描述非线性关系,模型形式为:
其中,P是抗压强度,D是煤岩的密度,PT是孔隙率,a,b,c是模型参数,n,m是描述非线性关系的指数,ε是误差项,Dn、表示密度和孔隙率的非线性影响;
参数估计:使用最大似然估计方法来估计模型参数a,b,c,n,m;
模型验证:通过留一交叉验证技术验证模型的能力;
S2:基于S1的煤岩的初始物理特性测量,获得切割路径横截面的煤岩基本结构,基于煤岩基本结构-强度模型,分析切割路径上的抗压强度变化,从而预测裂纹风险;
S3:在磨料水射流切割过程中,实时监测沿切割路径的声发射特征,以捕捉切割中局部切割位置的煤岩微裂纹发生情况,基于微裂纹发生情况验证基本结构-强度模型;
S4:根据验证后基本结构-强度模型,以及S2中预测的裂纹风险,实时调整切割路径和/或磨料水射流的工作参数,包括改变切割速度、调整压力或流量,以优化切割过程并最小化裂纹风险。
2.根据权利要求1所述的基于磨料水射流截割的煤岩裂纹风险预测方法,其特征在于,所述最大似然估计方法具体包括:
基于多项式回归模型中的误差项ε遵循正态分布N(0,σ2),构建似然函数,似然函数是给定参数时,观测到的数据出现概率的函数,似然函数L表示为所有独立观测数据抗压强度的联合概率密度函数:
其中Pi,Di,分别是第i个观测样本的抗压强度、密度和孔隙率,σ2表示误差项的方差,反映模型中随机误差的大小;
最大化似然函数:参数a,b,c,m,m,σ2的估计值是使似然函数最大化的值,通过对似然函数取对数并应用数值优化方法来实现,因为对数转换后的乘积形式变为求和形式,表示为:
计算和优化:利用梯度下降法数值优化技术来求解对数似然函数的最大值,计算对数似然函数关于每个参数的导数,并迭代调整参数值以找到使对数似然函数最大化的参数组合。
3.根据权利要求2所述的基于磨料水射流截割的煤岩裂纹风险预测方法,其特征在于,所述留一交叉验证技术表示为LOOCV,用于验证多项式回归模型能力,验证步骤如下:
准备数据集,数据集中包括若干数据点,每个数据点包括煤岩的密度、孔隙率和抗压强度,数据集中的每个样本都将轮流用作测试数据;
实施LOOCV:对于数据集中的每一个样本:
A.留出一个样本:从数据集中选取一个样本作为测试数据,其余的样本作为训练数据;
B.建立模型:使用训练数据建立多项式回归模型;
C.模型预测:使用留出的测试样本来评估模型的预测能力,计算预测值与实际值之间的差异;
计算平均误差:重复实施LOOCV,直到每个样本都被用作测试数据一次,记录每一次的预测误差,然后计算误差的平均值,平均值表示模型的总体预测能力;
分析结果:分析LOOCV的结果来评估多项式回归模型的稳定性和预测能力,低平均误差值表示模型具有良好预测性能。
4.根据权利要求3所述的基于磨料水射流截割的煤岩裂纹风险预测方法,其特征在于,所述S2具体包括:
重建横截面图像:基于X射线计算机断层扫描技术扫描结果,通过图像重建技术生成煤岩的横截面图像,以显示预定切割路径区域的内部结构;
通过横截面图像,计算预定切割路径区域内煤岩密度,根据横截面图像中孔隙和裂缝的分布和大小,评估切割路径区域的孔隙率;
映射切割路径,将密度和孔隙率数据映射到预定的切割路径上,将具体的图像数据点与切割路径的几何位置对应起来,从而获得沿着切割路径的密度和孔隙率分布。
5.根据权利要求4所述的基于磨料水射流截割的煤岩裂纹风险预测方法,其特征在于,所述S2中的预测裂纹风险包括对于切割路径上的每一点或某区域,将相应的密度和孔隙率数据代入基本结构-强度模型中计算预测的抗压强度,获得切割过程中各个点或区域的抗压强度估计;
分析沿切割路径的抗压强度预测值,关注低于临界抗压强度的区域,为高风险裂纹区域。
6.根据权利要求5所述的基于磨料水射流截割的煤岩裂纹风险预测方法,其特征在于,所述S3具体包括:
安装声发射传感器在切割路径附近,用于实时捕捉煤岩在切割点的声发射信号,指示煤岩结构变化和裂纹形成活动;
收集传感器数据,并利用快速傅里叶变换FFT对振动数据进行频率分析,以及对声发射信号进行时频分析,确定切割中煤岩的动态响应特性;
结合实时监测数据和预定切割路径信息,使用支持向量机SVM模型识别煤岩微裂纹形成的早期迹象。
7.根据权利要求6所述的基于磨料水射流截割的煤岩裂纹风险预测方法,其特征在于,所述快速傅里叶变换FFT将时域信号转换为频域信号,表示为:
其中,F(k)是频域中的第k个频率分量,f(n)是时域中的第n个样本,N是样本总数;
所述时频分析用于理解声发射信号随时间的变化,通过短时傅里叶变换进行时频分析,以确定在切割过程中何时何地发生变化。
8.根据权利要求7所述的基于磨料水射流截割的煤岩裂纹风险预测方法,其特征在于,所述支持向量机SVM模型具体包括:
基于频率分析以及时频分析结果中提取特征,包括峰值幅度、能量、频率内容和持续时间,形成特征向量xi用于SVM训练和预测;
SVM模型建立:通过一个超平面分隔不同类别的数据点,最大化不同类别之间的边界,对于微裂纹识别问题,将其视为二分类问题,即微裂纹存在,即正类,或不存在,即负类;
通过最小化目标函数来训练SVM模型:
其中,w是超平面的法向量,C是正则化参数,ξi是松弛变量,N是训练样本数量,目标是找到超平面的法向量w和偏差项b使得目标函数最小化;
决策函数为:f(x)=wTx+b,当f(x)的值大于零时,x被分类为正类,即微裂纹存在,否则为负类,即微裂纹不存在;
使用训练好的SVM模型对新收集的声发射数据进行分类,以判断煤岩是否在某个切割位置形成微裂纹。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410508380.3A CN118094359B (zh) | 2024-04-26 | 2024-04-26 | 基于磨料水射流截割的煤岩裂纹风险预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410508380.3A CN118094359B (zh) | 2024-04-26 | 2024-04-26 | 基于磨料水射流截割的煤岩裂纹风险预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN118094359A CN118094359A (zh) | 2024-05-28 |
CN118094359B true CN118094359B (zh) | 2024-07-05 |
Family
ID=91159694
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202410508380.3A Active CN118094359B (zh) | 2024-04-26 | 2024-04-26 | 基于磨料水射流截割的煤岩裂纹风险预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN118094359B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN118376505B (zh) * | 2024-06-24 | 2024-08-20 | 内蒙古工业大学 | 一种用于浮石混凝土养护期的发育程度评估方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105510142A (zh) * | 2016-01-15 | 2016-04-20 | 太原理工大学 | 一种煤岩多相不同流体三轴压裂试验装置与试验方法 |
CN113982687A (zh) * | 2021-09-26 | 2022-01-28 | 安徽理工大学 | 冲击地压负反馈云模拟监测预警系统构建方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103940690B (zh) * | 2014-04-16 | 2016-01-13 | 中国矿业大学 | 一种多功能煤岩单齿截割试验装置及方法 |
CN111272883A (zh) * | 2020-03-04 | 2020-06-12 | 成都理工大学 | 一种基于声发射模型的岩石破裂模式智能探测识别方法 |
CN111859781A (zh) * | 2020-06-16 | 2020-10-30 | 重庆大学 | 一种采动煤岩多场响应快速获取方法 |
CN112949202A (zh) * | 2021-03-19 | 2021-06-11 | 交通运输部科学研究院 | 一种基于贝叶斯网络的岩爆概率预测方法 |
CN114861489B (zh) * | 2022-04-18 | 2023-09-19 | 河海大学 | rPCK代理模型辅助的结构动参数辨识方法 |
-
2024
- 2024-04-26 CN CN202410508380.3A patent/CN118094359B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105510142A (zh) * | 2016-01-15 | 2016-04-20 | 太原理工大学 | 一种煤岩多相不同流体三轴压裂试验装置与试验方法 |
CN113982687A (zh) * | 2021-09-26 | 2022-01-28 | 安徽理工大学 | 冲击地压负反馈云模拟监测预警系统构建方法 |
Also Published As
Publication number | Publication date |
---|---|
CN118094359A (zh) | 2024-05-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN118094359B (zh) | 基于磨料水射流截割的煤岩裂纹风险预测方法 | |
WO2019201176A1 (zh) | 列车部件裂纹损伤预测方法和装置 | |
US9909953B2 (en) | Method and apparatus for non-destructive detection of tire anomalies | |
CN108984893B (zh) | 一种基于梯度提升方法的趋势预测方法 | |
JPWO2018150616A1 (ja) | 異常音検出装置、異常度計算装置、異常音生成装置、異常音検出学習装置、異常信号検出装置、異常信号検出学習装置、これらの方法及びプログラム | |
WO2024041233A1 (zh) | 桥梁结构多因素耦合作用疲劳损伤与寿命的评估方法 | |
CN116728783B (zh) | 一种基于3d打印机的仿真方法及系统 | |
CN111967140B (zh) | 一种考虑混合不确定性的性能退化实验建模与分析方法 | |
CN117782957B (zh) | 一种船用电缆老化性能的测试方法及系统 | |
CN117233347B (zh) | 一种碳钢球化等级测量方法、系统及设备 | |
JP4299508B2 (ja) | 製造プロセスにおける操業と品質の関連分析装置、関連分析方法及びコンピュータ読み取り可能な記憶媒体 | |
JP2011145764A (ja) | 操業分析装置、操業分析方法、コンピュータプログラム及びコンピュータ読み取り可能な記憶媒体 | |
Gupta et al. | Statistical pattern analysis of ultrasonic signals for fatigue damage detection in mechanical structures | |
CN114382490A (zh) | 一种盾构机刀盘磨损评估预测方法 | |
Kundu et al. | PCA-ANN based approach for remaining useful life prediction for roller ball bearings | |
Wu et al. | Failure mode division and remaining useful life prognostics of multi-indicator systems with multi-fault | |
CN115392019A (zh) | 一种热轧工作辊实时剩余寿命的预测方法及其系统 | |
CN118246128B (zh) | 一种钢筋混凝土框架结构质量检测方法、介质及系统 | |
CN117893100B (zh) | 基于卷积神经网络的质量评价数据更新模型的构建方法 | |
Keshtgar et al. | Probabilistic approach for nondestructive detection of fatigue crack initiation and sizing | |
CN118094951B (zh) | 一种高压气瓶的内胆分析方法、设备及介质 | |
CN110929425B (zh) | 一种基于Lasso算法的纸张纵向伸长率的软测量方法 | |
US20230095525A1 (en) | System and method for evaluating residual life of components made of composite materials | |
Song et al. | Rail Fatigue Crack Classification Based on Improved DPC Algorithm Using Acoustic Emission Technology | |
CN115856099A (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 |