Method and device for calculating area of cross-zone oil-gas-bearing region of oil-gas reserve calculation unit

Abstract

The invention relates to a method and device for calculating the area of a cross-zone oil-gas-containing region of an oil-gas reserve calculation unit, and the method comprises the steps: converting the geodetic coordinates of the oil-gas-containing region into latitude and longitude coordinates when the oil-gas-containing region crosses a zone; determining the center of the oil-gas-containing area according to the minimum longitude and the maximum longitude of the oil-gas-containing area; constructing a region coordinate system for the oil-gas-containing region by using a center, and converting the latitude and longitude coordinates into region coordinates of the oil-gas-containing region; according to the area coordinates of the oil-gas-containing area, the oil-gas-containing area is divided into a plurality of sub-areas; and accumulating and calculating the area of each sub-region to obtain the area of the oil-gas-containing region. According to the method, cross-zone automatic processing and calculation are achieved, the cross-zone removing effect of reserve parameter calculation is achieved, overall reserve parameter calculation is made to return to a calculation mode similar to the mode without cross-zone, the precision and working time efficiency of high-reserve parameter calculation are guaranteed, and automation, intelligence and reliability of reserve calculation are achieved.

Description

Technical Field

[0001] The present specification relates to the field of oil and gas field exploration and development, and in particular to a method and device for calculating the area of ​​cross-zone oil and gas containing regions in an oil and gas reserve calculation unit. Background Art

[0002] Oil and gas reserves are the core assets of oil companies and an important indicator of their healthy development potential. There are various types of oil and gas reserves, including new (predicted reserves, controlled reserves, proven reserves) reserve calculation, reserve recalculation, reserve accounting, reserve calibration and settlement, and SEC reserve estimation. The technical requirements and emphases are different. The calculation accuracy control nodes are numerous and wide-ranging (including calculation models, calculation parameters, and technical methods, etc.). If any of the nodes is not well controlled, the calculation accuracy may cause error transmission, resulting in poor reserve calculation accuracy. Therefore, it is difficult to accurately calculate reserves.

[0003] The delineation and calculation of oil and gas bearing areas are extremely important parameters in reserve calculations, and are also the most sensitive parameters, which are related to the accuracy and reliability of reserve calculations. According to the provisions of industry standards and enterprise standards, the accurate calculation of oil and gas bearing areas involves relatively complex coordinate projection transformations, etc. If the oil and gas bearing area is across zones, the current calculation is very complex considering the calculation accuracy, and is mainly done manually. The usual processing method involves zoning the oil and gas bearing area, and then applying the reserve calculation formula to each area.

[0004] The existing technology uses the zone calculation method to first split the oil and gas area according to the zone, and then carry out the reserve calculation work in each zone. This method is difficult to apply cross-zone parameters when calculating the reserve parameters of the same set of oil and gas reservoirs. The trend of the reserve parameters of the oil and gas reservoirs needs to be carried out in conjunction with the existing data of the entire reservoir. Although the data crosses the zone, it must participate in the overall reserve parameter calculation as part of the overall trend of the oil and gas reservoir; in addition, it is difficult to visualize after crossing the zone, and the entire workflow is cumbersome and the workload is huge. Moreover, the manual process cannot automatically transmit the changes in 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. Summary of the invention

[0005] In order to solve the problems of the prior art, the embodiments of this specification provide a method and device for calculating the area of ​​cross-zone oil and gas containing regions in an oil and gas reserve calculation unit.

[0006] An embodiment of the present specification provides a method for calculating the area of ​​cross-zone oil and gas containing regions in an oil and gas reserve calculation unit, the method comprising: when the oil and gas containing region crosses zones, converting the geodetic coordinates of the oil and gas containing region into longitude and latitude coordinates; determining the center of the oil and gas containing region based on the minimum longitude and maximum longitude of the oil and gas containing region; constructing a regional coordinate system for the oil and gas containing region using the center, and converting the longitude and latitude coordinates into regional coordinates of the oil and gas containing region; dividing the oil and gas containing region into a plurality of sub-regions based on the regional coordinates of the oil and gas containing region; and accumulating and calculating the area of ​​each sub-region to obtain the area of ​​the oil and gas containing region.

[0007] According to one aspect of an embodiment of the present specification, a regional coordinate system is constructed for the oil and gas containing region using the center, and converting the longitude and latitude coordinates into regional coordinates for the oil and gas containing region includes: using the vertical line where the center is located as the central meridian of the regional coordinate system; and converting the longitude and latitude coordinates of the oil and gas containing region into regional coordinates.

[0008] According to one aspect of an embodiment of the present specification, dividing the oil and gas containing area into a plurality of sub-areas according to the regional coordinates of the oil and gas containing area includes: selecting longitude and latitude points corresponding to the minimum angle on the boundary of the oil and gas containing area, and using the longitude and latitude points as corner points; dividing the boundary of the oil and gas containing area into a plurality of line segments according to the corner points; and dividing the line segments according to the number of proportional subdivisions of the space in the longitude and latitude directions, thereby dividing the oil and gas containing area into a plurality of sub-areas.

[0009] According to one aspect of an embodiment of the present specification, the method further includes: determining the longitude and latitude coordinates of the vertices in the sub-region according to the subdivision grid coordinates in the proportional subdivision space; and calculating the area of ​​the sub-region according to the longitude and latitude coordinates of the vertices.

[0010] According to one aspect of an embodiment of the present specification, the sub-region is a triangle, and accumulating the area of ​​each sub-region includes: determining the distance between any two points of the triangle based on the longitude and latitude coordinates of the triangle to obtain the length of each side of the triangle; and determining the area of ​​each triangle based on the side length.

[0011] According to one aspect of an embodiment of the present specification, after calculating the area of ​​the oil and gas containing region, the method further includes: converting the regional coordinates of the oil and gas containing region into longitude and latitude coordinates; and converting the longitude and latitude coordinates into geodetic coordinates.

[0012] According to one aspect of the embodiments of this specification, whether the oil and gas containing area crosses the belt is determined in the following manner: based on the geodetic coordinates corresponding to all points on the boundary of the oil and gas containing area, the belt number corresponding to the geodetic coordinates of each point is determined; if the belt numbers corresponding to the geodetic coordinates of all points are inconsistent, it is determined that the oil and gas containing area crosses the belt.

[0013] An embodiment of the present specification provides an apparatus for calculating the area of ​​cross-zone oil and gas containing regions in an oil and gas reserve calculation unit, the apparatus comprising: a coordinate conversion unit, for converting the geodetic coordinates of the oil and gas containing region into longitude and latitude coordinates when the oil and gas containing region crosses zones; a center determination unit, for determining the center of the oil and gas containing region based on the minimum longitude and maximum longitude of the oil and gas containing region; a regional coordinate construction unit, for constructing a regional coordinate system for the oil and gas containing region using the center, and converting the longitude and latitude coordinates into regional coordinates of the oil and gas containing region; a segmentation unit, for dividing the oil and gas containing region into a plurality of sub-regions based on the regional coordinates of the oil and gas containing region; and an area calculation unit, for accumulating and calculating the area of ​​each sub-region to obtain the area of ​​the oil and gas containing region.

[0014] The embodiments of the present specification also provide a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the shale full-diameter core resistivity anisotropy test acquisition method is implemented.

[0015] This manual realizes automatic cross-band processing and calculation, realizes the "cross-band removal" effect of reserve parameter calculation, and makes the overall reserve parameter calculation return to the calculation method similar to the case without cross-band, ensuring the accuracy and timeliness of high reserve parameter calculation, which is conducive to realizing the automation, intelligence and reliability of reserve calculation. BRIEF DESCRIPTION OF THE DRAWINGS

[0016] In order to more clearly illustrate the embodiments of this specification or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of this specification. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0017] Figure 1 It is a schematic diagram of a method for calculating the area of ​​cross-zone oil and gas containing regions of an oil and gas reserve calculation unit according to an embodiment of this specification;

[0018] Figure 2 Shown is a flow chart of a method for determining the regional coordinates of an oil and gas-bearing area according to an embodiment of this specification;

[0019] Figure 3 The figure shows a flow chart of a method for dividing an oil and gas-bearing area according to an embodiment of the present specification;

[0020] Figure 4 The figure is a flow chart of a method for calculating the area of ​​a triangle according to an embodiment of the present specification;

[0021] Figure 5The figure is a flow chart of a method for determining the area of ​​a triangle according to an embodiment of the present specification;

[0022] Figure 6 The figure is a flow chart of a coordinate conversion method according to an embodiment of the present specification;

[0023] Figure 7 The figure shows a schematic diagram of the structure of a device for calculating the area of ​​oil and gas containing regions across a belt in an oil and gas reserve calculation unit according to an embodiment of the present specification;

[0024] Figure 8 The figure shows a schematic diagram of the coordinate space mapping between a proportional subdivision space and an oil and gas containing area in an embodiment of this specification;

[0025] Fig. 9 The figure shows a schematic diagram of the structure of a computer device according to an embodiment of the present specification.

[0026] Description of the accompanying symbols:

[0027] 701, coordinate conversion unit;

[0028] 702. Center determination unit;

[0029] 703. Regional coordinate construction unit;

[0030] 704, segmentation unit;

[0031] 705. Area calculation unit;

[0032] 902. Computer equipment;

[0033] 904, processor;

[0034] 906. Memory;

[0035] 908, driving mechanism;

[0036] 910, input / output module;

[0037] 912. Input device;

[0038] 914. Output device;

[0039] 916. Presentation equipment;

[0040] 918. Graphical user interface;

[0041] 920, network interface;

[0042] 922, communication link;

[0043] 924. Communication bus. DETAILED DESCRIPTION

[0044] In order to enable those skilled in the art to better understand the technical solutions in this specification, the technical solutions in the embodiments of this specification will be clearly and completely described below in conjunction with the drawings in the embodiments of this specification. Obviously, the described embodiments are only part of the embodiments of this specification, not all of the embodiments. Based on the embodiments in this specification, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of this specification.

[0045] It should be noted that the terms "first", "second", etc. in the description and claims of this specification and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequence. It should be understood that the data used in this way can be interchanged where appropriate, so that the embodiments of this specification described herein can be implemented in an order other than those illustrated or described herein. In addition, the terms "including" and "having" and any variations thereof are intended to cover non-exclusive inclusions, for example, a process, method, device, product or equipment that includes a series of steps or units is not necessarily limited to those steps or units clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or equipment.

[0046] This specification provides method operation steps as described in the embodiments or flow charts, but more or fewer operation steps may be included based on routine or non-creative work. The order of steps listed in the embodiments is only one way of executing the order of many steps and does not represent the only execution order. When the system or device product is executed in practice, it can be executed in the order of the method shown in the embodiments or the drawings or in parallel.

[0047] It should be noted that the method and device for calculating the area of ​​cross-zone oil and gas containing regions in the oil and gas reserve calculation unit of this specification can be used in the field of oil and gas field exploration and development. This specification does not limit the application field of the method and device for calculating the area of ​​cross-zone oil and gas containing regions in the oil and gas reserve calculation unit.

[0048] Figure 1 The flowchart of a method for calculating the area of ​​cross-zone oil and gas containing regions of an oil and gas reserve calculation unit according to an embodiment of the present specification is shown. The method includes:

[0049] Step 101, when the oil and gas-bearing region crosses the belt, the geodetic coordinates of the oil and gas-bearing region are converted into longitude and latitude coordinates. In the process of oil and gas exploration and development, due to the complexity of geological structure, the oil and gas-bearing region may cross different geological belts or structural units. Therefore, this step first determines whether the oil and gas-bearing region crosses the geological belt, that is, whether the cross-belt phenomenon occurs.

[0050] Specifically, according to the geodetic coordinates corresponding to all points on the boundary of the oil and gas bearing area, determine the belt number corresponding to the geodetic coordinates of each point. In this way, it is determined whether all points in the oil and gas bearing area are located in the same geological belt. Obtain all points in the oil and gas bearing area, each of which has a corresponding belt number. Determine whether all belt numbers of an oil and gas bearing area are consistent. If they are inconsistent, it can be determined that the oil and gas bearing area crosses the belt, and cross-belt calculation is required. Among them, the geodetic coordinate points of each point in the known oil and gas bearing area are point sequences composed of coordinates X and coordinates Y, respectively {(X i ,Y i )}, i = 1, 2, ..., N. The following formula can be used to calculate the zone number corresponding to each geodetic coordinate point: zone = x / 1000000, where x represents the geodetic coordinate of a point in the oil and gas-bearing area.

[0051] After determining the cross-zone of the oil and gas bearing area, the geodetic coordinates of the oil and gas bearing area are converted into longitude and latitude coordinates. Specifically, it is known that a represents the major semi-axis and f represents the flattening, and the latitude and longitude B and L of the oil and gas bearing area are calculated.

[0052] represents the minor semi-axis of the ellipse,

[0053] represents the first eccentricity,

[0054] represents the second eccentricity,

[0055] Noffset=0, indicating the north offset of domestic geodetic coordinates.

[0056] Eoffset=500000, indicating the easting offset of domestic geodetic coordinates.

[0057] L0 = zone*6-3, indicating the central meridian divided by 6 degrees.

[0058]

[0059] Where, x = X-Eoffset, y = Y-Noffset;

[0060] m0=a·(1-e 2 )

[0061]

[0062]

[0063]

[0064]

[0065]

[0066]

[0067]

[0068]

[0069]

[0070] Use the above iterative formula to calculate B f , the convergence condition is

[0071] T f =tan(B f )

[0072] ∈ 2 =e′ 2 cos(B f ) 2

[0073]

[0074]

[0075]

[0076]

[0077]

[0078] B represents the latitude of the calculation point, L represents the longitude of the calculation point, X represents the geodetic east coordinate, Y represents the geodetic north coordinate, and other parameters B f , T f ,∈,W f 、M f m0, m2, m4, m6, m8, a0, a2, a4, a6, a8 are all intermediate variables with no actual meaning. f Represents the radius of curvature of the camber circle.

[0079] Step 102, determining the center of the oil and gas containing area according to the minimum longitude and the maximum longitude of the oil and gas containing area.

[0080] Because the oil and gas-bearing area involves cross-zones, the geodetic coordinates of the oil and gas-bearing area do not have a unified coordinate system in the two zones. Therefore, in order to maintain the calculation accuracy, it is necessary to build a set of regional coordinate systems in which the coordinate points have a unified measurement system. Specifically, to establish a regional coordinate system, we must first calculate the central longitude of the current oil and gas-bearing area. The center of the oil and gas-bearing area is determined by the following formula:

[0081] L0 ′ =[min(L i )+max(L i )] / 2; where max(L i ) represents the maximum longitude of the oil and gas bearing area, min(L i ) represents the minimum longitude of the oil and gas bearing area, L i Represents the longitude of the i-th point in the oil and gas bearing area.

[0082] Step 103, using the center to construct a regional coordinate system for the oil and gas containing region, and converting the longitude and latitude coordinates into regional coordinates of the oil and gas containing region.

[0083] Specifically, the center of the oil and gas-bearing area determined in step 102 is used as the central meridian of the regional coordinate system to be constructed. Further, the longitude and latitude coordinates of the oil and gas-bearing area are converted into regional coordinates. For a detailed description of converting longitude and latitude coordinates into regional coordinates, see Figure 2 Description, this application will not be repeated here.

[0084] This step realizes the conversion of the cross-zone oil and gas-bearing areas into the same area, avoiding the calculation error of the oil and gas-bearing areas in different coordinate systems due to the cross-zone. Therefore, the oil and gas-bearing areas can be visualized with the same regional coordinates and displayed in one plane.

[0085] Step 104: divide the oil and gas containing region into a plurality of sub-regions according to the region coordinates of the oil and gas containing region.

[0086] In order to accurately calculate the area of ​​the oil and gas-bearing region, it is necessary to perform finite element segmentation on the oil and gas-bearing surface area and divide it into finer sub-areas to improve the area calculation accuracy of the oil and gas-bearing region.

[0087] In the embodiments of this specification, a uniformly proportioned subdivision space is used to map with the physical space of the oil and gas containing area, and the oil and gas containing area is firstly divided by dividing the uniformly proportioned subdivision space into a plurality of uniform grids.

[0088] In this step, multiple segments are evenly divided in the longitude and latitude directions of the uniformly proportioned subdivision space, thereby dividing the uniformly proportioned subdivision space into multiple grids. At the same time, multiple points are selected on the boundary of the oil and gas-bearing region, and the boundary of the oil and gas-bearing region is divided into multiple edges according to the multiple points. Each edge is divided according to the number of segments in the longitude and latitude directions of the uniformly proportioned subdivision space, and the oil and gas-bearing region is divided into multiple sub-areas.

[0089] For example, the uniformly proportioned subdivision space is divided into 100 by 100 grids, so that the coordinates of each grid can be determined. Correspondingly, the oil and gas-bearing area can also be divided into a corresponding number of grids, and the coordinates of the grids in the oil and gas-bearing area can be determined according to the coordinates of each grid in the uniformly proportioned subdivision space. For a detailed description of dividing the oil and gas-bearing area into multiple sub-areas, see Figure 8 .

[0090] Step 105, cumulatively calculating the area of ​​each sub-region to obtain the area of ​​the oil and gas containing region.

[0091] In this step, the area of ​​each sub-region is calculated according to the coordinates of each sub-region. Finally, the areas of all sub-regions in the oil and gas-bearing region are added together to obtain the area of ​​the oil and gas-bearing region. In this step, the oil and gas-bearing region is divided into finer sub-regions, which can improve the calculation accuracy of the area of ​​the oil and gas-bearing region.

[0092] This manual comprehensively considers the cross-band coordinate conversion of the three-level coordinate system fusion, effectively supports the realization of high-precision planar visualization and interactive processing of oil and gas-bearing areas, forms a complete automated processing flow, and realizes user-unaware cross-band calculation, relative position calculation and oil and gas area calculation.

[0093] Figure 2 The flowchart of a method for determining the regional coordinates of an oil and gas-bearing area according to an embodiment of this specification is shown, which specifically includes the following steps:

[0094] Step 201, taking the vertical line where the center is located as the central meridian of the regional coordinate system. Specifically, taking the center of the oil and gas bearing area determined in step 102 as the central meridian of the regional coordinate system to be constructed, further constructing the regional coordinate system according to the center of the oil and gas bearing area, thereby ensuring that each point of the oil and gas bearing area in the regional coordinate system is in the same zone without the phenomenon of crossing zones.

[0095] Step 202, converting the longitude and latitude coordinates of the oil and gas-bearing area into regional coordinates. In this step, the regional coordinates of each point in the oil and gas-bearing area are calculated according to the following formula:

[0096]

[0097] Among them, Y' represents the regional coordinates of a calculation point in the oil and gas area. R1 represents the intermediate quantity, a represents the major semi-axis, B represents the latitude of the calculation point, and e represents the first eccentricity;

[0098] represents the intermediate variable, where R2 represents the intermediate quantity, B represents the latitude of the calculation point, and a0, a2, a4, a6, and a8 are calculated by the following formulas:

[0099]

[0100]

[0101]

[0102]

[0103] Among them, m0, m2, m4, m6, and m8 are obtained by the following formulas:

[0104] m0=a·(1-e 2 )

[0105]

[0106]

[0107]

[0108]

[0109] Furthermore, the regional coordinates of the calculation points in the oil and gas-bearing area are calculated according to the following formula:

[0110]

[0111] Among them, R1 represents the intermediate quantity, B represents the latitude of the calculation point, τ 2 =e ′2 cos(B) 2 , τ represents an intermediate quantity and has no specific meaning.

[0112] After the regional coordinate system is constructed in this step, the regional coordinates of the oil and gas-bearing area can be directly calculated accurately because of the unified measurement system. The user can pick up the regional coordinate system at a certain point in the visualization system, and can also obtain the geodetic coordinates and longitude and latitude coordinates at the same time.

[0113] Figure 3 The flowchart of a method for dividing oil and gas-bearing regions in an embodiment of this specification is shown, which specifically includes the following steps:

[0114] Step 301, selecting the longitude and latitude points corresponding to the minimum angle on the boundary of the oil and gas containing area, and taking the longitude and latitude points as corner points.

[0115] In the embodiments of this specification, the boundary of the oil and gas region is an ordered closed curve, which is composed of a sequence of multiple points. The longitude and latitude coordinates of each point are expressed as follows: {(B1, L1), (B2, L2), ... (BN , L N )}. In this step, the i-th point is randomly selected on the closed curve of the boundary of the oil and gas region, and two consecutive points before and after the point are taken to obtain three consecutive points, which are expressed as: (B i-1 , L i-1 ), (B i , L i ), (B i+1 , L i+1 ). Calculate the ith point as the starting point and construct two vectors with the other two points in the three points, which are expressed as: vector C = (B i-1 -B i , L i-1 -L i ), D=(B i+1 -B i , L i+1 -L i ). And further calculate the angle θ between the two vectors. The calculation formula is as follows:

[0116]

[0117] Furthermore, the angles formed by other multiple points on the closed curve of the oil and gas area boundary are calculated in the same way, and the four longitude and latitude points with the smallest angles are selected as the corner points of the closed curve of the oil and gas area boundary in the order of angle size, which are expressed as: (B c0 , L c0 ), (B c1 , L c1 ), (B c2 , L c2 ), (B c3 , L c3 ) to automatically obtain corner points.

[0118] Before executing this step, the method further includes: converting the regional coordinates of the oil and gas-bearing area into longitude and latitude coordinates. The conversion process converts the geodetic coordinates into the inverse of longitude and latitude coordinates, and this step will not be described in detail here.

[0119] Step 302: Divide the boundary of the oil and gas-bearing region into a plurality of line segments according to the corner points.

[0120] In this step, the four points are taken as the intersection points of the boundaries of the oil and gas-bearing area, and the closed curve is divided into four line segments, which can be recorded as I0, I1, I2 and I3. Figure 8 shown.

[0121] Step 303, dividing the line segments respectively according to the number of subdivided grids in the longitude and latitude directions of the proportional subdivided space, and dividing the oil and gas bearing area into a plurality of sub-areas.

[0122] In this step, the proportional subdivision space is Figure 8 As shown, it is a mathematical space that has a mapping relationship with the physical space of the oil and gas-bearing area. In the proportional subdivision space, the number of subdivision grids in the longitude and latitude directions can be set to N1 and N2 respectively. If the length of the proportional subdivision space in the longitude and latitude directions is regarded as 1, the coordinates of the subdivision grid can be expressed as: The unit length of the subdivision grid in the longitude direction of the proportional subdivision space is The unit length of the subdivision grid in the longitude direction of the proportional subdivision space is

[0123] The fragments after the boundary division of the oil and gas-bearing area are also divided according to the number of divisions of the proportional subdivision space in the longitude and latitude directions, forming N1×N2 subdivision grids, and the subdivision grids are quadrilaterals. Connecting the two diagonal points of each subdivision grid, each quadrilateral subdivision grid can be split into two triangular face elements, and the triangular face element is the sub-area of ​​the oil and gas-bearing area, then the proportional subdivision space as a whole forms 2×N1×N2 triangular face elements.

[0124] In this step, the coordinates of a certain point on the subdivided grid of the proportional subdivided space and the mapping relationship between the proportional subdivided space and the physical space of the oil and gas region can be determined. For example, the coordinates (γ, ε) of a certain point on the subdivided grid can be converted into the coordinates corresponding to the oil and gas region. Specifically, the coordinates of the two points corresponding to γ ​​are calculated proportionally on the I0 and I2 line segments, and the coordinates of the two points corresponding to ε are calculated proportionally on the I1 and I3 line segments, which can be expressed as: (B l0 , L l0 ), (B l1 , L l1 ), (B l2 , L l2 ), (B l3 , L l3 ).

[0125] Figure 4 The flowchart of a method for calculating the area of ​​a triangle according to an embodiment of the present specification is shown, which specifically includes the following steps:

[0126] Step 401, determining the longitude and latitude coordinates of the vertices in the sub-region according to the subdivision grid coordinates in the proportional subdivision space.

[0127] Specifically, first determine the coordinates of a point in the subdivided grid, such as (γ, ε), and further use the following formula to calculate the latitude B and longitude coordinates L of the subdivided grid point and the corresponding sub-region in the oil and gas-bearing region.

[0128] The specific calculation formula is as follows:

[0129] B=(1-γ)·B l1 +γ·B l0 +(1-ε)·B l2 +ε·B l3 -(γ·ε·B c1 +ε·(1-γ)·B c3 +γ·(1-ε)·B c0 +(1-ε)(1-γ)·B c2

[0130] L=(1-γ)·L l1 +γ·L l0 +(1-ε)·L l2 +ε·L l3 -(γ·ε·L c1 +ε·(1-γ)·L c3 +γ·(1-ε)·L c0 +(1-ε)(1-γ)·L c2

[0131] Among them, γ represents the horizontal coordinate of a point in the proportional subdivision space, ε represents the vertical coordinate of the point in the proportional subdivision space, B represents latitude, L represents longitude, and B l1 , B l0 , B l3 Indicates different edges, (B l0 , L l0 ), (B l1 , L l1 ),,(B l2 , L l2 ), (B l3 , L l3 ) represents the four corners of the sub-region, or represents the four vertices of the sub-region. Figure 8 As shown, the four vertices of the sub-area are: (B c1 , L c1 ), (B c2 , L c2 ), (B c3 , L c3 )(B c4 , L c4 ).

[0132] Step 402, calculating the area of ​​the sub-region according to the longitude and latitude coordinates of the vertex.

[0133] In this step, the area of ​​the triangular face elements is calculated one by one according to the area calculation formula, and the sum of all the triangular face elements is accumulated, so that the high-precision oil and gas-bearing area can be calculated. For a detailed description of the area calculation formula, see Figure 5 .

[0134] Figure 5 The flowchart of a method for determining the area of ​​a triangle according to an embodiment of the present specification is shown, which specifically includes the following steps:

[0135] Step 501, determine the distance between any two points of the triangle according to the longitude and latitude coordinates of the triangle, and obtain the length of each side of the triangle. Calculate the distance between two points on the earth ellipsoid according to the following formula:

[0136] Assume that the longitude and latitude coordinates of two points on the Earth ellipsoid are (B1, L1) and (B2, L2).

[0137] λ = ΔL = L2-L1; λ represents the latitude difference between two points;

[0138] tanU1=(1-f)tan(B1);

[0139] tanU2=(1-f)tan(B2);

[0140]

[0141] sinU1=tanU1*cosU1;

[0142]

[0143] sinU2=tanU2*cosU2; the above tanU1, tanU2, cosU1, sinU1, cosU2, sinU2 all represent trigonometric functions related to the longitudes of two points.

[0144] The iterative method is further used to calculate the following trigonometric functions related to the longitude and latitude of the two points:

[0145] sinσ 2 =(cosU2·sinλ) 2 +(cosU1·sinu2-sinU1·cosU2·cosλ) 2 ;

[0146]

[0147] Cosσ=sinU1·sinU2+cosU1·cosU2·cosλ;

[0148] σ=atan2(sino,cosσ); where σ represents the angular distance between two points on the sphere;

[0149]

[0150] cosα 2 =1-sinα 2 ;

[0151] Among them, sinσ, cosU2, sinλ all represent trigonometric functions related to the longitude of two points.

[0152] Among them, c represents the intermediate quantity, σ M It represents the spherical angular distance from the second finger to the midpoint of the two geodes;

[0153] λ i+1 =ΔL+(1-c)·f·sinα·{σ+c·sinσ·[cos2σ M +c·cosσ·(-1+2cos2σ M 2 )]};

[0154] Until the λ value of the two calculations is less than 10 -12 Until convergence.

[0155] u 2 =cosα 2 ·e′ 2 ;

[0156] Among them, A represents the intermediate amount;

[0157] Among them, B represents the intermediate amount;

[0158]

[0159] Where B represents latitude, Δσ represents the correction term for the length of the geodesic line, which is used to correct the difference between the length of the geodesic line and the simple straight-line distance caused by the curvature and flattening of the earth, and σ M represents the spherical angular distance from the equator to the midpoint of the two geodesic points, and σ represents the angular distance between two points on the sphere;

[0160] s=b·A·(σ-Δσ), where s represents the precise distance between two points on the Earth ellipsoid, σ represents the angular distance between two points on the spherical surface, and Δσ represents the correction term for the geoid length, which is used to correct the difference between the geoid length and the simple straight-line distance caused by the curvature and flattening of the Earth.

[0161] Using the above formula, we can calculate the distance between the three vertices of a triangle and the distance between any two vertices, thereby calculating the length of the three sides of the triangle.

[0162] Step 502: Determine the area of ​​each triangle according to the side length.

[0163] Use the following formula to calculate the area of ​​a triangle:

[0164]

[0165] Among them, AΔ represents the area of ​​the ground triangle, s1, s2, s3 represent the lengths of the sides of the triangle respectively,

[0166] After calculating the area of ​​each cut triangle, the total area of ​​the oil and gas-bearing region can be calculated.

[0167] Figure 6 The figure is a flow chart of a coordinate conversion method according to an embodiment of the present specification, which specifically includes the following steps:

[0168] Step 601, convert the regional coordinates of the oil and gas-bearing area into longitude and latitude coordinates. The specific conversion formula is as follows:

[0169]

[0170]

[0171] Among them, X' represents the regional abscissa of the point on the boundary of the oil and gas area, L0' represents the central meridian of the center point of the oil and gas area, and B f , T f ,∈ is an intermediate variable, N f Represents the radius of curvature of the camber circle.

[0172] Step 602 , converting the longitude and latitude coordinates into geodetic coordinates. In this step, the geodetic coordinates corresponding to the calculation point in the oil and gas-bearing region are calculated according to the following formula, which is similar to the formula in step 202 .

[0173]

[0174] Among them, Noffset = 0: domestic geodetic coordinate northing offset, Eoffset = 500000: domestic geodetic coordinate easting offset. R1 represents the intermediate quantity, B represents the latitude of the calculation point, and τ represents the intermediate quantity.

[0175] Some embodiments of this specification also include: when the boundary of the oil and gas area changes or is adjusted, it is necessary to convert the boundary area coordinates of the oil and gas area back to longitude and latitude coordinates. Specifically, the following steps are included:

[0176] Converting the regional coordinates of the oil and gas-bearing area into longitude and latitude coordinates;

[0177] Convert latitude and longitude coordinates to geodetic coordinates.

[0178] The conversion formula for converting the regional coordinates of the oil and gas region into longitude and latitude coordinates is the same as that in step 601; the formula for converting the longitude and latitude coordinates into the geodetic coordinates corresponding to the calculation point in the oil and gas region is similar to the formula in step 602. The embodiments of this specification are not described in detail here.

[0179] like Figure 7 The figure shows a schematic diagram of the structure of a device for calculating the area of ​​oil and gas containing regions across zones in an oil and gas reserve calculation unit according to an embodiment of the present specification. The figure describes the basic structure of the device for calculating the area of ​​oil and gas containing regions across zones in an oil and gas reserve calculation unit. The functional units and modules therein can be implemented in software or in general chips or specific chips. The device specifically includes:

[0180] A coordinate conversion unit 701 is used to convert the geodetic coordinates of the oil and gas-bearing area into longitude and latitude coordinates when the oil and gas-bearing area crosses the belt;

[0181] A center determination unit 702, configured to determine the center of the oil and gas containing area according to the minimum longitude and the maximum longitude of the oil and gas containing area;

[0182] A regional coordinate construction unit 703, configured to construct a regional coordinate system for the oil and gas containing region using the center, and convert the longitude and latitude coordinates into regional coordinates of the oil and gas containing region;

[0183] A segmentation unit 704 is used to segment the oil and gas containing region into a plurality of sub-regions according to the regional coordinates of the oil and gas containing region;

[0184] The area calculation unit 705 is used to cumulatively calculate the area of ​​each sub-region to obtain the area of ​​the oil and gas containing region.

[0185] Figure 8 The figure shows a schematic diagram of the coordinate space mapping between a proportional subdivision space and an oil and gas containing area according to an embodiment of the present specification.

[0186] The left side of the figure is the digital space, and the right side is the oil and gas area space. The proportional subdivision space has a mapping relationship with the physical space of the oil and gas area. In the longitude and latitude directions of the proportional subdivision space, the number of subdivision grids is set respectively to determine the coordinates of each subdivision point in the digital space.

[0187] According to the longitude and latitude coordinates of each point on the boundary of the oil and gas area, four points on the boundary of the oil and gas area are determined, which are: (B c1 ,L c1 ),(B c2 ,L c2 ),(B c3 ,L c3 )(B c4 ,L c4), and divide the closed curve into 4 line segments, which can be recorded as I0, I1, I2 and I3. According to the number of subdivided grids in longitude and latitude in the proportional subdivision space, the line segments of the oil and gas bearing area are also subdivided to obtain the oil and gas bearing area space as shown in the figure.

[0188] like Fig. 9 As shown, it is a schematic diagram of a computer device provided by an embodiment of this specification. The oil and gas reserve calculation unit cross-zone oil and gas area calculation method described in this application can be applied to the computer device. The computer device 902 may include one or more processors 904, such as one or more central processing units (CPUs), and each processing unit may implement one or more hardware threads. The computer device 902 may also include any memory 906, which is used to store any kind of information such as code, settings, data, etc. Non-limiting, for example, the memory 906 may include any one or more combinations of the following: any type of RAM, any type of ROM, flash memory device, hard disk, optical disk, etc. More generally, any memory can use any technology to store information. Further, any memory can provide volatile or non-volatile retention of information. Further, any memory can represent a fixed or removable component of the computer device 902. In one case, when the processor 904 executes an associated instruction stored in any memory or a combination of memories, the computer device 902 can perform any operation of the associated instruction. The computer device 902 also includes one or more drive mechanisms 908 for interacting with any storage, such as a hard disk drive mechanism, an optical disk drive mechanism, etc.

[0189] The computer device 902 may also include an input / output module 910 (I / O) for receiving various inputs (via input devices 912) and for providing various outputs (via output devices 914). A specific output mechanism may include a presentation device 916 and an associated graphical user interface (GUI) 918. In other embodiments, the input / output module 910 (I / O), the input device 912, and the output device 914 may not be included, and the computer device 902 may be used as a computer device in a network. The computer device 902 may also include one or more network interfaces 920 for exchanging data with other devices via one or more communication links 922. One or more communication buses 924 couple the components described above together.

[0190] The communication link 922 may be implemented in any manner, for example, through a local area network, a wide area network (e.g., the Internet), a point-to-point connection, etc., or any combination thereof. The communication link 922 may include any combination of hardwired links, wireless links, routers, gateway functions, name servers, etc., governed by any protocol or combination of protocols.

[0191] Corresponds to Figures 1 to 6 The method in the embodiment of the present specification also provides a computer-readable storage medium, on which a computer program is stored, and when the computer program is executed by a processor, the steps of the above method are executed.

[0192] The embodiment of the present specification also provides a computer-readable instruction, wherein when the processor executes the instruction, the program therein causes the processor to execute the following Figures 1 to 6 The method shown.

[0193] It should be understood that in the various embodiments of this specification, the size of the serial numbers of the above-mentioned processes 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 this specification.

[0194] It should also be understood that in the embodiments of this specification, the term "and / or" is only a description of the association relationship of the associated objects, indicating that three relationships can exist. For example, A and / or B can represent: A exists alone, A and B exist at the same time, and B exists alone. In addition, the character " / " in this specification generally indicates that the associated objects before and after are in an "or" relationship.

[0195] Those of ordinary skill in the art will appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed in this specification can be implemented in electronic hardware, computer software, or a combination of the two. In order to clearly illustrate the interchangeability of hardware and software, the composition and steps of each example have been generally described in the above description according to function. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professional and technical personnel can use different methods to implement the described functions for each specific application, but such implementation should not be considered to be beyond the scope of this specification.

[0196] Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working processes of the systems, devices and units described above can refer to the corresponding processes in the aforementioned method embodiments and will not be repeated here.

[0197] In the several embodiments provided in this specification, it should be understood that the disclosed systems, devices and methods can be implemented in other ways. For example, the device embodiments described above are only schematic. For example, the division of the units is only a logical function division. There may be other division methods in actual implementation, such as multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. In addition, the mutual coupling or direct coupling or communication connection shown or discussed can be an indirect coupling or communication connection through some interfaces, devices or units, or it can be an electrical, mechanical or other form of connection.

[0198] The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed on multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the embodiments of this specification.

[0199] In addition, each functional unit in each embodiment of this specification may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit. The above integrated unit may be implemented in the form of hardware or in the form of software functional units.

[0200] If the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of this specification is essentially or the part that contributes to the prior art, or all or part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including a number of instructions to enable a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the method described in each embodiment of this specification. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), disk or optical disk and other media that can store program codes.

[0201] Specific embodiments are used in this specification to illustrate the principles and implementation methods of this specification. The description of the above embodiments is only used to help understand the methods and core ideas of this specification. At the same time, for those skilled in the art, according to the ideas of this specification, there will be changes in the specific implementation methods and application scope. In summary, the content of this specification should not be understood as a limitation on this specification.

Claims

1. A method for calculating the area of ​​cross-zone oil and gas containing regions in an oil and gas reserve calculation unit, characterized in that: The method comprises: When the oil and gas-bearing area crosses the belt, the geodetic coordinates of the oil and gas-bearing area are converted into longitude and latitude coordinates; Determining the center of the oil and gas containing area according to the minimum longitude and the maximum longitude of the oil and gas containing area; Using the center to construct a regional coordinate system for the oil and gas containing area, converting the longitude and latitude coordinates into regional coordinates of the oil and gas containing area; Dividing the oil and gas-bearing region into a plurality of sub-regions according to the regional coordinates of the oil and gas-bearing region; The area of ​​each sub-region is cumulatively calculated to obtain the area of ​​the oil and gas containing region.

2. The method according to claim 1, characterized in that Using the center to construct a regional coordinate system for the oil and gas-bearing region, converting the latitude and longitude coordinates into regional coordinates of the oil and gas-bearing region includes: Taking the vertical line where the center is located as the central meridian of the regional coordinate system; Convert the longitude and latitude coordinates of the oil and gas bearing area into regional coordinates.

3. The method according to claim 1, characterized in that According to the regional coordinates of the oil and gas-bearing region, the oil and gas-bearing region is divided into a plurality of sub-regions including: Selecting the longitude and latitude points corresponding to the minimum angle on the boundary of the oil and gas-bearing area, and using the longitude and latitude points as corner points; Dividing the boundary of the oil and gas-bearing region into a plurality of line segments according to the corner points; According to the number of subdivided grids in the longitude and latitude directions of the proportional subdivided space, the line segments are divided respectively to divide the oil and gas bearing area into a plurality of sub-areas.

4. The method according to claim 3, characterized in that The method further comprises: Determine the longitude and latitude coordinates of the vertices in the sub-region according to the subdivision grid coordinates in the proportional subdivision space; Calculate the area of ​​the sub-region based on the latitude and longitude coordinates of the vertex.

5. The method according to claim 4, characterized in that The sub-region is a triangle, and the cumulative calculation of the area of ​​each sub-region includes: According to the longitude and latitude coordinates of the triangle, determine the distance between any two points of the triangle and obtain the length of each side of the triangle; Based on the side lengths, determine the area of ​​each triangle.

6. The method according to claim 1, characterized in that After calculating the area of ​​the oil and gas-bearing region, the method further includes: Converting the regional coordinates of the oil and gas-bearing area into longitude and latitude coordinates; Convert latitude and longitude coordinates to geodetic coordinates.

7. The method according to claim 6, characterized in that The method comprises: According to the geodetic coordinates corresponding to all points on the boundary of the oil and gas-bearing area, determine the belt number corresponding to the geodetic coordinates of each point; If the zone numbers corresponding to the geodetic coordinates of all points are inconsistent, it is determined that the oil and gas-bearing area crosses the zone.

8. A device for calculating the area of ​​cross-zone oil and gas containing regions in an oil and gas reserve calculation unit, characterized in that: The device comprises: A coordinate conversion unit is used to convert the geodetic coordinates of the oil and gas-bearing area into longitude and latitude coordinates when the oil and gas-bearing area crosses the belt; A center determination unit, used to determine the center of the oil and gas containing area according to the minimum longitude and the maximum longitude of the oil and gas containing area; A regional coordinate construction unit, used to construct a regional coordinate system for the oil and gas containing region using the center, and convert the latitude and longitude coordinates into regional coordinates of the oil and gas containing region; A segmentation unit, used for segmenting the oil and gas-bearing region into a plurality of sub-regions according to the regional coordinates of the oil and gas-bearing region; The area calculation unit is used to cumulatively calculate the area of ​​each sub-area to obtain the area of ​​the oil and gas containing area.

9. A computer device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the computer program, the method according to any one of claims 1 to 7 is implemented.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the method according to any one of claims 1 to 7 is implemented.

Images