Rock-soil body porosity distribution in-situ estimation method based on borehole endoscopic image

The stereoscopic image of the drilled wall surface is obtained through binocular imaging technology, and the concaveness and porosity are calculated, which solves the problem that the porosity distribution of the entire pore segment of the rock and soil body cannot be obtained in real time in the existing technology, and realizes efficient porosity measurement and distribution analyze.

CN119290706APending Publication Date: 2025-01-10INST OF ROCK & SOIL MECHANICS CHINESE ACAD OF SCI
View PDF 0 Cites 0 Cited by

Patent Information

Application Number
CN202411529215.2
Authority / Receiving Office
CN · China
Patent Type
Applications(China)
Current Assignee / Owner
Filing Date
2024-10-30
Publication Date
2025-01-10

AI Technical Summary

Technical Problem

The existing methods of porosity measurement of rock and soil bodies require centering or relying on laboratory tests, and the porosity distribution of the entire pore segment cannot be obtained in real time, and the results are not intuitive enough.

Method used

Binocular imaging technology is used to obtain continuous three-dimensional images of the drilled wall surface, calculate the concaveness of any pixel point on the hole wall, analyze the concaveness distribution characteristics, determine the pore calculation range, calculate the pore volume and total volume, and then estimate the pores of the rock and soil body Spend.

Benefits of technology

The three-dimensional volume porosity of the rock and soil body is achieved in situ, and the porosity at any depth and orientation can be calculated, which simplifies the continuous analysis of porosity distribution characteristics and improves the simplicity and efficiency of measurement.

✦ Generated by Eureka AI based on patent content.

Smart Images

  • Figure CN119290706A_ABST
    Figure CN119290706A_ABST
Patent Text Reader

Abstract

The invention provides a rock-soil body porosity distribution in-situ estimation method based on a borehole endoscopic image. According to the method, a hole wall three-dimensional image is obtained through a binocular imaging technology, the distance between any pixel point of the wall face of a drilled hole and a central axis of the drilled hole is calculated based on the three-dimensional image, the concavity of any pixel point of the wall face of the drilled hole is calculated, the calculation range of the porosity of the rock-soil body is determined, the porosity of the rock-soil body is calculated, and the distribution characteristics of the porosity of the rock-soil body can be statistically analyzed. According to the method, the three-dimensional volume porosity of the hole wall rock mass is obtained in an in-situ mode, the porosity at any depth and in any direction can be calculated, continuous analysis of the pore characteristics of the whole hole can be achieved, and the method does not depend on a rock core and is simple and efficient.
Need to check novelty before this filing date? Find Prior Art

Description

Technical Field

[0001] The present invention relates to the field of engineering technology, and in particular to an in-situ estimation method for porosity distribution of rock and soil bodies based on borehole endoscopic images. Background Art

[0002] Pores are the unfilled spaces between solid mineral particles in rock and soil materials, and are also the main storage and transportation channels for groundwater. In rock and soil landslide engineering, pore water becomes an important factor affecting landslide disasters. Therefore, the measurement of pore characteristics of rock and soil in landslide engineering is crucial.

[0003] Porosity is a common pore characterization parameter and one of the important indicators for measuring the water storage capacity of rock and soil. Conventional rock and soil porosity measurement methods are mainly through laboratory measurement and in-situ measurement.

[0004] Laboratory rock and soil porosity measurement requires in-situ coring, and then the standard cylindrical rock and soil samples are subjected to water saturation, mercury injection and other tests to obtain their porosity. The laboratory measurement method can accurately measure the pore content in the sample, but it is a huge workload to obtain the porosity distribution of the entire hole section, and the porosity cannot be obtained in the hole section where it is difficult to coring.

[0005] In-situ measurement is an important means to obtain the continuous distribution of porosity. At present, conventional in-situ porosity testing methods include electrical logging and radioactive logging, which usually require interpretation to obtain porosity. The results are not intuitive and relatively qualitative. Summary of the invention

[0006] Based on this, it is necessary to provide an in-situ estimation method for the porosity distribution of rock and soil based on borehole endoscopic images, which can obtain the three-dimensional volume porosity of the rock and soil on the hole wall in an in-situ manner without relying on the core.

[0007] The present invention adopts the following technical solution:

[0008] The present invention provides an in-situ estimation method for porosity distribution of rock and soil bodies based on borehole endoscopic images, comprising the following steps:

[0009] A continuous stereoscopic image of the borehole wall is obtained by using binocular imaging technology, and the stereoscopic image of the borehole wall contains depth and azimuth information; the concavity of any pixel point on the borehole wall is calculated based on the stereoscopic image of the borehole wall; the concavity distribution characteristics are analyzed to determine the calculation range of the rock mass porosity, and the pore volume and the total volume of the rock mass are calculated; the porosity within the determined range is calculated according to the definition of the porosity of the rock mass.

[0010] Furthermore, the porosity distribution characteristics are statistically analyzed based on the porosity calculation results and the orientation and depth information in the pore wall image.

[0011] Specifically, the in-situ estimation method of the porosity distribution of rock and soil based on borehole endoscopic images includes the following steps:

[0012] The binocular imaging probe in the binocular imaging technology is placed at the axis position of the circular borehole in the rock and soil body;

[0013] Acquire a continuous stereoscopic image of the borehole wall of the rock and soil body by using binocular imaging, wherein the stereoscopic image of the borehole wall contains depth and orientation information;

[0014] A polar coordinate system is established with the center of the borehole as the origin. The following definitions are used: r is the distance between the borehole radial direction and the central axis, r = R is the borehole wall, θ is the azimuth, θ = 0° is the due north direction, z is the borehole depth, z = 0 at the borehole mouth, and the coordinates of any pixel point on the borehole wall are expressed as P rθz ;

[0015] The distance between any pixel point on the hole wall and the drilling axis is calculated using the binocular ranging algorithm and recorded as L rθz , the distance between any pixel point on the hole wall and the hole wall is recorded as D rθz , calculate the concavity D of any pixel point on the hole wall rθz =L rθz -R;

[0016] The concavity D of any pixel point on the hole wall rθz Perform statistical analysis to determine the concavity D rθz The average value D' is used to determine the calculation range of the rock and soil porosity: r = R ~ D', z = h1 ~ h2, h1 and h2 are the actual values ​​of the drilling depth;

[0017] The total pore volume is calculated based on the sum of the spatial volumes represented by all pixel points in the porosity calculation range of the rock and soil mass:

[0018] V P =∑V P rθz

[0019] Among them, V P is the total pore volume, V P rθz is the spatial volume represented by any pixel point, V P rθz The calculation formula is:

[0020] V P rθz =S×D rθz'

[0021] Where S is the actual area of ​​the pixel; when the concavity D rθz When it is greater than or equal to the average concavity D', the characterization concavity is D rθz' =D', when the pixel concavity D rθzWhen the concavity is less than the average concavity D', the concavity is represented by D rθz' =D rθz ; When the pixel concavity D rθz = 0, the pixel point is located on the inner wall of the borehole, and the corresponding pore volume at this position is 0;

[0022] Calculate the total volume V of the rock and soil mass A :

[0023] V A =N×S×D'

[0024] Where N is the number of image pixels, S is the actual area of ​​the pixel, and D' is the average concavity;

[0025] The porosity of rock and soil is calculated according to the definition of rock and soil porosity.

[0026] The in-situ estimation method of the porosity distribution of rock and soil bodies based on borehole endoscopic images may also include the step of calculating the depth porosity or the position porosity.

[0027] The depth porosity is defined as the ratio of the pore volume to the total volume on any depth line, and the calculation formula is:

[0028]

[0029] in, is the depth porosity, h1 is the image start depth, h2 is the image end depth, Δh=h2-h1, V is the number of pixels in each row, H is the number of pixels in each column, and V×H=N.

[0030] The azimuthal porosity is defined as the ratio of the pore volume on any azimuthal line to the total volume, and the calculation formula is:

[0031]

[0032] in, is the depth porosity, H is the number of pixels in each column, and D' is the average concavity.

[0033] The present invention can also provide a rock and soil body porosity in-situ measurement device, comprising:

[0034] An image acquisition module, used to acquire a stereoscopic image of the borehole wall of the rock and soil body;

[0035] a data analysis module for obtaining depth and orientation information from the hole wall stereo image using the hole wall stereo image;

[0036] The calculation module is used to calculate the concavity of any pixel point on the hole wall according to the depth and orientation information, analyze the concavity distribution characteristics, determine the calculation range of rock mass porosity, and calculate the porosity.

[0037] Furthermore, the calculation module is also used to calculate the depth porosity or position porosity and output the distribution characteristics of the porosity.

[0038] The present invention can also provide a rock and soil body porosity in-situ measurement device, comprising:

[0039] Memory, for storing computer programs;

[0040] A processor is used to implement the steps of the above-mentioned method for in-situ estimation of porosity distribution of rock and soil bodies based on borehole endoscopic images when executing the computer program.

[0041] The present invention may also provide a computer-readable storage medium having a computer program stored thereon, and when the computer program is executed by a processor, the steps of the above-mentioned in-situ estimation method of the porosity distribution of rock and soil bodies based on borehole endoscopic images are implemented.

[0042] Compared with the prior art, the core advantages of the present invention are:

[0043] The present invention obtains the three-dimensional volume porosity of the hole wall rock and soil in situ, does not rely on the rock core, can calculate the porosity at any depth and orientation, and then analyze the distribution characteristics of the porosity, and can realize continuous analysis of the pore characteristics of the entire hole, which is simpler and more efficient. BRIEF DESCRIPTION OF THE DRAWINGS

[0044] Figure 1 Schematic diagram of binocular imaging technology measurement.

[0045] Figure 2 Schematic diagram of the hole wall image obtained by binocular imaging technology.

[0046] Figure 3 Schematic diagram of drilling coordinate system.

[0047] Figure 4 Schematic diagram of hole wall concavity.

[0048] Figure 5 Schematic diagram of the calculation range of rock and soil volume.

[0049] Figure 6 This is a top-down diagram of the scope of rock and soil volume calculation.

[0050] Figure 7 Schematic diagram of pore volume and pixel points representing spatial volume.

[0051] Figure 8 Schematic diagram of porosity depth distribution characteristics.

[0052] Fig. 9 Schematic diagram of porosity distribution characteristics.

[0053] Fig.10 This is the hole wall image obtained in Example 1.

[0054] Fig.11 The results are compared between the measurement results of the method of the present invention and the core water saturation method in Example 1.

[0055] Fig.12 This is the hole wall image obtained in Example 2.

[0056] Fig.13 The following is a comparison of the measurement results using the method of the present invention and the core water saturation method in Example 2.

[0057] Fig.14 This is the hole wall image obtained in Example 3.

[0058] Fig.15 The following is a comparison of the measurement results using the method of the present invention and the core water saturation method in Example 3. DETAILED DESCRIPTION

[0059] The present invention is further described in detail below in conjunction with specific embodiments so that those skilled in the art can more clearly understand the present invention. The following embodiments are only used to illustrate the present invention, but are not only used to limit the scope of the present invention. Based on the specific embodiments in the present invention, all other embodiments obtained by those of ordinary skill in the art without making creative work belong to the protection scope of the present invention. In the embodiments of the present invention, if not specifically specified, the technical means used are conventional means well known to those skilled in the art.

[0060] like Figure 1 The equipment for in-situ measurement of rock and soil porosity includes a binocular camera, a binocular imaging probe, an electronic compass (for determining orientation information) and a depth measurement mechanism. With the help of binocular imaging technology, 360° image acquisition of the borehole wall is performed along the depth direction of the borehole. At the same time, with the help of binocular imaging technology, the depth and orientation information of the borehole wall can be collected, and the synchronization of image, depth and orientation information can be ensured, which can form the following Figure 2 The hole wall expansion and curling diagrams are shown.

[0061] like Figures 3 to 9 The present invention provides an in-situ estimation method for porosity distribution of rock and soil bodies based on borehole endoscopic images, comprising the following steps:

[0062] S1, using binocular imaging technology to obtain continuous high-resolution stereo images of the borehole wall (referred to as "high-resolution stereo images of the borehole wall"), the image contains depth and orientation information.

[0063] During the image acquisition process, it is necessary to ensure that the binocular imaging probe used for image acquisition is located at the drilling axis.

[0064] S2, based on the high-resolution stereo image of the hole wall obtained by binocular imaging technology, calculates the concavity of any pixel point on the hole wall.

[0065] like Figure 3 As shown in the figure, after obtaining the high-resolution stereoscopic image of the borehole wall, a polar coordinate system is established with the center of the borehole mouth as the coordinate origin, where R is the borehole radius, r is the distance between the borehole radial direction and the central axis, r = R is the borehole wall, θ is the azimuth, θ = 0° is the north direction, z is the borehole depth, and z = 0 at the borehole mouth. Therefore, the coordinates of each pixel point can be expressed as P rθz .

[0066] like Figure 4 As shown, the concavity of any pixel point on the hole wall is further calculated. The concavity is defined as the distance between any pixel point on the hole wall and the borehole wall, denoted as D rθz The binocular ranging algorithm is used to calculate the distance between any pixel point on the hole wall and the drilling axis, which is recorded as L rθz .

[0067] Calculate the concavity D of any pixel point on the hole wall rθz , the calculation formula is:

[0068] D rθz =L rθz -R

[0069] Where R is the drilling radius.

[0070] S3, analysis of concavity D xy Distribution characteristics, determine the appropriate range of pore characteristics analysis of rock and soil.

[0071] When analyzing the pore distribution characteristics of rock and soil, it is necessary to clarify the scope of the rock and soil. The volume calculation range of the rock and soil should be based on the concavity D rθz The characteristics of the soil body are reasonably set to avoid the pore characteristics within the calculation range of the soil body being unrepresentative. The specific steps are:

[0072] The concavity D of any pixel point on the hole wall rθz Perform statistical analysis to determine the concavity D rθz The average of the two is denoted as D'.

[0073] The calculation range of rock and soil volume is determined by four boundaries, namely r = R, r = D', z = h1 and z = h2, such as Figure 5 and Figure 6 As shown in the figure, r=R and r=D' are fixed boundaries of the volume calculation range of the rock and soil mass, while z=h1 and z=h2 can be set according to actual needs.

[0074] For example, when analyzing the pore characteristics of the entire borehole, h1 = 0 and h2 is the actual depth of the borehole. When analyzing the pore characteristics of the rock mass at a specific depth, h1 and h2 are the specific depth ranges.

[0075] S4, calculate the pore volume and the total volume of the rock and soil within the determined range;

[0076] After determining the volume calculation range of the rock and soil mass, the pore volume and the total volume of the rock and soil mass within the determined range must first be calculated. The pore volume is the sum of the volumes of all pores within the determined volume calculation range, and the rock and soil mass volume is the range defined by the four boundaries. The specific steps are:

[0077] like Figure 7 As shown in Figure 2, the pore volume can be determined by the sum of the spatial volumes represented by all pixels, and the calculation formula is:

[0078] V P =∑V P rθz

[0079] Among them, V P is the total pore volume, V P rθz is the spatial volume represented by any pixel point, V P rθz It can be determined by the product of the actual area represented by the pixel and the concavity it represents;

[0080] V P rθz =S×D rθz'

[0081] Among them, S is the actual area of ​​the pixel point, which is determined by the number of pixels of the binocular camera and the size of the photographed object. In the hole wall image of the same borehole obtained by the same device, the actual area S of the pixel point is a fixed value. rθz' The pixel represents the concavity. When the pixel concavity D rθz When it is greater than or equal to the average concavity D', the characterization concavity is D rθz' =D', when the pixel concavity D rθz When the concavity is less than the average concavity D', the concavity is represented by D rθz' =D rθz When the pixel concavity D rθz =0, the pixel point is located on the inner wall of the borehole, and the corresponding pore volume at this position is 0.

[0082] Calculate the total volume of the rock and soil. The total volume of the rock and soil is the volume of the space enclosed by the boundaries r = R, r = D', z = h1 and z = h2, which can be calculated as:

[0083] V A =N×S×D'

[0084] Among them, V A is the total volume of the rock and soil mass, and N is the number of image pixels.

[0085] S5, calculate the porosity within a certain range according to the definition of rock and soil porosity.

[0086] The porosity of rock and soil is defined as the ratio of the pore volume in the rock and soil to the volume of the rock and soil. The calculation formula is:

[0087]

[0088] In the formula, The porosity of rock mass.

[0089] It can be seen from the formula that when calculating the porosity of rock and soil, this method only needs to obtain the concavity of the pixel points of the hole wall through binocular imaging technology to calculate the porosity of the rock and soil, which is relatively simple compared with other methods.

[0090] S6, statistically analyzing the distribution characteristics of the porosity based on the porosity calculation results and the orientation and depth information in the pore wall image.

[0091] A porosity calculation method is given in S5. After obtaining the concavity of the hole wall pixel points, this method can calculate the porosity of the rock and soil within any depth range.

[0092] In order to statistically analyze the distribution characteristics of rock and soil porosity, it is necessary to refine the calculation range of porosity. Combining binocular imaging technology, the characteristics of pore wall depth and orientation can be obtained, and the concepts of depth porosity and azimuth porosity can be defined.

[0093] Through the calculation and analysis of depth porosity and azimuthal porosity, the distribution characteristics of rock and soil porosity are statistically analyzed. The specific steps are:

[0094] 1) Calculate the depth porosity

[0095] The depth porosity is defined as the ratio of the pore volume to the total volume on any depth line (specifically any row of pixels). The calculation formula is:

[0096]

[0097] in, is the depth porosity, h1 is the image start depth, h2 is the image end depth, Δh=h2-h1, V is the number of pixels in each row, H is the number of pixels in each column, and V×H=N.

[0098] According to the above formula, H depth porosities are calculated.

[0099] 2) Calculate azimuthal porosity

[0100] Azimuthal porosity is defined as the ratio of the pore volume on any azimuthal line (specifically any column of pixels) to the total volume. The calculation formula is:

[0101]

[0102] in, is the depth porosity.

[0103] According to the above formula, V azimuthal porosities are calculated.

[0104] Combining the depth and azimuth information in the image, as well as the calculated H depth porosities and V azimuth porosities, the depth distribution map and azimuth distribution map of porosity are drawn, as shown in Figure 8 and Fig. 9 shown.

[0105] The following is verified with actual application cases.

[0106] Example 1

[0107] Using the in-situ estimation method of the porosity distribution of rock and soil based on borehole endoscopic images of the present invention, the in-situ test site is located in Xifeng County, Guizhou Province. With the help of a digital borehole stereo imaging device, the borehole wall image is obtained, such as Fig.10 As shown; the porosity calculation method of the present invention is used to calculate the porosity between 11 and 13 m of the pore wall, and the calculation result is: the average porosity

[0108] In order to verify the reliability of the above estimation method, the core samples at the depth of 11 to 13 m in the corresponding boreholes were uniformly sampled and processed into three standard cylindrical specimens with a diameter of 50 mm and a height of 100 mm. Conventional core saturation method porosity test was carried out.

[0109] like Fig.11 As shown, the porosity test results of the three rock samples by the saturated water method are 3.86%, 3.53% and 2.78% respectively, and the average porosity is 3.39%. The test results of the present invention differ from the test results of the core saturated water method by 0.83%. Considering that different measurement methods have certain operating errors, the measurement results of the method described in the present invention are basically reliable.

[0110] Example 2

[0111] The in-situ estimation method of the porosity distribution of rock and soil based on borehole endoscopic images of the present invention is used. The in-situ test site is located in Yongwen County, Guizhou Province. With the help of a digital borehole stereo imaging device, the borehole wall image is obtained, such as Fig.12 As shown; the porosity calculation method of the present invention is used to calculate the porosity between 23 and 25 m of the pore wall, and the calculation result is: the average porosity

[0112] In order to verify the reliability of the above estimation method, the cores at the depth of 23 to 25 m in the corresponding boreholes were uniformly sampled and processed into three standard cylindrical specimens with a diameter of 50 mm and a height of 100 mm. Conventional core saturation method porosity test was carried out.

[0113] like Fig.13 As shown, the porosity test results of the three rock samples by the saturated water method are 6.43%, 8.18% and 5.85% respectively, and the average porosity is 6.82%. The test results of the present invention differ from the test results of the core saturated water method by 0.40%. Considering that different measurement methods have certain operating errors, the measurement results of the method described in the present invention are basically reliable.

[0114] Example 3

[0115] The in-situ estimation method of the porosity distribution of rock and soil based on borehole endoscopic images of the present invention is used. The in-situ test site is located in Kunming City, Yunnan Province. With the help of a digital borehole stereo imaging device, the borehole wall image is obtained, such as Fig.14 As shown; the porosity calculation method of the present invention is used to calculate the porosity between 1 and 3 m of the pore wall, and the calculation result is: the average porosity

[0116] In order to verify the reliability of the above estimation method, the cores at the depth of 1 to 3 m in the corresponding boreholes were uniformly sampled and processed into 4 standard cylindrical specimens with a diameter of 50 mm and a height of 100 mm. Conventional core saturation method porosity test was carried out.

[0117] like Fig.13 As shown, the porosity test results of the 4 rock samples by the saturated water method are 5.51%, 3.05%, 3.97% and 4.46% respectively, and the average porosity is 4.25%. The test results of the present invention differ from the test results of the core saturated water method by 0.58%. Considering that different measurement methods have certain operating errors, the measurement results of the method described in the present invention are basically reliable.

[0118] The present invention can further provide an in-situ measurement device for rock and soil porosity, comprising: an image acquisition module for acquiring a stereoscopic image of the borehole wall of the rock and soil body; a data analysis module for acquiring depth and orientation information from the stereoscopic image of the borehole wall using the stereoscopic image of the borehole wall; a calculation module for calculating the concavity of any pixel point on the borehole wall according to the depth and orientation information, analyzing the distribution characteristics of the concavity, determining the calculation range of the rock mass porosity, and calculating the porosity. The calculation module is also used to calculate the depth porosity or the position porosity, and output the distribution characteristics of the porosity.

[0119] The present invention can further provide an in-situ measurement device for rock and soil porosity, comprising: a memory for storing a computer program; a processor for implementing the steps of an in-situ estimation method for rock and soil porosity distribution based on borehole endoscopic images when executing the computer program.

[0120] The present invention can further provide a computer-readable storage medium having a computer program stored thereon, and when the computer program is executed by a processor, the steps of the in-situ estimation method of the porosity distribution of rock and soil bodies based on borehole endoscopic images are implemented.

[0121] It is necessary to point out that the above embodiments are limited to further elaboration and explanation of the technical solution of the present invention, and are not further limitations of the technical solution of the present invention. The method of the present invention is only a preferred implementation scheme, and is not used to limit the protection scope of the present invention. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included in the protection scope of the present invention.

Claims

1. A method for in-situ estimation of porosity distribution of rock and soil based on borehole endoscopic images, characterized in that: The following steps are involved: Using binocular imaging technology to obtain a continuous stereoscopic image of the borehole wall, wherein the stereoscopic image of the borehole wall contains depth and orientation information; Calculating the concavity of any pixel point on the hole wall based on the hole wall stereo image; Analyze the concavity distribution characteristics, determine the calculation range of rock mass porosity, and calculate the pore volume and total rock mass volume; The porosity within a certain range is calculated based on the definition of rock and soil porosity.

2. The in-situ estimation method for porosity distribution of rock and soil bodies based on borehole endoscopic images according to claim 1 is characterized in that: The following steps are involved: Place the binocular imaging probe at the axis of the circular borehole in the rock mass; Acquire a continuous stereoscopic image of the borehole wall of the rock and soil body by using binocular imaging, wherein the stereoscopic image of the borehole wall contains depth and orientation information; A polar coordinate system is established with the center of the borehole as the origin. The following definitions are used: r is the distance between the borehole radial direction and the central axis, r = R is the borehole wall, θ is the azimuth, θ = 0° is the due north direction, z is the borehole depth, z = 0 at the borehole mouth, and the coordinates of any pixel point on the borehole wall are expressed as P rθz ; The distance between any pixel point on the hole wall and the drilling axis is calculated using the binocular ranging algorithm and recorded as L rθz , the distance between any pixel point on the hole wall and the hole wall is recorded as D rθz , calculate the concavity D of any pixel point on the hole wall rθz =L rθz -R; The concavity D of any pixel point on the hole wall rθz Perform statistical analysis to determine the concavity D rθz The average value D' is used to determine the calculation range of the rock and soil porosity: r = R ~ D', z = h1 ~ h2, h1 and h2 are the actual values ​​of the drilling depth; The total pore volume is calculated based on the sum of the spatial volumes represented by all pixel points in the porosity calculation range of the rock and soil mass: In P =∑V P rθz Among them, V P is the total pore volume, V P rθz is the spatial volume represented by any pixel point, V P rθz The calculation formula is: In P rθz =S×D rθz' Where S is the actual area of ​​the pixel; when the concavity D rθz When it is greater than or equal to the average concavity D', the concavity is represented by D rθz' =D', when the pixel concavity D rθz When the concavity is less than the average concavity D', the concavity is represented by D rθz' =D rθz ; When the pixel concavity D rθz = 0, the pixel point is located on the inner wall of the borehole, and the corresponding pore volume at this position is 0; Calculate the total volume V of the rock and soil mass A : In A =N×S×D' Where N is the number of image pixels, S is the actual area of ​​the pixel, and D' is the average concavity; The porosity of rock and soil is calculated according to the definition of rock and soil porosity.

3. The in-situ estimation method for porosity distribution of rock and soil bodies based on borehole endoscopic images according to claim 2 is characterized in that: The method may further include the step of calculating the depth porosity or the position porosity; The depth porosity is defined as the ratio of the pore volume to the total volume on any depth line, and the calculation formula is: in, is the depth porosity, h1 is the image start depth, h2 is the image end depth, Δh = h2 - h1, V is the number of pixels in each row, H is the number of pixels in each column, and V × H = N; The azimuthal porosity is defined as the ratio of the pore volume on any azimuthal line to the total volume, and the calculation formula is: in, is the depth porosity, H is the number of pixels in each column, and D' is the average concavity.

4. The in-situ estimation method for porosity distribution of rock and soil bodies based on borehole endoscopic images according to claim 3 is characterized in that: Based on the porosity calculation results of the rock and soil mass and the orientation and depth information in the pore wall stereo image, the porosity distribution characteristics are statistically drawn and analyzed.

5. An in-situ measurement device for rock and soil porosity, characterized in that: include: An image acquisition module, used to acquire a stereoscopic image of the borehole wall of the rock and soil body; a data analysis module for obtaining depth and orientation information from the hole wall stereo image using the hole wall stereo image; The calculation module is used to calculate the concavity of any pixel point on the hole wall according to the depth and orientation information, analyze the concavity distribution characteristics, determine the calculation range of rock mass porosity, and calculate the porosity.

6. The in-situ measurement device for rock and soil porosity according to claim 5, characterized in that: The calculation module is also used to calculate the depth porosity or position porosity and output the distribution characteristics of the porosity.

7. An in-situ measurement device for porosity of rock and soil, characterized in that: include: Memory, for storing computer programs; A processor is used to implement the steps of the in-situ estimation method of the porosity distribution of rock and soil bodies based on borehole endoscopic images as described in any one of claims 1 to 4 when executing the computer program.

8. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, which, when executed by a processor, implements the steps of the in-situ estimation method for porosity distribution of rock and soil bodies based on borehole endoscopic images as described in any one of claims 1 to 4.