Real-time visualization method for seismic crack prediction result data

Abstract

The invention provides a real-time visualization method for seismic fracture prediction result data, and the method comprises the steps: taking the fracture direction and fracture strength of a seismic fracture prediction result as a basis, and carrying out the block construction of fracture visualization results of different scales according to the fracture graph visualization density; according to the method, a strategy of image drawing of multi-scale block images can be rapidly loaded through proportional scale mapping, real-time display of large-scale crack prediction result data is realized, and common browsing operations such as real-time scaling and translation of a large-scale earthquake prediction crack image are realized. According to the method, the requirements of crack screening and visualization modes are comprehensively considered, and technical support is provided for real-time interactive display of large-scale crack prediction results.

Description

Technical Field

[0001] The invention relates to the field of geophysical exploration seismic interpretation and comprehensive research, and in particular to a real-time visualization method for earthquake crack prediction result data. Background Art

[0002] In the exploration of various types of reservoirs such as tight sandstone, oil shale, igneous rock and buried hill in oil fields, fractures are the main migration and storage space of oil and gas, and fracture prediction is an important link in complex oil and gas exploration. The traditional seismic fracture prediction results have a small amount of data and have little impact on the display and real-time interaction of fractures. In the case of a small scale, the conventional real-time drawing method can meet the real-time interaction requirements of fracture prediction results.

[0003] With the continuous expansion of seismic exploration scale and the demand for comprehensive display of fracture prediction results across blocks, the scale of seismic fracture display continues to expand. The conventional real-time drawing method has a large workload and low efficiency, and it is difficult to meet the needs of real-time interaction. There is a lack of real-time interactive display methods for large-scale fracture prediction results. Summary of the invention

[0004] In view of the above problems, the present invention is proposed to provide a real-time visualization method for earthquake fracture prediction result data that overcomes the above problems or at least partially solves the above problems.

[0005] According to one aspect of the present invention, a real-time visualization method for earthquake crack prediction result data is provided, the visualization method comprising:

[0006] Step 100: Read and analyze the crack prediction result data to obtain the seismic information of each seismic channel in the data;

[0007] Step 200: Given a target gather and a fracture visualization target, pre-draw fracture visualization results of different scales in blocks according to fracture graphic visualization density and projection ratio;

[0008] Step 300: Given a logical center coordinate and a projection scale coefficient, the logical coordinate range is calculated using the viewport size, and the crack pre-drawn blocks within the range are loaded according to the scale mapping and the logical coordinate range, and the selected pre-drawn results are scaled and displayed in the viewport;

[0009] Step 400: Change the logical center coordinates and repeat step 300;

[0010] Step 500: Change the mapping scale factor and repeat step 300;

[0011] Step 600: Change the drawing target gather and repeat step 200;

[0012] Step 700: Change the visualization target and repeat step 200.

[0013] Optionally, the seismic information includes: line number, track number, crack direction and crack strength.

[0014] Optionally, the crack prediction result data is data with a line scale exceeding 10°.

[0015] Optionally, in step 200, a target gather is given, including part of the gather data in the entire gather data; and a fracture visualization target is given, including a display method of the fracture.

[0016] Optionally, the display methods of the cracks specifically include: directional ellipse, rose diagram, and line diagram.

[0017] Optionally, in step 200, pre-drawing the fracture visualization results of different scales in blocks according to the fracture graphic visualization density and projection ratio specifically includes:

[0018] Construct a pre-drawn multi-scale image model and a quadtree pyramid model. Given the size of an initial image, construct the corresponding image and draw the target gather image on the image.

[0019] Construct four images of the same size, and draw the target gather in blocks on the four images respectively, with one quarter of the target gather area drawn on each image;

[0020] Construct eight images of the same size, and draw the target gather in blocks on each of the eight images. Draw one-fourth of the target gather area in each image, and build them in sequence until the crack drawing scale of all images reaches a certain threshold, completing the block pre-drawing of the image.

[0021] Optionally, in step 200, the generation of the block pre-drawn image is sequentially generated and cached in a background parallel manner, or is generated and cached in a ready-to-use manner.

[0022] Optionally, in step 200, the selection of crack lines drawn on different blocks is not only related to whether the center point of the line is in the block area, but also to whether the drawing area of ​​the crack symbol intersects with the block image.

[0023] Optionally, in step 300, a logical center coordinate is given, the logical center of the display gather is set, and the logical center coordinate is mapped to the center of the fracture display area;

[0024] Given a projection scale factor, set the display scale and calculate the logical coordinate range through the projection affine transformation.

[0025] Optionally, in step 300, loading the pre-drawn crack blocks within the range according to the scale mapping relationship and the logical coordinate range includes:

[0026] Establishing the mapping relationship between the projection scale coefficient and different scale levels of the pre-drawn multi-scale image model;

[0027] The scale level is calculated according to the projection scale coefficient, and the pre-drawn image at the corresponding scale level is loaded through the logical coordinate range.

[0028] Optionally, in step 300, scaling and displaying the selected pre-drawing result in the viewport includes:

[0029] The selected blocks are scaled and drawn at the corresponding position of the viewport according to the precise position ratio;

[0030] If there is no corresponding scale data in the multi-scale fracture model, draw it directly.

[0031] The present invention provides a real-time visualization method for earthquake crack prediction result data. Based on the crack direction and crack strength of earthquake crack prediction results, the method adopts a block construction of crack visualization results of different scales according to the density of crack graphic visualization, and can quickly load multi-scale block images through proportional scale mapping. The image drawing strategy realizes the real-time display of large-scale crack prediction result data and realizes common browsing operations such as real-time zooming and translation of large-scale earthquake prediction crack maps. The method comprehensively considers the needs of crack screening and visualization methods, and provides technical support for the real-time interactive display of large-scale crack prediction results.

[0032] The above description is only an overview of the technical solution of the present invention. In order to more clearly understand the technical means of the present invention, it can be implemented according to the contents of the specification. In order to make the above and other purposes, features and advantages of the present invention more obvious and easy to understand, the specific implementation methods of the present invention are listed below. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0034] Figure 1 A real-time visualization method for earthquake crack prediction result data provided by an embodiment of the present invention;

[0035] Figure 2 A quadtree-based multi-scale crack pre-rendering model provided in an embodiment of the present invention;

[0036] Figure 3 Demonstration of different visualization target display modes provided by embodiments of the present invention;

[0037] Figure 4 This is a diagram of the screening results of cracks of different strengths and directions provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0038] The exemplary embodiments of the present disclosure will be described in more detail below with reference to the accompanying drawings. Although the exemplary embodiments of the present disclosure are shown in the accompanying drawings, it should be understood that the present disclosure can be implemented in various forms and should not be limited by the embodiments set forth herein. On the contrary, these embodiments are provided to enable a more thorough understanding of the present disclosure and to fully convey the scope of the present disclosure to those skilled in the art.

[0039] The terms "comprises" and "having" and any variations thereof in the description embodiments, claims and drawings of the present invention are intended to cover non-exclusive inclusions, for example, including a series of steps or units.

[0040] The technical solution of the present invention is further described in detail below in conjunction with the accompanying drawings and embodiments.

[0041] A real-time visualization method for large-scale earthquake crack prediction result data is provided. The method constructs crack visualization results of different scales in blocks according to the density of crack graphic visualization. Multi-scale block images can be quickly loaded through proportional scale mapping to achieve real-time display of large-scale crack prediction result data. It also supports common crack selection result browsing operations such as real-time zooming and panning.

[0042] The present invention adopts the following technical solution:

[0043] Step 1: Read and analyze the crack prediction result data to obtain the line number, channel number, crack direction and crack strength of each seismic channel in the data;

[0044] It can read the crack prediction text format data in the following format, parse the data, and obtain the location, crack direction, and crack strength information of each crack;

[0045] Crack prediction result text format data

[0046]

[0047] The reading and parsing results are in the same format.

[0048] The crack direction ranges from 0 to 360, where 0 degrees is horizontal to the right, indicating the east direction, and the crack intensity ranges from 0 to 1, where 0 indicates the weakest crack intensity and 1 indicates the strongest crack intensity.

[0049] Step 2: Given a target gather and a fracture visualization target, pre-draw fracture visualization results of different scales in blocks according to the fracture graphic visualization density and projection ratio;

[0050] Given a target gather, the desired target gather can be obtained through given conditions, such as the gather corresponding to the direction angle ∈[θ0,θ1] and the crack strength ∈[s0,s1].

[0051] Given a crack visualization target, determine the visualization method of the crack, such as filled ellipse, rose diagram, color separation line diagram, etc. Figure 3 As shown, different visualization target display methods are shown.

[0052] like Figure 2 As shown, a multi-scale crack image model and a pyramid model are constructed based on the quadtree method.

[0053] (1) Given an initial image of size (w0, h0), construct a corresponding image with a width of w0 and a height of h0, obtain the image drawing object, use the drawing object to draw the crack images of all traces of the target gather according to the set visualization target, and cache them with code 0.

[0054] (2) Construct four images of image size (w0, h0), which are encoded as {00, 01, 10, 11} respectively, and calculate the logical range of the target gather {(Inline min ,Inline max ),(Cdp min ,Cdp max )}, calculate the Inline,Cdp range according to the following formula: d Inline =Inline max -Inline min , d Cdp =Cdp max -Cdp min . Divide the logical coordinates into four equal parts to obtain four logical ranges Four logical coordinate regions are obtained, corresponding to images coded as {00, 01, 10, 11} respectively. The crack images of all traces intersecting with each logical range are drawn according to the set visualization target, and the crack image drawing of one scale is completed. The drawing results are cached with the code {00, 01, 10, 11}.

[0055] (3) According to the method shown in (2), the crack visualization results of multiple scales such as 8 images and 16 images are constructed in sequence, and the rendering results are cached in a quadtree hierarchical encoding.

[0056] (4) Construct each scale according to the method shown in (2) and (3) until the crack mapping scale of all images reaches a certain density threshold Th stop ,in

[0057] Th stop = Number of images to be drawn / image area, that is, Th≤Th stop , stop building the multi-scale crack image model, and record the maximum scale L max .

[0058] A mapping relationship between the projection scale coefficient and different scale levels of the pre-drawn multi-scale image model is established, that is, a corresponding relationship between the projection scale coefficient S and the scale level L is established.

[0059] S=aL+b

[0060] Changing a and b can achieve different mapping relationships between projection scale coefficient and scale level.

[0061] Step 3: Given a logical center coordinate and a projection scale coefficient, the logical coordinate range is calculated using the viewport size, and the crack pre-drawn blocks within the range are loaded according to the scale mapping and the logical coordinate range, and the selected pre-drawn results are scaled and displayed in the viewport;

[0062] Given a logical center coordinate (Inline0, Cdp0).

[0063] Map the logical center coordinates to the center of the crack display area (viewport); if the display area (viewport) is a rectangular area (x, y, w, h), the center point coordinates are Map the logical coordinate point (Inline0, Cdp0) to (x center ,y center ).

[0064] Given a projection scale coefficient S0, the logical coordinate range {(Inline1, Cdp1), (Inline2, Cdp2)} is calculated by projection affine transformation. The affine transformation is as follows:

[0065]

[0066] The scale level L0 is calculated according to the scale coefficient S0, and the pre-drawn image at the corresponding scale level can be selected through the logical coordinate range {(Inline1, Cdp1), (Inline2, Cdp2)}.

[0067] Calculate the pixel range (w p ,h p), and calculates its scale ratio with the viewport extents Calculate the position of each image in the viewport:

[0068]

[0069] Draw the selected pre-drawn result proportionally at (x,y) in the viewport.

[0070] If there is no corresponding scale data in the multi-scale fracture model, that is, L0>L max Draw directly.

[0071] like Figure 4 As shown in the figure, the screening results of cracks with different intensities and directions.

[0072] Step 4: Change the logical center coordinates and repeat step 3;

[0073] Change the logical center coordinates to (Inline1, Cdp1), repeat step 3, complete the drawing of the crack image with the changed center, and achieve the effect of view translation.

[0074] Step 5: Change the mapping scale factor and repeat step 3;

[0075] Change the radial transformation scale factor of the view and logical coordinates and repeat the above calculation to achieve the scaling effect of the crack map.

[0076] Repeat step 4 for each offset range.

[0077] Step 6: Change the target gather and repeat step 2;

[0078] Change the target gather, for example, select a gather that satisfies the direction angle ∈[θ2,θ3] and the crack strength ∈[s2,s3].

[0079] Step 7: Change the visualization target and repeat step 2.

[0080] Change the visualization target, such as changing the visualization target from a rose diagram to a line diagram, repeat step 2 to build a multi-scale fracture model, and then select the display.

[0081] Implementation Example 1

[0082] In a specific embodiment 1 of the present invention, Figure 1 As shown, Figure 1 This is a flow chart of a method for real-time display of large-scale crack prediction results of the present invention. The real-time display method includes the following steps:

[0083] Step 101: Read and analyze the crack prediction result data to obtain information such as the line number, track number, crack direction and crack strength of each seismic track in the data.

[0084] The crack results after reading are as follows:

[0085]

[0086] Step 102: Given a target gather and a crack visualization target, pre-draw the crack visualization results of different scales in blocks according to the crack graphic visualization density and projection ratio. Select the gather corresponding to the direction angle ∈ [0,360] and the crack intensity ∈ [0,1], the logical range {(1,100), (1,70)}, and visualize the cracks in the form of a color-coded line graph. Construct the corresponding image with a width of 600 and a height of 420 (the aspect ratio is the same as the aspect ratio of the gather range), draw the crack image on the image, and cache it with code 0, and set the density threshold Th stop =50 / (600*420), construct four images of size (600,420), respectively encoded as {00,01,10,11}, divide the logical coordinates into four equal parts, and obtain four logical ranges {(1,50),(1,35)}, {(51,100),(1,35)}, {(1,50),(36,70)}, {(51,100),(36,70)} corresponding to the images encoded as {00,01,10,11}, calculate Th for each image, if Th>Th stop Draw the crack images of all the paths intersecting each logical range according to the set visualization target, complete the crack image drawing of one scale, and cache the drawing results with the code {00,01,10,11}. Repeat the above process, visualize the crack results of all scales and cache the drawing results, until Th≤Th stop , stop building the multi-scale crack image model and record the maximum scale of 6.

[0087] A corresponding relationship between the projection scale coefficient S and the scale level L is established: S=0.05L.

[0088] Step 103: Given a logical center coordinate, given a projection scale factor, the logical coordinate range is calculated using the viewport size, and the crack pre-drawn blocks within the range are loaded according to the scale mapping and the logical coordinate range, and the selected pre-drawn results are scaled and displayed in the viewport. Given a logical center coordinate (50,35), the range of the display area (viewport) is (0,0,1024,768), the center of the viewport is (512,384), and given a projection scale factor of 0.05, the calculated logical coordinate range is {(24.4,15.8), (75.6,54.2)}. According to S=0.05, the scale level S=0.05, L=1 is calculated. The logical coordinate range {(24.4,15.8),(75.6,54.2)} and the L=1 level cache images {(1,50),(1,35)}, {(51,100),(1,35)}, {(1,50),(36,70)}, {(51,100),(36,70)} all have intersections, so the images encoded as {00,01,10,11} are all selected. The logical coordinate range {(24.4,15.8),(75.6,54.2)} corresponds to the pixel range (614,460) at scale level L=1, and the scaling ratio with the viewport range is calculated. Calculate the drawing positions of the {00,01,10,11} image in the viewport: (-468,-296), (532,-296), (-468,404), (532,404) and draw the image at the corresponding positions in proportion.

[0089] Step 104: Change the logical center coordinates to (45, 20), and repeat step 103 to achieve the effect of view translation.

[0090] Step 105: Change the mapping scale factor S=0.01, repeat step 103; achieve the view zoom-in effect, change the mapping scale factor S=0.1, repeat step 103; achieve the view zoom-out effect.

[0091] Step 106: Change the target gather, for example, select a gather that satisfies the direction angle ∈ [60, 275] and the crack strength ∈ [0.6, 0.8], and repeat step 102.

[0092] Step 107: Change the visualization target from a line graph to a rose graph, and repeat step 102.

[0093] Beneficial effects: The present invention is based on the crack direction and crack strength of earthquake crack prediction results, and adopts a block construction of crack visualization results of different scales according to the density of crack graphic visualization. It can quickly load multi-scale block images through proportional scale mapping, realize the real-time display of large-scale crack prediction result data, and realize the common browsing operations such as real-time zooming and translation of large-scale earthquake prediction crack maps. It comprehensively considers the needs of crack screening and visualization methods, and provides technical support for the real-time interactive display of large-scale crack prediction results.

[0094] The above specific implementation methods further illustrate the objectives, technical solutions and beneficial effects of the present invention in detail. It should be understood that the above are only specific implementation methods of the present invention and are not intended to limit the scope of protection of the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention should be included in the scope of protection of the present invention.

Claims

1. A real-time visualization method for earthquake crack prediction result data, characterized in that: The visualization method comprises: Step 100: Read and analyze the crack prediction result data to obtain the seismic information of each seismic channel in the data; Step 200: Given a target gather and a fracture visualization target, pre-draw fracture visualization results of different scales in blocks according to fracture graphic visualization density and projection ratio; Step 300: Given a logical center coordinate and a projection scale coefficient, the logical coordinate range is calculated using the viewport size, and the crack pre-drawn blocks within the range are loaded according to the scale mapping and the logical coordinate range, and the selected pre-drawn results are scaled and displayed in the viewport; Step 400: Change the logical center coordinates and repeat step 300; Step 500: Change the mapping scale factor and repeat step 300; Step 600: Change the drawing target gather and repeat step 200; Step 700: Change the visualization target and repeat step 200.

2. A real-time visualization method for earthquake crack prediction result data according to claim 1, characterized in that: The seismic information includes: line number, track number, crack direction and crack strength.

3. A real-time visualization method for earthquake crack prediction result data according to claim 1, characterized in that: The crack prediction result data is data with a line scale exceeding 10°.

4. A real-time visualization method for earthquake crack prediction result data according to claim 1, characterized in that: In step 200, a target gather is given, including part of the gather data in the whole gather data; and a fracture visualization target is given, including the display method of the fracture.

5. The real-time visualization method for earthquake crack prediction result data according to claim 1, characterized in that: The display methods of the cracks specifically include: directional ellipse, rose diagram, and line diagram.

6. A real-time visualization method for earthquake crack prediction result data according to claim 1, characterized in that: In step 200, pre-drawing the fracture visualization results in blocks of different scales according to the fracture graphic visualization density and projection ratio specifically includes: Construct a pre-drawn multi-scale image model and a quadtree pyramid model. Given the size of an initial image, construct the corresponding image and draw the target gather image on the image. Construct four images of the same size, and draw the target gather in blocks on the four images respectively, with one quarter of the target gather area drawn on each image; Construct eight images of the same size, and draw the target gather in blocks on each of the eight images. Draw one-fourth of the target gather area on each image, and build them in sequence until the crack drawing scale of all images reaches a certain threshold, completing the block pre-drawing of the image.

7. The real-time visualization method for earthquake crack prediction result data according to claim 1, characterized in that: In step 200, the generation of the block pre-drawn image is sequentially generated and cached in a background parallel manner, or is generated and cached in a ready-to-use manner.

8. The real-time visualization method for earthquake crack prediction result data according to claim 1, characterized in that: In step 200, the selection of crack lines drawn on different blocks is not only related to whether the center point of the line is in the block area, but also related to whether the drawing area of ​​the crack symbol intersects with the block image.

9. The real-time visualization method for earthquake crack prediction result data according to claim 1, characterized in that: In the step 300, a logical center coordinate is given, the logical center of the display gather is set, and the logical center coordinate is mapped to the center of the fracture display area; Given a projection scale factor, set the display scale and calculate the logical coordinate range through the projection affine transformation.

10. The real-time visualization method for earthquake crack prediction result data according to claim 1, characterized in that: In step 300, loading the pre-drawn crack blocks within the range according to the scale mapping relationship and the logical coordinate range includes: Establishing the mapping relationship between the projection scale coefficient and different scale levels of the pre-drawn multi-scale image model; The scale level is calculated according to the projection scale coefficient, and the pre-drawn image at the corresponding scale level is loaded through the logical coordinate range.

11. A real-time visualization method for earthquake crack prediction result data according to claim 1, characterized in that: In step 300, scaling and displaying the selected pre-drawing result in the viewport includes: The selected blocks are scaled and drawn at the corresponding position of the viewport according to the precise position ratio; If there is no corresponding scale data in the multi-scale fracture model, draw it directly.

Images