CN104036538B - Soil-rock mixture three-dimensional microscopical structure is rebuild and analysis method and system - Google Patents

Soil-rock mixture three-dimensional microscopical structure is rebuild and analysis method and system Download PDF

Info

Publication number
CN104036538B
CN104036538B CN201410293183.0A CN201410293183A CN104036538B CN 104036538 B CN104036538 B CN 104036538B CN 201410293183 A CN201410293183 A CN 201410293183A CN 104036538 B CN104036538 B CN 104036538B
Authority
CN
China
Prior art keywords
soil
block
stone
block stone
dimensional
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.)
Expired - Fee Related
Application number
CN201410293183.0A
Other languages
Chinese (zh)
Other versions
CN104036538A (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.)
Tsinghua University
Original Assignee
Tsinghua University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tsinghua University filed Critical Tsinghua University
Priority to CN201410293183.0A priority Critical patent/CN104036538B/en
Publication of CN104036538A publication Critical patent/CN104036538A/en
Application granted granted Critical
Publication of CN104036538B publication Critical patent/CN104036538B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

本发明提出一种土石混合体三维细观结构重建与分析方法及系统,该方法包括:制备土石混合体试样;对土石混合体试样进行CT扫描以得到土石混合体试样的CT断层图像序列;获取CT断层图像序列中每一张CT断层图像中的块石断面区域,并对块石断面区域进行分割以得到多个独立的块石断面;对每一张CT断层图像中的多个独立的块石断面进行配准,并根据配准结果生成块石的三维细观结构模型;对块石的三维细观结构模型进行分析以确定土石混合体试样的含石量、粒度分布、块石空间分布、块石几何特性和块石形态特性。根据本发明的实施例可方便地生成三维细观结构模型,并具有土石混合体分析结果精确、可靠的优点,同时也可直接用于数值计算分析。

The present invention proposes a method and system for reconstruction and analysis of a three-dimensional mesoscopic structure of a soil-rock mixture, the method comprising: preparing a soil-rock mixture sample; performing CT scanning on the soil-rock mixture sample to obtain a CT tomographic image of the soil-rock mixture sample Sequence; Obtain the block section area in each CT tomographic image in the CT tomographic image sequence, and segment the block section area to obtain multiple independent block stone sections; The independent block section is registered, and the three-dimensional mesoscopic structure model of the block is generated according to the registration result; the three-dimensional mesoscopic structure model of the block is analyzed to determine the stone content, particle size distribution, Block spatial distribution, block geometry and block shape characteristics. According to the embodiment of the present invention, the three-dimensional mesoscopic structural model can be conveniently generated, and has the advantages of accurate and reliable analysis results of the soil-rock mixture, and can also be directly used for numerical calculation and analysis.

Description

土石混合体三维细观结构重建与分析方法及系统Method and system for reconstruction and analysis of 3D mesoscopic structure of soil-rock mixture

技术领域technical field

本发明涉及岩土介质分析技术领域,尤其涉及一种土石混合体三维细观结构重建与分析方法及系统。The invention relates to the technical field of rock-soil medium analysis, in particular to a method and system for reconstructing and analyzing a three-dimensional mesoscopic structure of an earth-rock mixture.

背景技术Background technique

土石混合体是自然界中广泛存在的一类复杂的岩土介质,也是地质工程、岩土工程中经常遇到的一类地质体。由于土石混合体内部超径块石的存在,使得其物理力学性质很难通过传统的室内或野外试验进行获取。土石混合体物理力学性质的研究一直以来是岩土工程界面临的难题。Soil-rock mixture is a kind of complex rock-soil medium that widely exists in nature, and it is also a kind of geological body that is often encountered in geological engineering and geotechnical engineering. Due to the existence of super-diameter boulders inside the soil-rock mixture, its physical and mechanical properties are difficult to obtain through traditional indoor or field tests. The research on the physical and mechanical properties of soil-rock mixtures has always been a difficult problem in the field of geotechnical engineering.

随着数值计算技术的发展,数值试验为土石混合体物理力学性质的研究并揭示这类岩土体的变形破裂机理提供了新的技术支持。在开展土石混合体的数值试验研究时,建立与物理力学试验一致的计算分析模型,是确保数值计算可靠性的重要基础。然而,目前在土石混合体的三维细观结构模型生成方面,基本是采用将块石简化为球体、凸多面体等简单的几何形体,然后采用计算机随机生成的方式,产生相应的细观结构模型。这种方法难于保证所建立的模型与物理试验模型的一致性,从而使得计算分析结果难以与物理试验很好的对比,参考价值不高。With the development of numerical calculation technology, numerical experiments provide new technical support for the study of physical and mechanical properties of soil-rock mixtures and reveal the deformation and fracture mechanism of such rock-soil masses. When carrying out numerical test research on soil-rock mixture, establishing a calculation and analysis model consistent with physical and mechanical tests is an important basis for ensuring the reliability of numerical calculations. However, at present, in terms of generating the three-dimensional mesostructure model of soil-rock mixture, the basic method is to simplify the block rocks into simple geometric shapes such as spheres and convex polyhedrons, and then use the method of computer random generation to generate the corresponding mesostructure model. This method is difficult to ensure the consistency between the established model and the physical test model, which makes it difficult to compare the calculation and analysis results with the physical test, and the reference value is not high.

发明内容Contents of the invention

本发明旨在至少在一定程度上解决相关技术中的技术问题之一。The present invention aims to solve one of the technical problems in the related art at least to a certain extent.

为此,本发明的第一个目的在于提出一种土石混合体三维细观结构重建与分析方法。该方法可方便地生成土石混合体的三维细观结构模型,并可直接用于土石混合体的细观结构分析,具有分析结果精确可靠的优点。Therefore, the first object of the present invention is to propose a method for reconstruction and analysis of the three-dimensional mesoscopic structure of soil-rock mixture. This method can conveniently generate a three-dimensional mesostructure model of soil-rock mixture, and can be directly used in the analysis of mesostructure of soil-rock mixture, which has the advantage of accurate and reliable analysis results.

本发明的第二个目的在于提出一种土石混合体三维细观结构重建与分析系统。The second object of the present invention is to propose a three-dimensional mesoscopic structure reconstruction and analysis system of soil-rock mixture.

为了实现上述目的,本发明的第一方面的实施例公开了一种土石混合体三维细观结构重建与分析方法,包括以下步骤:制备土石混合体试样;对所述土石混合体试样进行CT扫描以得到所述土石混合体试样的CT断层图像序列;获取所述CT断层图像序列中每一张CT断层图像中的块石断面区域,并对所述块石断面区域进行分割以得到多个独立的块石断面;对所述每一张CT断层图像中的多个独立的块石断面进行配准,并根据配准结果生成块石的三维细观结构模型;以及对所述块石的三维细观结构模型进行分析以确定所述土石混合体试样的含石量、粒度分布、块石空间分布、块石几何特性和块石形态特性。In order to achieve the above object, the embodiment of the first aspect of the present invention discloses a three-dimensional mesostructure reconstruction and analysis method of soil-rock mixture, comprising the following steps: preparing a sample of soil-rock mixture; CT scan to obtain the CT tomographic image sequence of the soil-rock mixture sample; obtain the block section area in each CT tomographic image in the CT tomographic image sequence, and segment the block section area to obtain A plurality of independent block stone sections; registering the multiple independent block stone sections in each CT tomographic image, and generating a three-dimensional mesoscopic structure model of the block stone according to the registration result; and The three-dimensional mesoscopic structure model of rocks is analyzed to determine the rock content, particle size distribution, spatial distribution of rocks, geometrical characteristics of rocks and morphological characteristics of rocks in the soil-rock mixture sample.

另外,根据本发明上述实施例的土石混合体三维细观结构重建与分析方法还可以具有如下附加的技术特征:In addition, the method for reconstructing and analyzing the three-dimensional mesoscopic structure of soil-rock mixture according to the above-mentioned embodiments of the present invention may also have the following additional technical features:

在一些示例中,通过有机玻璃制样筒制备所述土石混合体试样,通过CT扫描机对所述土石混合体试样进行扫描。In some examples, the soil-rock mixture sample is prepared by a plexiglass sample cylinder, and the soil-rock mixture sample is scanned by a CT scanner.

在一些示例中,所述获取CT断层图像序列中每一张CT断层图像中的块石断面区域,并对所述块石断面区域进行分割以得到多个独立的块石断面,进一步包括:对所述CT断层图像序列中每一张CT断层图像进行滤波去噪;对滤波去噪后的每一张CT断层图像进行二值化处理以确定所述块石断面区域和土体区域;通过腐蚀和膨胀的图像形态学运算对所述块石断面区域进行分割以得到多个独立的块石断面。In some examples, the acquiring the block section area in each CT tomographic image in the CT tomographic image sequence, and segmenting the block section area to obtain multiple independent block section sections further includes: Perform filtering and denoising on each CT tomographic image in the CT tomographic image sequence; perform binarization processing on each CT tomographic image after filtering and denoising to determine the block section area and soil area; and the image morphology operation of expansion to segment the section area of the block to obtain multiple independent sections of the block.

在一些示例中,所述对每一张CT断层图像中的多个独立的块石断面进行配准,并根据配准结果生成块石的三维细观结构模型,进一步包括:依次对每一张CT断层图像上的所有块石断面的边界进行识别与编号,并计算每个块石断面的几何特征指标;对所有CT断层图像上的块石断面的边界进行配准以得到所述配准结果;根据所述块石断面的几何特征指标和所述配准结果生成所述块石的三维细观结构模型。In some examples, the registration of multiple independent block sections in each CT tomographic image, and generating a three-dimensional mesoscopic structure model of the block according to the registration results further includes: Identify and number the boundaries of all block sections on the CT tomogram, and calculate the geometric feature index of each block section; register the boundaries of all block sections on the CT tomogram to obtain the registration result ; Generate a three-dimensional mesoscopic structure model of the block according to the geometric feature index of the block section and the registration result.

在一些示例中,所述对块石的三维细观结构模型进行分析以确定土石混合体试样的含石量、粒度分布、块石空间分布、块石几何特性和块石形态特性,进一步包括:对块石的三维细观结构模型进行渲染和展示,以根据渲染和展示结果确定所述土石混合体试样的含石量、粒度分布、块石空间分布、块石几何特性和块石形态特性。In some examples, analyzing the three-dimensional mesoscopic structure model of the block to determine the rock content, particle size distribution, spatial distribution of the block, geometric properties of the block, and morphological properties of the block, further comprising : Render and display the three-dimensional mesoscopic structure model of the block, so as to determine the rock content, particle size distribution, spatial distribution of the block, block geometry and block shape of the soil-rock mixture sample according to the rendering and display results characteristic.

本发明的第二方面的实施例公开了一种土石混合体三维细观结构重建与分析系统,包括:土石混合体试样制备装置,用于制备土石混合体试样;CT断层图像扫描装置,用于对所述土石混合体试样进行CT扫描以得到所述土石混合体试样的CT断层图像序列;三维模型重建模块,用于获取所述CT断层图像序列中每一张CT断层图像中的块石断面区域,并对所述块石断面区域进行分割以得到多个独立的块石断面,以及对所述每一张CT断层图像中的多个独立的块石断面进行配准,并根据配准结果生成块石的三维细观结构模型;以及三维细观结构信息统计分析模块,用于对所述块石的三维细观结构模型进行分析以确定所述土石混合体试样的含石量、粒度分布、块石空间分布、块石几何特性和块石形态特性。The embodiment of the second aspect of the present invention discloses a three-dimensional mesostructure reconstruction and analysis system of soil-rock mixture, including: a soil-rock mixture sample preparation device for preparing soil-rock mixture samples; a CT tomographic image scanning device, It is used to perform CT scanning on the soil-rock mixture sample to obtain a CT tomographic image sequence of the soil-rock mixture sample; a three-dimensional model reconstruction module is used to obtain each CT tomographic image in the CT tomographic image sequence block section area, and segment the block section area to obtain multiple independent block section sections, and register the multiple independent block section sections in each CT tomographic image, and A three-dimensional mesoscopic structure model of the block is generated according to the registration result; and a statistical analysis module of the three-dimensional mesostructure information is used to analyze the three-dimensional mesostructure model of the block to determine the content of the soil-rock mixture sample Stone quantity, particle size distribution, block spatial distribution, block geometric characteristics and block shape characteristics.

另外,根据本发明上述实施例的土石混合体三维细观结构重建与分析系统还可以具有如下附加的技术特征:In addition, the three-dimensional mesoscopic structure reconstruction and analysis system of soil-rock mixture according to the above-mentioned embodiments of the present invention may also have the following additional technical features:

在一些示例中,所述土石混合体试样制备装置为有机玻璃制样筒,所述CT断层图像扫描装置为CT扫描机。In some examples, the soil-rock mixture sample preparation device is a plexiglass sample cylinder, and the CT tomographic image scanning device is a CT scanner.

在一些示例中,所述三维模型重建模块用于:对所述CT断层图像序列中每一张CT断层图像进行滤波去噪;对滤波去噪后的每一张CT断层图像进行二值化处理以确定所述块石断面区域和土体区域;通过腐蚀和膨胀的图像形态学运算对所述块石断面区域进行分割以得到多个独立的块石断面。In some examples, the three-dimensional model reconstruction module is configured to: perform filtering and denoising on each CT tomographic image in the CT tomographic image sequence; perform binarization processing on each CT tomographic image after filtering and denoising To determine the block section area and soil area; segment the block section area through image morphology operations of erosion and expansion to obtain multiple independent block section areas.

在一些示例中,所述三维模型重建模块还用于:依次对每一张CT断层图像上的所有块石断面的边界进行识别与编号,并计算每个块石断面的几何特征指标;对所有CT断层图像上的块石断面的边界进行配准以得到所述配准结果;根据所述块石断面的几何特征指标和所述配准结果生成所述块石的三维细观结构模型。In some examples, the three-dimensional model reconstruction module is also used to: sequentially identify and number the boundaries of all block sections on each CT tomographic image, and calculate the geometric feature index of each block section; The boundary of the block section on the CT tomographic image is registered to obtain the registration result; a three-dimensional mesoscopic structure model of the block is generated according to the geometric feature index of the block section and the registration result.

在一些示例中,还包括:三维显示计算机装置,用于对块石的三维细观结构模型进行渲染和展示。In some examples, it also includes: a three-dimensional display computer device for rendering and displaying the three-dimensional mesoscopic structure model of the block stone.

根据本发明的实施例,能够根据土石混合体试样方便地建立其对应的三维结构模型,并对试样内部块石的粒度组成、形态等特性进行定量分析,同时实现模型的三维可视化展示与存储,为土石混合体的物理力学数值试验提供支持。According to the embodiment of the present invention, the corresponding three-dimensional structural model can be conveniently established according to the soil-rock mixture sample, and the characteristics such as the particle size composition and shape of the rock inside the sample can be quantitatively analyzed, and the three-dimensional visual display and display of the model can be realized at the same time. Storage, providing support for physical and mechanical numerical tests of soil-rock mixtures.

相较于现有的方式,本发明的实施例能在对土石混合体三维细观结构建模时,利用土石混合体CT断层序列图像,在不破坏试样的情况下实现土石混合体的三维细观结构模型重建,扫描后的试样可继续进行室内物理力学试验研究。根据重建后的土石混合体试样模型,可对其内部块石几何形态特征、粒度组成、空间分布等信息进行统计分析,并可对模型进行渲染、剖切等查询操作。同时,重建后的模型可直接存储为通用数据格式,以用于数值试验研究。本发明的实施例,大大提高了所生成模型的精细化程度,在数值试验与物理力学试验之间建立了桥梁纽带,从而使得在此基础上开展的物理力学性质数值试验研究更加合理与可靠。Compared with the existing methods, the embodiment of the present invention can realize the three-dimensional structure of the soil-rock mixture without destroying the sample by using CT tomographic images of the soil-rock mixture when modeling the three-dimensional mesoscopic structure of the soil-rock mixture. The mesostructure model is reconstructed, and the scanned samples can continue to be studied in the physical and mechanical experiments in the laboratory. According to the reconstructed soil-rock mixture sample model, statistical analysis can be carried out on the geometric characteristics, particle size composition, spatial distribution and other information of the internal block stones, and query operations such as rendering and sectioning can be performed on the model. At the same time, the reconstructed model can be directly stored in a common data format for numerical experiment research. The embodiment of the present invention greatly improves the degree of refinement of the generated model, and establishes a bridge link between the numerical test and the physical-mechanical test, so that the numerical test research of physical-mechanical properties carried out on this basis is more reasonable and reliable.

本发明附加的方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。Additional aspects and advantages of the invention will be set forth in part in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention.

附图说明Description of drawings

图1是本发明一个实施例中的有机玻璃制样筒的构造图;Fig. 1 is the structural diagram of the plexiglass sample cylinder in one embodiment of the present invention;

图2是本发明一个实施例中的CT扫描机的示意图;Figure 2 is a schematic diagram of a CT scanner in one embodiment of the present invention;

图3是根据本发明一个实施例的土石混合体三维细观结构重建与分析方法的块石分割的流程图;Fig. 3 is a flow chart of the block-stone segmentation of the three-dimensional mesoscopic structure reconstruction and analysis method of soil-rock mixture according to an embodiment of the present invention;

图4是根据本发明一个实施例的土石混合体三维细观结构重建与分析方法的生成块石的三维细观结构模型的流程图;Fig. 4 is a flow chart of generating a three-dimensional mesoscopic structure model of a rock in a method for reconstructing and analyzing a three-dimensional mesoscopic structure of a soil-rock mixture according to an embodiment of the present invention;

图5是本发明一个实施例中的复杂岩土介质细观结构三维重建系统的结构图;Fig. 5 is a structural diagram of a three-dimensional reconstruction system of complex rock-soil medium mesostructure in an embodiment of the present invention;

图6是根据本发明一个实施例的土石混合体三维细观结构重建与分析方法的流程图;以及Fig. 6 is a flow chart of a method for reconstructing and analyzing a three-dimensional mesoscopic structure of a soil-rock mixture according to an embodiment of the present invention; and

图7是根据本发明一个实施例的土石混合体三维细观结构重建与分析系统的结构图。Fig. 7 is a structural diagram of a three-dimensional mesostructure reconstruction and analysis system for soil-rock mixtures according to an embodiment of the present invention.

具体实施方式detailed description

下面详细描述本发明的实施例,所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。Embodiments of the present invention are described in detail below, examples of which are shown in the drawings, wherein the same or similar reference numerals designate the same or similar elements or elements having the same or similar functions throughout. The embodiments described below by referring to the figures are exemplary and are intended to explain the present invention and should not be construed as limiting the present invention.

下文的公开提供了许多不同的实施例或例子用来实现本发明的不同结构。为了简化本发明的公开,下文中对特定例子的部件和设置进行描述。当然,它们仅仅为示例,并且目的不在于限制本发明。此外,本发明可以在不同例子中重复参考数字和/或字母。这种重复是为了简化和清楚的目的,其本身不指示所讨论各种实施例和/或设置之间的关系。此外,本发明提供了各种特定的工艺和材料的例子,但是本领域普通技术人员可以意识到其他工艺的可应用于性和/或其他材料的使用。另外,以下描述的第一特征在第二特征之“上”的结构可以包括第一和第二特征形成为直接接触的实施例,也可以包括另外的特征形成在第一和第二特征之间的实施例,这样第一和第二特征可能不是直接接触。The following disclosure provides many different embodiments or examples for implementing different structures of the present invention. To simplify the disclosure of the present invention, components and arrangements of specific examples are described below. Of course, they are only examples and are not intended to limit the invention. Furthermore, the present invention may repeat reference numerals and/or letters in different instances. This repetition is for the purpose of simplicity and clarity and does not in itself indicate a relationship between the various embodiments and/or arrangements discussed. In addition, various specific examples of processes and materials are provided herein, but one of ordinary skill in the art will recognize the applicability of other processes and/or the use of other materials. Additionally, configurations described below in which a first feature is "on" a second feature may include embodiments where the first and second features are formed in direct contact, and may include additional features formed between the first and second features. For example, such that the first and second features may not be in direct contact.

在本发明的描述中,需要说明的是,除非另有规定和限定,术语“安装”、“相连”、“连接”应做广义理解,例如,可以是机械连接或电连接,也可以是两个元件内部的连通,可以是直接相连,也可以通过中间媒介间接相连,对于本领域的普通技术人员而言,可以根据具体情况理解上述术语的具体含义。In the description of the present invention, it should be noted that unless otherwise specified and limited, the terms "installation", "connection" and "connection" should be understood in a broad sense, for example, it can be a mechanical connection or an electrical connection, or it can be two The internal communication of each element may be directly connected or indirectly connected through an intermediary. Those skilled in the art can understand the specific meanings of the above terms according to specific situations.

以下结合附图描述根据本发明实施例的土石混合体三维细观结构重建与分析方法及系统。The method and system for three-dimensional mesoscopic structure reconstruction and analysis of soil-rock mixture according to the embodiment of the present invention will be described below with reference to the accompanying drawings.

该方法应用于土石混合体的三维细观结构建模与分析,包括:This method is applied to the modeling and analysis of 3D mesostructure of soil-rock mixture, including:

土石混合体试样CT三维扫描步骤,采用CT扫描机对土石混合体试样进行扫描,获取土石混合体试样的CT断层序列图像(即CT断层图像序列);CT断层图像预处理步骤,采用数字图像处理技术对获得的土石混合体CT断面层图像进行处理,采用阈值分割技术将块石与周围的土体区分开,采用边缘检测技术得到CT断层图像内每个块石断面的边界点信息,并对其进行分块标注;土石混合体内部块石配准步骤,将获得的相邻CT断层图像切片内的分块标注的每个块石断面边界信息,根据其形态特性及空间位置进行匹配识别,得到构成每个块石的所有CT切片上的块石断面边界;块石三维模型重建步骤,根据匹配后的块石断面信息,构建表征每个块石表面的三角网格,建立土石混合体试样内部所有块石的三维模型;块石几何形态统计分析步骤,根据重建后的土石混合体内部块石模型,对块石的几何形态特征、粒度组成、空间分布等信息进行统计分析;土石混合体三维细观结构可视化及输出步骤,将重建的三维细观结构模型根据块石的粒径大小对块石进行着色渲染,允许用户进行旋转、切片等查询操作,同时允许用户将模型输出为通用的数据格式用于数值计算。The CT three-dimensional scanning step of the soil-rock mixture sample uses a CT scanner to scan the soil-rock mixture sample to obtain a CT tomographic sequence image (ie, a CT tomographic image sequence) of the soil-rock mixture sample; the CT tomographic image preprocessing step adopts Digital image processing technology processes the obtained CT cross-sectional image of the soil-rock mixture, uses threshold segmentation technology to distinguish the rock from the surrounding soil, and uses edge detection technology to obtain the boundary point information of each rock section in the CT tomographic image , and mark it in blocks; the internal block rock registration step of the soil-rock mixture is to carry out the boundary information of each block block marked in blocks in the obtained adjacent CT tomographic slices according to its morphological characteristics and spatial position. Matching and identification to obtain the boundary of the block section on all CT slices that constitute each block; the step of reconstruction of the 3D model of the block is to construct a triangular mesh representing the surface of each block according to the information of the section of the block after matching, and establish a soil-rock The three-dimensional model of all the blocks inside the mixture sample; the step of statistical analysis of block geometry, based on the reconstructed internal block model of the soil-rock mixture, statistical analysis of the geometric characteristics, particle size composition, spatial distribution and other information of the blocks ; The 3D mesoscopic structure visualization and output steps of the soil-rock mixture, the reconstructed 3D mesoscopic structure model is used to color and render the rocks according to the particle size of the rocks, allowing users to perform query operations such as rotation and slicing. The output is a common data format for numerical calculations.

具体而言,图6是根据本发明一个实施例的土石混合体三维细观结构重建与分析方法的流程图。Specifically, FIG. 6 is a flowchart of a method for reconstructing and analyzing a three-dimensional mesostructure of a soil-rock mixture according to an embodiment of the present invention.

如图6所示,根据本发明一个实施例的土石混合体三维细观结构重建与分析方法,包括如下步骤:As shown in Fig. 6, the reconstruction and analysis method for the three-dimensional mesoscopic structure of soil-rock mixture according to an embodiment of the present invention includes the following steps:

步骤S601:制备土石混合体试样。例如可通过有机玻璃制样筒制备土石混合体试样。Step S601: Prepare soil-rock mixture samples. For example, the soil-rock mixture sample can be prepared through a plexiglass sample cylinder.

步骤S602:对土石混合体试样进行CT扫描以得到土石混合体试样的CT断层图像序列。例如通过但不限于CT扫描机对土石混合体试样进行扫描。Step S602: Perform CT scanning on the soil-rock mixture sample to obtain a CT tomographic image sequence of the soil-rock mixture sample. For example, a soil-rock mixture sample is scanned by but not limited to a CT scanner.

步骤S603:获取CT断层图像序列中每一张CT断层图像中的块石断面区域,并对块石断面区域进行分割以得到多个独立的块石断面。Step S603: Obtain the section area of the block in each CT tomographic image in the CT tomographic image sequence, and segment the section area of the block to obtain multiple independent sections of the block.

具体地说,首先对CT断层图像序列中每一张CT断层图像进行滤波去噪,然后对滤波去噪后的每一张CT断层图像进行二值化处理以确定所述块石断面区域和土体区域,最后可通过腐蚀和膨胀的图像形态学运算对所述块石断面区域进行分割以得到多个独立的块石断面。Specifically, firstly, filter and denoise each CT tomographic image in the CT tomographic image sequence, and then perform binarization processing on each CT tomographic image after filtering and denoising to determine the section area of the rock and soil Finally, the block section area can be segmented by image morphology operations of erosion and expansion to obtain multiple independent block sections.

步骤S604:对每一张CT断层图像中的多个独立的块石断面进行配准,并根据配准结果生成块石的三维细观结构模型。Step S604: Register multiple independent block sections in each CT tomographic image, and generate a three-dimensional mesoscopic structural model of the block according to the registration results.

作为一个具体的示例,例如:依次对每一张CT断层图像上的所有块石断面的边界进行识别与编号,并计算每个块石断面的几何特征指标,并对所有CT断层图像上的块石断面的边界进行配准以得到所述配准结果,最后可根据块石断面的几何特征指标和配准结果生成块石的三维细观结构模型。As a specific example, for example: sequentially identify and number the boundaries of all block sections on each CT tomographic image, and calculate the geometric feature index of each block section, and calculate the block on all CT tomographic images The boundary of the stone section is registered to obtain the registration result, and finally a three-dimensional mesoscopic structural model of the block can be generated according to the geometric feature index of the block section and the registration result.

步骤S605:对块石的三维细观结构模型进行分析以确定土石混合体试样的含石量、粒度分布、块石空间分布、块石几何特性和块石形态特性。Step S605: Analyzing the three-dimensional mesoscopic structure model of the rock to determine the rock content, particle size distribution, spatial distribution of the rock, geometric characteristics of the rock, and morphological characteristics of the rock.

具体地说,首先可对块石的三维细观结构模型进行渲染和展示,这样,可根据渲染和展示结果确定土石混合体试样的含石量、粒度分布、块石空间分布、块石几何特性和块石形态特性等。Specifically, firstly, the three-dimensional mesoscopic structural model of the block can be rendered and displayed. In this way, the rock content, particle size distribution, spatial distribution of the block, and block geometry of the soil-rock mixture sample can be determined according to the rendering and display results. characteristics and morphological characteristics of stones, etc.

下面结合附图1-5对本发明实施例的方法进行更为详细的说明。The method of the embodiment of the present invention will be described in more detail below with reference to the accompanying drawings 1-5.

如图1所示,有机玻璃制样筒40(以下简称制样筒)主要包括试样帽1、对开侧筒2、底盘3。沿对开侧筒2高度方向设置有环形加强肋4,对开侧筒2两翼上的螺母孔5用于制样时对开侧筒2的固定以及完成扫描后开始试验前对开侧筒2的拆除。试样帽1上的阀门6通过试样帽1内的通气孔与试样连通,用于控制制样的过程中的抽真空操作。凹槽7用于放置并固定试验中的加载装置。试样帽1与底盘3上设置有橡皮筋凹槽8,用于制样的过程中橡皮膜的捆扎与固定。对开侧筒2底部与底盘3上的螺母孔9用于将两者固定在一起并共同承载试样重量。As shown in FIG. 1 , the plexiglass sample preparation cylinder 40 (hereinafter referred to as the sample preparation cylinder) mainly includes a sample cap 1 , a split side cylinder 2 , and a chassis 3 . Annular reinforcing ribs 4 are arranged along the height direction of the split side cylinder 2, and the nut holes 5 on both wings of the split side cylinder 2 are used for fixing the split side cylinder 2 during sample preparation and for split side cylinder 2 before starting the test after scanning demolition. The valve 6 on the sample cap 1 communicates with the sample through the vent hole in the sample cap 1, and is used to control the vacuuming operation during the sample preparation process. Groove 7 is used to place and fix the loading device in the test. The sample cap 1 and the chassis 3 are provided with a rubber band groove 8, which is used for binding and fixing the rubber film during the sample preparation process. The nut holes 9 on the bottom of the split side cylinder 2 and the chassis 3 are used to fix the two together and jointly carry the weight of the sample.

在本实施中,安装固定好底盘3、橡皮膜、对开侧筒2后,将土石混合体分层装入制样筒40中,并分层击实至所需密实度,待试样(即土石混合体试样)装满后安装试样帽1,抽真空后连同有机玻璃制样筒40一起进行试样的CT扫描,扫描完成后拆掉对开侧筒2将试样安装到试验仪器上开展试验。In this implementation, after installing and fixing the chassis 3, the rubber film, and the split side cylinder 2, the soil-rock mixture is loaded into the sample preparation cylinder 40 in layers, and compacted in layers to the required compactness, and the sample ( That is, the soil-rock mixture sample) is filled and the sample cap 1 is installed. After vacuuming, the CT scan of the sample is carried out together with the plexiglass sample preparation cylinder 40. After the scanning is completed, the split side cylinder 2 is removed to install the sample in the test Test on the instrument.

如图2所示,CT扫描机主要包括扫描架10、扫描床11、计算机系统12、CT机控制操作盘13。其中计算机系统12和CT机控制操作盘13是执行扫描操作及数据处理存储的控制中心,扫描架10装有X线球管、滤线器、准直器、探测器等,用于X射线的发射与接收,扫描床11是运载扫描试样的工具,能够根据控制要求垂直和水平纵向运动。As shown in FIG. 2 , the CT scanner mainly includes a scanning frame 10 , a scanning bed 11 , a computer system 12 , and a CT machine control panel 13 . Wherein the computer system 12 and the CT machine control panel 13 are the control centers for performing scanning operations and data processing and storage, and the scanning frame 10 is equipped with X-ray tubes, line filters, collimators, detectors, etc. Transmitting and receiving, the scanning bed 11 is a tool for carrying the scanning sample, and can move vertically and horizontally according to control requirements.

在本实施中,将土石混合体试样14连同有机玻璃制样筒40平放于扫描床11之上,试样14轴向与扫描床11的水平纵向运动方向一致,调整扫描床11的高度,调整土石混合体试样横截面中心与扫描架10发射的十字对中线中心重合后,将扫描床11向扫描架10水平推送,使试样14前端面恰好位于扫描架球管所在平面内,位置归零以标记扫描起点。In this implementation, the soil-rock mixture sample 14 is placed flat on the scanning bed 11 together with the plexiglass sample cylinder 40, the axial direction of the sample 14 is consistent with the horizontal and vertical movement direction of the scanning bed 11, and the height of the scanning bed 11 is adjusted. After adjusting the center of the cross-section of the soil-rock mixture sample to coincide with the center of the cross centerline emitted by the scanning frame 10, push the scanning bed 11 horizontally to the scanning frame 10, so that the front end of the sample 14 is just located in the plane where the scanning frame ball tube is located. The position is zeroed to mark the start of the scan.

通过计算机系统12建立该次扫描记录,设置扫描对象类型、球管X线的发射能量、CT断层间距等参数并选定扫描范围,之后通过CT机控制操作盘13启动CT机,完成扫描,获得的CT切片序列图像(即CT断层图像序列)可保存至外部存储设备中。The scanning record is established by the computer system 12, parameters such as the type of the scanning object, the emission energy of the tube X-ray, and the distance between CT slices are set, and the scanning range is selected, and then the CT machine is started through the CT machine control panel 13 to complete the scan and obtain The CT slice sequence images (ie, CT tomographic image sequence) can be saved to an external storage device.

如图3所示,块石分割的过程如下所示:As shown in Figure 3, the process of block stone segmentation is as follows:

步骤S10,导入上述CT扫描得到的土石混合体试样CT切片序列图像的数据,并将其存储到数据库中。Step S10, importing the data of the CT slice sequence image of the soil-rock mixture sample obtained by the above-mentioned CT scanning, and storing it in the database.

步骤S11,对扫描得到的二维CT切片序列图像进行批量化的滤波去噪处理,调整滤波器模板至合适尺寸,通过均值滤波或中值滤波等方法将二维切片图像中的部分噪点去除,使图像平滑化。Step S11, perform batch filtering and denoising processing on the scanned two-dimensional CT slice sequence images, adjust the filter template to an appropriate size, and remove part of the noise in the two-dimensional slice images by mean filtering or median filtering, etc. Smoothes the image.

步骤S12,分析经初步滤波处理后二维切片图像的灰度直方图,确定区分“土体”与“块石”的灰度阈值,据此对二维CT切片序列图像按照需求进行批量二值化处理(块石为研究对象,土体及孔隙为背景),当然,本发明的实施例并不限于此,例如也可进行三值化处理(即:块石与孔隙两种研究对象,土体作为背景)。Step S12, analyze the grayscale histogram of the two-dimensional slice image after the preliminary filtering process, determine the grayscale threshold value for distinguishing "soil" and "block stone", and perform batch binary value on the two-dimensional CT slice sequence image according to the requirements (the rock is the research object, and the soil and the pores are the background), of course, the embodiments of the present invention are not limited thereto, for example, the ternary processing (that is: the rock and the hole are two research objects, and the soil body as the background).

步骤S13,通过腐蚀、膨胀等简单图像形态学运算的组合使用对二值化的二维CT切片序列图像进行批量分割,以将相互接触的两个或多个独立块石断面分离,为后续块石断面边界的标识、配准及三维重建提供基础。Step S13, using a combination of simple image morphology operations such as erosion and expansion to perform batch segmentation on the binarized two-dimensional CT slice sequence images, so as to separate two or more independent stone sections that are in contact with each other, and form subsequent blocks The identification, registration and 3D reconstruction of the boundary of the stone section provide the basis.

如图4所示,生成块石的三维细观结构模型可通过如下方式进行:As shown in Figure 4, the three-dimensional mesoscopic structure model of block stone can be generated in the following ways:

步骤S20,导入对块石断面分割后的CT断层图像序列(即:土石混合体二维CT切片序列图像)。Step S20, importing the CT tomographic image sequence (ie, the two-dimensional CT slice sequence image of the soil-rock mixture) after segmenting the section of the rock.

步骤S21,依次对每张CT切片图像上的所有块石断面边界进行识别与编号,并计算其周长、面积、形心坐标等几何特征指标,将其存储在二维块石断面边界信息数据库中。Step S21, sequentially identify and number all block section boundaries on each CT slice image, calculate their perimeter, area, centroid coordinates and other geometric feature indicators, and store them in the two-dimensional block section boundary information database middle.

步骤S22,选取土石混合体试样顶面(或底面)的CT切片图像为起始基准图像,按顺序(或逆序)将所有二维CT切片序列图像上的全部块石断面边界与其前一张基准CT图像进行边界配准。根据形态学相似原理,同一个块石在相邻两切片上的断面边界具有相近的几何特征指标,遍历第i张(1<i≤n,n为一个土石混合体试样的二维CT序列中的切片图像总数)CT切片图像上的所有块石断面边界,检测出其中与第i-1张CT切片图像上第j个块石断面边界Si-1,j配准的块石断面边界Si,k,并赋予两者相同的块石编号,将得到的三维块石边界数据存储于三维块石信息数据库中。对于某切片图像中旧块石边界消失与新块石边界出现的情况,及时更新整体块石编号与数据库信息。Step S22, select the CT slice image of the top surface (or bottom surface) of the soil-rock mixture sample as the initial reference image, and sequentially (or reversely) compare all the boundaries of all rock sections on all two-dimensional CT slice sequence images with the previous one Baseline CT images for boundary registration. According to the principle of morphological similarity, the cross-sectional boundaries of the same block on two adjacent slices have similar geometric feature indicators, and traverse the i-th sheet (1<i≤n, n is a two-dimensional CT sequence of a soil-rock mixture sample The total number of slice images in) all block section boundaries on the CT slice image, and detect the block section boundary that is registered with the jth block section boundary S i-1, j on the i-1th CT slice image S i, k , and give them the same block number, and store the obtained three-dimensional block boundary data in the three-dimensional block information database. For the situation that the old block boundary disappears and the new block boundary appears in a certain slice image, the overall block number and database information are updated in time.

步骤S23,利用S21步得到的二维断面边界信息数据库及S22步得到的三维块石信息数据库中所存储的三维块石边界数据,构建每个块石的表面三角网格,实现土石混合体试样内部块石的三维重建,从而获得能反映土石混合体细观结构的三维模型。Step S23, using the two-dimensional cross-section boundary information database obtained in step S21 and the three-dimensional block boundary data stored in the three-dimensional block information database obtained in step S22, to construct the triangular mesh on the surface of each block to realize the soil-rock mixture test. The 3D reconstruction of the internal rocks of the sample can be used to obtain a 3D model that can reflect the mesoscopic structure of the soil-rock mixture.

步骤S24,将建立的土石混合体试样三维模型保存输出为图片或通用的数据格式。通过将三维重建土石混合体试样中的块石模型保存为*.stl或*.gts等通用数据格式,可生成数值模型开展模拟计算并与相应试样的物理试验结果相互对比验证,同时也便于用其它CAD、CAM等软件进行查看和展示。Step S24, saving and outputting the established three-dimensional model of the soil-rock mixture sample as a picture or a common data format. By saving the block model in the 3D reconstructed soil-rock mixture sample as a general data format such as *.stl or *.gts, a numerical model can be generated for simulation calculation and compared with the physical test results of the corresponding sample. It is convenient to use other CAD, CAM and other software to view and display.

如图5所示,是本发明一个实施例中复杂岩土介质细观结构三维重建系统30,主要由三维模型重建模块31、三维细观结构信息统计分析模块32及三维试样可视化模块33三部分构成。该系统主要应用于计算机装置15中,该计算机装置包括通过数据总线连接的存储器16、处理器17、输入设备18和输出设备19。As shown in Fig. 5, it is a three-dimensional reconstruction system 30 of complex rock-soil medium mesostructure in an embodiment of the present invention, mainly composed of a three-dimensional model reconstruction module 31, a three-dimensional mesostructure information statistical analysis module 32 and a three-dimensional sample visualization module 33 Partial composition. The system is mainly implemented in a computer device 15 comprising a memory 16, a processor 17, an input device 18 and an output device 19 connected by a data bus.

所述三维模型重建模块31,包括:图像预处理、块石断面边界标识与配准、三维块石表面重建及模型输出等功能,根据CT扫描得到的土石混合体试样断层序列图像,基于上述方法进行预处理及切片间的块石断面边界匹配,最终实现土石混合体内部块石模型的三维重建,并将三维细观结构模型以通用的数据格式输出,以便数值计算。The three-dimensional model reconstruction module 31 includes functions such as image preprocessing, block section boundary identification and registration, three-dimensional block surface reconstruction, and model output. The method preprocesses and matches the boundaries of block and rock sections between slices, and finally realizes the 3D reconstruction of the block model inside the soil-rock mixture, and outputs the 3D mesostructure model in a common data format for numerical calculation.

所述三维细观结构信息统计分析模块32,用于对土石混合体三维重建试样细观结构信息进行统计分析。根据试样内部块石的三维重建边界计算单个块石的三轴尺寸、长细比、表面积、体积等几何特性指标,球度、棱角性等形态特性指标,及块石长轴定向性等空间排列特性,并存入三维块石信息数据库中。据此统计土石混合体含石量、粒度分布组成及空间分布等结构信息。The three-dimensional mesostructure information statistical analysis module 32 is used for performing statistical analysis on the mesostructure information of the soil-rock mixture three-dimensional reconstruction sample. According to the three-dimensional reconstruction boundary of the block stone inside the sample, calculate the geometric characteristic indexes such as triaxial size, slenderness ratio, surface area, volume, etc. Arrange the characteristics and store them in the three-dimensional block stone information database. Based on this, the structural information such as stone content, particle size distribution composition and spatial distribution of the soil-rock mixture are counted.

所述三维试样可视化模块33,用于对重建后的土石混合体试样三维模型进行可视化显示。根据试样内部块石的三维重建数据,采用三维可视化技术实现三维模型及模型剖切面的渲染、展示,并能通过人机交互实现对三维模型内部结构剖切面的查看及内部块石几何特性、形态特性信息的实时查询。The three-dimensional sample visualization module 33 is used for visually displaying the reconstructed three-dimensional model of the soil-rock mixture sample. According to the 3D reconstruction data of the internal blocks of the sample, the 3D visualization technology is used to realize the rendering and display of the 3D model and the model section plane, and the viewing of the internal structure section plane of the 3D model and the geometric characteristics of the internal block stones can be realized through human-computer interaction. Real-time query of morphological property information.

根据本发明实施例的方法,能够根据土石混合体试样方便地建立其对应的三维结构模型,并对试样内部块石的粒度组成、形态等特性进行定量分析,同时实现模型的三维可视化展示与存储,为土石混合体的物理力学数值试验提供支持。According to the method of the embodiment of the present invention, the corresponding three-dimensional structure model can be conveniently established according to the soil-rock mixture sample, and the characteristics such as the particle size composition and shape of the block stone inside the sample can be quantitatively analyzed, and the three-dimensional visual display of the model can be realized at the same time and storage, providing support for physical and mechanical numerical tests of soil-rock mixtures.

相较于现有的方式,本发明的方法能在对土石混合体三维细观结构建模时,利用土石混合体CT断层序列图像,在不破坏试样的情况下实现土石混合体的三维细观结构模型重建,扫描后的试样可继续进行室内物理力学试验研究。根据重建后的土石混合体试样模型,可对其内部块石几何形态特征、粒度组成、空间分布等信息进行统计分析,并可对模型进行渲染、剖切等查询操作。同时,重建后的模型可直接存储为通用数据格式,以用于数值试验研究。本发明的方法,大大提高了所生成模型的精细化程度,在数值试验与物理力学试验之间建立了桥梁纽带,从而使得在此基础上开展的物理力学性质数值试验研究更加合理与可靠。Compared with the existing methods, the method of the present invention can realize the three-dimensional detailed structure of the soil-rock mixture by using CT tomographic images of the soil-rock mixture when modeling the three-dimensional mesoscopic structure of the soil-rock mixture without destroying the sample. The structural model is reconstructed, and the scanned samples can continue to be studied in the physical and mechanical experiments in the laboratory. According to the reconstructed soil-rock mixture sample model, statistical analysis can be carried out on the geometric characteristics, particle size composition, spatial distribution and other information of the internal block stones, and query operations such as rendering and sectioning can be performed on the model. At the same time, the reconstructed model can be directly stored in a common data format for numerical experiment research. The method of the invention greatly improves the refinement degree of the generated model, and establishes a bridge link between the numerical test and the physical-mechanical test, thereby making the numerical test research of physical-mechanical properties carried out on this basis more reasonable and reliable.

本发明的进一步实施例提供了一种土石混合体三维细观结构重建与分析系统。如图7所示,结合图1、2和5,该系统700包括:土石混合体试样制备装置40(如有机玻璃制样筒)、CT断层图像扫描装置50(如CT扫描机)、三维模型重建模块31和三维细观结构信息统计分析模块32。A further embodiment of the present invention provides a three-dimensional mesoscopic structure reconstruction and analysis system of soil-rock mixture. As shown in Figure 7, in conjunction with Figures 1, 2 and 5, the system 700 includes: a soil-rock mixture sample preparation device 40 (such as a plexiglass sample cylinder), a CT tomographic image scanning device 50 (such as a CT scanner), a three-dimensional Model reconstruction module 31 and three-dimensional mesostructure information statistical analysis module 32.

其中,土石混合体试样制备装置40用于制备土石混合体试样。CT断层图像扫描装置50用于对所述土石混合体试样进行扫描以得到所述土石混合体试样的CT断层图像序列。三维模型重建模块31用于获取所述CT断层图像序列中每一张CT断层图像中的块石断面区域,并对所述块石断面区域进行分割以得到多个独立的块石断面,以及对所述每一张CT断层图像中的多个独立的块石断面进行配准,并根据配准结果生成块石的三维细观结构模型。三维细观结构信息统计分析模块32用于对所述块石的三维细观结构模型进行分析以确定土石混合体试样的含石量、粒度分布、块石空间分布、块石几何特性和块石形态特性。Wherein, the earth-rock mixture sample preparation device 40 is used for preparing the earth-rock mixture sample. The CT tomographic image scanning device 50 is used to scan the soil-rock mixture sample to obtain a CT tomographic image sequence of the soil-rock mixture sample. The three-dimensional model reconstruction module 31 is used to obtain the block section area in each CT tomogram in the CT tomogram sequence, and segment the block section area to obtain multiple independent block sections, and A plurality of independent block sections in each CT tomographic image are registered, and a three-dimensional mesoscopic structure model of the block is generated according to the registration result. The three-dimensional mesostructure information statistical analysis module 32 is used to analyze the three-dimensional mesostructure model of the block to determine the stone content, particle size distribution, spatial distribution of the block, geometric characteristics of the block and block size of the soil-rock mixture sample. stone morphological properties.

进一步地,三维模型重建模块31用于:对所述CT断层图像序列中每一张CT断层图像进行滤波去噪;对滤波去噪后的每一张CT断层图像进行二值化处理以确定所述块石断面区域和土体区域;通过腐蚀和膨胀的图像形态学运算对所述块石断面区域进行分割以得到多个独立的块石断面。Further, the three-dimensional model reconstruction module 31 is configured to: perform filtering and denoising on each CT tomographic image in the CT tomographic image sequence; perform binarization processing on each CT tomographic image after filtering and denoising to determine the The block section area and the soil area; the block section area is segmented by image morphology operation of erosion and expansion to obtain multiple independent block sections.

此外,三维模型重建模块31还用于:依次对每一张CT断层图像上的所有块石断面的边界进行识别与编号,并计算每个块石断面的几何特征指标;对所有CT断层图像上的块石的边界进行配准以得到所述配准结果;根据所述块石断面的几何特征指标和所述配准结果生成所述块石的三维细观结构模型。In addition, the three-dimensional model reconstruction module 31 is also used to: sequentially identify and number the boundaries of all block sections on each CT tomographic image, and calculate the geometric feature index of each block section; The boundary of the block is registered to obtain the registration result; a three-dimensional mesoscopic structure model of the block is generated according to the geometric feature index of the section of the block and the registration result.

结合图5的所示,还包括:三维显示计算机装置33(即三维试样可视化模块),用于对块石的三维细观结构模型进行渲染和展示。In combination with what is shown in FIG. 5 , it also includes: a three-dimensional display computer device 33 (ie, a three-dimensional sample visualization module), which is used to render and display the three-dimensional mesoscopic structure model of the block stone.

需要说明的是,本发明实施例的系统的具体实现方式与方法部分的具体实现方式类似,为了减少冗余,不做赘述。It should be noted that, the specific implementation manner of the system in the embodiment of the present invention is similar to the specific implementation manner of the method part, and will not be repeated in order to reduce redundancy.

根据本发明实施例的系统,能够根据土石混合体试样方便地建立其对应的三维结构模型,并对试样内部块石的粒度组成、形态等特性进行定量分析,同时实现模型的三维可视化展示与存储,为土石混合体的物理力学数值试验提供支持。According to the system of the embodiment of the present invention, the corresponding three-dimensional structure model can be conveniently established according to the soil-rock mixture sample, and the particle size composition, shape and other characteristics of the block stone inside the sample can be quantitatively analyzed, and the three-dimensional visual display of the model can be realized at the same time and storage, providing support for physical and mechanical numerical tests of soil-rock mixtures.

相较于现有的方式,本发明的系统能在对土石混合体三维细观结构建模时,利用土石混合体CT断层序列图像,在不破坏试样的情况下实现土石混合体的三维细观结构模型重建,扫描后的试样可继续进行室内物理力学试验研究。根据重建后的土石混合体试样模型,可对其内部块石几何形态特征、粒度组成、空间分布等信息进行统计分析,并可对模型进行渲染、剖切等查询操作。同时,重建后的模型可直接存储为通用数据格式,以用于数值试验研究。本发明的系统,大大提高了所生成模型的精细化程度,在数值试验与物理力学试验之间建立了桥梁纽带,从而使得在此基础上开展的物理力学性质数值试验研究更加合理与可靠。Compared with the existing methods, the system of the present invention can use the CT tomographic sequence images of the soil-rock mixture when modeling the three-dimensional mesoscopic structure of the soil-rock mixture to realize the three-dimensional detailed structure of the soil-rock mixture without destroying the sample. The structural model is reconstructed, and the scanned samples can continue to be studied in the physical and mechanical experiments in the laboratory. According to the reconstructed soil-rock mixture sample model, statistical analysis can be carried out on the geometric characteristics, particle size composition, spatial distribution and other information of the internal block stones, and query operations such as rendering and sectioning can be performed on the model. At the same time, the reconstructed model can be directly stored in a common data format for numerical experiment research. The system of the present invention greatly improves the degree of refinement of the generated model, and establishes a bridge between the numerical test and the physical-mechanical test, so that the numerical test research of physical-mechanical properties carried out on this basis is more reasonable and reliable.

在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外,在不相互矛盾的情况下,本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。In the description of this specification, descriptions referring to the terms "one embodiment", "some embodiments", "example", "specific examples", or "some examples" mean that specific features described in connection with the embodiment or example , structure, material or characteristic is included in at least one embodiment or example of the present invention. In this specification, the schematic representations of the above terms are not necessarily directed to the same embodiment or example. Furthermore, the described specific features, structures, materials or characteristics may be combined in any suitable manner in any one or more embodiments or examples. In addition, those skilled in the art can combine and combine different embodiments or examples and features of different embodiments or examples described in this specification without conflicting with each other.

尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。Although the embodiments of the present invention have been shown and described above, it can be understood that the above embodiments are exemplary and should not be construed as limiting the present invention, those skilled in the art can make the above-mentioned The embodiments are subject to changes, modifications, substitutions and variations.

Claims (8)

1. a kind of soil-rock mixture three-dimensional microscopical structure is rebuild and analysis method, it is characterised in that comprise the following steps:
Prepare soil-rock mixture sample;
CT scan is carried out to the soil-rock mixture sample to obtain the CT tomographic sequences of the soil-rock mixture sample;
The block stone cross sectional area in each CT faultage image in the CT tomographic sequences is obtained, and to described piece of stone section Split to obtain multiple independent block stone sections in region;
Multiple independent block stone section in each CT faultage image carries out registration, and generates block according to registration result The three-dimensional microscopical structure model of stone, further includes:
The border of all pieces of stone sections on each CT faultage image is identified successively with numbering, and calculate each block stone The geometric properties index of section;
Border to the block stone section on all CT faultage images carries out registration to obtain the registration result, the step of the registration Suddenly include:Each block stone section boundary information for the piecemeal mark in the section of adjacent C T faultage images that will be obtained, according to its shape Step response and locus carry out match cognization, obtain constituting the block stone section border in all CT sections of each block stone;
Geometric properties index and the registration result according to described piece of stone section generate the three-dimensional microscopical structure model of block stone, raw The step of into the three-dimensional microscopical structure model, includes:According to the block stone section information after matching, build and characterize each block stone table The triangle gridding in face, sets up all pieces of threedimensional models of stone in soil-rock mixture sample inside;
And
Three-dimensional microscopical structure model to described block of stone is analyzed to determine rock-soil ratio, the granularity of the soil-rock mixture sample Distribution, block stone spatial distribution, block stone geometrical property and block stone morphological character.
2. method according to claim 1, it is characterised in that wherein, the native stone is prepared by lucite sample preparation cylinder Mixture sample, is scanned by CT Scanner to the soil-rock mixture sample.
3. method according to claim 1, it is characterised in that each CT tomography in the acquisition CT tomographic sequences Block stone cross sectional area in image, and described piece of stone cross sectional area is split to obtain multiple independent block stone sections, enter One step includes:
Denoising is filtered to each CT faultage image in the CT tomographic sequences;
Each CT faultage image after to filtering and noise reduction carries out binary conversion treatment to determine described piece of stone cross sectional area and the soil body Region;
Described piece of stone cross sectional area is split by the morphological image computing corroded and expand multiple independent to obtain Block stone section.
4. the method according to claim any one of 1-3, it is characterised in that the three-dimensional microscopical structure model to block stone It is analyzed to determine rock-soil ratio, size distribution, block stone spatial distribution, block stone geometrical property and the block stone of soil-rock mixture sample Morphological character, further includes:
The three-dimensional microscopical structure model of block stone is rendered and shown, is rendered and shown that result determines that the native stone is mixed with basis The rock-soil ratio of fit sample, size distribution, block stone spatial distribution, block stone geometrical property and block stone morphological character.
5. a kind of soil-rock mixture three-dimensional microscopical structure is rebuild and analysis system, it is characterised in that including:
Soil-rock mixture sample preparation device, for preparing soil-rock mixture sample;
CT faultage image scanning means, for carrying out CT scan to obtain the soil-rock mixture to the soil-rock mixture sample The CT tomographic sequences of sample;
Reconstructing three-dimensional model module, for obtaining the block stone section in the CT tomographic sequences in each CT faultage image Region, and described piece of stone cross sectional area is split to obtain multiple independent block stone sections, the reconstructing three-dimensional model mould Block is additionally operable to:
The border of all pieces of stone sections on each CT faultage image is identified successively with numbering, and calculate each block stone The geometric properties index of section;
Border to the block stone section on all CT faultage images carries out registration to obtain the registration result, the step of the registration Suddenly include:Each block stone section boundary information for the piecemeal mark in the section of adjacent C T faultage images that will be obtained, according to its shape Step response and locus carry out match cognization, obtain constituting the block stone section border in all CT sections of each block stone;
Geometric properties index and the registration result according to described piece of stone section generate the three-dimensional microscopical structure model of block stone, raw The step of into the three-dimensional microscopical structure model, includes:According to the block stone section information after matching, build and characterize each block stone table The triangle gridding in face, sets up all pieces of threedimensional models of stone in soil-rock mixture sample inside;
And
Three-dimensional microscopical structure Information Statistics analysis module, is analyzed with true for the three-dimensional microscopical structure model to described block of stone Rock-soil ratio, size distribution, block stone spatial distribution, block stone geometrical property and the block stone form for determining the soil-rock mixture sample are special Property.
6. system according to claim 5, it is characterised in that wherein, the soil-rock mixture sample preparation device is to have Machine glass sample preparation cylinder, the CT faultage images scanning means is CT Scanner.
7. system according to claim 5, it is characterised in that the reconstructing three-dimensional model module is used for:
Denoising is filtered to each CT faultage image in the CT tomographic sequences;
Each CT faultage image after to filtering and noise reduction carries out binary conversion treatment to determine described piece of stone cross sectional area and the soil body Region;
Described piece of stone cross sectional area is split by the morphological image computing corroded and expand multiple independent to obtain Block stone section.
8. the system according to claim any one of 5-7, it is characterised in that also include:Three-dimensional Display computer installation, uses Rendered and shown in the three-dimensional microscopical structure model to block stone.
CN201410293183.0A 2014-06-25 2014-06-25 Soil-rock mixture three-dimensional microscopical structure is rebuild and analysis method and system Expired - Fee Related CN104036538B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410293183.0A CN104036538B (en) 2014-06-25 2014-06-25 Soil-rock mixture three-dimensional microscopical structure is rebuild and analysis method and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410293183.0A CN104036538B (en) 2014-06-25 2014-06-25 Soil-rock mixture three-dimensional microscopical structure is rebuild and analysis method and system

Publications (2)

Publication Number Publication Date
CN104036538A CN104036538A (en) 2014-09-10
CN104036538B true CN104036538B (en) 2017-05-31

Family

ID=51467296

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410293183.0A Expired - Fee Related CN104036538B (en) 2014-06-25 2014-06-25 Soil-rock mixture three-dimensional microscopical structure is rebuild and analysis method and system

Country Status (1)

Country Link
CN (1) CN104036538B (en)

Families Citing this family (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104268830A (en) * 2014-09-11 2015-01-07 河海大学 Digital-image-based method for determining permeability coefficient of nonuniform geotechnical material
CN104331922B (en) * 2014-11-01 2017-03-15 西安科技大学 A kind of Soil Microstructure three-dimensional rebuilding method
CN105424430A (en) * 2015-11-05 2016-03-23 河海大学 Manufacturing method of cemented soil and stone mixed sample
CN105956314A (en) * 2016-05-16 2016-09-21 河海大学 Numerical test method capable of creating different earth-rock mixing ratios
CN107829477B (en) * 2017-02-20 2020-05-22 浙江贝乐卫浴科技有限公司 Intelligent sanitary appliance controlled through Internet of things
CN107680058A (en) * 2017-09-28 2018-02-09 哈尔滨工业大学深圳研究生院 Permanent delay smooth filtering method based on green space in CT image denoisings
CN108376422A (en) * 2018-01-05 2018-08-07 中山大学 A kind of continuous across the scale measurement characterizing method of rock-soil material
CN108961388B (en) * 2018-06-05 2023-03-21 哈尔滨工业大学深圳研究生院 Microstructure three-dimensional modeling method, microstructure three-dimensional modeling device, microstructure three-dimensional modeling equipment and computer storage medium
CN110907317A (en) * 2018-09-14 2020-03-24 中石化石油工程技术服务有限公司 Rock particle analysis method, apparatus, device and computer readable storage device
CN109556905A (en) * 2018-09-28 2019-04-02 中山大学 Resin particle syncretizing effect analyzer
CN109145520B (en) * 2018-10-22 2023-03-24 重庆大学 Earth-rock mixture tunnel design method based on digital image and big data
CN109509251B (en) * 2018-11-08 2023-02-28 重庆大学 Multi-factor three-dimensional soil-rock mixture generation method
CN109671056B (en) * 2018-12-03 2020-10-27 西安交通大学 A method for detecting porosity defects in composite sleepers based on X-ray images
CN109632429B (en) * 2019-01-30 2020-06-16 中国水利水电科学研究院 A kind of sample preparation method of soil-rock mixture biaxial compression test
CN110069844B (en) * 2019-04-15 2020-10-30 武汉大学 Mesoscopic model generation method considering rock structure characteristics and mineral composition
CN110363855B (en) * 2019-07-22 2020-12-11 四川大学 Transparency modeling method of rockfill dam
CN115644951B (en) * 2022-11-20 2025-01-28 成都真实维度科技有限公司 A precise puncture sampling system based on real-time modeling
CN118658161B (en) * 2024-08-21 2024-11-08 中国航空工业集团公司金城南京机电液压工程研究中心 Method and device for identifying microstructure of three-dimensional braided composite materials using machine vision

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102609980A (en) * 2012-01-18 2012-07-25 西安建筑科技大学 Three-dimensional reconstruction method for concrete CT (computed tomography) image
CN103679810A (en) * 2013-12-26 2014-03-26 海信集团有限公司 Method for three-dimensional reconstruction of liver computed tomography (CT) image

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102609980A (en) * 2012-01-18 2012-07-25 西安建筑科技大学 Three-dimensional reconstruction method for concrete CT (computed tomography) image
CN103679810A (en) * 2013-12-26 2014-03-26 海信集团有限公司 Method for three-dimensional reconstruction of liver computed tomography (CT) image

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
土石混合体研究现状及发展趋势;徐文杰,张海洋;《水利水电科技进展》;20130131;第33卷(第1期);第80-88页 *
土石混合体细观结构及力学特性数值模拟研究;徐文杰 等;《岩石力学与工程学报》;20070228;第26卷(第2期);第300-311页 *
基于 CT 图像的混凝土细观三维有限元模型重建;郝书亮;《中国优秀硕士学位论文全文数据库信息科技辑》;20111215(第S1期);摘要,第2.4节,第3.1节,第3.4-3.5节 *
基于数字图像处理的土石混合体细观结构;徐文杰 等;《辽宁工程技术大学学报》;20080228;第27卷(第1期);第51-53页 *

Also Published As

Publication number Publication date
CN104036538A (en) 2014-09-10

Similar Documents

Publication Publication Date Title
CN104036538B (en) Soil-rock mixture three-dimensional microscopical structure is rebuild and analysis method and system
CN109697752A (en) Based on rock core CT image hole information extraction quantitatively characterizing rock core heterogeneity method
CN106053168B (en) Concrete based on 3D printing technique carefully sees the method for visualizing of three-phase structure
CN103698803B (en) Rock pore structure characterization method and device
CN102867302B (en) Core fracture identification method based on three-dimensional image information processing
CN102768699B (en) Method based on CT image Accurate Reconstruction dissimilar materials microcosmic FEM mesh
CN110136249A (en) A method for 3D visualization of reservoir rock pores and fractures and simulation of gas flow
CN110146525A (en) A Prediction Method of Coal Porosity and Permeability Parameters Based on Fractal Theory and CT Scanning
CN104809756B (en) Asphalt mixture gap spatial construction method based on X-ray CT images
CN105445160B (en) The void characteristics and its extracting method of a kind of asphalt
CN104008573B (en) Grading analysis method and particle based on three-dimensional configuration place stent
CN111524229A (en) Three-dimensional geometric morphology information extraction system and method for rock particles
CN109490165A (en) The method for characterizing the non-fabric selectivity reservoir space of carbonate rock
CN105913432A (en) Aorta extracting method and aorta extracting device based on CT sequence image
CN106324002A (en) Carbonatite pore structure characterization method based on rock classification and multi-scale digital cores
CN104766336B (en) A kind of solid engines three dimensional CT defect extraction and labeling method
Laycock et al. Combining x-ray micro-CT technology and 3D printing for the digital preservation and study of a 19th Century Cantonese chess piece with intricate internal structure
Friese et al. Analysis of tomographic mineralogical data using YaDiV—Overview and practical case study
CN104217437A (en) Lesion area partitioning method for prostate KVCT image
CN107945278A (en) The three-dimensionalreconstruction and finite element method of a kind of organism structure
CN114565722A (en) Three-dimensional model monomer realization method
CN108399354A (en) The method and apparatus of Computer Vision Recognition tumour
CN108362716A (en) A kind of historical relic material detection determination method and detection device based on Medical CT
CN102512140A (en) Locating optical projection coherence tomography imaging rotation center method
CN103955943A (en) Non-supervision change detection method based on fuse change detection operators and dimension driving

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170531

Termination date: 20200625

CF01 Termination of patent right due to non-payment of annual fee