Contour map-based calculation unit reserve parameter calculation method, device and equipment

Abstract

The invention discloses a calculation unit reserve parameter calculation method, device and equipment based on a contour map, and the method comprises the steps: carrying out the closing processing of a contour map which is obtained in advance and contains oil and gas area constraints, and obtaining a closed contour map; determining the parameter value of the boundary line according to the parameter value of the intersection point of the contour line and the boundary line in the closed contour map; performing triangular topology identification on all point information in the closed contour map to obtain all initial triangular patch elements; refining the initial triangular patch elements meeting a preset rule in each initial triangular patch element to obtain all processed triangular patch elements; according to the parameter values of the contour line and the boundary line, determining the area and the value in the area of each processed triangular patch element; and according to the area of each processed triangular patch element and the value in the area, determining reserve parameters constrained by the oil-gas-containing area. According to the method, the reserve parameters can be automatically, rapidly and more accurately calculated under the condition that human intervention is not needed.

Description

Technical Field

[0001] The invention relates to a method, a device and equipment for calculating the reserves parameters of a calculation unit based on a contour map. Background Art

[0002] Oil and gas reserves are the core assets of oil companies and an important indicator of their healthy development potential. Domestic and foreign oil companies attach great importance to reserve calculations and invest a lot of manpower, material and financial resources in it every year. There are various types of oil and gas reserves, including new reserve calculations, reserve recalculations, reserve accounting, reserve calibration and settlement, and SEC reserve estimation. Among them, new reserve calculations also include predicted reserves, controlled reserves, and proven reserves. The technical requirements and emphases of these oil and gas reserve calculations are different. The calculation accuracy control nodes are numerous and wide-ranging, including calculation models, calculation parameters, and technical methods. If any node fails to control the calculation accuracy, error transmission will occur, resulting in poor reserve calculation accuracy. Therefore, it is difficult to accurately calculate reserves.

[0003] The effective thickness, effective porosity and oil and gas saturation of single wells and calculation units are extremely important parameters in reserve calculation, which are related to the accuracy and reliability of reserve calculation. According to the provisions of industry standards and enterprise standards, there are three main methods for calculating the effective thickness, effective porosity and oil and gas saturation of calculation units: well point arithmetic average method, well point area trade-off method and contour area trade-off method. Among them, the contour area trade-off method is the most effective method for finely characterizing the distribution characteristics of parameters, but this method requires accurate calculation of the area between each contour line, and then the parameter value calculated by the weighted method of the area of ​​the average value of the two contour lines (the area is the area between the two contour lines) is the reserve parameter value calculated by the contour area trade-off method. Summary of the invention

[0004] In order to better realize the calculation of the reserve parameters of the calculation unit using the contour area trade-off method, the embodiment of the present application provides a method, device and equipment for calculating the reserve parameters of the calculation unit based on the contour map.

[0005] In a first aspect, an embodiment of the present application provides a method for calculating reserve parameters based on a contour map, the method comprising:

[0006] Perform closing processing on the pre-obtained contour map of oil and gas area constraints to obtain a closed contour map;

[0007] Determining the parameter value of the boundary line according to the parameter value of the intersection point between the contour line and the boundary line in the closed contour map;

[0008] Performing triangular topology recognition on all point information in the closed contour map to obtain all initial triangular facets;

[0009] Refining the initial triangular surface fragments that meet the preset rules in each of the initial triangular surface fragments to obtain all processed triangular surface fragments;

[0010] Determine the area and the value within the area of ​​each processed triangular face element according to the parameter values ​​of the contour lines and the boundary lines;

[0011] According to the area of ​​each processed triangular facet and the value within the area, the reserve parameters constrained by the oil and gas area are determined.

[0012] In one or some optional implementations of the embodiment of the present application, the closing process is performed on the pre-obtained contour map of oil and gas area constraints to obtain a closed contour map, including:

[0013] Determine whether each contour line in the contour line map of the oil and gas area constraint obtained in advance has an intersection with any boundary line;

[0014] If so, delete the remaining line segments after the intersection, determine all intersection points as control points, and obtain a new contour map;

[0015] If not, extend the contour line until it intersects with the boundary line, determine all intersection points as control points, and obtain a new contour map;

[0016] Determine whether all curves in the new contour map are closed;

[0017] If so, the closed contour map is obtained.

[0018] In one or some optional implementations of the embodiments of the present application, obtaining the parameter value of the boundary line according to the parameter value of the intersection of the contour line and the boundary line in the closed contour map includes:

[0019] The parameter value of the boundary line is obtained by interpolating or extrapolating the parameter value of the intersection point of the contour line and the boundary line in the closed contour map.

[0020] In one or some optional implementations of the embodiment of the present application, the initial triangular face element that meets the preset rules includes an initial triangular face element that crosses the contour line and an initial triangular face element whose center is not within the oil and gas containing area;

[0021] The step of thinning the initial triangular surface fragments that meet the preset rules in each of the initial triangular surface fragments to obtain all processed triangular surface fragments includes:

[0022] Add the intersection points of the initial triangles that cross the contour lines and the contour lines to the control points;

[0023] According to the control points, constraining the initial triangular surface elements that cross the contour lines to obtain the processed triangular surface elements;

[0024] as well as,

[0025] Deleting the initial triangular facets whose inner parts are not within the oil and gas containing area;

[0026] Among them, whether the center of the initial triangular surface element is within the oil and gas containing area is determined in the following way: according to the parameter values ​​of the contour lines and boundary lines, the center of the initial triangular surface element is determined to judge whether the center of the initial triangular surface element is within the oil and gas containing area.

[0027] In one or some optional implementations of the embodiment of the present application, constraining the initial triangular surface fragments that cross the contour line according to the control points to obtain the processed triangular surface fragments includes:

[0028] According to the control points, the area of ​​the triangular surface fragment composed of the three most adjacent control points in the initial triangular surface fragment crossing the contour line is converged to be close to zero, so as to obtain the processed triangular surface fragment.

[0029] In one or some optional implementations of the embodiment of the present application, determining the area and the value within the area of ​​each processed triangular face fragment according to the parameter values ​​of the contour lines and the boundary lines includes:

[0030] Determining the parameter values ​​of each vertex of each processed triangular face fragment according to the parameter values ​​of the contour lines and the boundary lines;

[0031] Determining the area of ​​all the processed triangular surface fragments according to the parameter values ​​of each vertex of each processed triangular surface fragment;

[0032] Determine whether the positions of the vertices of the processed triangular face fragments are on the same contour line;

[0033] If yes, then the value within the area of ​​the processed triangle face element is the value of the same contour line;

[0034] If not, the area value of the processed triangular face element is determined according to the weighted quadrilateral formed by the center point of the processed triangular face element and the midpoints of each side.

[0035] In one or some optional implementations of the embodiment of the present application, the determining the value within the area of ​​the processed triangular face element according to the weighted quadrilateral formed by the center point of the processed triangular face element and the midpoints of each side includes:

[0036] According to the weighted quadrilateral formed by the center point of the processed triangular face element and the midpoints of each side, the value within the area of ​​the processed triangular face element is determined based on the following formula 1:

[0037]

[0038] Wherein, ΔV is the value within the area of ​​the processed triangle element; S A , S B and S C are the areas of the quadrilaterals formed by the center points of the processed triangle face fragments and the midpoints of each side; V A 、V B and V C They are respectively the parameter values ​​of the contour lines or boundary lines where the vertices of the processed triangle face fragments are located.

[0039] In one or some optional implementations of the embodiment of the present application, the determination of the reserve parameters constrained by the oil and gas area according to the area of ​​each processed triangular facet and the value within the area includes:

[0040] According to the area of ​​each processed triangular facet and the value within the area, the reserve parameter constrained by the oil and gas area is determined based on the following formula 2:

[0041]

[0042] Wherein, V is the reserve parameter constrained by the oil and gas area; ΔS is the area of ​​the processed triangular face element; ΔV is the value within the area of ​​the processed triangular face element; and N=3.

[0043] In a second aspect, an embodiment of the present application provides a device for calculating reserve parameters based on a contour map, the device comprising:

[0044] A closing processing module is used to perform closing processing on the pre-obtained contour map of oil and gas area constraints to obtain a closed contour map;

[0045] A first determination module is used to determine the parameter value of the boundary line according to the parameter value of the intersection point of the contour line and the boundary line in the closed contour map;

[0046] A topology recognition module is used to perform triangular topology recognition on all point information in the closed contour map to obtain all initial triangular facets;

[0047] A refinement module, used for refining the initial triangular surface fragments that meet the preset rules in each of the initial triangular surface fragments to obtain all processed triangular surface fragments;

[0048] A second determination module is used to determine the area and the value within the area of ​​each processed triangular face element according to the parameter values ​​of the contour lines and the boundary lines;

[0049] The third determination module is used to determine the reserve parameters constrained by the oil and gas area according to the area of ​​each processed triangular facet and the value within the area.

[0050] In a third aspect, an embodiment of the present application provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the method for calculating reserve parameters based on a contour map as described above.

[0051] In a fourth aspect, an embodiment of the present application provides a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the computer program, the method for calculating reserve parameters based on contour maps as described above is implemented.

[0052] In a fifth aspect, an embodiment of the present application provides a computer program product comprising instructions, which, when executed on a computer device, enables the computer device to execute the method for calculating reserve parameters based on contour maps as described above.

[0053] In a sixth aspect, an embodiment of the present application provides a chip, the chip including a processor and a communication interface, the communication interface and the processor are coupled, and the processor is used to run a computer program or instruction to implement the above-mentioned method for calculating reserve parameters based on contour maps.

[0054] The beneficial effects of the above technical solution provided by the embodiment of the present application include at least:

[0055] The method for calculating reserve parameters based on contour maps provided in the embodiment of the present application performs triangular topology recognition on the contour maps constrained by oil and gas area, splits the contour maps into several triangular facets, and obtains the reserve parameters constrained by oil and gas area by calculating the area of ​​each triangular facet and the value within the area. The method can automatically realize fast and more accurate calculation of reserve parameters without human intervention, and can achieve the purpose of ensuring data accuracy while obtaining extremely high work efficiency, and achieve the technical effect of improving quality and greatly improving timeliness in reserve calculation work, and realize automatic, high-timeliness and high-quality reserve parameter calculation based on contour maps.

[0056] Other features and advantages of the present invention will be described in the following description, and partly become apparent from the description, or understood by practicing the present invention. The purpose and other advantages of the present invention can be realized and obtained by the structures particularly pointed out in the written description, claims, and drawings.

[0057] The technical solution of the present invention is further described in detail below through the accompanying drawings and embodiments. BRIEF DESCRIPTION OF THE DRAWINGS

[0058] The accompanying drawings are used to provide a further understanding of the present invention and constitute a part of the specification. Together with the embodiments of the present invention, they are used to explain the present invention and do not constitute a limitation of the present invention. In the accompanying drawings:

[0059] Figure 1 A schematic diagram of the steps of a method for calculating reserve parameters based on a contour map provided in an embodiment of the present application;

[0060] Figure 2 A pre-obtained contour map provided for an embodiment of the present application;

[0061] Figure 3 The contour map after closing processing provided in the embodiment of the present application;

[0062] Figure 4 A schematic diagram of all initial triangular facets in the contour map provided in the embodiment of the present application;

[0063] Figure 5 A schematic diagram with labels of all initial triangular facets in the contour map provided in the embodiment of the present application;

[0064] Figure 6 A schematic diagram of all triangular facets in the contour map provided in the embodiment of the present application;

[0065] Figure 7 A schematic diagram of a quadrilateral whose vertices are formed by the center of a triangular facet element and the midpoints of each side and the vertices provided in an embodiment of the present application;

[0066] Figure 8 A flowchart of a method for calculating reserve parameters based on contour maps provided in an embodiment of the present application;

[0067] Fig. 9 A schematic diagram of the structure of a device for calculating reserve parameters based on a contour map provided in an embodiment of the present application. DETAILED DESCRIPTION

[0068] In the following description, specific details such as specific system structures, technologies, etc. are provided for the purpose of illustration rather than limitation, so as to provide a thorough understanding of the embodiments of the present application. However, it should be clear to those skilled in the art that the present application may also be implemented in other embodiments without these specific details. In other cases, detailed descriptions of well-known systems, devices, circuits, and methods are omitted to prevent unnecessary details from obstructing the description of the present application.

[0069] It should be understood that when used in the present specification and the appended claims, the term "comprising" indicates the presence of described features, wholes, steps, operations, elements and / or components, but does not exclude the presence or addition of one or more other features, wholes, steps, operations, elements, components and / or combinations thereof.

[0070] It should also be understood that the term “and / or” used in the specification and appended claims refers to any and all possible combinations of one or more of the associated listed items, and includes these combinations.

[0071] As used in the specification and appended claims of this application, the term "if" can be interpreted as "when" or "uponce" or "in response to determining" or "in response to detecting", depending on the context. Similarly, the phrase "if it is determined" or "if [described condition or event] is detected" can be interpreted as meaning "uponce it is determined" or "in response to determining" or "uponce [described condition or event] is detected" or "in response to detecting [described condition or event]", depending on the context.

[0072] In addition, in the description of the present application specification and the appended claims, the terms "first", "second", "third", etc. are only used to distinguish the descriptions and cannot be understood as indicating or implying relative importance.

[0073] References to "one embodiment" or "some embodiments" etc. described in the specification of this application mean that one or more embodiments of the present application include specific features, structures or characteristics described in conjunction with the embodiment. Therefore, the statements "in one embodiment", "in some embodiments", "in some other embodiments", "in some other embodiments", etc. that appear in different places in this specification do not necessarily refer to the same embodiment, but mean "one or more but not all embodiments", unless otherwise specifically emphasized in other ways. The terms "including", "comprising", "having" and their variations all mean "including but not limited to", unless otherwise specifically emphasized in other ways.

[0074] It should be understood that the size of the serial numbers of the steps in the following embodiments does not mean the order of execution. The execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiments of the present application.

[0075] In order to illustrate the technical solution of the present application, a specific embodiment is provided below for illustration.

[0076] The inventors found that in the prior art, the most difficult part of calculating effective thickness and other reserve parameters using the contour area trade-off method is how to accurately determine the area between two contour lines. To determine the area between contour lines, it is necessary to manually circle the area POLYGON one by one, or to automatically generate POLYGON by selecting contour line objects one by one through related software. Then, it is necessary to determine the parameter value applied to the area between contour lines. If the area mainly surrounded by contour lines can be basically determined as the average value of the two contour line values, in many cases, such as fault block reservoirs, it is often difficult to determine the parameter value of the area surrounded by contour lines and the boundary line of the oil and gas reservoir. This will bring three problems. First, it is difficult to determine the area between contour lines. To ensure accuracy, because the accuracy of manual delineation is poor, the resulting errors directly affect the calculation accuracy of relevant reserve parameters, and the problem of selecting values ​​in the contour area involving boundary lines (fault lines) is even greater; second, the workload is huge. Usually, in order to better control the trend of reserve parameter changes, the number of contour lines will be large, which brings a lot of picking and selecting workload. If there are more calculation units, the picking and selecting workload will be even greater, and the reserve work will usually be changed and adjusted many times, changed once, and repeated once, this part of the workload will be huge; third, manual delineation cannot automatically transmit the change of reserves to the entire reserve calculation process, resulting in a broken chain of data transmission, which makes it impossible to realize the automation and intelligence of reserve calculation. Therefore, it is urgent to develop a contour map based on the effective thickness, effective porosity and oil and gas saturation determined by the logging interpretation results, automatically and accurately determine the reserve parameters, improve the accuracy and timeliness of the calculation of relevant reserve parameters, and ensure the reliability of reserve calculation. Based on this, the inventors have made the present invention after further research and development, providing a method, device and equipment for calculating unit reserve parameters based on contour maps.

[0077] Embodiment 1

[0078] The present application embodiment provides a method for calculating reserve parameters based on contour maps, referring to Figure 1 As shown, the method includes:

[0079] S101: Closing the previously obtained contour map with oil and gas area constraints to obtain a closed contour map.

[0080] In the above step S101, the closing process is performed on the pre-obtained contour map of oil and gas area constraints to obtain a closed contour map, which specifically includes:

[0081] Determine whether each contour line in the contour line map of the oil and gas area constraint obtained in advance has an intersection with any boundary line;

[0082] If so, delete the remaining line segments after the intersection, determine all intersection points as control points, and obtain a new contour map;

[0083] If not, extend the contour line until it intersects with the boundary line, determine all intersection points as control points, and obtain a new contour map;

[0084] Determine whether all curves in the new contour map are closed;

[0085] If so, the closed contour map is obtained.

[0086] In the embodiment of the present application, for the pre-obtained contour map constrained by the oil and gas area, it is necessary to perform closing processing, that is, to ensure that all the curves in the contour map are intersecting closed curves. First, define the boundary line and contour line object types of the pre-obtained contour map constrained by the oil and gas area, obtain all the boundary lines and contour lines in the pre-obtained contour map constrained by the oil and gas area, and then calculate the intersection of the contour line and the boundary line. Since both ends of the contour line have intersections with the boundary line under normal circumstances, the closure of the pre-obtained contour map constrained by the oil and gas area can be judged by judging whether there is an intersection between the contour line and the boundary line in the pre-obtained contour map constrained by the oil and gas area. If there is an intersection, delete the remaining line segments after the intersection, and determine all intersections as control points; if there is no intersection, extend the contour line to intersect with the boundary line, and determine all intersections as control points. It should be noted that if the contour line has an intersection with other contour lines, the endpoints of the contour line at the intersection are deleted, and no calculation and judgment are performed. After the above closing process, a new contour map is obtained. It is necessary to determine whether all the curves in the new contour map are closed. If so, a closed contour map is obtained; if not, the closing process continues.

[0087] In the embodiment of the present application, the specific implementation method of determining the boundary lines and contour line object types of the contour map of the oil and gas area constraint can refer to the detailed description in the prior art and is not specifically limited here.

[0088] In a specific embodiment, Figure 2 As shown in the figure, it is a pre-obtained contour map of oil and gas area constraints. Some contour lines in the figure are not completely closed with the boundary line, that is, there is no intersection. It is necessary to extend the contour line to intersect with the boundary line, such as Figure 3 As shown in the figure, it is a closed contour map obtained after closing processing. Each contour line has an intersection with the boundary line, and all curves are closed.

[0089] S102: Determine a parameter value of the boundary line according to the parameter value of the intersection point between the contour line and the boundary line in the closed contour map.

[0090] In the above step S102, the step of determining the parameter value of the boundary line according to the parameter value of the intersection point between the contour line and the boundary line in the closed contour map specifically includes:

[0091] The parameter value of the boundary line is obtained by interpolating or extrapolating the parameter value of the intersection point of the contour line and the boundary line in the closed contour map.

[0092] In the embodiment of the present application, it is necessary to assign parameter values ​​to the boundary line. Since the boundary line and the contour line are different, there is no parameter value information. After the closing process of step S101, since the contour line and the boundary line have an intersection, there is parameter value information of at least two points on the boundary line, that is, there is a parameter value trend. By interpolating or extrapolating the parameter value information of these two points, parameter value information can be assigned to all points of a boundary line, that is, parameter values ​​are assigned to the boundary line.

[0093] S103: Perform triangular topology recognition on all point information in the closed contour map to obtain all initial triangular facets.

[0094] In the above step S103, after obtaining the closed contour map, it is necessary to perform triangular topological structure recognition on the point information involved in all the primitives in the closed contour map, and obtain all the initial triangular facets through topological structure recognition and establishment. The specific implementation method of triangular topological recognition can refer to the detailed description in the prior art and is not specifically limited here.

[0095] In a specific embodiment, Figure 4 As shown in FIG. 1 , all initial triangular face elements are obtained by triangulating all point information in a closed contour map through triangular topology recognition.

[0096] S104: Thinning the initial triangular surface fragments that meet the preset rules in each of the initial triangular surface fragments to obtain all processed triangular surface fragments.

[0097] In the above step S104, the initial triangular surface elements that meet the preset rules include initial triangular surface elements that cross the contour line and initial triangular surface elements whose center points are not within the oil and gas bearing area;

[0098] The step of thinning the initial triangular surface fragments that meet the preset rules in each of the initial triangular surface fragments to obtain all processed triangular surface fragments specifically includes:

[0099] Add the intersection points of the initial triangles that cross the contour lines and the contour lines to the control points;

[0100] According to the control points, constraining the initial triangular surface elements that cross the contour lines to obtain the processed triangular surface elements;

[0101] as well as,

[0102] Deleting the initial triangular facets whose inner parts are not within the oil and gas containing area;

[0103] Among them, whether the center of the initial triangular surface element is within the oil and gas containing area is determined in the following way: according to the parameter values ​​of the contour lines and boundary lines, the center of the initial triangular surface element is determined to judge whether the center of the initial triangular surface element is within the oil and gas containing area.

[0104] In the embodiment of the present application, triangular topology identification based on all control points in the closed contour map cannot directly establish the correct contour area topology structure. The initial triangular face elements may cross the contour lines or fault lines due to insufficient control over the positions of the control points. In addition, due to changes in the concave and convex properties of the boundary lines of the oil and gas containing areas, some triangular face elements may not be within the oil and gas containing areas.

[0105] In a specific embodiment, Figure 5 As shown in FIG. 1 , all initial triangular face elements in a closed contour map are shown. Regions B, C, and D in the figure are initial triangular face elements that cross the contour lines, and region A in the figure is an initial triangular face element whose center is not within the oil and gas bearing area.

[0106] In the embodiment of the present application, it is necessary to refine the initial triangular facets that meet the preset rules. Assume that the oil and gas area is domain S, and each contour line control domain is S i , then S=∑S i , assuming that the initial triangular face element is the basic unit, denoted by Δ j , then S i =∑Δ j To ensure the recovery of the correct oil and gas area and the establishment of the topological structure of the contour area, it is necessary to ensure that Δ j ∈S i Therefore, the preset rule can be set as a triangle patch that crosses the contour line or the center is not in the oil and gas area. The initial triangle patch that meets the preset rule, that is, the initial triangle patch that crosses the contour line and the initial triangle patch that the center is not in the oil and gas area, is refined so that the minimum number of triangle patches whose center is in the oil and gas area and does not cross the contour line is the optimal convergence condition. The gradient method can be used to find the optimal convergence condition. The specific implementation process can refer to the detailed description of the prior art and will not be repeated here.

[0107] In the embodiment of the present application, constraining the initial triangular face fragments that cross the contour line according to the control points to obtain the processed triangular face fragments specifically includes:

[0108] According to the control points, the area of ​​the triangular surface fragment composed of the three most adjacent control points in the initial triangular surface fragment crossing the contour line is converged to be close to zero, so as to obtain the processed triangular surface fragment.

[0109] In the embodiment of the present application, the initial triangular facets crossing the contour lines are refined in the following manner:

[0110] Assuming that the initial triangular face element that meets the preset rules crosses the contour control domain, that is, crosses the contour line, there must be two intersection points on the crossed contour line to split the initial triangular face element, and form areas of and The intersection of the contour line and the initial triangular face element is added as a new control point to the original oil and gas area constraint, so that the area of ​​the new triangular face element composed of the three nearest control points in the initial triangular face element crossing the contour line is gradually reduced. That is, when the area of ​​the triangle composed of the three nearest control points in the initial triangle across the contour line approaches zero, the convergence condition is reached, and the processed triangle is obtained. Among them, k is the number of the area operation between two contour lines. If there are multiple contour lines, it is necessary to automatically take multiple times according to the value of the contour line.

[0111] In the embodiment of the present application, the initial triangular facets whose inner parts are not within the oil and gas-bearing area are refined in the following manner:

[0112] First, the incenters of all initial triangular face elements are calculated to determine whether the incenter of each initial triangular face element is within the oil and gas bearing area. If not, the initial triangular face element is deleted.

[0113] In an embodiment of the present application, when refining the initial triangular surface patch, it is possible to first determine whether the initial triangular surface patch crosses the contour line, and then determine whether the center of the initial triangular surface patch is within the oil and gas containing area. Alternatively, it is possible to first determine whether the center of the initial triangular surface patch is within the oil and gas containing area, and then determine whether the initial triangular surface patch crosses the contour line. Of course, it is also possible to simultaneously and in parallel determine whether the initial triangular surface patch crosses the contour line and determine whether the center of the initial triangular surface patch is within the oil and gas containing area.

[0114] For example, if the method for refining the initial triangular surface element is to first determine whether the center of the initial triangular surface element is within the oil and gas containing area, and then determine whether the initial triangular surface element crosses the contour line, at this time, if the center of the initial triangular surface element is within the oil and gas containing area, it is further determined whether the initial triangular surface element crosses the contour line.

[0115] In a specific embodiment, Figure 6 As shown, after all the initial triangular face elements are refined, all the processed triangular face elements are obtained to achieve correct area segmentation between contour lines.

[0116] S105: Determine the area and the value within the area of ​​each processed triangular face element according to the parameter values ​​of the contour lines and the boundary lines.

[0117] In the above step S105, the area and the value within the area of ​​each processed triangular face element are determined according to the parameter values ​​of the contour lines and the boundary lines, specifically including:

[0118] Determining the parameter values ​​of each vertex of each processed triangular face fragment according to the parameter values ​​of the contour lines and the boundary lines;

[0119] Determining the area of ​​all the processed triangular surface fragments according to the parameter values ​​of each vertex of each processed triangular surface fragment;

[0120] Determine whether the positions of the vertices of the processed triangular face fragments are on the same contour line;

[0121] If yes, then the value within the area of ​​the processed triangle face element is the value of the same contour line;

[0122] If not, the area value of the processed triangular face element is determined according to the weighted quadrilateral formed by the center point of the processed triangular face element and the midpoints of each side.

[0123] In an embodiment of the present application, after finding all the processed triangular face fragments that contain the oil and gas area constraint, it is necessary to determine the area of ​​all the processed triangular face fragments. Based on the parameter values ​​of the contour lines or boundary lines where the vertices of each processed triangular face fragment are located, the parameter values ​​of each vertex of each processed triangular face fragment can be obtained, and then the area of ​​each processed triangular face fragment can be obtained.

[0124] In the embodiment of the present application, when determining the value within the area of ​​each processed triangular face fragment, it is necessary to first determine whether the positions of the vertices of the processed triangular face fragment are on the same contour line, and determine the value according to the determination result according to the value determination principle shown in Table 1 below:

[0125]

[0126] Table 1 Triangle face element value table

[0127] Since all the vertices of the processed triangles will fall on the contour lines or boundary lines, the contour lines are the lines that describe the trend of the parameter values, and the vertex values ​​refer to the parameter values ​​at the vertex positions. In the processed triangles, if the three vertices are on the same contour line, the three points are the same, and the area of ​​the processed triangle is ΔV = V A =V A =V C; If two vertices are on the same contour line or three vertices are not on the same contour line, then the two points are the same or different. It is necessary to determine the value within the area of ​​the processed triangle surface element according to the weight of the quadrilateral formed by the center point of the processed triangle surface element and the midpoints of each side. Figure 7 As shown in the figure, it is a schematic diagram of the center point of the triangle face element, the midpoints of each side, and the quadrilateral formed by the vertices after processing. A , S B and S C are the areas of the quadrilaterals formed by the incenter of the processed triangle and the midpoints of each side, V A 、V B and V C They are the parameter values ​​of the contour lines or boundary lines where the vertices of the processed triangle face fragments are located.

[0128] In the embodiment of the present application, the weighting of the quadrilateral formed by the center point of the processed triangular face element and the midpoints of each side to determine the value within the area of ​​the processed triangular face element specifically includes:

[0129] According to the weighted quadrilateral formed by the center point of the processed triangular face element and the midpoints of each side, the value within the area of ​​the processed triangular face element is determined based on the following formula 1:

[0130]

[0131] Wherein, ΔV is the value within the area of ​​the processed triangle element; S A , S B and S C are the areas of the quadrilaterals formed by the center points of the processed triangle face fragments and the midpoints of each side; V A 、V B and V C They are respectively the parameter values ​​of the contour lines or boundary lines where the vertices of the processed triangle face fragments are located.

[0132] S106: Determine the reserve parameters constrained by the oil and gas area according to the area of ​​each processed triangular facet and the value within the area.

[0133] In the above step S106, the reserve parameters constrained by the oil and gas area are determined according to the area of ​​each processed triangular facet and the value within the area, which specifically includes:

[0134] According to the area of ​​each processed triangular facet and the value within the area, the reserve parameter constrained by the oil and gas area is determined based on the following formula 2:

[0135]

[0136] Wherein, V is the reserve parameter constrained by the oil and gas area; ΔS is the area of ​​the processed triangular face element; ΔV is the value within the area of ​​the processed triangular face element; and N=3.

[0137] In the embodiment of the present application, after obtaining the area and the value within the area of ​​each processed triangular face element, based on the above formula 2, the reserve parameters of the oil and gas area constraint can be obtained, including the effective thickness, effective porosity and oil and gas saturation of the oil and gas area constraint and other reserve parameters.

[0138] In the present application embodiment, Figure 8 As shown in the figure, it is a workflow diagram of the method for calculating reserve parameters based on contour map, wherein the oil and gas area constraint is an input, corresponding to the range POLYGON of the reserve parameter contour trend map; the reserve parameter contour trend map is the contour map of the oil and gas area constraint obtained in advance; the triangular mesh pre-division is the triangular topology identification, and all the initial triangular surface elements are obtained; the objective function convergence is the refinement process of the initial triangular surface elements, and when it does not converge, the trend control point is calculated, that is, the intersection of the contour line and the initial triangular surface element is added as a new control point to the original oil and gas area constraint; the topological relationship processing is to constrain the initial triangular surface elements that cross the contour line; the singular value processing is to delete the initial triangular surface elements whose inner core is not in the oil and gas area.

[0139] In the past, when calculating reserve parameters by the contour area trade-off method, the method of manually drawing POLYGON or manually picking contour lines to close and calculate POLYGON was usually adopted, which has a large error risk. This risk is particularly prominent in oil and gas area areas where boundary lines are involved. More importantly, the workflow is complicated. Because the amount of contour line data is large, the picking workload is huge. If there are many calculation units, the workload caused by any adjustment of related parameters during reserve calculation is huge, and the work efficiency is very low. The method for calculating reserve parameters based on contour maps provided in this application records the complete topological information of the contour area boundary lines within the oil and gas area, so it can realize the correct identification of the contour area boundary lines within the oil and gas area, and provide the contour area distribution parameters in the geological reserve calculation. Compared with the previous contour area calculation, it not only improves the efficiency by several times or even dozens of times, but also respects the input data in terms of quality, achieves the same quality as the input data, and improves the data accuracy.

[0140] The method for calculating reserve parameters based on contour maps provided in the embodiment of the present application performs triangular topology recognition on the contour maps constrained by oil and gas area, splits the contour maps into several triangular facets, and obtains the reserve parameters constrained by oil and gas area by calculating the area of ​​each triangular facet and the value within the area. The method can automatically realize fast and more accurate calculation of reserve parameters without human intervention, and can achieve the purpose of ensuring data accuracy while obtaining extremely high work efficiency, and achieve the technical effect of improving quality and greatly improving timeliness in reserve calculation work, and realize automatic, high-timeliness and high-quality reserve parameter calculation based on contour maps.

[0141] Embodiment 2

[0142] Based on the same inventive concept, the embodiment of the present application also provides a device for calculating reserve parameters based on contour maps, referring to Fig. 9 As shown, the device comprises:

[0143] A closing processing module 101 is used to perform closing processing on the pre-obtained contour map of oil and gas area constraints to obtain a closed contour map;

[0144] A first determination module 102, configured to determine a parameter value of a boundary line according to a parameter value of an intersection point between a contour line and a boundary line in the closed contour map;

[0145] A topology recognition module 103 is used to perform triangular topology recognition on all point information in the closed contour map to obtain all initial triangular facets;

[0146] A thinning module 104 is used to thin each of the initial triangular surface fragments that meets the preset rules to obtain all processed triangular surface fragments;

[0147] A second determination module 105 is used to determine the area and the value within the area of ​​each processed triangular facet according to the parameter values ​​of the contour lines and the boundary lines;

[0148] The third determination module 106 is used to determine the reserve parameters constrained by the oil and gas area according to the area of ​​each processed triangular facet and the value within the area.

[0149] Embodiment 3

[0150] Based on the same inventive concept, an embodiment of the present application further provides a computer-readable storage medium on which a computer program is stored. When the program is executed by a processor, the method for calculating reserve parameters based on a contour map as described in the above-mentioned embodiment 1 is implemented.

[0151] Embodiment 4

[0152] Based on the same inventive concept, an embodiment of the present application also provides a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, the method for calculating reserve parameters based on the contour map as described in the above-mentioned embodiment 1 is implemented.

[0153] Embodiment 5

[0154] Based on the same inventive concept, an embodiment of the present application also provides a computer program product comprising instructions. When the computer program product is run on a computer device, the computer device executes the method for calculating reserve parameters based on contour maps as described in the above-mentioned embodiment 1.

[0155] Embodiment 6

[0156] Based on the same inventive concept, an embodiment of the present application also provides a chip, the chip includes a processor and a communication interface, the communication interface and the processor are coupled, and the processor is used to run a computer program or instruction to implement the method for calculating reserve parameters based on the contour map as described in the above embodiment 1.

[0157] Those skilled in the art will appreciate that embodiments of the present invention may be provided as methods, systems, or computer program products. Therefore, the present invention may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage and optical storage, etc.) containing computer-usable program code.

[0158] The present invention is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of the processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0159] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0160] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process in the computer or other programmable device. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0161] Obviously, those skilled in the art can make various changes and modifications to the present invention without departing from the spirit and scope of the present invention. Thus, if these modifications and variations of the present invention fall within the scope of the claims of the present invention and their equivalents, the present invention is also intended to include these modifications and variations.

Claims

1. A method for calculating reserve parameters based on contour maps, characterized in that: include: Performing closing processing on the pre-obtained contour map of oil and gas area constraints to obtain a closed contour map; Determining the parameter value of the boundary line according to the parameter value of the intersection point of the contour line and the boundary line in the closed contour map; Performing triangular topology recognition on all point information in the closed contour map to obtain all initial triangular facets; Refining the initial triangular surface fragments that meet the preset rules in each of the initial triangular surface fragments to obtain all processed triangular surface fragments; Determine the area and the value within the area of ​​each processed triangular face element according to the parameter values ​​of the contour lines and the boundary lines; According to the area of ​​each processed triangular facet and the value within the area, the reserve parameters constrained by the oil and gas area are determined.

2. The method according to claim 1, characterized in that The method of closing the pre-obtained contour map of oil and gas area constraints to obtain a closed contour map includes: Determine whether each contour line in the contour line map of the oil and gas area constraint obtained in advance has an intersection with any boundary line; If so, delete the remaining line segments after the intersection, determine all intersection points as control points, and obtain a new contour map; If not, extend the contour line until it intersects with the boundary line, determine all intersection points as control points, and obtain a new contour map; Determine whether all curves in the new contour map are closed; If so, the closed contour map is obtained.

3. The method according to claim 1, characterized in that The step of obtaining the parameter value of the boundary line according to the parameter value of the intersection point between the contour line and the boundary line in the closed contour map comprises: The parameter value of the boundary line is obtained by interpolating or extrapolating the parameter value of the intersection point of the contour line and the boundary line in the closed contour map.

4. The method according to claim 1, characterized in that The initial triangular surface elements that meet the preset rules include initial triangular surface elements that cross the contour line and initial triangular surface elements whose center points are not within the oil and gas bearing area; The step of thinning the initial triangular surface fragments that meet the preset rules in each of the initial triangular surface fragments to obtain all processed triangular surface fragments includes: Add the intersection points of the initial triangles that cross the contour lines and the contour lines to the control points; According to the control points, constraining the initial triangular surface elements that cross the contour lines to obtain the processed triangular surface elements; as well as, Deleting the initial triangular facets whose inner parts are not within the oil and gas containing area; Wherein, whether the center of the initial triangular face element is within the oil and gas containing area is determined by the following method: According to the parameter values ​​of the contour lines and the boundary lines, the center point of the initial triangular surface element is determined to judge whether the center point of the initial triangular surface element is within the oil and gas containing area.

5. The method according to claim 4, characterized in that The constraining the initial triangular surface element that crosses the isoline according to the control point to obtain the processed triangular surface element comprises: According to the control points, the area of ​​the triangular surface fragment composed of the three most adjacent control points in the initial triangular surface fragment crossing the contour line is converged to be close to zero, so as to obtain the processed triangular surface fragment.

6. The method according to claim 1, characterized in that Determining the area and the value within the area of ​​each processed triangular face element according to the parameter values ​​of the contour lines and the boundary lines includes: Determining the parameter values ​​of each vertex of each processed triangular face fragment according to the parameter values ​​of the contour lines and the boundary lines; Determining the area of ​​all the processed triangular surface fragments according to the parameter values ​​of each vertex of each processed triangular surface fragment; Determine whether the positions of the vertices of the processed triangular face fragments are on the same contour line; If yes, then the value within the area of ​​the processed triangle face element is the value of the same contour line; If not, the area value of the processed triangular face element is determined according to the weighted quadrilateral formed by the center point of the processed triangular face element and the midpoints of each side.

7. The method according to claim 6, characterized in that The step of determining the value within the area of ​​the processed triangular face element according to the weighted quadrilateral formed by the center point of the processed triangular face element and the midpoints of each side includes: According to the weighted quadrilateral formed by the center point of the processed triangular face element and the midpoints of each side, the value within the area of ​​the processed triangular face element is determined based on the following formula 1: Wherein, ΔV is the value within the area of ​​the processed triangle element; S A , S B and S C are the areas of the quadrilaterals formed by the center points of the processed triangle face fragments and the midpoints of each side; V A 、V B and V C They are respectively the parameter values ​​of the contour lines or boundary lines where the vertices of the processed triangle face fragments are located.

8. The method according to claim 1, characterized in that Determining the reserve parameters constrained by the oil and gas area according to the area of ​​each processed triangular facet and the value within the area includes: According to the area of ​​each processed triangular facet and the value within the area, the reserve parameter constrained by the oil and gas area is determined based on the following formula 2: Wherein, V is the reserve parameter constrained by the oil and gas area; ΔS is the area of ​​the processed triangular face element; ΔV is the value within the area of ​​the processed triangular face element; and N=3.

9. A device for calculating reserve parameters based on contour maps, characterized in that: include: A closing processing module is used to perform closing processing on the pre-obtained contour map of oil and gas area constraints to obtain a closed contour map; A first determination module is used to determine the parameter value of the boundary line according to the parameter value of the intersection point of the contour line and the boundary line in the closed contour map; A topology recognition module is used to perform triangular topology recognition on all point information in the closed contour map to obtain all initial triangular facets; A refinement module, used for refining the initial triangular surface fragments that meet the preset rules in each of the initial triangular surface fragments to obtain all processed triangular surface fragments; A second determination module is used to determine the area and the value within the area of ​​each processed triangular face element according to the parameter values ​​of the contour lines and the boundary lines; The third determination module is used to determine the reserve parameters constrained by the oil and gas area according to the area of ​​each processed triangular facet and the value within the area.

10. A computer-readable storage medium, wherein instructions are stored in the computer-readable storage medium, and when the instructions are executed on a terminal, the terminal executes the method for calculating reserve parameters based on a contour map as described in any one of claims 1 to 8.

11. A computer device, characterized in that: The method comprises a memory, a processor and a computer program stored in the memory and executable on the processor, wherein when the processor executes the computer program, the method for calculating reserve parameters based on contour maps as described in any one of claims 1 to 8 is implemented.

12. A computer program product comprising instructions, which, when executed on a computer device, enables the computer device to execute the method for calculating reserve parameters based on contour maps as claimed in any one of claims 1 to 8.

13. A chip, comprising a processor and a communication interface, wherein the communication interface and the processor are coupled, and the processor is used to run a computer program or instruction to implement the method for calculating reserve parameters based on a contour map as described in any one of claims 1 to 8.

Images