CN108256258A - Cemented fill mechanical response characteristic Forecasting Methodology based on SEM image - Google Patents

Cemented fill mechanical response characteristic Forecasting Methodology based on SEM image Download PDF

Info

Publication number
CN108256258A
CN108256258A CN201810098527.0A CN201810098527A CN108256258A CN 108256258 A CN108256258 A CN 108256258A CN 201810098527 A CN201810098527 A CN 201810098527A CN 108256258 A CN108256258 A CN 108256258A
Authority
CN
China
Prior art keywords
cemented filling
sample
sem
image
filling body
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.)
Granted
Application number
CN201810098527.0A
Other languages
Chinese (zh)
Other versions
CN108256258B (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.)
Xian University of Science and Technology
Original Assignee
Xian University of Science and Technology
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 Xian University of Science and Technology filed Critical Xian University of Science and Technology
Priority to CN201810098527.0A priority Critical patent/CN108256258B/en
Publication of CN108256258A publication Critical patent/CN108256258A/en
Application granted granted Critical
Publication of CN108256258B publication Critical patent/CN108256258B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00Investigating 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/22Investigating 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 measuring secondary emission from the material
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/23Clustering techniques
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06Power analysis or power optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Analytical Chemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Computer Hardware Design (AREA)
  • Health & Medical Sciences (AREA)
  • Biochemistry (AREA)
  • Pathology (AREA)
  • Chemical & Material Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Biology (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

本发明公开了一种基于SEM图像的胶结充填体力学响应特性预测方法,包括步骤:一、制成SEM扫描电镜样品;二、扫描形成SEM电镜扫描图像并存储到计算机中;三、对SEM电镜扫描图像进行高斯滤波处理;四、得到多个胶结充填体聚类图像;五、确定胶结充填体微观孔隙图,得到胶结充填体微观孔隙二值图;六、将高斯滤波处理后的SEM电镜扫描图像与胶结充填体微观孔隙二值图合并,得到测试样本图像;七、对测试样本图像进行正规化处理;八、将正规化测试样本图像输入预先构建的Tensorflow深度学习力学响应预测网络中,得到单轴力学响应预测结果。本发明预测效率高,耗费人力物力少,对研究胶结充填体强度和稳定性具有重要意义。

The invention discloses a method for predicting the mechanical response characteristics of a cemented filling body based on an SEM image, which comprises the steps of: 1. making a SEM scanning electron microscope sample; 2. scanning to form a SEM scanning electron microscope scanning image and storing it in a computer; The scanned image is processed by Gaussian filtering; 4. Obtain multiple cluster images of cemented filling bodies; 5. Determine the microscopic pore map of the cemented filling body, and obtain the binary image of the microscopic pores of the cemented filling body; 6. Scan the SEM electron microscope after Gaussian filtering The image is merged with the microscopic pore binary image of the cemented filling body to obtain the test sample image; 7. Normalize the test sample image; 8. Input the normalized test sample image into the pre-built Tensorflow deep learning mechanical response prediction network to obtain Uniaxial mechanical response prediction results. The invention has high prediction efficiency, consumes less manpower and material resources, and has great significance for studying the strength and stability of cemented filling bodies.

Description

基于SEM图像的胶结充填体力学响应特性预测方法Prediction method of mechanical response characteristics of cemented filling body based on SEM image

技术领域technical field

本发明属于胶结充填采矿技术领域,具体涉及一种基于SEM图像的胶结充填体力学响应特性预测方法。The invention belongs to the technical field of cemented filling mining, and in particular relates to a method for predicting the mechanical response characteristics of cemented filling bodies based on SEM images.

背景技术Background technique

随着国家科学技术的发展,对节能环保技术的要求也越来越高,传统的胶结充填采矿使用水泥作为胶凝材料,水泥的成本高达充填总成本的75%。通过研究发展,尾砂中含有活性氧化硅和氧化铝,使用尾砂代替部分水泥作为胶结材料,不仅能够降低尾砂的排放量,有效降低充填采矿的成本,还能够提高充填体强度,减少地面坍塌面积,对环境的保护也起着积极推动的作用。因此,选矿厂排出的尾砂逐渐成为矿山胶结充填的主要骨料。胶结充填体作为胶结充填采矿法的核心内容,它涉及到矿山安全和矿山经济效益。以“充填体作用机理、充填体强度、充填体的合理匹配及充填体的力学响应特性”为内容的充填体力学,几年来受到了采矿届的高度重视。近年来已召开了八届国际充填学术会议,在充填体力学的许多方面有了很大的进展,学术界普遍认为充填体力学性能是严重影响和制约胶结充填采矿法的关键因素。With the development of national science and technology, the requirements for energy saving and environmental protection technology are getting higher and higher. Traditional cemented backfill mining uses cement as the cementitious material, and the cost of cement is as high as 75% of the total backfill cost. Through research and development, tailings contain active silica and alumina, and using tailings instead of part of cement as cementing material can not only reduce the discharge of tailings, effectively reduce the cost of filling mining, but also improve the strength of the filling body and reduce the ground The collapsed area also plays an active role in promoting the protection of the environment. Therefore, the tailings discharged from the concentrator gradually become the main aggregate of cemented filling in mines. The cemented filling body is the core content of the cemented filling mining method, which involves mine safety and mine economic benefits. Filling body mechanics, with the content of "filling body mechanism, filling body strength, reasonable matching of filling body and mechanical response characteristics of filling body", has been highly valued by the mining industry for several years. In recent years, eight international academic conferences on filling have been held, and great progress has been made in many aspects of filling body mechanics. The academic circles generally believe that the mechanical properties of filling body are the key factors that seriously affect and restrict the cemented filling mining method.

采用尾砂实现不同水灰配比、不同养护龄期等对胶结充填体的力学性质具有直接的影响关系。尾砂作为充填采空区最常用的充填骨料之一,在解决充填骨料不足的同时,为极厚矿体矿柱回采时贫化率低、损失率大、“三下”资源开采安全性低以及深部岩体地压控制难等问题的解决提供了有效途径。许多研究者对尾砂膏体充填的成分配比、稳定过程及机械强度做了深入的研究。例如,Kesimal A等人研究了脱泥铜铅锌尾砂与膏体强度的关系,发现尾砂颗粒大小分布对胶结充填体强度有较大的影响;在2003年第16期第10卷的期刊《Minerals Engineering》(矿物工程)上发表了文章The effect of desliming bysedimentation on paste backfill performance(脱矿泥的矿体充填体沉积性能影响);Fall等人研究了养护温度对尾砂胶结充填体的强度的影响;在2010年第4期第10卷的期刊《Engineering Geology》(工程地质)上发表了文章A Contribution to understandingthe effects of curing temperature on the mechanical properties of minecemented tailings backfill(温度对尾砂胶结充填体力学性能影响的贡献);G Xiu等人采用不同比例尾砂胶结强度与不同浓度下在实验室进行了实验,揭示尾砂在微观方面化学反应机理,对充填体稳定性的宏观尺寸的影响研究;在2012年第6期第14卷的期刊《International Journal of Digital Content Technology&Its Applications》(数字内容技术及其应用)上发表了文章Microstructure Test and Macro Size Effect on theStability of Cemented Tailings Backfill(微观结构试验及宏观尺寸对胶结尾砂充填体的稳定性影响);李文臣等人通过胶砂实验,研究了硫酸盐对胶结充填体单轴抗压强度与弹性模量关系的影响;在2016年第1期第42卷的期刊《中国煤炭》上发表了文章《硫酸盐对胶结充填体单轴抗压强度与弹性模量关系影响研究》。但是,现有技术中,对胶结充填体力学响应特性预测多采用实验测试的方法,测试周期长、效率低,耗费的人力物力高,影响了新的胶结充填体的快速推广应用,容易造成采矿工期的拖延。The use of tailings to achieve different water-cement ratios and different curing ages have a direct impact on the mechanical properties of the cemented backfill. As one of the most commonly used filling aggregates for filling gobs, tailings not only solve the shortage of filling aggregates, but also have a low dilution rate and a large loss rate when mining pillars of extremely thick ore bodies, and the "three down" resource mining is safe It provides an effective way to solve the problems of low ground pressure and difficult ground pressure control in deep rock mass. Many researchers have done in-depth research on the composition ratio, stabilization process and mechanical strength of tailings paste filling. For example, Kesimal A et al. studied the relationship between deslimed copper-lead-zinc tailings and paste strength, and found that the particle size distribution of tailings has a great influence on the strength of cemented filling bodies; The effect of desliming bysedimentation on paste backfill performance was published on "Minerals Engineering"; Fall et al. studied the effect of curing temperature on the strength of tailings cemented backfill In 2010, the journal "Engineering Geology" (Engineering Geology), Issue 4, Volume 10, published an article A Contribution to understanding the effects of curing temperature on the mechanical properties of minecemented tailings backfill. The contribution of mechanical properties); G Xiu et al. used different proportions of tailings cementation strength and different concentrations to conduct experiments in the laboratory to reveal the chemical reaction mechanism of tailings in the microscopic aspect and the influence of the macroscopic size on the stability of the filling body. ; Published the article Microstructure Test and Macro Size Effect on the Stability of Cemented Tailings Backfill (Microstructure Test and Macro Size Effect on the Stability of Cemented Tailings Backfill) in the journal "International Journal of Digital Content Technology & Its Applications" (Digital Content Technology and Its Applications) in the 6th Volume 14 of 2012 The influence of macroscopic size on the stability of cement-tailed sand filling); Li Wenchen et al. studied the influence of sulfate on the relationship between uniaxial compressive strength and elastic modulus of cemented filling through mortar experiments; in the first issue of 2016 The 42nd volume of the journal "China Coal" published an article "Research on the Effect of Sulfate on the Relationship between Uniaxial Compressive Strength and Elastic Modulus of Cemented Filling Body". However, in the existing technology, the prediction of the mechanical response characteristics of cemented filling bodies is mostly based on experimental testing methods. The test period is long, the efficiency is low, and the cost of manpower and material resources is high, which affects the rapid promotion and application of new cemented filling bodies. Delay in construction period.

发明内容Contents of the invention

本发明所要解决的技术问题在于针对上述现有技术中的不足,提供一种基于SEM图像的胶结充填体力学响应特性预测方法,其方法步骤简单,设计新颖合理,实现方便快捷,预测效率高,周期短,耗费的人力物力少,对于研究胶结充填体的强度和稳定性具有重要意义,实用性强,应用范围广,推广应用价值高。The technical problem to be solved by the present invention is to provide a method for predicting the mechanical response characteristics of cemented filling bodies based on SEM images, which has simple steps, novel and reasonable design, convenient and fast implementation, and high prediction efficiency. The cycle is short, and it consumes less manpower and material resources. It is of great significance for the study of the strength and stability of the cemented filling body. It has strong practicability, a wide range of applications, and high promotion and application value.

为解决上述技术问题,本发明采用的技术方案是:一种基于SEM图像的胶结充填体力学响应特性预测方法,其特征在于,该方法包括以下步骤:In order to solve the above-mentioned technical problems, the technical solution adopted in the present invention is: a method for predicting the mechanical response characteristics of cemented filling bodies based on SEM images, characterized in that the method includes the following steps:

步骤一、从胶结充填体试样上取一部分制成SEM扫描电镜样品;Step 1, taking a part of the cemented filling body sample to make a SEM scanning electron microscope sample;

步骤二、采用SEM扫描电镜对SEM扫描电镜样品进行扫描,形成SEM电镜扫描图像并存储到计算机中;Step 2, using SEM scanning electron microscope to scan the SEM scanning electron microscope sample, forming a SEM scanning electron microscope scanning image and storing it in the computer;

步骤三、所述计算机调用高斯滤波处理模块对SEM电镜扫描图像进行高斯滤波处理,得到高斯滤波处理后的SEM电镜扫描图像;Step 3, the computer calls the Gaussian filter processing module to perform Gaussian filter processing on the SEM electron microscope scanned image, and obtains the SEM electron microscope scanned image after the Gaussian filter process;

步骤四、所述计算机调用FCM模糊聚类处理模块对进行高斯滤波处理后的SEM电镜扫描图像进行孔隙图像提取,得到与聚类中心数目相等的多个胶结充填体聚类图像;Step 4, the computer invokes the FCM fuzzy clustering processing module to extract the pore image from the Gaussian filter-processed SEM electron microscope scanning image, and obtain a plurality of cemented filling body clustering images equal to the number of clustering centers;

步骤五、所述计算机将灰度值最小一类的胶结充填体聚类图像确定为胶结充填体微观孔隙图,并对胶结充填体微观孔隙图进行二值化处理,得到胶结充填体微观孔隙二值图;Step 5. The computer determines the cluster image of the cemented filling body with the smallest gray value as the microscopic pore map of the cemented filling body, and performs binary processing on the microscopic pore map of the cemented filling body to obtain the microscopic pore size of the cemented filling body. value map;

步骤六、所述计算机将步骤三中得到的高斯滤波处理后的SEM电镜扫描图像与步骤五中得到的胶结充填体微观孔隙二值图进行合并,得到测试样本图像;Step 6, the computer merges the Gaussian filter-processed SEM scanning image obtained in step 3 with the binary image of the microscopic pores of the cemented filling body obtained in step 5 to obtain a test sample image;

步骤七、所述计算机对测试样本图像进行正规化处理,形成像素为960×960的正规化测试样本图像;Step 7, the computer normalizes the test sample image to form a normalized test sample image with pixels of 960×960;

步骤八、所述计算机将步骤七中得到的正规化测试样本图像输入预先构建的Tensorflow深度学习力学响应预测网络中,得到单轴力学响应预测结果。Step 8: The computer inputs the normalized test sample image obtained in step 7 into the pre-built Tensorflow deep learning mechanical response prediction network to obtain a uniaxial mechanical response prediction result.

上述的基于SEM图像的胶结充填体力学响应特性预测方法,其特征在于:步骤一中所述SEM扫描电镜样品的长度、宽度和高度均为10mm。The above method for predicting the mechanical response characteristics of cemented filling bodies based on SEM images is characterized in that: the length, width and height of the SEM scanning electron microscope sample in step 1 are all 10 mm.

上述的基于SEM图像的胶结充填体力学响应特性预测方法,其特征在于:步骤三中所述计算机调用高斯滤波处理模块对SEM电镜扫描图像进行高斯滤波处理采用的公式为L(x,y)=I(x,y)*G(x,y),其中,I(x,y)表示SEM电镜扫描图像,G(x,y)为高斯滤波函数,L(x,y)为高斯滤波处理后的SEM电镜扫描图像,x为图像的横坐标,y为图像的纵坐标。The above method for predicting the mechanical response characteristics of cemented filling bodies based on SEM images is characterized in that: the computer described in step 3 calls the Gaussian filter processing module to perform Gaussian filter processing on the SEM electron microscope scanning image. The formula used is L(x, y)= I(x,y)*G(x,y), where I(x,y) represents the SEM scanning image, G(x,y) is the Gaussian filter function, and L(x,y) is the Gaussian filter The SEM electron microscope scanning image, x is the abscissa of the image, and y is the ordinate of the image.

上述的基于SEM图像的胶结充填体力学响应特性预测方法,其特征在于:步骤四中所述计算机调用FCM模糊聚类处理模块对进行高斯滤波处理后的SEM电镜扫描图像进行孔隙图像提取,得到与聚类中心数目相等的多个胶结充填体聚类图像的具体过程为:The above method for predicting the mechanical response characteristics of cemented filling bodies based on SEM images is characterized in that: the computer calls the FCM fuzzy clustering processing module in step 4 to extract the pore images from the SEM electron microscope scanning images after the Gaussian filter processing, and obtain the same as The specific process of clustering images of multiple cemented fillings with the same number of cluster centers is as follows:

步骤401、定义采用基于样本加权的FCM模糊聚类算法,目标函数为满足极值的约束条件为其中,U为模糊矩阵且U=[u11,u22,…,ucn],uik为矩阵U的元素且uik表示第k个样本点属于第i类的隶属度,n为样本点总数,c为聚类中心数目;V={v1,v2,...vc}是c个类的聚类中心,wk为样本点xk的权值,dik为样本点xk到中心点vi的欧式距离,vi为V的元素,xk为样本集X的第k个样本点且X={x1,x2,...xn},m为隶属度uik的权重指数且m>1;Step 401, define the FCM fuzzy clustering algorithm based on sample weighting, and the objective function is The constraint condition for satisfying the extremum value is Among them, U is a fuzzy matrix and U=[u 11 ,u 22 ,…,u cn ], u ik is an element of matrix U and u ik represents the membership degree of the kth sample point belonging to the i-th class, and n is the sample point The total number, c is the number of cluster centers; V={v 1 ,v 2 ,...v c } is the cluster centers of c classes, w k is the weight of sample point x k , d ik is the sample point x The Euclidean distance from k to the center point v i , v i is the element of V, x k is the kth sample point of the sample set X and X={x 1 ,x 2 ,...x n }, m is the degree of membership The weight index of u ik and m>1;

步骤402、设置聚类中心数目c的值、隶属度uik的权重指数m的值和最小迭代误差ε的值;Step 402, setting the value of the number of cluster centers c, the value of the weight index m of the degree of membership u ik and the value of the minimum iteration error ε;

步骤403、用公式更新样本点xk的权值wk;uτj为矩阵U的元素且uτj表示第j个样本点属于第τ类的隶属度,1≤τ≤c,1≤j≤n;vτ为V的元素;uij为矩阵U的元素且uij表示第j个样本点属于第i类的隶属度;Step 403, use the formula Update the weight w k of the sample point x k ; u τj is an element of the matrix U and u τj represents the membership degree of the jth sample point belonging to the τth class, 1≤τ≤c, 1≤j≤n; v τj is The element of V; u ij is the element of matrix U and u ij represents the membership degree of the j-th sample point belonging to the i-th class;

步骤404、用公式更新uik;其中,drk为样本点xk到中心点vr的欧式距离,1≤r≤c;Step 404, use the formula Update u ik ; among them, d rk is the Euclidean distance from the sample point x k to the center point v r , 1≤r≤c;

步骤405、用公式更新viStep 405, use the formula update v i ;

步骤406、判断是否满足||J(t+1)-J(t)||<ε,当满足||J(t+1)-J(t)||<ε时,聚类停止,提取得到与聚类中心数目相等的多个胶结充填体聚类图像;否则,返回步骤403;其中,t为时间。Step 406. Judging whether ||J (t+1) -J(t)||<ε is satisfied, and when ||J (t+1) -J(t)||<ε is satisfied, the clustering stops and extraction Obtain multiple cluster images of cemented filling bodies equal to the number of cluster centers; otherwise, return to step 403; wherein, t is time.

上述的基于SEM图像的胶结充填体力学响应特性预测方法,其特征在于:步骤402中设置聚类中心数目c的值为4,设置隶属度uik的权重指数m的值为2,设置最小迭代误差ε的取值为0.3。The above method for predicting the mechanical response characteristics of cemented filling bodies based on SEM images is characterized in that: in step 402, the value of the number c of cluster centers is set to 4, the value of the weight index m of the degree of membership u ik is set to 2, and the minimum iteration The value of the error ε is 0.3.

上述的基于SEM图像的胶结充填体力学响应特性预测方法,其特征在于:步骤八中所述Tensorflow深度学习力学响应特性预测网络的构建方法为:The above method for predicting the mechanical response characteristics of cemented filling bodies based on SEM images is characterized in that: the construction method of the Tensorflow deep learning mechanical response characteristic prediction network described in step 8 is:

步骤801、从多个编号后的各个胶结充填体试样上取一部分制成SEM扫描电镜样品,剩余部分作为单轴抗压强度测试样品;并对多个SEM扫描电镜样品和多个抗压强度测试样品一一对应编号;Step 801, take a part of each cemented filling body sample after multiple numbers to make a SEM scanning electron microscope sample, and the remaining part is used as a uniaxial compressive strength test sample; and for multiple SEM scanning electron microscope samples and multiple compressive strength samples The test samples are numbered one by one;

步骤802、采用胶结充填体单轴抗压强度测试装置分别对多个抗压强度测试样品进行单轴抗压强度测试,并对测得的多个抗压强度测试样品的单轴抗压强度取均值,得到胶结充填体试样的单轴抗压强度;Step 802: Use the uniaxial compressive strength test device for the cemented filling body to perform uniaxial compressive strength tests on multiple compressive strength test samples respectively, and take the measured uniaxial compressive strength of the multiple compressive strength test samples as The average value is obtained to obtain the uniaxial compressive strength of the cemented filling body sample;

步骤803、获取Tensorflow深度学习力学响应特性预测网络的训练样本图像,具体过程为:Step 803, obtaining the training sample image of the Tensorflow deep learning mechanical response characteristic prediction network, the specific process is:

步骤8031、采用SEM扫描电镜分别对多个SEM扫描电镜样品进行多次扫描,形成多个SEM电镜扫描图像并存储到计算机中;所述SEM电镜扫描图像的数量至少为150个;Step 8031, using the SEM scanning electron microscope to scan multiple SEM scanning electron microscope samples respectively to form multiple SEM scanning electron microscope scanning images and store them in the computer; the number of the SEM scanning electron microscope scanning images is at least 150;

步骤8032、所述计算机调用高斯滤波处理模块分别对多个SEM电镜扫描图像进行高斯滤波处理,得到多个高斯滤波处理后的SEM电镜扫描图像;Step 8032, the computer invokes the Gaussian filter processing module to perform Gaussian filter processing on a plurality of SEM electron microscope scanned images respectively, to obtain a plurality of Gaussian filtered SEM electron microscope scanned images;

步骤8033、所述计算机调用FCM模糊聚类处理模块分别对多个进行高斯滤波处理后的SEM电镜扫描图像进行孔隙图像提取,得到多组胶结充填体聚类图像,每组胶结充填体聚类图像中胶结充填体聚类图像的数量与聚类中心数目相等;Step 8033, the computer invokes the FCM fuzzy clustering processing module to extract pore images from the Gaussian-filtered SEM electron microscope scanning images to obtain multiple groups of cemented filling cluster images, each group of cemented filling cluster images The number of cluster images of medium cemented filling body is equal to the number of cluster centers;

步骤8034、所述计算机将每组胶结充填体聚类图像中灰度值最小一类的胶结充填体聚类图像确定为胶结充填体微观孔隙图,并对多个胶结充填体微观孔隙图进行二值化处理,得到多个胶结充填体微观孔隙二值图;Step 8034, the computer determines the cemented backfill cluster image with the smallest gray value in each group of cemented backfill cluster images as the microscopic pore map of the cemented backfill, and double-checks the multiple cemented backfill microscopic pore maps Value-based processing to obtain multiple microscopic pore binary maps of cemented filling bodies;

步骤8035、所述计算机将步骤8032中得到的多个高斯滤波处理后的SEM电镜扫描图像与步骤8034中得到多个的胶结充填体微观孔隙二值图按照编号对应进行合并,得到多个训练样本图像;Step 8035, the computer combines the multiple Gaussian-filtered SEM scanning images obtained in step 8032 with the multiple binary images of microscopic pores of the cemented filling body obtained in step 8034 according to the corresponding numbers to obtain multiple training samples image;

步骤804、所述计算机分别对多个训练样本图像进行正规化处理,形成多个像素为960×960的正规化训练样本图像;Step 804, the computer performs normalization processing on a plurality of training sample images respectively to form a plurality of normalized training sample images with a pixel size of 960×960;

步骤805、所述计算机构建一个卷积网络核的层数为五层、输入层为正规化训练样本图像、输出层为正规化训练样本图像对应的单轴抗压强度的Tensorflow深度学习网络,将其存储的多个正规化训练样本图像作为训练样本,对Tensorflow深度学习网络进行训练,得到Tensorflow深度学习力学响应特性预测网络;所述Tensorflow深度学习力学特性响应预测网络五层卷积网络核的大小从一层到第五层分别为3x3,2x2,3x3,2x2,2x2。Step 805, the computer constructs a Tensorflow deep learning network with five layers of convolutional network cores, an input layer of normalized training sample images, and an output layer of uniaxial compressive strength corresponding to the normalized training sample images. A plurality of normalized training sample images stored in it are used as training samples to train the Tensorflow deep learning network to obtain the Tensorflow deep learning mechanical response characteristic prediction network; the Tensorflow deep learning mechanical characteristic response prediction network is the size of the five-layer convolutional network kernel From the first floor to the fifth floor are 3x3, 2x2, 3x3, 2x2, 2x2 respectively.

上述的基于SEM图像的胶结充填体力学响应特性预测方法,其特征在于:步骤802中所述胶结充填体单轴抗压强度测试装置包括座垫和固定连接在座垫顶部的多根拉杆,以及用于给胶结充填体试样施加轴向压力的轴向加压传力机构和用于给轴向加压传力机构提供动力的轴向加压动力系统;所述座垫的底部固定连接有多个底座,所述座垫的顶部设置有用于放置胶结充填体试样的试样放置槽,所述座垫上位于所述试样放置槽的中心位置处设置有排水阀;多根拉杆的中部设置有用于固定多根拉杆的固定架,多根拉杆的上部固定连接有顶部装载板;所述轴向加压传力机构包括安装在顶部装载板上的气缸,所述气缸的活塞杆向下设置,所述气缸的活塞杆底部连接有压力传递板;所述轴向加压动力系统包括压缩空气气源和加压控制器,以及一端与压缩空气气源连接、另一端与气缸连接的气体输送管;所述气体输送管上从连接压缩空气气源到连接气缸的位置依次设置有气动三联件、压力传感器和气缸控制电磁阀,所述压力传感器与加压控制器的输入端连接,所述气缸控制电磁阀与加压控制器的输出端连接,所述加压控制器通过通信模块与计算机连接。The above method for predicting mechanical response characteristics of cemented filling body based on SEM images is characterized in that: the uniaxial compressive strength test device for cemented filling body described in step 802 includes a seat cushion and a plurality of tie rods fixedly connected to the top of the seat cushion, and An axial pressure force transmission mechanism for applying axial pressure to the cemented filling body sample and an axial pressure power system for providing power to the axial pressure force transmission mechanism; the bottom of the seat cushion is fixedly connected with A plurality of bases, the top of the seat cushion is provided with a sample placement groove for placing the cemented filling body sample, and the seat cushion is provided with a drain valve at the center of the sample placement groove; the plurality of pull rods The middle part is provided with a fixed frame for fixing multiple pull rods, and the upper part of the multiple pull rods is fixedly connected with a top loading plate; the axial pressure force transmission mechanism includes a cylinder installed on the top loading plate, and the piston rod of the cylinder is directed toward the The bottom of the piston rod of the cylinder is connected with a pressure transmission plate; the axial pressurization power system includes a compressed air source and a pressurization controller, and a device with one end connected to the compressed air source and the other end connected to the cylinder Gas delivery pipe; the gas delivery pipe is sequentially provided with a pneumatic triple piece, a pressure sensor and a cylinder control solenoid valve from the position connecting the compressed air source to the connection cylinder, and the pressure sensor is connected to the input end of the pressurization controller, The cylinder control solenoid valve is connected with the output end of the pressurization controller, and the pressurization controller is connected with the computer through the communication module.

上述的基于SEM图像的胶结充填体力学响应特性预测方法,其特征在于:所述试样放置槽内设置有套装在胶结充填体试样底部的O型密封圈,所述座垫上设置有位于所述试样放置槽周围的多孔石。The method for predicting the mechanical response characteristics of the cemented filling body based on the SEM image is characterized in that: the sample placement groove is provided with an O-ring set on the bottom of the cemented filling body sample, and the seat cushion is provided with a The sample is placed in a porous stone around the groove.

上述的基于SEM图像的胶结充填体力学响应特性预测方法,其特征在于:所述加压控制器为可编程逻辑控制器,所述通信模块为RS-485通信模块。The above method for predicting mechanical response characteristics of cemented filling bodies based on SEM images is characterized in that: the pressurization controller is a programmable logic controller, and the communication module is an RS-485 communication module.

上述的基于SEM图像的胶结充填体力学响应特性预测方法,其特征在于:步骤802中所述采用胶结充填体单轴抗压强度测试装置分别对多个抗压强度测试样品进行单轴抗压强度测试,其中对每个抗压强度测试样品进行单轴抗压强度测试的具体过程为:The above method for predicting mechanical response characteristics of cemented filling body based on SEM images is characterized in that: in step 802, the uniaxial compressive strength test device for cemented filling body is used to test the uniaxial compressive strength of multiple compressive strength test samples respectively. Test, wherein the specific process of uniaxial compressive strength test for each compressive strength test sample is:

步骤8021、将O型密封圈放入所述试样放置槽内后,将胶结充填体试样放入所述试样放置槽内,使胶结充填体试样的中心与气缸的活塞杆和压力传递板的中心相对应;并在座垫上放入位于所述试样放置槽周围的多孔石;Step 8021. After putting the O-ring into the sample placement groove, put the cemented filling body sample into the sample placement groove, so that the center of the cemented filling body sample is in contact with the piston rod and pressure of the cylinder. Corresponding to the center of the transfer plate; and put porous stones located around the sample placement groove on the seat cushion;

步骤8022、打开压缩空气气源,通过调节气动三联件调节压缩空气气源输出的压缩空气的气压,加压控制器通过控制气缸控制电磁阀换向,控制气缸的活塞杆向下或向上运动,对胶结充填体试样施加压力或卸载压力,将胶结充填体试样破裂时加压控制器采集到的压力传感器检测的压力值记为F,加压控制器将压力值F传输给计算机,计算机根据公式计算得到抗压强度测试样品的单轴抗压强度P;其中,S为抗压强度测试样品的顶面面积;当气缸的活塞杆向下运动时,带动压力传递板向下运动,通过压力传递板给胶结充填体试样施加压力,当气缸的活塞杆向上运动时,带动压力传递板向上运动,压力传递板离开胶结充填体试样的上表面,卸载压力。Step 8022, turn on the compressed air source, adjust the air pressure of the compressed air output by the compressed air source by adjusting the pneumatic triple unit, the pressurization controller controls the reversing of the solenoid valve by controlling the cylinder, and controls the piston rod of the cylinder to move downward or upward, Apply pressure or unload pressure to the cemented filling body sample, record the pressure value detected by the pressure sensor collected by the pressure controller when the cemented filling body sample ruptures as F, and the pressure controller transmits the pressure value F to the computer, and the computer According to the formula Calculate the uniaxial compressive strength P of the compressive strength test sample; where, S is the top surface area of the compressive strength test sample; when the piston rod of the cylinder moves downward, it drives the pressure transmission plate to move downward, through the pressure transmission The plate exerts pressure on the cemented filling body sample. When the piston rod of the cylinder moves upward, it drives the pressure transmission plate to move upward, and the pressure transmission plate leaves the upper surface of the cemented filling body sample to unload the pressure.

本发明与现有技术相比具有以下优点:Compared with the prior art, the present invention has the following advantages:

1、本发明采用扫描电子显微镜(Scanning Electronic Microscopy,SEM)扫描采集样品图像,采用FCM模糊聚类处理方法提取胶结充填体微观孔隙图,再采用Tensorflow深度学习网络建立图像到力学响应特性之间的端到端的预测模型,从而预测胶结充填体的力学响应特性,方法步骤简单,设计新颖合理,实现方便快捷,预测效率高,周期短,耗费的人力物力少。1. The present invention uses scanning electron microscope (Scanning Electronic Microscopy, SEM) to scan and collect sample images, uses FCM fuzzy clustering processing method to extract the microscopic pore map of cemented filling body, and then uses Tensorflow deep learning network to establish the relationship between the image and the mechanical response characteristics The end-to-end prediction model is used to predict the mechanical response characteristics of the cemented filling body. The method has simple steps, novel and reasonable design, convenient and quick implementation, high prediction efficiency, short cycle time, and less manpower and material resources.

2、本发明在采用FCM模糊聚类处理方法提取胶结充填体微观孔隙图像前,还采用了高斯滤波方法对SEM电镜扫描图像进行高斯滤波处理,有助于得到更准确的预测结果。2. Before using the FCM fuzzy clustering processing method to extract the microscopic pore image of the cemented filling body, the present invention also uses the Gaussian filter method to perform Gaussian filter processing on the SEM electron microscope scanning image, which helps to obtain more accurate prediction results.

3、本发明采用FCM模糊聚类处理方法提取胶结充填体微观孔隙图像,能够避免样本空间中不同样本矢量对聚类结果的不同影响。3. The present invention adopts the FCM fuzzy clustering processing method to extract the microscopic pore image of the cemented filling body, which can avoid different influences of different sample vectors in the sample space on the clustering results.

4、本发明采用将高斯滤波处理后的SEM电镜扫描图像与胶结充填体微观孔隙二值图合并的方法得到测试样本图像,作为Tensorflow深度学习力学响应预测网络的输入,能够避免原始SEM电镜扫描图像中冗余信息的影响,提高了力学响应特性预测精度。4. The present invention adopts the method of merging the Gaussian filter-processed SEM electron microscope scanning image with the cemented filling body microscopic pore binary image to obtain the test sample image as the input of the Tensorflow deep learning mechanical response prediction network, which can avoid the original SEM electron microscope scanning image The impact of redundant information in the medium improves the prediction accuracy of mechanical response characteristics.

5、本发明采用Tensorflow深度学习网络构建了SEM电镜扫描图像与胶结充填体力学响应特性之间的关系,一次进行构建Tensorflow深度学习力学响应预测网络,能够多次方便快捷的使用,使得进行胶结充填体力学响应特性预测时,无需再多次做实验测试,只需将SEM电镜扫描图像采集到计算机中,即可自动完成力学响应特性预测的整个过程,方便快捷。5. The present invention uses the Tensorflow deep learning network to construct the relationship between the SEM electron microscope scanning image and the mechanical response characteristics of the cemented filling body, and constructs the Tensorflow deep learning mechanical response prediction network once, which can be used conveniently and quickly for many times, so that the cemented filling can be carried out When predicting the mechanical response characteristics of the body, there is no need to do experimental tests many times, just collect the scanning images of the SEM electron microscope into the computer, and the whole process of predicting the mechanical response characteristics can be automatically completed, which is convenient and fast.

6、单轴抗压强度是反映充填体力学性能的一个重要参数,它能在一定程度上反应充填体的强度和稳定性;本发明构建Tensorflow深度学习力学响应预测网络时,输出为单轴抗压强度,并将单轴抗压强度作为单轴力学响应预测结果,对于研究胶结充填体的强度和稳定性具有重要意义。6. The uniaxial compressive strength is an important parameter reflecting the mechanical properties of the filling body, and it can reflect the strength and stability of the filling body to a certain extent; when the present invention constructs the Tensorflow deep learning mechanical response prediction network, the output is the uniaxial compressive strength It is of great significance to study the strength and stability of cemented filling bodies by using the uniaxial compressive strength as the prediction result of uniaxial mechanical response.

7、本发明采用自主研发制造的胶结充填体单轴抗压强度测试装置对抗压强度测试样品进行抗压强度测试,胶结充填体单轴抗压强度测试装置的结构简单,实现及使用操作方便,且能够测得准确的单轴抗压强度。7. The present invention adopts the self-developed and manufactured uniaxial compressive strength test device for the cemented filling body to test the compressive strength test sample. The uniaxial compressive strength test device for the cemented filling body has a simple structure and is easy to implement and operate , and can measure the accurate uniaxial compressive strength.

8、本发明研究的力学响应特性是胶结充填体的重要特性,胶结充填体又是胶结充填采矿法的核心内容,因此本发明的方法不仅能够为研究新型胶结充填体贡献力量,还能够为降低尾砂的排放量、降低充填采矿的成本、保护环境、提高矿石回采率、缓解深井高温、优化矿区环境和控制地表沉降等贡献力量;本发明的实用性强,应用范围广,推广应用价值高。8. The mechanical response characteristics studied by the present invention are important characteristics of the cemented backfill, and the cemented backfill is the core content of the cemented backfill mining method. Therefore, the method of the present invention can not only contribute to the study of new cemented backfills, but also contribute to the reduction of Discharge of tailings, reduce the cost of filling mining, protect the environment, increase the ore recovery rate, alleviate the high temperature of deep wells, optimize the environment of the mining area, and control the surface subsidence; the present invention has strong practicability, wide application range, and high application value .

综上所述,本发明方法步骤简单,设计新颖合理,实现方便快捷,预测效率高,周期短,耗费的人力物力少,对于研究胶结充填体的强度和稳定性具有重要意义,实用性强,应用范围广,推广应用价值高。In summary, the method of the present invention has simple steps, novel and reasonable design, convenient and quick implementation, high prediction efficiency, short cycle, and less manpower and material resources. The application range is wide, and the promotion and application value is high.

下面通过附图和实施例,对本发明的技术方案做进一步的详细描述。The technical solutions of the present invention will be described in further detail below with reference to the accompanying drawings and embodiments.

附图说明Description of drawings

图1为本发明的方法流程框图。Fig. 1 is a flow chart of the method of the present invention.

图2为本发明具体实施例中高斯滤波处理后的SEM电镜扫描图像。Fig. 2 is a SEM electron microscope scanning image after Gaussian filter processing in a specific embodiment of the present invention.

图3A为本发明具体实施例中对进行高斯滤波处理后的SEM电镜扫描图像进行FCM模糊聚类处理得到的亮的胶结充填体聚类图像。Fig. 3A is a bright cemented filling body clustering image obtained by performing FCM fuzzy clustering processing on the SEM electron microscope scanning image after Gaussian filtering processing in a specific embodiment of the present invention.

图3B为本发明具体实施例中对进行高斯滤波处理后的SEM电镜扫描图像进行FCM模糊聚类处理得到的较亮的胶结充填体聚类图像。Fig. 3B is a brighter cemented filling cluster image obtained by performing FCM fuzzy clustering on the Gaussian-filtered SEM scanning image in a specific embodiment of the present invention.

图3C为本发明具体实施例中对进行高斯滤波处理后的SEM电镜扫描图像进行FCM模糊聚类处理得到的暗的胶结充填体聚类图像。Fig. 3C is a cluster image of dark cemented filling bodies obtained by performing FCM fuzzy clustering on the SEM scanning image after Gaussian filtering in a specific embodiment of the present invention.

图3D为本发明具体实施例中对进行高斯滤波处理后的SEM电镜扫描图像进行FCM模糊聚类处理得到的最暗的胶结充填体聚类图像。Fig. 3D is the darkest cemented filling body clustering image obtained by performing FCM fuzzy clustering processing on the SEM electron microscope scanning image after Gaussian filtering processing in a specific embodiment of the present invention.

图4为本发明具体实施例中的胶结充填体微观孔隙二值图。Fig. 4 is a binary map of microscopic pores of a cemented filling body in a specific embodiment of the present invention.

图5为本发明具体实施例中的测试样本图像。Fig. 5 is a test sample image in a specific embodiment of the present invention.

图6为本发明胶结充填体单轴抗压强度测试装置的结构示意图。Fig. 6 is a schematic structural view of the uniaxial compressive strength testing device of the cemented filling body of the present invention.

附图标记说明:Explanation of reference signs:

1—气体输送管; 2—气缸; 3—压力传递板;1—gas delivery pipe; 2—cylinder; 3—pressure transmission plate;

4—压缩空气气源; 5—排水阀; 6—O型密封圈;4—compressed air source; 5—drain valve; 6—O-ring;

7—多孔石; 8—拉杆; 9—顶部装载板;7—porous stone; 8—tie rod; 9—top loading plate;

10—座垫; 11—试样放置槽; 12—气动三联件;10—seat cushion; 11—sample placement slot; 12—pneumatic triple piece;

13—压力传感器; 14—气缸控制电磁阀; 15—底座;13—pressure sensor; 14—cylinder control solenoid valve; 15—base;

16—通信模块; 17—计算机; 18—加压控制器;16—communication module; 17—computer; 18—pressurization controller;

19—胶结充填体试样。19—Consolidated filling body sample.

具体实施方式Detailed ways

如图1所示,本发明的基于SEM图像的胶结充填体力学响应特性预测方法,包括以下步骤:As shown in Figure 1, the method for predicting the mechanical response characteristics of cemented filling bodies based on SEM images of the present invention includes the following steps:

步骤一、从胶结充填体试样19上取一部分制成SEM扫描电镜样品;Step 1, taking a part from the cemented filling body sample 19 to make a SEM scanning electron microscope sample;

具体实施时,制成SEM扫描电镜样品还进行了多次喷碳处理。During specific implementation, the SEM scanning electron microscope sample was also subjected to multiple carbon spraying treatments.

本实施例中,步骤一中所述SEM扫描电镜样品的长度、宽度和高度均为10mm。In this embodiment, the length, width and height of the SEM scanning electron microscope sample in step 1 are all 10 mm.

步骤二、采用SEM扫描电镜对SEM扫描电镜样品进行扫描,形成SEM电镜扫描图像并存储到计算机17中;Step 2, using SEM scanning electron microscope to scan the SEM scanning electron microscope sample, forming a SEM scanning electron microscope scanning image and storing it in the computer 17;

步骤三、所述计算机17调用高斯滤波处理模块对SEM电镜扫描图像进行高斯滤波处理,得到高斯滤波处理后的SEM电镜扫描图像;Step 3, the computer 17 invokes the Gaussian filter processing module to perform Gaussian filter processing on the SEM electron microscope scanned image, and obtains the SEM electron microscope scanned image after the Gaussian filter process;

本实施例中,步骤三中所述计算机17调用高斯滤波处理模块对SEM电镜扫描图像进行高斯滤波处理采用的公式为L(x,y)=I(x,y)*G(x,y),其中,I(x,y)表示SEM电镜扫描图像,G(x,y)为高斯滤波函数,L(x,y)为高斯滤波处理后的SEM电镜扫描图像,x为图像的横坐标,y为图像的纵坐标。In this embodiment, the computer 17 in step 3 invokes the Gaussian filter processing module to perform Gaussian filter processing on the SEM electron microscope scanning image using the formula L(x, y)=I(x, y)*G(x, y) , wherein, I(x, y) represents the SEM electron microscope scanning image, G(x, y) is the Gaussian filter function, L(x, y) is the SEM electron microscope scanning image after the Gaussian filter processing, and x is the abscissa of the image, y is the vertical coordinate of the image.

本实施例中,高斯滤波处理后的SEM电镜扫描图像如图2所示。In this embodiment, the SEM electron microscope scanning image after Gaussian filter processing is shown in FIG. 2 .

步骤四、所述计算机17调用FCM模糊聚类处理模块对进行高斯滤波处理后的SEM电镜扫描图像进行孔隙图像提取,得到与聚类中心数目相等的多个胶结充填体聚类图像;Step 4, the computer 17 invokes the FCM fuzzy clustering processing module to extract pore images from the Gaussian filter-processed SEM electron microscope scanning image, and obtain a plurality of cemented filling body cluster images equal to the number of cluster centers;

本实施例中,步骤四中所述计算机17调用FCM模糊聚类处理模块对进行高斯滤波处理后的SEM电镜扫描图像进行孔隙图像提取,得到与聚类中心数目相等的多个胶结充填体聚类图像的具体过程为:In this embodiment, the computer 17 in step 4 invokes the FCM fuzzy clustering processing module to extract pore images from the Gaussian filter-processed SEM electron microscope scanning image, and obtain multiple cemented filling body clusters equal to the number of cluster centers The specific process of the image is:

步骤401、定义采用基于样本加权的FCM模糊聚类算法,目标函数为满足极值的约束条件为其中,U为模糊矩阵且U=[u11,u22,…,ucn],uik为矩阵U的元素且uik表示第k个样本点属于第i类的隶属度,n为样本点总数(在图像中对应每个坐标点的灰度值),c为聚类中心数目(在图像中根据图像亮度进行分类);V={v1,v2,...vc}是c个类的聚类中心,wk为样本点xk的权值,dik为样本点xk到中心点vi的欧式距离,vi为V的元素,xk为样本集X的第k个样本点且X={x1,x2,...xn},m为隶属度uik的权重指数且m>1;采用基于样本加权的FCM模糊聚类算法,能够避免样本空间中不同样本矢量对聚类结果的不同影响;Step 401, define the FCM fuzzy clustering algorithm based on sample weighting, and the objective function is The constraint condition for satisfying the extremum value is Among them, U is a fuzzy matrix and U=[u 11 ,u 22 ,…,u cn ], u ik is an element of matrix U and u ik represents the membership degree of the kth sample point belonging to the i-th class, and n is the sample point The total number (the gray value corresponding to each coordinate point in the image), c is the number of cluster centers (classified according to the brightness of the image in the image); V={v 1 ,v 2 ,...v c } is c The cluster center of each class, w k is the weight of the sample point x k , d ik is the Euclidean distance from the sample point x k to the center point v i , v i is the element of V, x k is the kth of the sample set X sample points and X={x 1 ,x 2 ,...x n }, m is the weight index of the membership degree u ik and m>1; the FCM fuzzy clustering algorithm based on sample weighting can avoid the The different influences of different sample vectors on the clustering results;

步骤402、设置聚类中心数目c的值、隶属度uik的权重指数m的值和最小迭代误差ε的值;Step 402, setting the value of the number of cluster centers c, the value of the weight index m of the degree of membership u ik and the value of the minimum iteration error ε;

本实施例中,步骤402中设置聚类中心数目c的值为4,设置隶属度uik的权重指数m的值为2,设置最小迭代误差ε的取值为0.3。In this embodiment, in step 402, the value of the number c of cluster centers is set to 4, the value of the weight index m of the degree of membership u ik is set to 2, and the value of the minimum iteration error ε is set to 0.3.

步骤403、用公式更新样本点xk的权值wk;uτj为矩阵U的元素且uτj表示第j个样本点属于第τ类的隶属度,1≤τ≤c,1≤j≤n;vτ为V的元素;uij为矩阵U的元素且uij表示第j个样本点属于第i类的隶属度;Step 403, use the formula Update the weight w k of the sample point x k ; u τj is an element of the matrix U and u τj represents the membership degree of the jth sample point belonging to the τth class, 1≤τ≤c, 1≤j≤n; v τj is The element of V; u ij is the element of matrix U and u ij represents the membership degree of the j-th sample point belonging to the i-th class;

步骤404、用公式更新uik;其中,drk为样本点xk到中心点vr的欧式距离,1≤r≤c;Step 404, use the formula Update u ik ; among them, d rk is the Euclidean distance from the sample point x k to the center point v r , 1≤r≤c;

步骤405、用公式更新viStep 405, use the formula update v i ;

步骤406、判断是否满足||J(t+1)-J(t)||<ε,当满足||J(t+1)-J(t)||<ε时,聚类停止,提取得到与聚类中心数目相等的多个胶结充填体聚类图像;否则,返回步骤403;其中,t为时间。Step 406. Judging whether ||J (t+1) -J(t)||<ε is satisfied, and when ||J (t+1) -J(t)||<ε is satisfied, the clustering stops and extraction Obtain multiple cluster images of cemented filling bodies equal to the number of cluster centers; otherwise, return to step 403; wherein, t is time.

本实施例中,聚类中心数目c的值为4,得到四个胶结充填体聚类图像,如图3A~图3D所示,分别为亮、较亮、暗和最暗四个胶结充填体聚类图像。In this embodiment, the value of the number c of cluster centers is 4, and four clustered images of cemented backfills are obtained, as shown in Fig. 3A to Fig. 3D, which are bright, brighter, dark and darkest four cemented backfills respectively Cluster images.

步骤五、所述计算机17将灰度值最小一类的胶结充填体聚类图像(最暗的胶结充填体聚类图像)确定为胶结充填体微观孔隙图,并对胶结充填体微观孔隙图进行二值化处理,得到胶结充填体微观孔隙二值图;Step 5, the computer 17 determines the cemented filling body cluster image with the smallest gray value (the darkest cemented filling body cluster image) as the cemented filling body microscopic pore map, and performs the microscopic pore map of the cemented filling body Binary processing to obtain the binary map of the microscopic pores of the cemented filling body;

本实施例中,灰度值最小的胶结充填体聚类图像为图3D,即图3D为胶结充填体微观孔隙图,得到的胶结充填体微观孔隙二值图如图4所示。In this embodiment, the cluster image of the cemented filling body with the smallest gray value is shown in Figure 3D, that is, Figure 3D is the microscopic pore map of the cemented filling body, and the obtained binary image of the microscopic pores of the cemented filling body is shown in Figure 4 .

步骤六、所述计算机17将步骤三中得到的高斯滤波处理后的SEM电镜扫描图像与步骤五中得到的胶结充填体微观孔隙二值图进行合并,得到测试样本图像;由于原始SEM电镜扫描图像中的冗余信息较多,即一般认为黑色区域为孔隙部分,由于受扫描时不同因素的影响,可能存在扫描的明暗程度不同一等情况,因此,本发明将步骤三中得到的高斯滤波处理后的SEM电镜扫描图像与步骤五中得到的胶结充填体微观孔隙二值图进行合并,得到测试样本图像,能够避免冗余信息的影响,提高单轴力学响应特性预测精度。Step 6, the computer 17 merges the Gaussian filter processed SEM scanning image obtained in step 3 with the binary image of the microscopic pores of the cemented filling body obtained in step 5 to obtain a test sample image; since the original SEM scanning image There is a lot of redundant information in , that is, it is generally considered that the black area is a pore part. Due to the influence of different factors during scanning, there may be situations where the degree of lightness and darkness of the scan are different. Therefore, the present invention processes the Gaussian filter obtained in step 3 The final SEM electron microscope scanning image is merged with the microscopic pore binary image of the cemented filling body obtained in step 5 to obtain a test sample image, which can avoid the influence of redundant information and improve the prediction accuracy of uniaxial mechanical response characteristics.

例如,如图5所示,即为本实施例中得到的测试样本图像。For example, as shown in FIG. 5 , it is the test sample image obtained in this embodiment.

步骤七、所述计算机17对测试样本图像进行正规化处理,形成像素为960×960的正规化测试样本图像;Step 7, the computer 17 normalizes the test sample image to form a normalized test sample image with pixels of 960×960;

步骤八、所述计算机17将步骤七中得到的正规化测试样本图像输入预先构建的Tensorflow深度学习力学响应预测网络中,得到单轴力学响应预测结果。Step 8. The computer 17 inputs the normalized test sample image obtained in step 7 into the pre-built Tensorflow deep learning mechanical response prediction network to obtain a uniaxial mechanical response prediction result.

本实施例中,步骤八中所述Tensorflow深度学习力学响应特性预测网络的构建方法为:In this embodiment, the construction method of the Tensorflow deep learning mechanical response characteristic prediction network described in step eight is:

步骤801、从多个编号后的各个胶结充填体试样19上取一部分制成SEM扫描电镜样品,剩余部分作为单轴抗压强度测试样品;并对多个SEM扫描电镜样品和多个抗压强度测试样品一一对应编号;例如,多个胶结充填体试样19的编号分别为1、2、…、N,多个SEM扫描电镜样品的编号分别为A1、A2、…、AN,多个单轴抗压强度测试样品的编号分别为B1、B2、…、BN;Step 801, take a part of each cemented filling body sample 19 after multiple numbers to make a SEM scanning electron microscope sample, and the remaining part is used as a uniaxial compressive strength test sample; and multiple SEM scanning electron microscope samples and multiple compression The strength test samples are numbered one by one; for example, the numbers of multiple cemented filling body samples 19 are 1, 2, ..., N, and the numbers of multiple SEM scanning electron microscope samples are A1, A2, ..., AN, and multiple The numbers of the uniaxial compressive strength test samples are B1, B2, ..., BN;

具体实施时,制成SEM扫描电镜样品还进行了多次喷碳处理。During specific implementation, the SEM scanning electron microscope sample was also subjected to multiple carbon spraying treatments.

步骤802、采用胶结充填体单轴抗压强度测试装置分别对多个抗压强度测试样品进行单轴抗压强度测试,并对测得的多个抗压强度测试样品的单轴抗压强度取均值,得到胶结充填体试样19的单轴抗压强度;Step 802: Use the uniaxial compressive strength test device for the cemented filling body to perform uniaxial compressive strength tests on multiple compressive strength test samples respectively, and take the measured uniaxial compressive strength of the multiple compressive strength test samples as The mean value is obtained to obtain the uniaxial compressive strength of the cemented filling body sample 19;

本实施例中,如图6所示,步骤802中所述胶结充填体单轴抗压强度测试装置包括座垫10和固定连接在座垫10顶部的多根拉杆8,以及用于给胶结充填体试样19施加轴向压力的轴向加压传力机构和用于给轴向加压传力机构提供动力的轴向加压动力系统;所述座垫10的底部固定连接有多个底座15,所述座垫10的顶部设置有用于放置胶结充填体试样19的试样放置槽,所述座垫10上位于所述试样放置槽的中心位置处设置有排水阀5;多根拉杆8的中部设置有用于固定多根拉杆8的固定架11,多根拉杆8的上部固定连接有顶部装载板9;所述轴向加压传力机构包括安装在顶部装载板9上的气缸2,所述气缸2的活塞杆向下设置,所述气缸2的活塞杆底部连接有压力传递板3;所述轴向加压动力系统包括压缩空气气源4和加压控制器18,以及一端与压缩空气气源4连接、另一端与气缸2连接的气体输送管1;所述气体输送管1上从连接压缩空气气源4到连接气缸2的位置依次设置有气动三联件12、压力传感器13和气缸控制电磁阀14,所述压力传感器13与加压控制器18的输入端连接,所述气缸控制电磁阀14与加压控制器18的输出端连接,所述加压控制器18通过通信模块16与计算机17连接。具体实施时,所述压力传递板3由橡胶制成。采用橡胶制成压力传递板3,一方面,能够分配气缸2的活塞杆传递的压力,使压力更加均匀地施加在胶结充填体试样19顶部;另一方面,压力传递板3传递压力到胶结充填体试样19上时,不会对胶结充填体试样19的顶面造成损伤。In this embodiment, as shown in FIG. 6, the uniaxial compressive strength test device for the cemented filling body described in step 802 includes a seat cushion 10 and a plurality of pull rods 8 fixedly connected to the top of the seat cushion 10, and is used to give the cemented filling The body sample 19 applies an axial pressure force transmission mechanism for axial pressure and an axial pressure power system for providing power to the axial pressure force transmission mechanism; the bottom of the seat cushion 10 is fixedly connected with a plurality of bases 15. The top of the seat cushion 10 is provided with a sample placement groove for placing the cemented filling body sample 19, and the seat cushion 10 is provided with a drain valve 5 at the center of the sample placement groove; The middle part of the pull rods 8 is provided with a fixed frame 11 for fixing multiple pull rods 8, and the upper part of the multiple pull rods 8 is fixedly connected with a top loading plate 9; the axial pressure transmission mechanism includes a cylinder installed on the top loading plate 9 2. The piston rod of the cylinder 2 is set downwards, and the bottom of the piston rod of the cylinder 2 is connected with a pressure transmission plate 3; the axial pressurization power system includes a compressed air source 4 and a pressurization controller 18, and One end is connected to the compressed air source 4, and the other end is connected to the gas cylinder 2; the gas delivery pipe 1 is sequentially provided with a pneumatic triple piece 12, a pressure The sensor 13 and the cylinder control solenoid valve 14, the pressure sensor 13 is connected to the input end of the pressurization controller 18, the cylinder control solenoid valve 14 is connected to the output end of the pressurization controller 18, and the pressurization controller 18 It is connected with the computer 17 through the communication module 16 . During specific implementation, the pressure transmission plate 3 is made of rubber. The pressure transmission plate 3 is made of rubber. On the one hand, it can distribute the pressure transmitted by the piston rod of the cylinder 2, so that the pressure is more evenly applied to the top of the cemented filling body sample 19; on the other hand, the pressure transmission plate 3 transmits pressure to the cemented filling body. When the filling body sample 19 is put on, the top surface of the cemented filling body sample 19 will not be damaged.

本实施例中,所述试样放置槽内设置有套装在胶结充填体试样19底部的O型密封圈6,所述座垫10上设置有位于所述试样放置槽周围的多孔石7。通过设置O型密封圈6,能够防止给胶结充填体试样19施加轴向压力时胶结充填体试样19与座垫10硬接触造成胶结充填体试样19的损伤。通过设置多孔石7,能够吸收胶结充填体试样19渗出的水。In this embodiment, the O-shaped sealing ring 6 set on the bottom of the cemented filling body sample 19 is arranged in the sample placement groove, and the porous stone 7 located around the sample placement groove is arranged on the seat cushion 10. . By providing the O-ring 6 , it is possible to prevent damage to the cemented filling body sample 19 caused by hard contact between the cemented filling body sample 19 and the seat cushion 10 when axial pressure is applied to the cemented filling body sample 19 . By providing the porous stone 7, the water seeped from the cemented filling body sample 19 can be absorbed.

本实施例中,所述加压控制器18为可编程逻辑控制器,所述通信模块16为RS-485通信模块。In this embodiment, the pressurization controller 18 is a programmable logic controller, and the communication module 16 is an RS-485 communication module.

本实施例中,步骤802中所述采用胶结充填体单轴抗压强度测试装置分别对多个抗压强度测试样品进行单轴抗压强度测试,其中对每个抗压强度测试样品进行单轴抗压强度测试的具体过程为:In this embodiment, as described in step 802, the uniaxial compressive strength test device for the cemented filling body is used to perform the uniaxial compressive strength test on a plurality of compressive strength test samples respectively, wherein the uniaxial compressive strength test is performed on each compressive strength test sample. The specific process of compressive strength test is as follows:

步骤8021、将O型密封圈6放入所述试样放置槽内后,将胶结充填体试样19放入所述试样放置槽内,使胶结充填体试样19的中心与气缸2的活塞杆和压力传递板3的中心相对应;并在座垫10上放入位于所述试样放置槽周围的多孔石7;Step 8021, after putting the O-ring 6 into the sample placement groove, put the cemented filling body sample 19 into the sample placement groove, so that the center of the cemented filling body sample 19 is aligned with the center of the cylinder 2 The piston rod corresponds to the center of the pressure transmission plate 3; and the porous stone 7 positioned around the sample placement groove is placed on the seat cushion 10;

步骤8022、打开压缩空气气源4,通过调节气动三联件12调节压缩空气气源4输出的压缩空气的气压,加压控制器18通过控制气缸控制电磁阀14换向,控制气缸2的活塞杆向下或向上运动,对胶结充填体试样19施加压力或卸载压力,将胶结充填体试样19破裂时加压控制器18采集到的压力传感器13检测的压力值记为F,加压控制器18将压力值F传输给计算机17,计算机17根据公式计算得到抗压强度测试样品的单轴抗压强度P;其中,S为抗压强度测试样品的顶面面积;当气缸2的活塞杆向下运动时,带动压力传递板3向下运动,通过压力传递板3给胶结充填体试样19施加压力,当气缸2的活塞杆向上运动时,带动压力传递板3向上运动,压力传递板3离开胶结充填体试样19的上表面,卸载压力。Step 8022, open the compressed air source 4, adjust the air pressure of the compressed air output by the compressed air source 4 by adjusting the pneumatic triple unit 12, and the pressurization controller 18 controls the reversing of the solenoid valve 14 by controlling the cylinder to control the piston rod of the cylinder 2 Move downward or upward, apply pressure to the cemented filling body sample 19 or unload the pressure, record the pressure value detected by the pressure sensor 13 collected by the pressure controller 18 when the cemented filling body sample 19 ruptures as F, pressurization control The device 18 transmits the pressure value F to the computer 17, and the computer 17 according to the formula Calculate the uniaxial compressive strength P of the compressive strength test sample; where, S is the top surface area of the compressive strength test sample; when the piston rod of the cylinder 2 moves downward, it drives the pressure transmission plate 3 to move downward, through The pressure transmission plate 3 applies pressure to the cemented filling body sample 19. When the piston rod of the cylinder 2 moves upward, it drives the pressure transmission plate 3 to move upward, and the pressure transmission plate 3 leaves the upper surface of the cemented filling body sample 19 to unload the pressure.

步骤803、获取Tensorflow深度学习力学响应特性预测网络的训练样本图像,具体过程为:Step 803, obtaining the training sample image of the Tensorflow deep learning mechanical response characteristic prediction network, the specific process is:

步骤8031、采用SEM扫描电镜分别对多个SEM扫描电镜样品进行多次扫描,形成多个SEM电镜扫描图像并存储到计算机17中;所述SEM电镜扫描图像的数量至少为150个;Step 8031, using the SEM scanning electron microscope to scan multiple SEM scanning electron microscope samples respectively, forming multiple SEM scanning electron microscope scanning images and storing them in the computer 17; the number of the SEM scanning electron microscope scanning images is at least 150;

步骤8032、所述计算机17调用高斯滤波处理模块分别对多个SEM电镜扫描图像进行高斯滤波处理,得到多个高斯滤波处理后的SEM电镜扫描图像;Step 8032, the computer 17 invokes the Gaussian filter processing module to perform Gaussian filter processing on a plurality of SEM electron microscope scanned images respectively, to obtain a plurality of Gaussian filtered SEM electron microscope scanned images;

本实施例中,所述计算机17调用高斯滤波处理模块分别对多个SEM电镜扫描图像进行高斯滤波处理采用的公式为L(x,y)=I(x,y)*G(x,y),其中,I(x,y)表示SEM电镜扫描图像,G(x,y)为高斯滤波函数,L(x,y)为高斯滤波处理后的SEM电镜扫描图像,x为图像的横坐标,y为图像的纵坐标。In this embodiment, the computer 17 calls the Gaussian filter processing module to perform Gaussian filter processing on a plurality of SEM electron microscope scanning images respectively. The formula adopted is L(x, y)=I(x, y)*G(x, y) , wherein, I(x, y) represents the SEM electron microscope scanning image, G(x, y) is the Gaussian filter function, L(x, y) is the SEM electron microscope scanning image after the Gaussian filter processing, and x is the abscissa of the image, y is the vertical coordinate of the image.

步骤8033、所述计算机17调用FCM模糊聚类处理模块分别对多个进行高斯滤波处理后的SEM电镜扫描图像进行孔隙图像提取,得到多组胶结充填体聚类图像,每组胶结充填体聚类图像中胶结充填体聚类图像的数量与聚类中心数目相等;Step 8033, the computer 17 invokes the FCM fuzzy clustering processing module to extract pore images from a plurality of SEM electron microscope scanning images after Gaussian filter processing, and obtain multiple groups of cemented filling cluster images, each group of cemented filling clusters The number of cluster images of cemented filling bodies in the image is equal to the number of cluster centers;

本实施例中,所述计算机17调用FCM模糊聚类处理模块对进行高斯滤波处理后的SEM电镜扫描图像进行孔隙图像提取的具体过程与步骤四相同。In this embodiment, the specific process of the computer 17 invoking the FCM fuzzy clustering processing module to extract the pore image from the Gaussian filtered SEM image is the same as step 4.

步骤8034、所述计算机17将每组胶结充填体聚类图像中灰度值最小一类的胶结充填体聚类图像(最暗的胶结充填体聚类图像)确定为胶结充填体微观孔隙图,并对多个胶结充填体微观孔隙图进行二值化处理,得到多个胶结充填体微观孔隙二值图;Step 8034, the computer 17 determines the cemented backfill cluster image with the smallest gray value (the darkest cemented backfill cluster image) in each group of cemented backfill cluster images as the cemented backfill microscopic pore map, And carry out binarization processing on the microscopic pore maps of multiple cemented filling bodies, and obtain the microscopic pore binary maps of multiple cemented filling bodies;

步骤8035、所述计算机17将步骤8032中得到的多个高斯滤波处理后的SEM电镜扫描图像与步骤8034中得到多个的胶结充填体微观孔隙二值图按照编号对应进行合并,得到多个训练样本图像;Step 8035, the computer 17 merges the multiple Gaussian-filtered SEM scanning images obtained in step 8032 with the multiple binary images of the microscopic pores of the cemented filling body obtained in step 8034 according to the corresponding numbers to obtain multiple training sample image;

步骤804、所述计算机17分别对多个训练样本图像进行正规化处理,形成多个像素为960×960的正规化训练样本图像;Step 804, the computer 17 respectively performs normalization processing on a plurality of training sample images to form a plurality of normalized training sample images with a pixel size of 960×960;

具体实施时,当训练样本图像大于正规化要得到的像素960×960时,对图像进行等比例缩小,得到正规化训练样本图像;当训练样本图像小于正规化要得到的像素960×960时,对图像采用白色进行边沿扩充,得到正规化训练样本图像。During specific implementation, when the training sample image is larger than the pixel 960×960 to be obtained by normalization, the image is scaled down to obtain a normalized training sample image; when the training sample image is smaller than the pixel 960×960 to be obtained by normalization, The edge of the image is expanded with white, and the normalized training sample image is obtained.

步骤805、所述计算机17构建一个卷积网络核的层数为五层、输入层为正规化训练样本图像、输出层为正规化训练样本图像对应的单轴抗压强度的Tensorflow深度学习网络,将其存储的多个正规化训练样本图像作为训练样本,对Tensorflow深度学习网络进行训练,得到Tensorflow深度学习力学响应特性预测网络;所述Tensorflow深度学习力学特性响应预测网络五层卷积网络核的大小从一层到第五层分别为3x3,2x2,3x3,2x2,2x2。Step 805, the computer 17 constructs a Tensorflow deep learning network with five layers of convolutional network core layers, an input layer of normalized training sample images, and an output layer of uniaxial compressive strength corresponding to the normalized training sample images, A plurality of normalized training sample images stored in it are used as training samples, and the Tensorflow deep learning network is trained to obtain a Tensorflow deep learning mechanical response characteristic prediction network; the Tensorflow deep learning mechanical characteristic response prediction network has five layers of convolutional network cores The sizes from the first floor to the fifth floor are 3x3, 2x2, 3x3, 2x2, 2x2 respectively.

综上所述,本发明采用扫描电子显微镜扫描采集样品图像,采用高斯滤波方法对SEM电镜扫描图像进行高斯滤波处理,采用FCM模糊聚类处理方法提取胶结充填体微观孔隙图,将高斯滤波处理后的SEM电镜扫描图像与胶结充填体微观孔隙二值图合并的方法得到测试样本图像,作为Tensorflow深度学习力学响应预测网络的输入,再采用Tensorflow深度学习网络建立图像到力学响应特性之间的端到端的预测模型,从而预测胶结充填体的力学响应特性,方法步骤简单,设计新颖合理,实现方便快捷,预测效率高,周期短,耗费的人力物力少。In summary, the present invention uses a scanning electron microscope to scan and collect sample images, uses a Gaussian filter method to perform Gaussian filter processing on the SEM electron microscope scanned image, uses FCM fuzzy clustering processing method to extract the microscopic pore map of the cemented filling body, and performs Gaussian filter processing. The test sample image is obtained by merging the SEM electron microscope scanning image and the microscopic pore binary image of the cemented filling body, which is used as the input of the Tensorflow deep learning mechanical response prediction network, and then the Tensorflow deep learning network is used to establish the end-to-end relationship between the image and the mechanical response characteristics. The end prediction model is used to predict the mechanical response characteristics of the cemented filling body. The method has simple steps, novel and reasonable design, convenient and quick implementation, high prediction efficiency, short cycle time, and less manpower and material resources.

以上所述,仅是本发明的较佳实施例,并非对本发明作任何限制,凡是根据本发明技术实质对以上实施例所作的任何简单修改、变更以及等效结构变化,均仍属于本发明技术方案的保护范围内。The above are only preferred embodiments of the present invention, and do not limit the present invention in any way. All simple modifications, changes and equivalent structural changes made to the above embodiments according to the technical essence of the present invention still belong to the technical aspects of the present invention. within the scope of protection of the scheme.

Claims (10)

1.一种基于SEM图像的胶结充填体力学响应特性预测方法,其特征在于,该方法包括以下步骤:1. A method for predicting the mechanical response characteristics of cemented filling bodies based on SEM images, characterized in that the method may further comprise the steps: 步骤一、从胶结充填体试样(19)上取一部分制成SEM扫描电镜样品;Step 1, taking a part from the cemented filling body sample (19) to make a SEM scanning electron microscope sample; 步骤二、采用SEM扫描电镜对SEM扫描电镜样品进行扫描,形成SEM电镜扫描图像并存储到计算机(17)中;Step 2, using SEM scanning electron microscope to scan the SEM scanning electron microscope sample, forming a SEM scanning electron microscope scanning image and storing it in the computer (17); 步骤三、所述计算机(17)调用高斯滤波处理模块对SEM电镜扫描图像进行高斯滤波处理,得到高斯滤波处理后的SEM电镜扫描图像;Step 3, the computer (17) calls the Gaussian filter processing module to carry out Gaussian filter processing to the SEM electron microscope scanned image, and obtains the SEM electron microscope scanned image after the Gaussian filter process; 步骤四、所述计算机(17)调用FCM模糊聚类处理模块对进行高斯滤波处理后的SEM电镜扫描图像进行孔隙图像提取,得到与聚类中心数目相等的多个胶结充填体聚类图像;Step 4, the computer (17) invokes the FCM fuzzy clustering processing module to extract the pore image from the Gaussian filter-processed SEM electron microscope scanning image, and obtain a plurality of cemented filling body cluster images equal to the number of cluster centers; 步骤五、所述计算机(17)将灰度值最小一类的胶结充填体聚类图像确定为胶结充填体微观孔隙图,并对胶结充填体微观孔隙图进行二值化处理,得到胶结充填体微观孔隙二值图;Step 5. The computer (17) determines the cluster image of the cemented filling body with the smallest gray value as the microscopic pore map of the cemented filling body, and performs binary processing on the microscopic pore map of the cemented filling body to obtain the cemented filling body Microscopic pore binary map; 步骤六、所述计算机(17)将步骤三中得到的高斯滤波处理后的SEM电镜扫描图像与步骤五中得到的胶结充填体微观孔隙二值图进行合并,得到测试样本图像;Step 6, the computer (17) merges the Gaussian filter processed SEM scanning image obtained in step 3 with the binary image of the microscopic pores of the cemented filling body obtained in step 5 to obtain a test sample image; 步骤七、所述计算机(17)对测试样本图像进行正规化处理,形成像素为960×960的正规化测试样本图像;Step 7, the computer (17) normalizes the test sample image to form a normalized test sample image with pixels of 960×960; 步骤八、所述计算机(17)将步骤七中得到的正规化测试样本图像输入预先构建的Tensorflow深度学习力学响应预测网络中,得到单轴力学响应预测结果。Step 8. The computer (17) inputs the normalized test sample image obtained in step 7 into the pre-built Tensorflow deep learning mechanical response prediction network to obtain a uniaxial mechanical response prediction result. 2.按照权利要求1所述的基于SEM图像的胶结充填体力学响应特性预测方法,其特征在于:步骤一中所述SEM扫描电镜样品的长度、宽度和高度均为10mm。2. The method for predicting mechanical response characteristics of cemented filling bodies based on SEM images according to claim 1, wherein the length, width and height of the SEM scanning electron microscope sample in step 1 are all 10 mm. 3.按照权利要求1所述的基于SEM图像的胶结充填体力学响应特性预测方法,其特征在于:步骤三中所述计算机(17)调用高斯滤波处理模块对SEM电镜扫描图像进行高斯滤波处理采用的公式为L(x,y)=I(x,y)*G(x,y),其中,I(x,y)表示SEM电镜扫描图像,G(x,y)为高斯滤波函数,L(x,y)为高斯滤波处理后的SEM电镜扫描图像,x为图像的横坐标,y为图像的纵坐标。3. according to the method for predicting the mechanical response characteristic of cemented filling body based on SEM image according to claim 1, it is characterized in that: described in step 3, computer (17) transfers Gaussian filter processing module to carry out Gaussian filter process to SEM electron microscope scanning image using The formula is L(x,y)=I(x,y)*G(x,y), wherein, I(x,y) represents the SEM electron microscope scanning image, G(x,y) is the Gaussian filter function, L (x, y) is the SEM scanning image after Gaussian filtering, x is the abscissa of the image, and y is the ordinate of the image. 4.按照权利要求1所述的基于SEM图像的胶结充填体力学响应特性预测方法,其特征在于:步骤四中所述计算机(17)调用FCM模糊聚类处理模块对进行高斯滤波处理后的SEM电镜扫描图像进行孔隙图像提取,得到与聚类中心数目相等的多个胶结充填体聚类图像的具体过程为:4. according to the method for predicting the mechanical response characteristics of cemented filling bodies based on SEM images according to claim 1, it is characterized in that: the computer (17) described in step 4 calls the FCM fuzzy clustering processing module to carry out the SEM after Gaussian filter processing The specific process of extracting the pore image from the scanning electron microscope image to obtain multiple cemented filling body cluster images equal to the number of cluster centers is as follows: 步骤401、定义采用基于样本加权的FCM模糊聚类算法,目标函数为满足极值的约束条件为其中,U为模糊矩阵且U=[u11,u22,…,ucn],uik为矩阵U的元素且uik表示第k个样本点属于第i类的隶属度,n为样本点总数,c为聚类中心数目;V={v1,v2,...vc}是c个类的聚类中心,wk为样本点xk的权值,dik为样本点xk到中心点vi的欧式距离,vi为V的元素,xk为样本集X的第k个样本点且X={x1,x2,...xn},m为隶属度uik的权重指数且m>1;Step 401, define the FCM fuzzy clustering algorithm based on sample weighting, and the objective function is The constraint condition for satisfying the extremum value is Among them, U is a fuzzy matrix and U=[u 11 ,u 22 ,…,u cn ], u ik is an element of matrix U and u ik represents the membership degree of the kth sample point belonging to the i-th class, and n is the sample point The total number, c is the number of cluster centers; V={v 1 ,v 2 ,...v c } is the cluster centers of c classes, w k is the weight of sample point x k , d ik is the sample point x The Euclidean distance from k to the center point v i , v i is the element of V, x k is the kth sample point of the sample set X and X={x 1 ,x 2 ,...x n }, m is the degree of membership The weight index of u ik and m>1; 步骤402、设置聚类中心数目c的值、隶属度uik的权重指数m的值和最小迭代误差ε的值;Step 402, setting the value of the number of cluster centers c, the value of the weight index m of the degree of membership u ik and the value of the minimum iteration error ε; 步骤403、用公式更新样本点xk的权值wk;uτj为矩阵U的元素且uτj表示第j个样本点属于第τ类的隶属度,1≤τ≤c,1≤j≤n;vτ为V的元素;uij为矩阵U的元素且uij表示第j个样本点属于第i类的隶属度;Step 403, use the formula Update the weight w k of the sample point x k ; u τj is an element of the matrix U and u τj represents the membership degree of the jth sample point belonging to the τth class, 1≤τ≤c, 1≤j≤n; v τj is The element of V; u ij is the element of matrix U and u ij represents the membership degree of the j-th sample point belonging to the i-th class; 步骤404、用公式更新uik;其中,drk为样本点xk到中心点vr的欧式距离,1≤r≤c;Step 404, use the formula Update u ik ; among them, d rk is the Euclidean distance from the sample point x k to the center point v r , 1≤r≤c; 步骤405、用公式更新viStep 405, use the formula update v i ; 步骤406、判断是否满足||J(t+1)-J(t)||<ε,当满足||J(t+1)-J(t)||<ε时,聚类停止,提取得到与聚类中心数目相等的多个胶结充填体聚类图像;否则,返回步骤403;其中,t为时间。Step 406. Judging whether ||J (t+1) -J(t)||<ε is satisfied, and when ||J (t+1) -J(t)||<ε is satisfied, the clustering stops and extraction Obtain multiple cluster images of cemented filling bodies equal to the number of cluster centers; otherwise, return to step 403; wherein, t is time. 5.按照权利要求4所述的基于SEM图像的胶结充填体力学响应特性预测方法,其特征在于:步骤402中设置聚类中心数目c的值为4,设置隶属度uik的权重指数m的值为2,设置最小迭代误差ε的取值为0.3。5. according to the method for predicting the mechanical response characteristics of cemented filling bodies based on SEM images according to claim 4, it is characterized in that: in step 402, the value of the number c of cluster centers is set to 4, and the weight index m of the degree of membership u ik is set The value is 2, and the value of the minimum iteration error ε is set to 0.3. 6.按照权利要求1所述的基于SEM图像的胶结充填体力学响应特性预测方法,其特征在于:步骤八中所述Tensorflow深度学习力学响应特性预测网络的构建方法为:6. according to the method for predicting mechanical response characteristics of cemented filling body based on SEM image according to claim 1, it is characterized in that: the construction method of Tensorflow deep learning mechanical response characteristic prediction network described in step 8 is: 步骤801、从多个编号后的各个胶结充填体试样(19)上取一部分制成SEM扫描电镜样品,剩余部分作为单轴抗压强度测试样品;并对多个SEM扫描电镜样品和多个抗压强度测试样品一一对应编号;Step 801, take a part from each cemented filling body sample (19) after multiple numbers to make a SEM scanning electron microscope sample, and the remaining part is used as a uniaxial compressive strength test sample; and multiple SEM scanning electron microscope samples and multiple The compressive strength test samples correspond to the numbers one by one; 步骤802、采用胶结充填体单轴抗压强度测试装置分别对多个抗压强度测试样品进行单轴抗压强度测试,并对测得的多个抗压强度测试样品的单轴抗压强度取均值,得到胶结充填体试样(19)的单轴抗压强度;Step 802: Use the uniaxial compressive strength test device for the cemented filling body to perform uniaxial compressive strength tests on multiple compressive strength test samples respectively, and take the measured uniaxial compressive strength of the multiple compressive strength test samples as The mean value is obtained to obtain the uniaxial compressive strength of the cemented filling body sample (19); 步骤803、获取Tensorflow深度学习力学响应特性预测网络的训练样本图像,具体过程为:Step 803, obtaining the training sample image of the Tensorflow deep learning mechanical response characteristic prediction network, the specific process is: 步骤8031、采用SEM扫描电镜分别对多个SEM扫描电镜样品进行多次扫描,形成多个SEM电镜扫描图像并存储到计算机(17)中;所述SEM电镜扫描图像的数量至少为150个;Step 8031, using the SEM scanning electron microscope to scan multiple SEM scanning electron microscope samples respectively, forming multiple SEM scanning electron microscope scanning images and storing them in the computer (17); the number of the SEM scanning electron microscope scanning images is at least 150; 步骤8032、所述计算机(17)调用高斯滤波处理模块分别对多个SEM电镜扫描图像进行高斯滤波处理,得到多个高斯滤波处理后的SEM电镜扫描图像;Step 8032, the computer (17) calls the Gaussian filter processing module to perform Gaussian filter processing on a plurality of SEM electron microscope scanning images respectively, and obtains a plurality of SEM electron microscope scanning images after Gaussian filter processing; 步骤8033、所述计算机(17)调用FCM模糊聚类处理模块分别对多个进行高斯滤波处理后的SEM电镜扫描图像进行孔隙图像提取,得到多组胶结充填体聚类图像,每组胶结充填体聚类图像中胶结充填体聚类图像的数量与聚类中心数目相等;Step 8033, the computer (17) invokes the FCM fuzzy clustering processing module to extract pore images from a plurality of Gaussian-filtered SEM electron microscope scanning images to obtain clustered images of multiple groups of cemented filling bodies, each group of cemented filling bodies The number of cemented filling cluster images in the cluster image is equal to the number of cluster centers; 步骤8034、所述计算机(17)将每组胶结充填体聚类图像中灰度值最小一类的胶结充填体聚类图像确定为胶结充填体微观孔隙图,并对多个胶结充填体微观孔隙图进行二值化处理,得到多个胶结充填体微观孔隙二值图;Step 8034, the computer (17) determines the cemented filling body cluster image with the smallest gray value in each group of cemented filling body cluster images as the cemented filling body microscopic pore map, and compares the microscopic pores of multiple cemented filling bodies Binary processing is performed on the map to obtain binary maps of the microscopic pores of multiple cemented filling bodies; 步骤8035、所述计算机(17)将步骤8032中得到的多个高斯滤波处理后的SEM电镜扫描图像与步骤8034中得到多个的胶结充填体微观孔隙二值图按照编号对应进行合并,得到多个训练样本图像;Step 8035, the computer (17) combines the multiple Gaussian-filtered SEM electron microscope scanning images obtained in step 8032 with the multiple binary images of cemented filling bodies obtained in step 8034 according to the corresponding numbers to obtain multiple a training sample image; 步骤804、所述计算机(17)分别对多个训练样本图像进行正规化处理,形成多个像素为960×960的正规化训练样本图像;Step 804, the computer (17) performs normalization processing on a plurality of training sample images respectively, forming a plurality of normalized training sample images whose pixels are 960×960; 步骤805、所述计算机(17)构建一个卷积网络核的层数为五层、输入层为正规化训练样本图像、输出层为正规化训练样本图像对应的单轴抗压强度的Tensorflow深度学习网络,将其存储的多个正规化训练样本图像作为训练样本,对Tensorflow深度学习网络进行训练,得到Tensorflow深度学习力学响应特性预测网络;所述Tensorflow深度学习力学特性响应预测网络五层卷积网络核的大小从一层到第五层分别为3x3,2x2,3x3,2x2,2x2。Step 805, the computer (17) constructs a convolutional network with five layers of layers, the input layer is a normalized training sample image, and the output layer is the Tensorflow deep learning of the uniaxial compressive strength corresponding to the normalized training sample image The network uses a plurality of normalized training sample images stored in it as training samples to train the Tensorflow deep learning network to obtain a Tensorflow deep learning mechanical response characteristic prediction network; the Tensorflow deep learning mechanical characteristic response prediction network is a five-layer convolutional network The kernel sizes are 3x3, 2x2, 3x3, 2x2, 2x2 from the first layer to the fifth layer. 7.按照权利要求6所述的基于SEM图像的胶结充填体力学响应特性预测方法,其特征在于:步骤802中所述胶结充填体单轴抗压强度测试装置包括座垫(10)和固定连接在座垫(10)顶部的多根拉杆(8),以及用于给胶结充填体试样(19)施加轴向压力的轴向加压传力机构和用于给轴向加压传力机构提供动力的轴向加压动力系统;所述座垫(10)的底部固定连接有多个底座(15),所述座垫(10)的顶部设置有用于放置胶结充填体试样(19)的试样放置槽,所述座垫(10)上位于所述试样放置槽的中心位置处设置有排水阀(5);多根拉杆(8)的中部设置有用于固定多根拉杆(8)的固定架(11),多根拉杆(8)的上部固定连接有顶部装载板(9);所述轴向加压传力机构包括安装在顶部装载板(9)上的气缸(2),所述气缸(2)的活塞杆向下设置,所述气缸(2)的活塞杆底部连接有压力传递板(3);所述轴向加压动力系统包括压缩空气气源(4)和加压控制器(18),以及一端与压缩空气气源(4)连接、另一端与气缸(2)连接的气体输送管(1);所述气体输送管(1)上从连接压缩空气气源(4)到连接气缸(2)的位置依次设置有气动三联件(12)、压力传感器(13)和气缸控制电磁阀(14),所述压力传感器(13)与加压控制器(18)的输入端连接,所述气缸控制电磁阀(14)与加压控制器(18)的输出端连接,所述加压控制器(18)通过通信模块(16)与计算机(17)连接。7. The method for predicting the mechanical response characteristics of cemented filling bodies based on SEM images according to claim 6, characterized in that: the uniaxial compressive strength test device for cemented filling bodies described in step 802 includes seat cushions (10) and fixed connections A plurality of pull rods (8) on the top of the seat cushion (10), and an axial pressure force transmission mechanism for applying axial pressure to the cemented filling body sample (19) and a force transmission mechanism for axial pressure An axial pressurization power system that provides power; the bottom of the seat cushion (10) is fixedly connected with a plurality of bases (15), and the top of the seat cushion (10) is provided with a cemented filling body sample (19) The sample placement groove of the seat cushion (10) is provided with a drain valve (5) at the center of the sample placement groove; the middle part of the plurality of pull rods (8) is provided with a plurality of pull rods (8) ) fixed frame (11), the top of a plurality of pull rods (8) is fixedly connected with a top loading plate (9); the axial pressurization force transmission mechanism includes a cylinder (2) mounted on the top loading plate (9) , the piston rod of the cylinder (2) is set downwards, and the bottom of the piston rod of the cylinder (2) is connected with a pressure transmission plate (3); the axial pressurized power system includes a compressed air source (4) and Pressurization controller (18), and the gas delivery pipe (1) that one end is connected with compressed air source (4), and the other end is connected with cylinder (2); The position connecting the source (4) to the cylinder (2) is provided with a pneumatic triple piece (12), a pressure sensor (13) and a cylinder control solenoid valve (14) in sequence, and the pressure sensor (13) is connected with the pressurization controller (18 ), the cylinder control solenoid valve (14) is connected with the output end of the pressurization controller (18), and the pressurization controller (18) is connected with the computer (17) through the communication module (16). 8.按照权利要求7所述的基于SEM图像的胶结充填体力学响应特性预测方法,其特征在于:所述试样放置槽内设置有套装在胶结充填体试样(19)底部的O型密封圈(6),所述座垫(10)上设置有位于所述试样放置槽周围的多孔石(7)。8. The method for predicting the mechanical response characteristics of cemented filling body based on SEM images according to claim 7, characterized in that: an O-shaped seal set at the bottom of the cemented filling body sample (19) is arranged in the said sample placement groove ring (6), the seat cushion (10) is provided with porous stones (7) located around the sample placement groove. 9.按照权利要求8所述的基于SEM图像的胶结充填体力学响应特性预测方法,其特征在于:所述加压控制器(18)为可编程逻辑控制器,所述通信模块(16)为RS-485通信模块。9. The method for predicting mechanical response characteristics of cemented filling bodies based on SEM images according to claim 8, characterized in that: the pressurization controller (18) is a programmable logic controller, and the communication module (16) is RS-485 communication module. 10.按照权利要求8所述的基于SEM图像的胶结充填体力学响应特性预测方法,其特征在于:步骤802中所述采用胶结充填体单轴抗压强度测试装置分别对多个抗压强度测试样品进行单轴抗压强度测试,其中对每个抗压强度测试样品进行单轴抗压强度测试的具体过程为:10. The method for predicting the mechanical response characteristics of cemented filling bodies based on SEM images according to claim 8, characterized in that: in step 802, the uniaxial compressive strength testing device for cemented filling bodies is used to test the compressive strength of multiple The sample is subjected to a uniaxial compressive strength test, wherein the specific process of performing a uniaxial compressive strength test on each compressive strength test sample is: 步骤8021、将O型密封圈(6)放入所述试样放置槽内后,将胶结充填体试样(19)放入所述试样放置槽内,使胶结充填体试样(19)的中心与气缸(2)的活塞杆和压力传递板(3)的中心相对应;并在座垫(10)上放入位于所述试样放置槽周围的多孔石(7);Step 8021, after putting the O-ring (6) into the sample placement groove, put the cemented filling body sample (19) into the sample placement groove, so that the cemented filling body sample (19) Corresponding to the center of the piston rod of the cylinder (2) and the center of the pressure transmission plate (3); and on the seat cushion (10), put the porous stone (7) that is positioned around the sample placement groove; 步骤8022、打开压缩空气气源(4),通过调节气动三联件(12)调节压缩空气气源(4)输出的压缩空气的气压,加压控制器(18)通过控制气缸控制电磁阀(14)换向,控制气缸(2)的活塞杆向下或向上运动,对胶结充填体试样(19)施加压力或卸载压力,将胶结充填体试样(19)破裂时加压控制器(18)采集到的压力传感器(13)检测的压力值记为F,加压控制器(18)将压力值F传输给计算机(17),计算机(17)根据公式计算得到抗压强度测试样品的单轴抗压强度P;其中,S为抗压强度测试样品的顶面面积;当气缸(2)的活塞杆向下运动时,带动压力传递板(3)向下运动,通过压力传递板(3)给胶结充填体试样(19)施加压力,当气缸(2)的活塞杆向上运动时,带动压力传递板(3)向上运动,压力传递板(3)离开胶结充填体试样(19)的上表面,卸载压力。Step 8022, open the compressed air source (4), adjust the air pressure of the compressed air output by the compressed air source (4) by adjusting the pneumatic triple unit (12), and the pressurization controller (18) controls the solenoid valve (14) by controlling the cylinder ) reversing, control the downward or upward movement of the piston rod of the cylinder (2), apply pressure to the cemented filling body sample (19) or unload the pressure, and pressurize the controller (18) when the cemented filling body sample (19) ruptures ) The pressure value detected by the pressure sensor (13) collected is denoted as F, and the pressurization controller (18) transmits the pressure value F to the computer (17), and the computer (17) according to the formula Calculate the uniaxial compressive strength P of the compressive strength test sample; wherein, S is the top surface area of the compressive strength test sample; when the piston rod of the cylinder (2) moves downward, it drives the pressure transmission plate (3) to The downward movement applies pressure to the cemented filling body sample (19) through the pressure transmission plate (3). When the piston rod of the cylinder (2) moves upward, it drives the pressure transmission plate (3) to move upward, and the pressure transmission plate (3) Leave the upper surface of the cemented filling body sample (19) and unload the pressure.
CN201810098527.0A 2018-01-31 2018-01-31 Prediction method of mechanical response characteristics of cemented backfill based on SEM images Active CN108256258B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810098527.0A CN108256258B (en) 2018-01-31 2018-01-31 Prediction method of mechanical response characteristics of cemented backfill based on SEM images

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810098527.0A CN108256258B (en) 2018-01-31 2018-01-31 Prediction method of mechanical response characteristics of cemented backfill based on SEM images

Publications (2)

Publication Number Publication Date
CN108256258A true CN108256258A (en) 2018-07-06
CN108256258B CN108256258B (en) 2019-01-25

Family

ID=62743195

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810098527.0A Active CN108256258B (en) 2018-01-31 2018-01-31 Prediction method of mechanical response characteristics of cemented backfill based on SEM images

Country Status (1)

Country Link
CN (1) CN108256258B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109858551A (en) * 2019-01-31 2019-06-07 湖南大学 Method, apparatus, equipment and medium based on image recognition detection structure stress
CN113743424A (en) * 2021-09-15 2021-12-03 昆明理工大学 DBSCAN cluster analysis method for extracting characteristic points of full tailings dynamic flocculation sedimentation image

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2000073972A1 (en) * 1999-05-28 2000-12-07 Texas Tech University Adaptive vector quantization/quantizer
CN102622729A (en) * 2012-03-08 2012-08-01 北京邮电大学 Spatial self-adaptive block-matching image denoising method based on fuzzy set theory
CN103632381A (en) * 2013-12-08 2014-03-12 中国科学院光电技术研究所 Extended target tracking method for extracting feature points by using skeleton
CN103839274A (en) * 2014-03-25 2014-06-04 中国科学院光电技术研究所 Extended target tracking method based on geometric proportion relation
CN104750940A (en) * 2015-04-10 2015-07-01 中国矿业大学(北京) Dynamic design method for strength of cemented filling body of underground mining stope
CN106503397A (en) * 2016-11-16 2017-03-15 中国航空工业集团公司北京航空材料研究院 A kind of method for obtaining metal material recrystallized structure evolution crystal boundary mobility parameter
CN106934109A (en) * 2017-02-20 2017-07-07 南京航空航天大学 A kind of Ni3Al based alloys represent the construction method of elementary volume, volume element model

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2000073972A1 (en) * 1999-05-28 2000-12-07 Texas Tech University Adaptive vector quantization/quantizer
CN102622729A (en) * 2012-03-08 2012-08-01 北京邮电大学 Spatial self-adaptive block-matching image denoising method based on fuzzy set theory
CN103632381A (en) * 2013-12-08 2014-03-12 中国科学院光电技术研究所 Extended target tracking method for extracting feature points by using skeleton
CN103839274A (en) * 2014-03-25 2014-06-04 中国科学院光电技术研究所 Extended target tracking method based on geometric proportion relation
CN104750940A (en) * 2015-04-10 2015-07-01 中国矿业大学(北京) Dynamic design method for strength of cemented filling body of underground mining stope
CN106503397A (en) * 2016-11-16 2017-03-15 中国航空工业集团公司北京航空材料研究院 A kind of method for obtaining metal material recrystallized structure evolution crystal boundary mobility parameter
CN106934109A (en) * 2017-02-20 2017-07-07 南京航空航天大学 A kind of Ni3Al based alloys represent the construction method of elementary volume, volume element model

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
OKAN ONAL等: "Artificial neural network application on microstructure–compressive strength relationship of cement mortar", 《ADVANCES IN ENGINEERING SOFLTWARE》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109858551A (en) * 2019-01-31 2019-06-07 湖南大学 Method, apparatus, equipment and medium based on image recognition detection structure stress
CN113743424A (en) * 2021-09-15 2021-12-03 昆明理工大学 DBSCAN cluster analysis method for extracting characteristic points of full tailings dynamic flocculation sedimentation image

Also Published As

Publication number Publication date
CN108256258B (en) 2019-01-25

Similar Documents

Publication Publication Date Title
CN108229062B (en) Method based on sensibility micro-parameter prediction cemented fill mechanical response characteristic
CN106960070B (en) A Seepage Simulation Method for Coal Mass Reconstructed Based on Finite Element-Discrete Element CT
CN101539491B (en) Device for testing model with three-dimensional gradient nonuniform loading structure
Zhu et al. Network modeling of the evolution of permeability and dilatancy in compact rock
CN108763711A (en) A kind of Permeability Prediction method based on Core Scanning Image piecemeal numerical simulation
CN110186789B (en) Construction waste roadbed permanent deformation orthogonal estimation method based on grey system
CN108827774A (en) Coal petrography brittleness evaluation method
CN108256258A (en) Cemented fill mechanical response characteristic Forecasting Methodology based on SEM image
CN108287254B (en) A kind of cemented fill micro-parameter Sensitivity Analysis based on image procossing
CN116401966A (en) A Simulation Method for Time-varying Slurry Diffusion Process in Porous Fractured Rock Mass
CN111445965A (en) Design method of carbon fiber reinforced cement-based material based on deep learning
CN105604341B (en) A kind of high sand content earthen ruins reinforcement means
CN110006764A (en) Geotechnical and underground engineering model test traffic dynamic load simulation device and test method
Qi et al. Study on the Mechanical Properties of Red Clay under Drying‐Wetting Cycles
CN108445030A (en) Cemented fill mechanical response characteristic prediction technique based on image and micro-parameter
CN117333670A (en) Recycled concrete hydration process detection method based on Transformer and microscopic images
Maragh et al. SEM-EDS and microindentation-driven large-area high-resolution chemomechanical mapping and computational homogenization of cementitious materials
Li et al. Bond strength between recycled aggregate concrete and rebar: Interpretable machine learning modeling approach for performance estimation and engineering design
CN108240940A (en) A test system and method for mechanical response characteristics of cemented filling body
CN109033503B (en) Numerical simulation method for uniaxial compressive strength test of porous cement concrete
CN102799860A (en) Method for holographic recognition of microscopic image
CN109020338A (en) A kind of cement stabilized type iron tailings sand base material design method
CN209878515U (en) Traffic dynamic load simulation device for geotechnical and underground engineering model tests
CN204255818U (en) For the pilot system of the information acquisition of concrete surface Initial Cracking and identification
CN113343474B (en) Microscopic analysis method of broken stone reinforced polymer composite material

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