CN117192628A - Deep fracture water-bearing stratum distribution identification method - Google Patents
Deep fracture water-bearing stratum distribution identification method Download PDFInfo
- Publication number
- CN117192628A CN117192628A CN202311451190.4A CN202311451190A CN117192628A CN 117192628 A CN117192628 A CN 117192628A CN 202311451190 A CN202311451190 A CN 202311451190A CN 117192628 A CN117192628 A CN 117192628A
- Authority
- CN
- China
- Prior art keywords
- data
- stratum
- resistivity
- inversion
- deep
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 68
- 238000000034 method Methods 0.000 title claims abstract description 50
- 238000009826 distribution Methods 0.000 title claims abstract description 47
- 238000012360 testing method Methods 0.000 claims abstract description 31
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 30
- 238000010586 diagram Methods 0.000 claims abstract description 25
- 239000011148 porous material Substances 0.000 claims abstract description 25
- 238000005553 drilling Methods 0.000 claims abstract description 7
- 238000001514 detection method Methods 0.000 claims description 11
- 238000005192 partition Methods 0.000 claims description 7
- 239000000523 sample Substances 0.000 claims description 7
- 238000003325 tomography Methods 0.000 claims description 5
- 238000005516 engineering process Methods 0.000 claims description 4
- 238000003384 imaging method Methods 0.000 claims description 4
- 238000010276 construction Methods 0.000 abstract description 14
- 238000005755 formation reaction Methods 0.000 description 18
- 239000000463 material Substances 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000007405 data analysis Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The application discloses a deep fracture water-bearing stratum distribution identification method, which comprises the following steps: drilling a plurality of object exploratory holes; respectively lowering corresponding resistivity testing devices along the geophysical prospecting holes to obtain stratum data of a region to be tested; constructing an inter-pore stratum resistivity cloud chart based on two groups of stratum data of any two object exploration holes; acquiring two sets of inversion stratum data based on the inter-pore stratum resistivity cloud picture; obtaining two change rule curves based on the two sets of inversion stratum data; obtaining a resistivity baseline value based on the two sets of inversion formation data; constructing a change rule diagram based on the resistivity datum line value and the two change rule curves; acquiring distribution information of deep broken water-containing stratum based on a change rule diagram; the method can accurately identify the distribution condition of the deep broken water-containing stratum, can intuitively display the variation trend and rule of the deep broken water-containing stratum, and provides reliable guiding basis for the safe construction of deep underground engineering.
Description
Technical Field
The disclosure relates to the technical field of deep stratum detection, in particular to a deep fracture water-bearing stratum distribution identification method.
Background
The water inrush and water inrush accidents occur in the process of deep engineering construction, and research and engineering practice show that water damage becomes the primary risk hidden trouble of deep engineering construction, so that the accuracy and reliability of engineering geology and hydrogeology detection in a deep engineering area are improved, and the method is a prerequisite for scientifically evaluating the water inrush risk and effectively treating the water damage.
At present, most of deep projects observe water conditions by means of traditional water exploration investigation holes during construction, only water outlet conditions can be observed, acquired information is limited hydrologic information, the water conditions of stratum in front of the whole working face cannot be comprehensively estimated, and occurrence states of deep stratum cannot be quantitatively represented; the whole view of the whole stratum can not be obtained, so that the water burst risk is difficult to evaluate accurately, the water burst of the stratum in front of the working face can not be found and treated in time, the risk of safety accidents of engineering is greatly increased, and great threat is caused to the safety of staff.
Disclosure of Invention
In view of the above, the embodiments of the present disclosure provide a method for identifying distribution of deep broken water-containing formations, which at least partially solves the problem in the prior art that accurate distribution of deep broken water-containing formations cannot be obtained.
The embodiment of the disclosure provides a deep fracture water-bearing stratum distribution identification method, which comprises the following steps:
drilling Q object exploratory holes;
respectively lowering corresponding resistivity testing devices along the Q geophysical prospecting holes to obtain Q groups of stratum data of the region to be tested; each group of stratum data comprises a plurality of original resistivity data and a plurality of electrode coordinate data, and each original resistivity data corresponds to each electrode coordinate data one by one;
constructing an inter-pore stratum resistivity cloud chart based on two groups of stratum data of any two object exploration holes;
acquiring two sets of inversion stratum data based on the inter-pore stratum resistivity cloud picture; each group of inversion stratum data comprises a plurality of inversion resistivity data and a plurality of inversion coordinate point data, and the inversion resistivity data corresponds to the inversion coordinate point data one by one;
obtaining two change rule curves based on the two groups of inversion stratum data;
obtaining a resistivity datum value based on the two sets of inversion stratum data;
constructing a change rule diagram based on the resistivity datum line value and the two change rule curves;
acquiring distribution information of deep broken water-containing stratum based on the change rule diagram;
wherein Q is more than or equal to 2.
Optionally, the Q geophysical prospecting holes are drilled in parallel at equal intervals and equal depths, and the Q geophysical prospecting holes are all positioned in front of the deep engineering face;
the ratio of the depth of the object detection hole to the distance between the adjacent holes is 2.
Optionally, the resistivity testing device is a resistivity testing special cable comprising a preset number of electrodes;
the front end of the special cable for resistivity test is provided with a counter weight with preset weight.
Optionally, the constructing an inter-pore formation resistivity cloud chart based on two sets of formation data of any two of the object-exploration holes includes:
removing noise data in the two groups of stratum data to obtain a first data set; the first data set comprises a plurality of original resistivity data after noise data are removed and corresponding electrode coordinate data;
and carrying out distributed iterative inversion imaging on the first data set by adopting a tomography method to obtain the inter-pore stratum resistivity cloud picture.
Optionally, the acquiring two sets of inversion formation data based on the inter-pore formation resistivity cloud map includes:
dividing the inter-pore stratum resistivity cloud picture into a plurality of square grids with preset sizes;
and respectively extracting inversion stratum data of the two groups of object exploration holes according to the scale of the preset resistivity cloud map scale.
Optionally, the number of inverted resistivity data is consistent with the number of raw resistivity data.
Optionally, the obtaining a resistivity datum value based on the two sets of inversion stratum data includes:
respectively acquiring two reference values based on the two groups of inversion stratum data;
the average value of the two reference values is the resistivity reference value;
the reference value is，/>The method comprises the steps of carrying out a first treatment on the surface of the Wherein (1)>Is water content coefficient->Is->The number of measuring points in the individual probe holes, +.>Numbering test points->Is->Inversion resistivity data in the individual probe holes, +.>And the inversion coordinate point data corresponding to the inversion resistivity data.
Optionally, the method further comprises: judging whether the object detection hole contains water or not by adopting an ultrasonic detection technology;
if so, the first and second data are not identical,1 is shown in the specification;
if not, the method comprises the steps of,is-1.
Optionally, the obtaining deep broken water-containing stratum distribution information based on the change rule diagram includes:
based on the change rule diagram, a first curve area and a second curve area, of which the change rule curves are lower than the resistivity datum line value, are obtained;
and acquiring an intersection area of the first curve area and the second curve area, namely a distribution area of the deep broken water-containing stratum.
Optionally, the method further comprises: acquiring all test point coordinates based on the intersection area;
constructing a stratum partition fracture distribution map based on the test point coordinates;
and the abscissa of the stratum partition fracture distribution map is a preset span, and the ordinate is sounding.
According to the deep broken water-containing stratum distribution identification method disclosed by the application, a plurality of groups of stratum data are obtained through the drilling and resistivity testing device, and the inter-pore stratum resistivity cloud image and inversion stratum data are constructed, so that the distribution situation of the deep broken water-containing stratum can be accurately identified; according to the method, a plurality of groups of original resistivity data and electrode coordinate data are used as a basis, and resistivity information of a deep stratum is obtained through inversion stratum data and inversion coordinate point data, so that multidimensional analysis of a water-containing stratum is realized; based on a preset formula and inversion stratum data, the method can obtain a resistivity datum value, and a change rule diagram is constructed by combining two groups of inversion stratum data, so that the change trend and rule of the deep broken water-containing stratum can be intuitively displayed; the distribution information of the deep broken water-containing stratum can be obtained by analyzing the change rule diagram, which has important significance in the fields of engineering construction, geological exploration and the like, and can provide scientific basis and reference for related decisions.
The foregoing description is only an overview of the disclosed technology, and may be implemented in accordance with the disclosure of the present disclosure, so that the above-mentioned and other objects, features and advantages of the present disclosure can be more clearly understood, and the following detailed description of the preferred embodiments is given with reference to the accompanying drawings.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present disclosure, the drawings that are needed in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present disclosure, and other drawings may be obtained according to these drawings without inventive effort to a person of ordinary skill in the art.
FIG. 1 is a schematic flow chart of an embodiment of the present disclosure.
Fig. 2 is a schematic diagram of drilling two object holes according to an embodiment of the present disclosure.
FIG. 3 is a schematic flow chart of the method for constructing the resistivity cloud of the stratum between holes in FIG. 1.
FIG. 4 is a schematic diagram of an inter-pore formation resistivity cloud according to an embodiment of the present disclosure.
Fig. 5 is a flow chart of the method for obtaining the resistivity reference line value in fig. 1.
Fig. 6 is a schematic diagram of a change rule chart in an embodiment of the disclosure.
FIG. 7 is a flow chart of a method for obtaining information about the distribution of a deep fractured water-bearing formation of FIG. 1.
FIG. 8 is a schematic representation of the fracture water bearing formation profile information of FIG. 4.
FIG. 9 is a schematic representation of a zone fracture profile of a subterranean formation in accordance with one embodiment of the present disclosure.
Detailed Description
The present disclosure is described in further detail below with reference to the drawings and the embodiments. It is to be understood that the specific embodiments described herein are merely illustrative of the relevant content and not limiting of the present disclosure. It should be further noted that, for convenience of description, only a portion relevant to the present disclosure is shown in the drawings.
In addition, embodiments of the present disclosure and features of the embodiments may be combined with each other without conflict. The technical aspects of the present disclosure will be described in detail below with reference to the accompanying drawings in conjunction with embodiments.
Unless otherwise indicated, the exemplary implementations/embodiments shown are to be understood as providing exemplary features of various details of some ways in which the technical concepts of the present disclosure may be practiced. Thus, unless otherwise indicated, features of the various implementations/embodiments may be additionally combined, separated, interchanged, and/or rearranged without departing from the technical concepts of the present disclosure.
The use of cross-hatching and/or shading in the drawings is typically used to clarify the boundaries between adjacent components. As such, the presence or absence of cross-hatching or shading does not convey or represent any preference or requirement for a particular material, material property, dimension, proportion, commonality between illustrated components, and/or any other characteristic, attribute, property, etc. of a component, unless indicated. In addition, in the drawings, the size and relative sizes of elements may be exaggerated for clarity and/or descriptive purposes. While the exemplary embodiments may be variously implemented, the specific process sequences may be performed in a different order than that described. For example, two consecutively described processes may be performed substantially simultaneously or in reverse order from that described. Moreover, like reference numerals designate like parts.
When an element is referred to as being "on" or "over", "connected to" or "coupled to" another element, it can be directly on, connected or coupled to the other element or intervening elements may be present. However, when an element is referred to as being "directly on," "directly connected to," or "directly coupled to" another element, there are no intervening elements present. For this reason, the term "connected" may refer to physical connections, electrical connections, and the like, with or without intermediate components.
For descriptive purposes, the present disclosure may use spatially relative terms such as "under … …," under … …, "" under … …, "" lower, "" above … …, "" upper, "" above … …, "" higher "and" side (e.g., as in "sidewall"), etc., to describe one component's relationship to another (other) component as illustrated in the figures. In addition to the orientations depicted in the drawings, the spatially relative terms are intended to encompass different orientations of the device in use, operation, and/or manufacture. For example, if the device in the figures is turned over, elements described as "under" or "beneath" other elements or features would then be oriented "over" the other elements or features. Thus, the exemplary term "below" … … can encompass both an orientation of "above" and "below". Furthermore, the device may be otherwise positioned (e.g., rotated 90 degrees or at other orientations) and the spatially relative descriptors used herein interpreted accordingly.
The terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting. As used herein, the singular forms "a", "an" and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise. Furthermore, when the terms "comprises" and/or "comprising," and variations thereof, are used in the present specification, the presence of stated features, integers, steps, operations, elements, components, and/or groups thereof is described, but the presence or addition of one or more other features, integers, steps, operations, elements, components, and/or groups thereof is not precluded. It is also noted that, as used herein, the terms "substantially," "about," and other similar terms are used as approximation terms and not as degree terms, and as such, are used to explain the inherent deviations of measured, calculated, and/or provided values that would be recognized by one of ordinary skill in the art.
Referring to fig. 1, the present application provides a deep fracture water-bearing formation profile identification method comprising the steps of:
s100, drilling Q object exploratory holes; wherein Q is more than or equal to 2.
S200, respectively lowering corresponding resistivity testing devices along Q geophysical prospecting holes to obtain Q groups of stratum data of a region to be tested;
each group of stratum data comprises a plurality of original resistivity data and a plurality of electrode coordinate data, and each original resistivity data corresponds to each electrode coordinate data one by one.
S300, constructing an inter-pore stratum resistivity cloud chart based on two groups of stratum data of any two object exploration holes.
S400, acquiring two sets of inversion stratum data based on an inter-pore stratum resistivity cloud picture; each group of inversion stratum data comprises a plurality of inversion resistivity data and a plurality of inversion coordinate point data, and the inversion resistivity data corresponds to the inversion coordinate point data one by one.
In this embodiment, the number of inversion resistivity data is identical to the number of raw resistivity data.
S500, obtaining two change rule curves based on two sets of inversion stratum data;
obtaining a resistivity baseline value based on the two sets of inversion formation data;
constructing a change rule diagram based on the resistivity datum line value and the two change rule curves;
s600, obtaining deep broken water-containing stratum distribution information based on the change rule diagram.
According to the deep broken water-containing stratum distribution identification method disclosed by the application, a plurality of groups of stratum data are obtained through the drilling and resistivity testing device, and the inter-pore stratum resistivity cloud image and inversion stratum data are constructed, so that the distribution situation of the deep broken water-containing stratum can be accurately identified; according to the method, a plurality of groups of original resistivity data and electrode coordinate data are used as a basis, and resistivity information of a deep stratum is obtained through inversion stratum data and inversion coordinate point data, so that multidimensional analysis of a water-containing stratum is realized; based on a preset formula and inversion stratum data, the method can obtain a resistivity datum value, and a change rule diagram is constructed by combining two groups of inversion stratum data, so that the change trend and rule of the deep broken water-containing stratum can be intuitively displayed; the distribution information of the deep broken water-containing stratum can be obtained by analyzing the change rule diagram, which has important significance in the fields of engineering construction, geological exploration and the like, and can provide scientific basis and reference for related decisions.
In summary, the application can identify the occurrence state of Fang Deceng in front of the face of the deep engineering, can efficiently identify the distribution situation of the deep broken water-containing stratum, provides reliable basis for the safe construction of the deep underground engineering, provides important reference for decisions in the related fields, and has higher practicability and application value.
In this embodiment, two holes are drilled as an example.
Referring to fig. 2, two geophysical prospecting holes are drilled at equal intervals and equal depths in parallel, and the two geophysical prospecting holes are both positioned in front of the deep engineering face.
In this embodiment, the resistivity testing device is preferably a resistivity testing dedicated cable containing a predetermined number of electrodes; in addition, the front end of the special cable for resistivity test is provided with a counterweight hammer with preset weight, so that the special cable can be put down in place, and effective test data can be obtained.
Preferably, the weight is cylindrical, its diameter is not greater than the external diameter of the specific resistance test cable, and its weight is preferably 5kg-10kg.
And each object exploratory hole is internally provided with a specific resistance test special cable, and 48 electrodes are arranged on each specific resistance test special cable at equal intervals so as to acquire enough original specific resistance data.
In this embodiment, the signal collection, arrangement and storage of the specific resistance test dedicated cable are performed by the collector.
In this embodiment, the ratio of the depth of the object-detecting hole to the distance between adjacent holes is preferably 2; here, the adjacent hole pitch refers to the pitch between any two selected holes.
In this embodiment, the drill deflection of the borehole is no greater than 2% and the bore diameter is greater than the outer diameter of the resistivity test cable.
Referring to fig. 3 and 4, construction of an inter-pore formation resistivity cloud chart is performed based on two sets of formation data acquired in a geophysical prospecting hole 1 and a geophysical prospecting hole 2, and the construction method specifically includes:
s310, eliminating noise data in two groups of stratum data to obtain a first data set; the first data set comprises a plurality of original resistivity data after noise data are removed and corresponding electrode coordinate data;
s320, performing distributed iterative inversion imaging on the first data set by adopting a tomography method to obtain the inter-pore stratum resistivity cloud picture.
In this embodiment, since the electrode coordinate data corresponds to the original resistivity data one by one, only noise data of the original resistivity data in the formation data needs to be removed, only the original resistivity data of the remaining non-noise data and the corresponding electrode coordinate data are left, and the inter-pore formation resistivity cloud image is obtained by using these data as effective input data, and extracted into image quality.
By the construction method of the inter-pore stratum resistivity cloud picture, noise data in stratum data is removed, so that the interference and error of the data can be reduced, the accuracy and reliability of a result are improved, and the real stratum resistivity distribution condition can be better restored; after noise data are removed, the obtained first data set comprises clear original resistivity data and corresponding electrode coordinate data, and the data can be used as input of a tomography method for distributed iterative inversion imaging; the resistivity cloud image of the stratum among the holes can intuitively show the resistivity distribution condition of the stratum at the deep part by a tomography method, which is very helpful for analyzing and understanding the properties and characteristics of the stratum at the deep part and breaking the stratum with water, and provides important references for subsequent geological exploration and engineering construction.
In this embodiment, for two sets of inversion formation data, the specific method includes: dividing an inter-pore stratum resistivity cloud picture into a plurality of square grids with preset sizes; and respectively extracting inversion stratum data of the two groups of object exploratory holes according to the scale of the preset resistivity cloud chart.
Wherein, the preset size may be 1×1.
Referring to fig. 5, in this embodiment, the method for acquiring the resistivity reference line value specifically includes the following steps:
a510, respectively acquiring two reference values based on two groups of inversion stratum data and a preset formula;
and A520, calculating an average value of the two reference values, namely the resistivity reference value.
The preset formula is:。
wherein,for reference value, < >>Is water content coefficient->Is->The number of measuring points in the individual probe holes, +.>Numbering test points->Is->Inversion resistivity data in the individual probe holes, +.>Is inversion coordinate point data corresponding to the inversion resistivity data.
In this embodiment, the deep fracture water-bearing stratum distribution identification method disclosed by the application further comprises the following steps: judging whether the object detection hole contains water or not by adopting an ultrasonic detection technology, if so (namely, the object detection hole contains water), then1, if no (i.e. no water in the object probe hole), then +.>Is-1.
Referring to fig. 6, two obtained change rule curves and a resistivity datum line value are drawn in preset coordinates, and a change rule diagram of the resistivity between two adjacent geophysical prospecting holes and the datum line is obtained, wherein in the change rule diagram, the abscissa is sounding (namely, the position of a test point extracted from an inter-hole stratum resistivity cloud picture), and the ordinate is corresponding inversion resistivity data.
Referring to fig. 7 and 8, the method for obtaining the distribution information of the deep fractured water-containing stratum specifically comprises the following steps:
s610, based on a change rule diagram, acquiring a first curve area and a second curve area of which two change rule curves are lower than a resistivity datum line value;
s620, acquiring an intersection area of the first curve area and the second curve area, namely a distribution area of the deep broken water-containing stratum.
In the embodiment, the area lower than the resistivity datum value can be intuitively determined from the change rule diagram, the distribution quantity and specific distribution positions of deep broken water-containing stratum are obtained, the position to be grouting and the grouting quantity of the corresponding position can be accurately determined through analysis of the distribution quantity and the distribution area, fixed-point quantitative grouting is realized, accurate guidance is provided for actual construction, construction efficiency and grouting safety are improved, and the occurrence probability of water inrush and water inrush accidents in subsequent construction is greatly reduced.
In addition, compared with the traditional data processing method, the method simplifies the complex data analysis process into the curve area for checking the change rule diagram, can save time and labor cost, and reduces the possibility of human errors.
Referring to fig. 9, the deep fracture water-bearing formation distribution identification method disclosed in the present application further includes: acquiring coordinates of all test points based on an intersection area in the change rule diagram;
and drawing the obtained coordinates of a plurality of test points in a cross-section coordinate system to obtain a stratum partition fracture distribution map, wherein in the stratum partition fracture distribution map, the abscissa is a preset span, and the ordinate is sounding, namely, in the application, the distribution area of deep broken water-containing stratum can be reflected by the stratum partition fracture distribution map, so that the method is more visual and concise.
In the description of the present specification, reference to the terms "one embodiment/manner," "some embodiments/manner," "example," "a particular example," "some examples," etc., means that a particular feature, structure, material, or characteristic described in connection with the embodiment/manner or example is included in at least one embodiment/manner or example of the application. In this specification, the schematic representations of the above terms are not necessarily for the same embodiment/manner or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments/modes or examples. Furthermore, the various embodiments/modes or examples described in this specification and the features of the various embodiments/modes or examples can be combined and combined by persons skilled in the art without contradiction.
Furthermore, the terms "first," "second," and the like, are used for descriptive purposes only and are not to be construed as indicating or implying a relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defining "a first" or "a second" may explicitly or implicitly include at least one such feature. In the description of the present application, the meaning of "plurality" means at least two, for example, two, three, etc., unless specifically defined otherwise.
It will be appreciated by those skilled in the art that the above-described embodiments are merely for clarity of illustration of the disclosure, and are not intended to limit the scope of the disclosure. Other variations or modifications will be apparent to persons skilled in the art from the foregoing disclosure, and such variations or modifications are intended to be within the scope of the present disclosure.
Claims (10)
1. A method for identifying a distribution of a deep fracture water-bearing formation, comprising:
drilling Q object exploratory holes;
respectively lowering corresponding resistivity testing devices along the Q geophysical prospecting holes to obtain Q groups of stratum data of the region to be tested; each group of stratum data comprises a plurality of original resistivity data and a plurality of electrode coordinate data, and each original resistivity data corresponds to each electrode coordinate data one by one;
constructing an inter-pore stratum resistivity cloud chart based on two groups of stratum data of any two object exploration holes;
acquiring two sets of inversion stratum data based on the inter-pore stratum resistivity cloud picture; each group of inversion stratum data comprises a plurality of inversion resistivity data and a plurality of inversion coordinate point data, and the inversion resistivity data corresponds to the inversion coordinate point data one by one;
obtaining two change rule curves based on the two groups of inversion stratum data;
obtaining a resistivity datum value based on the two sets of inversion stratum data;
constructing a change rule diagram based on the resistivity datum line value and the two change rule curves;
acquiring distribution information of deep broken water-containing stratum based on the change rule diagram;
wherein Q is more than or equal to 2.
2. The method for identifying the distribution of the deep fractured water-bearing stratum according to claim 1, wherein Q geophysical prospecting holes are drilled at equal intervals and equal depths in parallel, and the Q geophysical prospecting holes are all positioned in front of a deep engineering face;
the ratio of the depth of the object detection hole to the distance between the adjacent holes is 2.
3. The method for identifying the distribution of a deep fractured water-bearing stratum according to claim 2, wherein the resistivity testing device is a resistivity testing special cable comprising a preset number of electrodes;
the front end of the special cable for resistivity test is provided with a counter weight with preset weight.
4. The method of claim 1, wherein the constructing an inter-pore formation resistivity cloud based on two sets of formation data for any two of the physical exploration pores comprises:
removing noise data in the two groups of stratum data to obtain a first data set; the first data set comprises a plurality of original resistivity data after noise data are removed and corresponding electrode coordinate data;
and carrying out distributed iterative inversion imaging on the first data set by adopting a tomography method to obtain the inter-pore stratum resistivity cloud picture.
5. The method of identifying a deep fractured water bearing formation profile of claim 4, wherein the obtaining two sets of inversion formation data based on the inter-pore formation resistivity cloud patterns comprises:
dividing the inter-pore stratum resistivity cloud picture into a plurality of square grids with preset sizes;
and respectively extracting inversion stratum data of the two groups of object exploration holes according to the scale of the preset resistivity cloud map scale.
6. The deep fractured water bearing formation profile identification method of claim 5, wherein the amount of inverted resistivity data corresponds to the amount of raw resistivity data.
7. The method of identifying a deep fractured water bearing formation profile according to claim 5, wherein said obtaining a resistivity baseline value based on two sets of said inverted formation data comprises:
respectively acquiring two reference values based on the two groups of inversion stratum data;
the average value of the two reference values is the resistivity reference value;
the reference value is，/>The method comprises the steps of carrying out a first treatment on the surface of the Wherein (1)>Is water content coefficient->Is->The number of measuring points in the individual probe holes, +.>Numbering test points->Is->Inversion resistivity data in the individual probe holes, +.>And the inversion coordinate point data corresponding to the inversion resistivity data.
8. The method of identifying a distribution of a deep fractured water bearing formation of claim 7, further comprising: judging whether the object detection hole contains water or not by adopting an ultrasonic detection technology;
if so, the first and second data are not identical,1 is shown in the specification;
if not, the method comprises the steps of,is-1.
9. The method for identifying the distribution of the deep fractured water-containing stratum according to claim 8, wherein the obtaining the distribution information of the deep fractured water-containing stratum based on the change rule diagram comprises:
based on the change rule diagram, a first curve area and a second curve area, of which the change rule curves are lower than the resistivity datum line value, are obtained;
and acquiring an intersection area of the first curve area and the second curve area, namely a distribution area of the deep broken water-containing stratum.
10. The method of identifying a deep fractured water bearing formation profile of claim 9, further comprising: acquiring all test point coordinates based on the intersection area;
constructing a stratum partition fracture distribution map based on the test point coordinates;
and the abscissa of the stratum partition fracture distribution map is a preset span, and the ordinate is sounding.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311451190.4A CN117192628B (en) | 2023-11-03 | 2023-11-03 | Deep fracture water-bearing stratum distribution identification method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311451190.4A CN117192628B (en) | 2023-11-03 | 2023-11-03 | Deep fracture water-bearing stratum distribution identification method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117192628A true CN117192628A (en) | 2023-12-08 |
| CN117192628B CN117192628B (en) | 2024-01-26 |
Family
ID=89000152
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311451190.4A Active CN117192628B (en) | 2023-11-03 | 2023-11-03 | Deep fracture water-bearing stratum distribution identification method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117192628B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117872494A (en) * | 2024-03-12 | 2024-04-12 | 北京科技大学 | Deep metal ore granite rock internal joint information acquisition method and prediction method |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2002073240A1 (en) * | 2001-03-13 | 2002-09-19 | Conoco Phillips Company | Method and process for prediction of subsurface fluid and rock pressures in the earth |
| KR20100007352A (en) * | 2008-07-14 | 2010-01-22 | (주)건설기술원 | A resistivity survey system for the saturation variation long period monitoring of the soft ground |
| CN104730585A (en) * | 2015-03-26 | 2015-06-24 | 山东科技大学 | Method for monitoring destroying depth of mining working face base plate in real time |
| CN110018519A (en) * | 2019-02-01 | 2019-07-16 | 福建省厦门环境监测中心站 | A kind of underground harmful engineering ground and hydrogeological survey system and method |
| CN112433253A (en) * | 2020-10-14 | 2021-03-02 | 高军 | Three-dimensional detection method for surrounding rock tunnel with weak broken flowing water |
| CN113156518A (en) * | 2021-04-20 | 2021-07-23 | 中国矿业大学 | Real-time advanced detection method for vector resistivity of water-containing disaster body |
| WO2021254394A1 (en) * | 2020-06-18 | 2021-12-23 | 山东大学 | Three-dimensional imaging method and system for earth surface comprehensive geophysical prospecting |
| CN113960695A (en) * | 2021-09-18 | 2022-01-21 | 山东大学 | Fine exploration method for water-rich karst in complex urban environment |
| CN114706124A (en) * | 2022-03-07 | 2022-07-05 | 湖南继善高科技有限公司 | Fluid identification method and system based on combination of electromagnetism and earthquake |
| CN115639612A (en) * | 2022-07-06 | 2023-01-24 | 西安石油大学 | Gas-water distribution analysis method and device for tight sandstone |
| CN116773099A (en) * | 2023-06-08 | 2023-09-19 | 河海大学 | Method for precisely positioning dam leakage channel by combining chromatography and tracing under complex conditions |
| CN116859478A (en) * | 2023-09-05 | 2023-10-10 | 中国地质大学（武汉） | Groundwater simulation method and system based on transient electromagnetic imaging |
-
2023
- 2023-11-03 CN CN202311451190.4A patent/CN117192628B/en active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2002073240A1 (en) * | 2001-03-13 | 2002-09-19 | Conoco Phillips Company | Method and process for prediction of subsurface fluid and rock pressures in the earth |
| KR20100007352A (en) * | 2008-07-14 | 2010-01-22 | (주)건설기술원 | A resistivity survey system for the saturation variation long period monitoring of the soft ground |
| CN104730585A (en) * | 2015-03-26 | 2015-06-24 | 山东科技大学 | Method for monitoring destroying depth of mining working face base plate in real time |
| CN110018519A (en) * | 2019-02-01 | 2019-07-16 | 福建省厦门环境监测中心站 | A kind of underground harmful engineering ground and hydrogeological survey system and method |
| WO2021254394A1 (en) * | 2020-06-18 | 2021-12-23 | 山东大学 | Three-dimensional imaging method and system for earth surface comprehensive geophysical prospecting |
| CN112433253A (en) * | 2020-10-14 | 2021-03-02 | 高军 | Three-dimensional detection method for surrounding rock tunnel with weak broken flowing water |
| CN113156518A (en) * | 2021-04-20 | 2021-07-23 | 中国矿业大学 | Real-time advanced detection method for vector resistivity of water-containing disaster body |
| CN113960695A (en) * | 2021-09-18 | 2022-01-21 | 山东大学 | Fine exploration method for water-rich karst in complex urban environment |
| CN114706124A (en) * | 2022-03-07 | 2022-07-05 | 湖南继善高科技有限公司 | Fluid identification method and system based on combination of electromagnetism and earthquake |
| CN115639612A (en) * | 2022-07-06 | 2023-01-24 | 西安石油大学 | Gas-water distribution analysis method and device for tight sandstone |
| CN116773099A (en) * | 2023-06-08 | 2023-09-19 | 河海大学 | Method for precisely positioning dam leakage channel by combining chromatography and tracing under complex conditions |
| CN116859478A (en) * | 2023-09-05 | 2023-10-10 | 中国地质大学（武汉） | Groundwater simulation method and system based on transient electromagnetic imaging |
Non-Patent Citations (2)
| Title |
|---|
| ODEYEMI EMMANUEL OLORUNTOBA等: "Geophysical Investigation of Aquifer Layer of Landmark University Community, Omu-Aran, Nigeria using Electrical Resistivity Method", 《JOURNAL OF ENGINEERING AND APPLIED SCIENCES》, vol. 14, no. 12, pages 4285 - 4289 * |
| 王嘉文等: "应用测井参数预测含水层渗透系数模型的研究", 《工程地球物理学报》, vol. 20, no. 4, pages 564 - 571 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117872494A (en) * | 2024-03-12 | 2024-04-12 | 北京科技大学 | Deep metal ore granite rock internal joint information acquisition method and prediction method |
| CN117872494B (en) * | 2024-03-12 | 2024-06-07 | 北京科技大学 | Deep metal ore granite rock internal joint information acquisition method and prediction method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117192628B (en) | 2024-01-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN117192628B (en) | Deep fracture water-bearing stratum distribution identification method | |
| US20090260883A1 (en) | Continuous measurement of heterogeneity of geomaterials | |
| CN103306671B (en) | Four-quadrant reservoir type identification method and system | |
| US20140188447A1 (en) | Estimating fracture dimensions from microseismic data | |
| US20180342073A1 (en) | Methods and systems for automatically analyzing an image representative of a formation | |
| WO2012085848A1 (en) | System and method for reconstructing microseismic event statistics from detection limited data | |
| US20220082517A1 (en) | Methods and systems for characterizing a porous rock sample employing combined capillary pressure and nmr measurements | |
| CN106706715A (en) | Polluted soil detection method based on three-dimensional high-density electrical resistivity method | |
| CN107861162B (en) | Microelectrode logging data-based natural crack identification method and system | |
| CN112363242A (en) | Reservoir fluid identification method and device based on logging fusion | |
| US20150023564A1 (en) | System and method for estimating porosity distribution in subterranean reservoirs | |
| CN107064454B (en) | A kind of rock core of reservoir fissure development degree determines method | |
| US20170363531A1 (en) | System and method for estimating porosity distribution in subterranean reservoirs | |
| CN111963162A (en) | Dual-parameter fluid property identification method based on lithology and water-based property | |
| US9921329B2 (en) | Automated method for selecting positions within formations from which to extract samples thereof | |
| CN110362840B (en) | Crack evaluation method and system for micro-spherical focused logging | |
| CN108457646B (en) | Method for determining reservoir fluid properties | |
| US10295691B2 (en) | Method and apparatus for processing log data | |
| CN113075748A (en) | Fracture effectiveness evaluation method based on imaging logging and acoustic remote detection logging data | |
| WO2021108444A1 (en) | Discrimination between subsurface formation natural fractures and stress induced tensile fractures based on borehole images | |
| CN117235428B (en) | Method for judging key seismic source in coal rock damage process | |
| CN117272213B (en) | Method, device, equipment and medium for sweeping underground pollutant by using ground physical comprehensive parameters | |
| CN118091786B (en) | Method and device for correcting linear density of cracks based on underground logging data | |
| CN114814979B (en) | Method, system, electronic equipment and medium for predicting deep mine goaf | |
| CN114185101B (en) | Uranium deposit ore finding method, device and equipment |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |