CN111523190A - Method for evaluating sealing performance between chambers of underground water sealed cave depot - Google Patents
Method for evaluating sealing performance between chambers of underground water sealed cave depot Download PDFInfo
- Publication number
- CN111523190A CN111523190A CN202010379730.2A CN202010379730A CN111523190A CN 111523190 A CN111523190 A CN 111523190A CN 202010379730 A CN202010379730 A CN 202010379730A CN 111523190 A CN111523190 A CN 111523190A
- Authority
- CN
- China
- Prior art keywords
- connectivity
- matrix
- fracture
- fracture network
- sealing performance
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computer Networks & Wireless Communication (AREA)
- Examining Or Testing Airtightness (AREA)
Abstract
The invention provides a method for evaluating the sealing performance between chambers of an underground water sealed cave depot, which comprises the following steps: establishing a discrete fracture network model according to fracture parameters between two adjacent caverns, and acquiring node information of the fracture network model, wherein the fracture parameters mainly comprise fracture length, inclination, dip angle and density; determining connectivity of the fracture network model from a starting boundary to an end boundary by combining a graph theory method according to the node information of the fracture network model; and obtaining the sealing evaluation result between the adjacent chambers according to the connectivity judgment result, wherein the connectivity and the sealing are inversely related. The invention can evaluate the sealing performance between adjacent caverns and ensure that oil leakage between the caverns cannot occur due to communication. The invention can also judge the connectivity between the caverns and adopt different measures to improve the sealing property according to the connectivity result.
Description
Technical Field
The invention relates to the field of civil engineering, in particular to a method for evaluating the sealing performance between underground water sealed caverns.
Background
According to the socioeconomic development and international energy situation analysis of China, the country has already stored the petroleum strategic reserve as the basic task of the national economy and social development. The underground water-sealed rock cavern oil depot has become a preferred mode of strategic oil storage due to the advantages of high safety, small occupied area, low operation cost, environmental friendliness and the like. The key point for evaluating the oil storage capacity of the underground water-sealed rock cavern oil storage is to evaluate the sealing performance of the cavern. The sealing performance between the cavities is particularly important for the condition that different oil products are stored between the adjacent cavities. The sealing performance between the chambers is good, so that the phenomenon of oil leakage between different oil products in adjacent chambers can be avoided. The tightness of the caverns is related to the cavern distance, the liquid level difference of oil products between adjacent caverns, the water curtain pressure and the fracture parameters of rock mass. Therefore, it is necessary to establish a method for evaluating the sealing performance between the caverns of the underground water-sealed cavern to evaluate the sealing performance between the caverns, and adopt a distinguishing measure to improve the sealing performance between the caverns according to different evaluation results.
Disclosure of Invention
According to the problem that an effective cavern sealing evaluation means is lacked in the prior art, the method for evaluating the sealing performance between the caverns of the underground water sealed cavern is provided, and theoretical analysis is carried out according to fracture parameters of cavern intervals to obtain the sealing performance effect between the adjacent caverns.
The technical means adopted by the invention are as follows:
a method for evaluating the sealing performance between chambers of an underground water sealed cave depot comprises the following steps:
s1, generating a discrete fracture network model by adopting a Monte Carlo simulation technical means according to fracture parameters between two adjacent caverns, and acquiring node information of the fracture network model, wherein the fracture parameters mainly comprise fracture length, inclination angle and density;
s2, determining connectivity of the fracture network model from a starting boundary to an end boundary by combining a graph theory method according to the node information of the fracture network model;
and S3, obtaining a sealing evaluation result between adjacent chambers according to the connectivity judgment result, wherein the connectivity and the sealing are inversely related.
Further, the node information of the fracture network model in step S1 includes:
the position relation of the nodes and the water heads of the piezometers at the nodes;
the pressure measuring pipe water head at the node is the sum of the position water head and the pressure water head at the node.
Further, step S2 specifically includes:
s21, obtaining an adjacency matrix A according to the fracture network model node information, wherein:
in the formula, AijIs an element of the matrix A, Hi、HjThe pressure measuring pipe water heads are respectively nodes i and j;
s22, determining a connectivity characterization matrix C from the adjacency matrix A based on the first transformation matrix W and the second transformation matrix D, wherein,
the first transformation matrix W is obtained from the following calculation:
in the formula, WijAre the elements of the matrix W which,
the second transformation matrix D is obtained from the following calculation:
in the formula, DijIs an element of the matrix D which,
the connectivity characterization matrix C is obtained according to the following calculation:
in the formula, CmnCharacterizing elements of a matrix C for connectivity, m representing nodes on a starting boundary and n representing nodes on an ending boundary;
s23, determining connectivity of the fracture model between the caverns according to the obtained connectivity characterization matrix C, specifically, if all elements in the connectivity characterization matrix C are 0, indicating that the fracture network between the caverns is not communicated, otherwise, indicating that the fracture network model between the caverns is communicated.
Further, the method further comprises:
and S4, improving the sealing performance between the chambers in a distinguishing mode according to the sealing performance evaluation result.
Further, step S4 specifically includes: calculating a communication parameter p according to the connectivity characterization matrix C, and selecting different tightness improving means according to the communication parameter p;
wherein the connected parameter p is obtained according to the following calculation:
in the formula, CmnAnd (3) representing the elements of the connectivity characterization matrix C, wherein M is the number of nodes on the initial boundary of the fracture network, and N is the number of nodes on the terminal boundary of the fracture network.
Further, selecting different tightness improving means according to the communication parameter p comprises:
when the communication parameter p belongs to the first interval, the sealing performance is improved by improving the water pressure of the water curtain;
when the communication parameter p belongs to a second interval, improving the sealing performance through grouting seepage prevention;
and the first interval and the second interval have no intersection.
Compared with the prior art, the invention has the following advantages:
the method combines a fracture network method and a graph theory method, can obtain the connectivity of the fracture network between the caverns, further obtain the sealing performance between the caverns, and adopt different measures to ensure the sealing performance between the caverns according to the sealing performance effect. The method theoretically analyzes the sealing effect between adjacent chambers according to the crack parameters of the cavern space, and has strong adaptability.
For the above reasons, the present invention can be widely applied to the field of civil engineering.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed to be used in the description of the embodiments or the prior art will be briefly introduced below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to these drawings without creative efforts.
FIG. 1 is a flowchart of a method for evaluating sealability according to the present invention.
FIG. 2 is a schematic diagram of a fracture network model in an embodiment.
Detailed Description
In order to make the technical solutions of the present invention better understood, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
As shown in fig. 1, the invention provides a method for evaluating the sealing performance between chambers of an underground water-sealed cave depot, which comprises the following steps:
s1, generating a discrete fracture network model by adopting technical means such as Monte Carlo simulation and the like according to fracture parameters between two adjacent caverns, and acquiring node information of the fracture network model, wherein the fracture parameters mainly comprise fracture length, inclination angle and density. The fracture network can be generated by adopting different means, the formed model is a random fracture system, the fracture parameters are different, and the generated models are different.
The node information of the fracture network model comprises the position relation of nodes and the piezometric tube water heads at the nodes, and the piezometric tube water heads at the nodes are the sum of the position water heads and the pressure water heads at the nodes. The node position water head and the pressure water head are mainly obtained by calculating a flow control equation (such as a cubic law) and boundary conditions, and the calculation is needed and can be carried out by adopting commercial software or other means.
S2, determining connectivity of the fracture network model from a starting boundary to an end boundary according to the node information of the fracture network model by combining a graph theory method, wherein the connectivity specifically comprises the following steps:
s21, obtaining an adjacency matrix A according to the fracture network model node information, wherein:
in the formula, AijIs an element of the matrix A, Hi、HjThe pressure measuring pipe water heads are respectively nodes i and j;
s22, determining a connectivity characterization matrix C from the adjacency matrix A based on the first transformation matrix W and the second transformation matrix D, wherein,
the first transformation matrix W is obtained from the following calculation:
in the formula, WijAre the elements of the matrix W which,
the second transformation matrix D is obtained from the following calculation:
in the formula, DijIs an element of the matrix D which,
the connectivity characterization matrix C is obtained according to the following calculation:
in the formula, CmnIs communicated withThe elements of the characteristic representation matrix C, m represents the nodes on the starting boundary, and n represents the nodes on the terminal boundary;
s23, determining connectivity of the fracture model between the caverns according to the obtained connectivity characterization matrix C, specifically, if all elements in the connectivity characterization matrix C are 0, indicating that the fracture network between the caverns is not communicated, otherwise, indicating that the fracture network model between the caverns is communicated.
And S3, obtaining the sealing evaluation result between the adjacent chambers according to the connectivity judgment result. As a preferred embodiment, the sealing performance between the caverns is judged mainly according to the connectivity of the fracture model between the caverns, and the non-communication of the fracture network model indicates that the sealing performance between the caverns is better, otherwise, the sealing performance between the caverns is poorer.
And S4, improving the sealing performance between the chambers in a distinguishing mode according to the sealing performance evaluation result. The method specifically comprises the following steps:
calculating a communication parameter p according to the connectivity characterization matrix C, and selecting different tightness improving means according to the communication parameter p;
wherein the connected parameter p is obtained according to the following calculation:
in the formula, CmnAnd (3) representing the elements of the connectivity characterization matrix C, wherein M is the number of nodes on the initial boundary of the fracture network, and N is the number of nodes on the terminal boundary of the fracture network.
When the communication parameter p belongs to the first interval, the sealing performance is improved by improving the water pressure of the water curtain; when the communication parameter p belongs to a second interval, improving the sealing performance through grouting seepage prevention; and the first interval and the second interval have no intersection. In a preferred embodiment of the present invention, the first interval is (0, 2% ], and the second interval is (2%, 100% ]).
The method for evaluating the sealing performance between chambers of the underground water-sealed cave depot is described in detail by combining the following examples, and the specific steps are as follows:
a discrete fracture network model is established, in order to simplify calculation, a horizontal fracture exists between chambers, as shown in FIG. 2, 9 nodes exist in the graph, and as can be seen from the graph, a node 1 is adjacent to nodes 2, 8 and 9, a node 2 is adjacent to nodes 1, 3 and 9, a node 3 is adjacent to nodes 2, 4 and 9, a node 4 is adjacent to nodes 3, 5 and 9, a node 5 is adjacent to nodes 4, 6 and 9, a node 6 is adjacent to nodes 5, 7 and 9, a node 7 is adjacent to nodes 6, 8 and 9, a node 8 is adjacent to nodes 1, 7 and 9, a node 3 is adjacent to nodes 1, 2, 3, 4, 5, 6, 7 and 8, and the piezometric tubes of each point are shown in parentheses in FIG. 2 in unit m.
And judging the communication of the discrete fracture network models among the caverns, obtaining a correlation matrix according to the node information, and calculating an adjacency matrix A, a first transformation matrix W, a second transformation matrix D and a connectivity characterization matrix C according to the formulas (1) to (4). And calculating that all elements in the connectivity characterization matrix C are 1, namely the fracture network model between the chambers is connected.
Therefore, the network models of the discrete cracks between the chambers are communicated, so that the sealing performance between the adjacent chambers is poor.
As can be seen from fig. 2, the number of nodes on the starting and ending boundaries is 3, that is, M is equal to N is equal to 3, and the connectivity parameter p is calculated by substituting equation (5) to be 100% > > 2%, that is, grouting is required to reduce connectivity of the fracture network between the caverns and ensure the sealing performance between the caverns.
Finally, it should be noted that: the above embodiments are only used to illustrate the technical solution of the present invention, and not to limit the same; while the invention has been described in detail and with reference to the foregoing embodiments, it will be understood by those skilled in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some or all of the technical features may be equivalently replaced; and the modifications or the substitutions do not make the essence of the corresponding technical solutions depart from the scope of the technical solutions of the embodiments of the present invention.
Claims (6)
1. The method for evaluating the sealing performance between the chambers of the underground water sealed cave depot is characterized by comprising the following steps of:
s1, generating a discrete fracture network model by adopting a Monte Carlo simulation technical means according to fracture parameters between two adjacent caverns, and acquiring node information of the fracture network model, wherein the fracture parameters mainly comprise fracture length, inclination angle and density;
s2, determining connectivity of the fracture network model from a starting boundary to an end boundary by combining a graph theory method according to the node information of the fracture network model;
and S3, obtaining a sealing evaluation result between adjacent chambers according to the connectivity judgment result, wherein the connectivity and the sealing are inversely related.
2. The method for evaluating the tightness between chambers in the underground water seal cavern as claimed in claim 1, wherein the node information of the fracture network model in the step S1 includes:
the position relation of the nodes and the water heads of the piezometers at the nodes;
the pressure measuring pipe water head at the node is the sum of the position water head and the pressure water head at the node, and the position water head and the pressure water head are obtained through calculation according to a flow control equation and boundary conditions.
3. The method for evaluating the sealing performance between the underground water seal cavern and the cavern as claimed in claim 1, wherein the step S2 specifically comprises:
s21, obtaining an adjacency matrix A according to the fracture network model node information, wherein:
in the formula, AijIs an element of the matrix A, Hi、HjThe pressure measuring pipe water heads are respectively nodes i and j;
s22, determining a connectivity characterization matrix C from the adjacency matrix A based on the first transformation matrix W and the second transformation matrix D, wherein,
the first transformation matrix W is obtained from the following calculation:
in the formula, WijAre the elements of the matrix W which,
the second transformation matrix D is obtained from the following calculation:
in the formula, DijIs an element of the matrix D which,
the connectivity characterization matrix C is obtained according to the following calculation:
in the formula, CmnCharacterizing elements of a matrix C for connectivity, m representing nodes on a starting boundary and n representing nodes on an ending boundary;
s23, determining connectivity of the fracture model between the caverns according to the obtained connectivity characterization matrix C, specifically, if all elements in the connectivity characterization matrix C are 0, indicating that the fracture network between the caverns is not communicated, otherwise, indicating that the fracture network model between the caverns is communicated.
4. The method for evaluating the tightness between chambers of the underground water seal cavern as claimed in claim 1, further comprising:
and S4, improving the sealing performance between the chambers in a distinguishing mode according to the sealing performance evaluation result.
5. The method for evaluating the tightness between the underground water seal cavern and the cavern as claimed in claim 4, wherein the step S4 specifically comprises: calculating a communication parameter p according to the connectivity characterization matrix C, and selecting different tightness improving means according to the communication parameter p;
wherein the connected parameter p is obtained according to the following calculation:
in the formula, CmnAnd (3) representing the elements of the connectivity characterization matrix C, wherein M is the number of nodes on the initial boundary of the fracture network, and N is the number of nodes on the terminal boundary of the fracture network.
6. The method for evaluating the tightness between chambers of the underground water seal cavern as claimed in claim 5, wherein the selection of different tightness improving means according to the communication parameter p comprises:
when the communication parameter p belongs to the first interval, the sealing performance is improved by improving the water pressure of the water curtain;
when the communication parameter p belongs to a second interval, improving the sealing performance through grouting seepage prevention;
and the first interval and the second interval have no intersection.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010379730.2A CN111523190B (en) | 2020-05-08 | 2020-05-08 | Method for evaluating tightness between underground water seal cave depot and cavern |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010379730.2A CN111523190B (en) | 2020-05-08 | 2020-05-08 | Method for evaluating tightness between underground water seal cave depot and cavern |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111523190A true CN111523190A (en) | 2020-08-11 |
| CN111523190B CN111523190B (en) | 2023-10-20 |
Family
ID=71908034
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010379730.2A Active CN111523190B (en) | 2020-05-08 | 2020-05-08 | Method for evaluating tightness between underground water seal cave depot and cavern |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111523190B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060062442A1 (en) * | 2004-09-16 | 2006-03-23 | Imaging Therapeutics, Inc. | System and method of predicting future fractures |
| CN106991244A (en) * | 2017-04-13 | 2017-07-28 | 河海大学 | A kind of method of Fracture Networks connectedness and seepage calculation based on graph theory |
| CN109376489A (en) * | 2018-12-09 | 2019-02-22 | 东北大学 | A kind of fine method for arranging in underground storage water curtain hole |
-
2020
- 2020-05-08 CN CN202010379730.2A patent/CN111523190B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060062442A1 (en) * | 2004-09-16 | 2006-03-23 | Imaging Therapeutics, Inc. | System and method of predicting future fractures |
| CN106991244A (en) * | 2017-04-13 | 2017-07-28 | 河海大学 | A kind of method of Fracture Networks connectedness and seepage calculation based on graph theory |
| CN109376489A (en) * | 2018-12-09 | 2019-02-22 | 东北大学 | A kind of fine method for arranging in underground storage water curtain hole |
Non-Patent Citations (5)
| Title |
|---|
| W.SDERSHOWITZ ET AL: "A Discrete Fracture Network Approach for Evaluation of Hydraulic Fracture Stimulation of Naturally Fractured Reservoris", 《44TH U.S ROCK MECHANICS SYMPOSIUM AND 5TH US CANADA ROCK MECHANICS SYMPOSIUM》, pages 10 - 475 * |
| 刘欢;王者超;范凌峰;翟玉斌;赵艺伟;: "某地下水封石油洞库围岩结构面统计规律与水幕试验初步结果分析", 化工与医药工程, no. 02, pages 6 - 12 * |
| 王晋丽;陈喜;黄远洋;张志才;: "岩体裂隙网络随机生成及连通性研究", 水文地质工程地质, no. 02, pages 30 - 35 * |
| 王者超;: "地下油气储存理论与技术", 国际学术动态, no. 02, pages 35 - 37 * |
| 许振浩 等: "地下石油储库节理裂隙岩体水封机理及水封有效性评价研究", 《科技成果》, vol. 40, no. 2, pages 30 - 35 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111523190B (en) | 2023-10-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107563027B (en) | For analyzing the WELL TEST INTERPRETATION MODEL and method of more fracture-cavity units cascaded structures | |
| Diao et al. | Global resilience analysis of water distribution systems | |
| CN102395200B (en) | Node positioning method in wireless sensor network and apparatus thereof | |
| US11391145B2 (en) | Dynamic fracture width calculation method for drilling fluid loss in fractured formation | |
| CN111577256B (en) | Quantitative evaluation method for permeability increasing effect of cross-layer drilling hydraulic punching | |
| CN105484724A (en) | Drilling downhole anomaly monitoring method | |
| CN114091287B (en) | Method for evaluating crack connectivity and optimizing crack parameters based on complex network theory | |
| CN111523190A (en) | Method for evaluating sealing performance between chambers of underground water sealed cave depot | |
| CN109763812B (en) | Hydrate exploitation mechanical stability quantitative analysis method based on mole-coulomb theory | |
| CN108386187A (en) | The method for analyzing reservoir structure based on oil well indicative curve | |
| CN111612908A (en) | Hydraulic engineering inundation range display method, device, equipment and storage medium | |
| Kamal et al. | Use of transient testing in reservoir management | |
| CN114218877B (en) | Fracture-cavity type oil reservoir numerical simulation method and system | |
| CN115419393B (en) | Plate method for evaluating interlayer packing performance of cement annular layer | |
| CN115994500A (en) | Method and system for evaluating dynamic change of fracture conductivity of shale gas well | |
| CN110552682A (en) | Descending analysis method for vertical well yield of carbonate rock triple medium gas reservoir | |
| CN103266870B (en) | A kind of determination method and system of oil-gas reservoir water blocking damage radius | |
| Zhao et al. | Prediction model for spontaneous combustion of coal around boreholes using bedding gas drainage | |
| CN113361039A (en) | Section optimization method and system for sealing gasket of shield tunnel segment joint | |
| CN111127236A (en) | Technical strategy selection method based on oil reservoir water drive recovery ratio constitution decomposition | |
| Diamond et al. | Practical aspects of gas material balance: theory and application | |
| Heitfeld et al. | The problem of water permeability in dam geology. | |
| CN111368446B (en) | Fracture-cavity type medium connectivity assessment method and system based on harmonic injection system | |
| Wang et al. | Connectivity and permeability for rock fracture network in both radial and unidirectional flow configurations | |
| CN112282741B (en) | Target reservoir well test analysis method, computer storage medium and computer 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 |