Detailed Description
Reference will now be made in detail to exemplary embodiments, examples of which are illustrated in the accompanying drawings. When the following description refers to the accompanying drawings, the same numbers in different drawings refer to the same or similar elements, unless otherwise indicated. The implementations described in the following exemplary embodiments do not represent all implementations consistent with embodiments of the application. Rather, they are merely examples of apparatus and methods consistent with aspects of embodiments of the application as detailed in the accompanying claims.
Important noun interpretation:
Reservoirs, porous rock layers that are capable of storing and transporting hydrocarbons (oil and gas) underground, reservoir rocks often have high porosity and permeability that enable the hydrocarbons to flow and aggregate therein, reservoirs being the core object of hydrocarbon exploration and development because they are the primary storage sites for hydrocarbon resources;
Reservoir model-a mathematical model that divides a complex three-dimensional space into a limited number of small cells, i.e. "grids" or "cells", so that a computer can process and analyze these spatial data, including grid points and the locations of the grid points;
Fault refers to the structural phenomenon that rock is broken due to stress and rock blocks on two sides of a broken surface have obvious displacement;
the fault model is used for describing geometrical characteristics and spatial positions of faults and comprises parameters such as strike, dip angle, inclination, break distance and the like of the faults;
attribute models describing physical and chemical attributes of the geologic body, such as porosity, permeability, rock type, temperature, pressure, etc.;
Connectivity in graph theory, connectivity refers to whether a path connection exists between any two vertexes in a graph;
Reservoir connectivity, which is the ability of fluids in a reservoir to flow, which describes the degree of connectivity of pores and cracks in the reservoir, affecting migration, aggregation and recovery of oil and gas;
the communication coefficient is an index for quantifying the connectivity between different communication bodies in the reservoir, reflects the capability of fluid (such as petroleum, natural gas or water) to flow between different areas inside the reservoir, and can be used for evaluating the connectivity inside the reservoir and the fluid transmission efficiency;
Minkowski and algorithm, which is the sum of two Euclidean space point sets, also called as expansion sets of the two spaces, and is visually represented by the union of the area swept by the set A along the marginal continuous motion of the set B for one circle and the set B, and can also be the union of the area swept by the set B along the boundary continuous motion of the set A and the set A;
polygon segmentation algorithm-complex polygons are decomposed into simpler geometric shapes, such as triangles or convex polygons;
The Bentley-Ottmann algorithm, an efficient computational geometry algorithm for finding all the intersections of a set of line segments on a plane, is implemented primarily by a scanline technique, in which a vertical line (called a scanline) is swept across the entire plane, while maintaining a dynamic set of active line segments, and reporting the intersections efficiently as they occur;
searching a spanning tree in the undirected communication graph with the weight so that the sum of the weights of all edges in the spanning tree is minimum;
Dijkstra algorithm, an algorithm for finding the shortest path from a single source point to all other vertices in the weighted graph;
the minimum spanning tree algorithm is used for finding a subset of the edges in the weighted undirected graph, the subset forms a tree, all vertexes in the graph are contained, and the total weight of the edges is minimum.
Under the condition that the reservoir structure has good continuity, geological properties (such as rock types, porosity, permeability and the like) show certain regularity and continuity along with the change of the space position, at the moment, connectivity analysis is carried out on the property matrix of the reservoir, adjacent points in the space can be correspondingly found, namely, the communicating bodies in the reservoir can be divided, and the flow path in the reservoir can be found, however, when the structure such as faults exist in the stratum, the property matrix of the reservoir and the coordinate matrix of the reservoir are not continuous at the fault position. A certain point on the fault plane is divided into two points after the fault, and the two points still have adjacency in the attribute matrix, but are no longer adjacency in the coordinate matrix, so that the adjacency points in the space can not be found correctly only by using the attribute matrix, and the problem of inaccurate connectivity analysis in the fault structure exists.
According to the method for confirming reservoir connectivity, due to the fact that element coordinates on the fault are discontinuous, adjacent elements cannot be found through continuous coordinates, so that elements which are possibly communicated cannot be found, two sides of the fault are marked as the first reservoir and the second reservoir respectively, the fault comprises the first reservoir section and the second reservoir section, the elements on the first reservoir section and the elements on the second reservoir section are screened out, the elements on the first reservoir section are adjacent to the nearest elements, the elements on the second reservoir section are screened out, multiple groups of elements on the fault are marked by the nearest space distance, the multiple groups of elements are the first elements and target second elements corresponding to the first elements respectively, then the element types of the multiple groups of elements are judged, when the element types of the first elements are matched with the element types of the target second elements, a second target communication body where the target second elements are located is determined, the first communication body where the first elements are located is located, the adjacent to the second communication body is located, the adjacent connection body on the second communication body is not calibrated, and the communication structure is analyzed accurately, and the communication problem is not analyzed.
Fig. 1 is a schematic view of a scenario of a method for confirming reservoir connectivity according to an embodiment of the present application, where, as shown in fig. 1, an execution body of the method may be a system for confirming reservoir connectivity, and the system may be a server. The server can be a computer, a notebook, a tablet or a mobile phone and other devices. The implementation manner of the execution body is not particularly limited in this embodiment, as long as the execution body can determine, according to a spatial distance between a first element of a section on a first reservoir and a second element of a section on a second reservoir in a reservoir model, a target second element corresponding to the first element in the second element, where the section on the first reservoir and the section on the second reservoir are the same section, and when the element type of the first element matches the element type of the target second element, it is determined that the second target communication body and the first communication body are the same communication body, where the target second element in the second reservoir is located, and the first communication body is the communication body where the first element in the first reservoir is located.
The following describes in detail the technical solutions of the embodiments of the present application and how the technical solutions of the embodiments of the present application solve the above technical problems with specific embodiments. The following embodiments may be combined with each other, and the same or similar concepts or processes may not be described in detail in some embodiments. Embodiments of the present application will be described below with reference to the accompanying drawings.
Fig. 2 is a flowchart illustrating a method for confirming reservoir connectivity according to an embodiment of the present application. The implementation body of the method may be a server or other servers, and the embodiment is not particularly limited herein, as shown in fig. 2, and the method may include:
S201, determining a target second element corresponding to the first element in the second element according to the space distance between the first element of the section on the first reservoir and the second element of the section on the second reservoir in the reservoir model, wherein the section on the first reservoir and the section on the second reservoir are the same section.
Wherein, the first reservoir may refer to a region of a fault at one side of a fault line, the second reservoir may refer to a region of the fault at the other side of the fault line, the first reservoir and the second reservoir are opposite, in this embodiment of the present application, the first reservoir represents a lifting disc in the fault structure, the second reservoir represents a descending disc in the fault structure, the lifting disc and the descending disc are used to describe reservoir structures at two sides of the fault, the lifting disc may refer to a reservoir above the fault plane, the descending disc may refer to a reservoir below the fault plane, for example, when one side of the fault plane is faced, the other side of the fault plane is observed, the lifting disc is a portion above the fault plane, and the descending disc is a portion below the fault plane.
The section on the first reservoir may refer to the plane in which the fault is located, and the section on the second reservoir is the same section for identifying the structure in the fault.
The first element may refer to an attribute element on the first reservoir section closest to the fault line, and the second element may refer to an attribute element on the second reservoir section closest to the fault line, where the first element and the second element have attribute values, and the geological conditions stored therein may represent sand, stone, or porosity. In the embodiment of the application, the first element refers to the attribute element closest to the fault line on the section of the ascending disc, and the second element refers to the attribute element closest to the fault line on the section of the descending disc, and the specific conditions of sand and stone can be stored by using the first element and the second element, for example, the type of the first element is sand, and the type of the second element is stone.
The spatial distance may refer to a linear distance of two points or objects in a three-dimensional space, and is used to represent a shortest path length, and a measurement method may be selected according to requirements of a use scene, for example, using spatial coordinate calculation, and for example, using longitude and latitude information and earth radius calculation.
The target second element may refer to an element selected from the second elements according to a requirement, and the specific requirement is set according to a use condition. In the embodiment of the application, the target second element is a target second element which is calibrated by taking a certain first element as a reference, calculating the space distance between all the second elements and the first element, selecting the second element corresponding to the smallest distance, and executing the process on all the first elements, wherein the execution sequence is not limited by the application, and can be set according to the size and the classification of the data quantity.
In an embodiment of the present application, determining, according to a spatial distance between a first element of a section on a first reservoir and a second element of a section on a second reservoir in a reservoir model, a target second element corresponding to the first element in the second element includes:
Determining a first element in the first attribute elements according to the distance between the first attribute elements of the fracture surface on the first reservoir and a target fault line on the fracture surface in the reservoir model, wherein the first attribute elements represent elements of a storage element type on the fracture surface of the first reservoir, and the target fault line represents a line intersecting a horizontal plane with the fracture surface;
Determining a second element in the second attribute elements according to the distance between the second attribute elements of the section on the second reservoir in the reservoir model and the target fault line on the section, wherein the second attribute elements represent elements of the storage element type on the section of the second reservoir;
And determining a target second element corresponding to the first element in the second element according to the space distance between the first element and the second element.
The reservoir model may refer to a model storing coordinate information of the whole reservoir structure, may store three-dimensional position coordinates, may store geographic coordinates including longitude, latitude and altitude, and is specifically selected according to requirements of data analysis, and unit selection of the reservoir model depends on a use scene, for example, regional geological investigation, may select kilometers (km) as units, for example, oilfield development, and may select meters (m) as units.
On the basis of a reservoir model, a fault model is established, wherein the fault model can refer to a model for storing a fault structure, and comprises a fault position, a fault surface and a fault trend, the unit of the fault model can be consistent with the unit of the reservoir model or inconsistent with the unit of the reservoir model, and coordinate information of the fault model and the fault model can be converted through the relation of a coordinate system.
An attribute model is built on the basis of a reservoir model, which may refer to a model that stores the attributes of various points in the reservoir structure, including physical, chemical, and fluid properties. In an embodiment of the application, the attribute model stores mainly rock types at various points in the reservoir structure, such as sand at a certain point in the reservoir, and mud at a certain point.
The section may refer to a section of a selected certain fault in a plurality of fault structures, and the section in the embodiment of the present application is obtained by extracting fault information from the fault structures, where the first reservoir and the second reservoir are respectively on two sides of the section.
The target fault line may refer to a fault line selected according to a preset requirement, the preset requirement is determined according to a use condition and a use requirement, for example, a linear fault line is selected, for example, a curve fault line is selected, in the embodiment of the present application, the target fault line is a fault line obtained by intersecting any horizontal plane through fault information in a fault model, the target fault line may have a plurality of pieces, and the target fault line has coordinate information of a plurality of points.
The attribute elements can refer to elements of storage element types in the reservoir model, the element types can refer to type information represented by the elements, the physical properties can be the physical properties or the chemical properties, wherein the first attribute elements can refer to attribute elements of sections on a first reservoir, the second attribute elements can refer to attribute elements of sections on a second reservoir, the first attribute elements and the second attribute elements are elements selected from the attribute elements according to position requirements, the first attribute elements and the second attribute elements can be multiple, and the structure requirements can refer to a certain position.
In the embodiment of the application, the first element is screened out according to the shortest distance by calculating the distance between the first attribute element and the target fault line, and the element which is closest to the fault line and possibly belongs to the fault structure is similarly screened out according to the shortest distance by calculating the distance between the second attribute element and the target fault line, so that the elements which are possibly belonging to the fault structure on the first reservoir section and the second reservoir section are respectively obtained.
In this embodiment of the present application, the target second element is an element with the smallest spatial distance between the second element and the first element.
The element with the smallest space distance between the second elements and the first elements can refer to the second element with the smallest space distance among a plurality of second elements corresponding to a certain first element, the second element is marked as a target second element, and the minimum distance can be a plurality of second elements, and at the moment, any second element corresponding to one of the minimum distances is selected.
In some embodiments, any first element is selected, the distances between the first element and all second elements are calculated, all distance results are compared, and the second element closest to the first element is selected as the target second element. For example, there is a first element M1, and distances r1, r2, and r3 from second elements N1, N2, and N3, where r2 is the smallest, N2 may be calibrated as a target second element corresponding to M1, the minimum distance relationship corresponding to M1 and N2 is stored, and for each first element having a first reservoir, there is a target second element corresponding to the first element, and each set of corresponding results is stored.
S202, when the element type of the first element is matched with the element type of the target second element, determining that the second target communication body and the first communication body are the same communication body, wherein the second target communication body is the communication body where the target second element is located in the second reservoir, and the first communication body is the communication body where the first element is located in the first reservoir.
The element type of the first element may refer to a geological attribute type of the first element, the element type of the target second element may refer to a geological attribute type of the target second element, in this embodiment of the present application, the element type of the first element is matched with the element type of the target second element, and may refer to that the element type of the first element is the same as the element type of the target second element, and when the element type of the first element is different from the element type of the target second element, the first element and the second element are marked as unmatched.
The second target continuum may refer to a continuum in which the target second element is located, for identifying element aggregates whose attributes are all target second elements, and the continuum may refer to element aggregates having the same property, the elements being spatially connected to each other. The first continuum may refer to a continuum in which the first element is located, and is used to identify element aggregates whose attributes are all the first element. In the embodiment of the application, the element attributes of the second target communication body and the first communication body may be consistent or inconsistent, and the attribute depends on the position of the element. For example, a certain first element at the first reservoir is sand, second elements at the second reservoir are sand and stone, the distances between all the second elements and the first elements are compared, the nearest second element is screened out and used as a target second element, and if the target second element and the first element are sand, the target second communication body and the first communication body are marked as the same communication body.
The same communication may refer to an aggregate having the same attribute element for indicating that there is flow capacity within the aggregate, and the same communication a may be indicated as a piece of sand of an area than if the element type of the same communication a is sand. When the element type of the same communicating body B is stone, the communicating body A and the communicating body B cannot flow.
In the embodiment of the application, firstly, the elements of the ascending disc section and the elements of the descending disc section are divided into the communicating bodies, then, the attribute elements of the ascending disc section are used as the reference, the attribute element closest to the descending disc section is found out, after the attribute types of the ascending disc section and the descending disc section are compared, and if the two attribute types belong to the same type, the communicating bodies where the ascending disc section and the descending disc section are positioned are marked as the same communicating body.
According to the method for confirming reservoir connectivity, provided by the embodiment of the application, through carrying out first-time connector division on element types on a first reservoir and a second reservoir of a reservoir model, preliminary screening of reservoir data is realized, the division range of connectors is reduced, on the basis, a target second element corresponding to the first element in the second element is determined according to the space distance between the first element of a section on the first reservoir and the second element of a section on the second reservoir in the reservoir model, and then judgment is carried out according to the element type of the first element and the element type of the target second element, and when the element types are matched, the second target connector and the first connector are calibrated to be the same connector, so that connectors on two sides of a fault are marked, and the connector division of the whole reservoir structure is accurate.
Fig. 3 is a second flowchart of a method for confirming reservoir connectivity according to an embodiment of the present application. The implementation body of the method may be a server or other servers, and the embodiment is not particularly limited herein, and as shown in fig. 3, before determining, according to a distance between a first attribute element of a first fracture on a first reservoir and a target fault line on the fracture in the reservoir model, the first element in the first attribute element, the method may include:
s301, determining fault line coordinates in the section coordinate information according to the section coordinate information of the section in the reservoir model.
The section coordinate information may refer to coordinate information of each point on the section, and is used to represent the position and the size of the section. The section coordinate information can be directly read from the fault model, can be stored in a matrix form, and can also be stored in an array form.
The fault line coordinate may refer to coordinate information marked as a fault line in the fault line coordinate information, and after the fault line coordinate information is obtained, the fault line coordinate may be manually read or directly read by a computer.
S302, obtaining a target fault line according to fault line coordinates in the section coordinate information;
s303, dividing the area where the target fault line is located in the reservoir model into a first area and a second area along a first preset direction.
The first preset direction may refer to a direction requirement set in advance according to a requirement, and the requirement may be a normal direction or a tangential direction, the first preset direction is the tangential direction of each horizontal plane of the fault structure, and a target fault line and a region where the target fault line is located are obtained by utilizing the tangential planes of the horizontal planes and the faults.
The region in which the target fault line is located may refer to a region including the target fault line, and a periphery of the target fault line for selecting attribute elements around the fault line. In the embodiment of the application, the peripheral area range of the target fault line is selected according to the unit of the reservoir model, so that attribute elements around the fault line are all divided into the area where the target fault line is located, and the method can be realized by combining a plurality of algorithms, and the selection of a specific algorithm is determined according to the size of data volume and the reservoir structure.
In some embodiments, the maximum grid spacing d is obtained according to the constructed reservoir model, so as to ensure unit setting of the area where the target fault line is located, a geometric extension distance is 4*d in the normal direction of the target fault line, a geometric extension distance in the tangential direction is 6*d, minkowski and an algorithm are selected to achieve division of the area where the target fault line is located, minkowski and a circle with a radius of 4d are obtained, and a banded area O is obtained as the area where the target fault line is located.
The first region may refer to a part of the target region after the target region is divided according to the target fault line, the second region may refer to another part of the divided region, and the first region and the second region are used for ensuring that the found attribute elements belong to either the first reservoir upper section or the second reservoir upper section, but not both of the first and second reservoir upper sections, or otherwise affect the division of the communication body on the fault.
In an embodiment of the present application, dividing a region where a target fault line is located in a reservoir model into a first region and a second region along a first preset direction includes:
acquiring a plane coordinate image of a target fault line in a reservoir model along a first preset direction;
determining the area where the target fault line is located according to the plane coordinate image of the target fault line in the reservoir model;
According to the region where the target fault line is located, carrying out extension treatment on the target fault line to obtain a treated target fault line;
and dividing the region where the target fault line is located into a first region and a second region according to the processed target fault line.
The plane coordinate image may refer to an image represented by two-dimensional coordinates, and is used for reflecting the position and the form of the target fault line on a certain horizontal plane.
The extension processing may refer to extending the specified line segment in a set direction, and may be in a tangential direction or in a normal direction, and in the embodiment of the present application, the target fault line is extended on the tangential line in the area where the target fault line is located, so that the extended target fault line intersects with the boundary of the area where the target fault line is located.
S304, determining attribute elements on the section in the first area according to the distance between the attribute elements of the first area and the target fault line, wherein the attribute elements on the section in the first area are characterized as first elements;
S305, determining attribute elements on the section in the second area according to the distance between the attribute elements of the second area and the target fault line, wherein the attribute elements on the section in the second area are characterized as second elements.
The distance between the attribute element of the first area and the target fault line may refer to the distance between each attribute element of the first area and the target fault line, the distance between the attribute element of the second area and the target fault line may refer to the distance between each attribute element of the second area and the target fault line, the attribute elements in the first area on the section may be multiple, and the attribute elements in the second area on the section may be multiple.
In some embodiments, a polygon segmentation algorithm is selected to realize the division of the first area and the second area, two ends of a fault line extend outwards to obtain a fault line area L, the fault line area L is subtracted by an area O where a target fault line is located, so that the area O where the target fault line is located is segmented into two polygons O1 and O2, the two polygons are respectively located at two sides of the fault line, thereby ensuring that the found attribute elements belong to a rising disk section or a falling disk section, wherein the extension of the fault line can be realized by using the algorithm, the form of the algorithm does not set requirements, in the embodiment of the application, the fault line is sampled by a two-point three-time Hermite interpolation algorithm to obtain a plane coordinate image of the target fault line in a reservoir model, and the two ends of the target fault line extend by 6*d distances along the fault tangent direction, so as to obtain the fault line area L.
For example, for different horizontal planes, different target fault lines are selected, for each target fault line, a region O where the corresponding target fault line is located and a fault line region L are obtained, a certain target fault line is taken as an example, the region O where the target fault line is located is divided into two polygons O1 and O2, attribute elements A1, A2 and A3 are obtained in the polygon O1, attribute elements B1, B2 and B3 are obtained in the polygon O2, distances between the attribute elements and the target fault line are d1, d2, d3, d4, d5 and d6 respectively, after comparing d1, d2 and d3, the nearest distance between A1 and the target fault line is obtained, then A1 is the attribute element on O1, and similarly, the attribute element on B1 is O2 can be obtained. Therefore, multiple groups of A1 and B1 can be obtained for multiple target fault lines, at the moment, attribute elements with higher spatial positions are selected and divided into elements of the section of the ascending disc, attribute elements with lower spatial positions are divided into elements of the section of the descending disc, and therefore the ascending disc and the descending disc can be marked.
According to the method for confirming reservoir connectivity, firstly, the fault line coordinates are determined through the section coordinate information of the section in the reservoir model, the target fault line is obtained, the target fault line is subjected to extension treatment, and is ensured to intersect with the regional boundary, so that the region where the target fault line is located is divided, a first region and a second region are obtained, then, according to the distances between the attribute elements of the first region and the second region and the target fault line, the attribute elements in the respective regions are determined, namely the first element and the second element, respectively, the attribute elements can be ensured to be correctly distributed to the first reservoir and the second reservoir, the accurate analysis of the spatial relationship and the connectivity between reservoirs is facilitated, and scientific basis is provided for geological modeling and resource development.
Fig. 4 is a flowchart illustrating a method for confirming reservoir connectivity according to an embodiment of the present application. The execution body of the method may be a server or other servers, and as shown in fig. 4, after determining that the second target communication body and the first communication body are the same communication body, when there are multiple same communication bodies and multiple fault planes between the same communication bodies, the method further includes:
s401, determining projection of the same communication body on a fault plane.
The fault plane may refer to a plane between one communicating body and another communicating body, and may refer to a plurality of fault planes between one communicating body and another communicating body, or may refer to a plurality of fault planes between different communicating bodies, for example, a plurality of fault planes between a communicating body a and a communicating body B, a fault plane between a communicating body B and a communicating body C, and a fault plane between a communicating body C and a communicating body D.
In some embodiments, the fault plane is divided into two cases of fault plane intersection and fault plane non-intersection, a plurality of fault line coordinates in the fault plane coordinate information are determined according to the fault plane after the same communication body is divided in the fault model, the intersection points of the plurality of fault lines are judged by using the fault line coordinates, and the plurality of fault lines in the fault plane coordinate information are divided by using the coordinates of the intersection points, so that the fault plane is divided into new sub-faults, and the fault plane intersection result is determined. For example, whether the fault planes are intersected with each other or not and the positions of n intersected points are determined by using a Bentley-Ottmann algorithm, so that the fault line is cut into n+1 sub fault lines, and n+1 sub fault planes are obtained and are respectively numbered as a1, a2, a.
The projection of the same communication body on the fault plane can refer to the mapping of the same communication body on the fault plane, which is used for reflecting the contact condition between the same communication body and the fault plane.
S402, when the same communication body has projection on the fault plane, the same communication body having projection on the fault plane is taken as a target communication body.
In the embodiment of the present application, whether the same communication body is on one side of the fault plane is determined by using the orthographic projection area, for example, when the orthographic projection area of the same communication body on the fault plane is greater than 0, the same communication body is indicated to be on one side of the fault plane, the same communication body is adjacent to the fault plane, the same communication body is determined to be the target communication body, when the orthographic projection area of the same communication body on the fault plane is equal to 0, the same communication body is indicated to be not on one side of the fault plane, the same communication body is not adjacent to the fault plane, and the same communication body is not the target communication body.
S403, determining the adjacent relation between the target communicating bodies according to the positions of the target communicating bodies in the reservoir model.
Wherein the location of the target communication in the reservoir model may refer to the relative location of the target communication in the reservoir model to the fault.
The adjacency between target communicating bodies may refer to one target communicating body being adjacent to another target communicating body for representing the relative positional relationship between the target communicating bodies.
In an embodiment of the present application, determining, according to a location of a target communication body in a reservoir model, an adjacency relationship between the target communication bodies includes:
Determining the position distance relation between the target communicating body and the fault plane according to the position of the target communicating body in the reservoir model;
And determining the adjacent relation between the target communicating bodies according to the position distance relation between the target communicating bodies and the fault plane.
The positional distance relationship between the target communication and the fault plane may refer to the distance between the target communication and the fault plane, and may be represented by the distance between the centroid of the target communication and the centroid of the fault plane, or may be represented by the distance between a point on the edge of the target communication and a point on the edge of the fault plane.
In some embodiments, the projection of the same communication on the fault plane is utilized to determine the faults adjacent to the same communication, and then the adjacent relation between the same communication is determined according to the faults adjacent to the same communication.
S404, generating and displaying a representation adjacency relation graph according to the adjacency relation among the target communicating bodies.
The adjacency graph can refer to the distribution of the object communicating bodies and the faults, the distribution can refer to the adjacency of one object communicating body with another object communicating body through the faults, and the distribution comprises the adjacency of all the object communicating bodies and all the faults. In the embodiment of the application, a minimum weighted undirected connected graph is used as an adjacency graph, wherein the weight is the Euclidean distance between a target communicating body and a fault plane.
In some embodiments, the construction of the adjacency graph is divided into two cases of an intersecting fault and a non-intersecting fault according to the intersecting fault, when the intersecting fault exists, the shortest path between all communicating bodies is searched through Dijkstra algorithm, only paths passing through one fault vertex are reserved, the minimum weighted undirected communication graph is constructed, then all sub fault vertexes are replaced by parent fault vertexes, and when the intersecting fault does not exist, the minimum weighted undirected communication graph is constructed based on the relation graph directly through the minimum spanning tree algorithm.
According to the method for confirming the connectivity of the reservoir, after the second target communication body and the first communication body are confirmed to be the same communication body, the target communication bodies adjacent to the fault plane are identified by determining projections of the same communication body on the fault planes, and the adjacent relation between the second target communication body and the first communication body is determined according to the positions of the target communication bodies in a reservoir model. Finally, a minimum weighted undirected graph is generated and displayed according to the adjacency relations, and is used as an adjacency relation graph to show the distribution situation of the target communicating body and the fault, wherein the weight is based on the Euclidean distance between the target communicating body and the fault plane, and the minimum weighted undirected graph is constructed by dividing the two different situations of the existence of the intersection fault and the nonexistence of the intersection fault, so that the connectivity between the communicating bodies is systematically evaluated and optimized.
Fig. 5 is a schematic diagram of the result of a method for confirming reservoir connectivity according to an embodiment of the present application, as shown in fig. 5, there is a region O around a target fault line where a corresponding target fault line is located, the region O where the target fault line is located is in a strip shape, a fault line region L can be obtained by extending the fault line, the fault line region L divides the region O where the target fault line is located, two polygons O1 and O2 are obtained, and respective attribute elements are respectively present in the O1 and O2.
Fig. 6a is a schematic diagram of a result of a method for confirming reservoir connectivity according to an embodiment of the present application, as shown in fig. 6a, where in the first embodiment, there are 5 faults in the first embodiment, and 5 connectors are obtained in total by determining the target second element, and by dividing the first connector and dividing the target second connector and marking the same connector, it can be seen that the fault lines do not intersect at this time.
Fig. 6b is a schematic diagram of the result of a method for confirming reservoir connectivity according to an embodiment of the present application, as shown in fig. 6b, in the first embodiment, a graph showing the adjacency relationship between a plurality of connected bodies and a plurality of fault lines is shown, in the first embodiment, the connected bodies 1, 2 and 3 are distributed around the fault 1, the connected bodies 3 and 4 are distributed on two sides of the fault 3, and the connected bodies 4 and 5 are distributed on two sides of the fault 4.
Fig. 6c is a schematic diagram showing the result of a method for confirming reservoir connectivity according to an embodiment of the present application, as shown in fig. 6c, in the first embodiment, a graph showing the communication coefficients of a plurality of communicating bodies is shown, it can be seen that the communication coefficient of the communicating body 1 and the communicating body 3 is 0.4, the communication coefficient of the communicating body 3 and the communicating body 4 is 0.56, it can be inferred that the connectivity between the communicating body 1 and the communicating body 3 is smaller than the connectivity between the communicating body 3 and the communicating body 4, which indicates that the flow capacity between the communicating body 1 and the communicating body 3 is weaker than the flow capacity between the communicating body 3 and the communicating body 4.
Fig. 7a is a schematic diagram showing the result of a method for confirming reservoir connectivity according to an embodiment of the present application, as shown in fig. 7a, in the second embodiment, the fault direction is different and there are intersecting faults, in the second embodiment, there are 4 faults, and through determination of the target second element, 5 communicating bodies are obtained in total by dividing the first communicating body and dividing the target second communicating body and marking the same communicating body, and it can be seen that the fault line 2 and the fault line 3 intersect at this time.
Fig. 7b is a schematic diagram of the result of a method for confirming reservoir connectivity according to an embodiment of the present application, as shown in fig. 7b, in a second embodiment, an adjacency graph of a plurality of connected bodies and a plurality of fault lines is shown, in a second embodiment, the connected body 1 includes a fault 1, the connected body 2, the connected body 3 and the connected body 4 are distributed around an intersection point of the fault line 2 and the fault line 3, the connected body 3 and the connected body 4 are distributed on two sides of the fault 3, and the connected body 4 and the connected body 5 are distributed on two sides of the fault 4.
Fig. 7c is a schematic diagram seventh of the result of a method for confirming reservoir connectivity according to an embodiment of the present application, as shown in fig. 7c, in a second embodiment, a communication coefficient diagram of a plurality of communicating bodies is shown, it can be seen that the communication coefficient of the communicating body 1 and the communicating body 3 is 0.46, the communication coefficient of the communicating body 2 and the communicating body 4 is 0.12, it can be inferred that the connectivity between the communicating body 1 and the communicating body 3 is larger than the connectivity between the communicating body 2 and the communicating body 4, which indicates that the flow capacity between the communicating body 1 and the communicating body 3 is stronger than the flow capacity between the communicating body 2 and the communicating body 4.
Fig. 8 is a schematic structural diagram of a device for confirming reservoir connectivity according to an embodiment of the present application, and as shown in fig. 8, the reservoir connectivity analysis device 80 includes a first determining module 801 and a second determining module 802. Wherein:
A first determining module 801, configured to determine, according to a spatial distance between a first element of a section on a first reservoir and a second element of a section on a second reservoir in the reservoir model, a target second element corresponding to the first element in the second element, where the section on the first reservoir and the section on the second reservoir are the same section;
In one possible implementation manner, according to a spatial distance between a first element of a section on a first reservoir and a second element of a section on a second reservoir in the reservoir model, determining a target second element corresponding to the first element in the second element, where the first determining module 801 is specifically configured to:
Determining a first element in the first attribute elements according to the distance between the first attribute elements of the fracture surface on the first reservoir and a target fault line on the fracture surface in the reservoir model, wherein the first attribute elements represent elements of a storage element type on the fracture surface of the first reservoir, and the target fault line represents a line intersecting a horizontal plane with the fracture surface;
Determining a second element in the second attribute elements according to the distance between the second attribute elements of the section on the second reservoir in the reservoir model and the target fault line on the section, wherein the second attribute elements represent elements of the storage element type on the section of the second reservoir;
And determining a target second element corresponding to the first element in the second element according to the space distance between the first element and the second element.
In one possible implementation, before determining the first element of the first attribute elements according to the distance between the first attribute element of the fracture on the first reservoir and the target fault line on the fracture in the reservoir model, the first determining module 801 may be further specifically configured to:
determining fault line coordinates in the section coordinate information according to the section coordinate information of the section in the reservoir model;
obtaining a target fault line according to the fault line coordinates in the section coordinate information;
Dividing a region where a target fault line is located in the reservoir model into a first region and a second region along a first preset direction;
Determining attribute elements on the section in the first area according to the distance between the attribute elements of the first area and the target fault line, wherein the attribute elements on the section in the first area are characterized as first elements;
and determining attribute elements on the section in the second area according to the distance between the attribute elements of the second area and the target fault line, wherein the attribute elements on the section in the second area are characterized as second elements.
In one possible implementation manner, along a first preset direction, the area where the target fault line is located in the reservoir model is divided into a first area and a second area, and the first determining module 801 may be further specifically configured to:
acquiring a plane coordinate image of a target fault line in a reservoir model along a first preset direction;
determining the area where the target fault line is located according to the plane coordinate image of the target fault line in the reservoir model;
According to the region where the target fault line is located, carrying out extension treatment on the target fault line to obtain a treated target fault line;
and dividing the region where the target fault line is located into a first region and a second region according to the processed target fault line.
The second determining module 802 is configured to determine that the second target communication is the same communication as the first communication when the element type of the first element matches the element type of the target second element, the second target communication is the communication where the target second element is located in the second reservoir, and the first communication is the communication where the first element is located in the first reservoir.
In one possible implementation manner, after determining that the second target communication body and the first communication body are the same communication body, when there are multiple identical communication bodies and multiple fault planes between the multiple identical communication bodies, the second determining module 802 may be further specifically configured to:
determining the projection of the same communication body on the fault plane;
When the same communicating body has projection on the fault plane, the same communicating body with projection on the fault plane is taken as a target communicating body;
determining the adjacent relation between the target communicating bodies according to the positions of the target communicating bodies in the reservoir model;
and generating and displaying a representation adjacency relation graph according to adjacency relations among the target communicating bodies.
In one possible implementation, the second determining module 802 may be further specifically configured to determine, according to a location of the target continuum in the reservoir model, an adjacency between the target continuum:
Determining the position distance relation between the target communicating body and the fault plane according to the position of the target communicating body in the reservoir model;
And determining the adjacent relation between the target communicating bodies according to the position distance relation between the target communicating bodies and the fault plane.
The embodiment of the application provides a device for confirming reservoir connectivity, which can be used for executing the reservoir connectivity analysis method in any of the embodiments, and the implementation principle and the technical effect are similar and are not repeated here.
It should be noted that, it should be understood that the division of the modules of the above apparatus is merely a division of a logic function, and may be fully or partially integrated into a physical entity or may be physically separated. And these modules may all be implemented in the form of software calls through the processing elements. Or may be implemented entirely in hardware. The method can also be realized in a form of calling software by a processing element, and the method can be realized in a form of hardware by a part of modules. In addition, all or part of the modules may be integrated together or may be implemented independently. The processing element here may be an integrated circuit with signal processing capabilities. In implementation, each step of the above method or each module above may be implemented by an integrated logic circuit of hardware in a processor element or an instruction in a software form.
Fig. 9 is a schematic structural diagram of a device for confirming reservoir connectivity according to an embodiment of the present application. As shown in fig. 9, the reservoir connectivity confirmation device 90 includes:
The reservoir connectivity verification device 90 may include one or more processor cores 901, one or more computer-readable storage medium memories 902, communication components 903, and the like. The processor 901, the memory 902, and the communication unit 903 are connected via a bus 904.
In a specific implementation, at least one processor 901 executes computer-executable instructions stored in memory 902, causing the at least one processor 901 to perform a reservoir communication analysis method as described above.
The specific implementation process of the processor 901 may refer to the above-mentioned method embodiment, and its implementation principle and technical effects are similar, and this embodiment will not be described herein again.
In the embodiment shown in fig. 9, it should be understood that the Processor may be a central processing unit (english: central Processing Unit, abbreviated as CPU), or may be other general purpose processors, digital signal processors (english: DIGITAL SIGNAL Processor, abbreviated as DSP), application-specific integrated circuits (english: application SPECIFIC INTEGRATED Circuit, abbreviated as ASIC), or the like. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like. The steps of a method disclosed in connection with the present invention may be embodied directly in a hardware processor for execution, or in a combination of hardware and software modules in a processor for execution.
The Memory may include high-speed Memory (Random Access Memory, RAM) or may further include Non-volatile Memory (NVM), such as at least one disk Memory.
The bus may be an industry standard architecture (Industry Standard Architecture, ISA) bus, an external device interconnect (PERIPHERAL COMPONENT, PCI) bus, or an extended industry standard architecture (Extended Industry Standard Architecture, EISA) bus, among others. The buses may be divided into address buses, data buses, control buses, etc. For ease of illustration, the buses in the drawings of the embodiments of the present application are not limited to only one bus or one type of bus.
The embodiments of the present application also provide a computer program product comprising a computer program which, when executed by a processor, implements the above-mentioned method.
The embodiment of the application also provides a computer readable storage medium, wherein computer execution instructions are stored in the computer readable storage medium, and when a processor executes the computer execution instructions, the method is realized.
The above-described readable storage medium may be implemented by any type or combination of volatile or non-volatile memory devices, such as Static Random Access Memory (SRAM), electrically erasable programmable read-only memory (EEPROM), erasable programmable read-only memory (EPROM), programmable read-only memory (PROM), read-only memory (ROM), magnetic memory, flash memory, magnetic or optical disk. A readable storage medium can be any available medium that can be accessed by a general purpose or special purpose computer.
An exemplary readable storage medium is coupled to the processor such the processor can read information from, and write information to, the readable storage medium. In the alternative, the readable storage medium may be integral to the processor. The processor and the readable storage medium may reside in an Application SPECIFIC INTEGRATED Circuits (ASIC). The processor and the readable storage medium may reside as discrete components in a device.
The division of units is merely a logical function division, and there may be another division manner in actual implementation, for example, a plurality of units or components may be combined or may be integrated into another system, or some features may be omitted or not performed. Alternatively, the coupling or direct coupling or communication connection shown or discussed with each other may be an indirect coupling or communication connection via some interfaces, devices or units, which may be in electrical, mechanical or other form.
The units described as separate units may or may not be physically separate, and units shown as units may or may not be physical units, may be located in one place, or may be distributed over a plurality of network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.
In addition, each functional unit in the embodiments of the present invention may be integrated in one processing unit, or each unit may exist alone physically, or two or more units may be integrated in one unit.
The functions, if implemented in the form of software functional units and sold or used as a stand-alone product, may be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention may be embodied essentially or in a part contributing to the prior art or in a part of the technical solution in the form of a software product stored in a storage medium, comprising several instructions for causing a computer device (which may be a personal computer, a server, a network device, etc.) to perform all or part of the steps of the method of the embodiments of the present invention. The storage medium includes a U disk, a removable hard disk, a Read-Only Memory (ROM), a random access Memory (RAM, random Access Memory), a magnetic disk, an optical disk, or other various media capable of storing program codes.
Those of ordinary skill in the art will appreciate that all or a portion of the steps of implementing the various method embodiments described above may be implemented by hardware associated with program instructions. The foregoing program may be stored in a computer readable storage medium. The program, when executed, performs the steps comprising the method embodiments described above, and the storage medium described above includes various media capable of storing program code, such as ROM, RAM, magnetic or optical disk.
Finally, it should be noted that other embodiments of the invention will be apparent to those skilled in the art from consideration of the specification and practice of the invention disclosed herein. This invention is intended to cover any adaptations, uses, or adaptations of the invention following, in general, the principles of the invention and including such departures from the present disclosure as come within known or customary practice within the art to which the invention pertains and as may be applied to the precise construction hereinbefore set forth and shown in the drawings and as follows in the scope of the appended claims. The scope of the invention is limited only by the appended claims.