Method and device for rapidly identifying hydrate of sand pore filling type in sea area

Abstract

The application discloses a kind of sea area sandy pore filling type hydrate on-site rapid identification method and device, the method includes inputting the gamma value and deep resistivity value of the well site logging data of on-site drilling, according to the low natural gamma cutoff value and resistivity baseline value determined, several reservoir intervals lower than the low natural gamma cutoff value and deep resistivity greater than the resistivity baseline value are divided out;Make chart one, the chart one is deep resistivity and longitudinal wave speed intersection chart;According to the several reservoir intervals divided out, the deep resistivity value RD and longitudinal wave speed value Vp corresponding to each reservoir interval are obtained, and RD value and Vp value are projected into the chart one, and the projection position is obtained;According to the area of projection position in chart one, the hydrate type of reservoir interval is judged.The application only needs to obtain deep resistivity RD, longitudinal wave speed Vp from the logging data of natural gas hydrate drilling well site, and can preliminarily judge the hydrate type.

Description

Technical Field

[0001] This invention relates to the field of natural gas hydrate technology, specifically to a method and apparatus for rapid on-site identification of marine sandy pore-filled hydrates. Background Technology

[0002] Natural gas hydrates are abundant, highly efficient, and clean energy resources, representing a strategic high ground for future global energy development. Studies have shown that sandy hydrate reservoirs possess characteristics such as high porosity, high saturation, and coarse grain size, making them crucial targets for efficient trial production under current technological conditions. Japan has selected sandy hydrate reservoirs as its targets in both of its hydrate trial productions. In conducting offshore hydrate drilling, how can sandy hydrate reservoirs be quickly and accurately identified on-site to provide a basis for formulating subsequent coring and development plans?

[0003] The conventional method first identifies sandy hydrate reservoirs by logging characteristics such as low natural gamma ray, high resistivity, and high P-wave velocity, then develops a coring plan, conducts coring trips, and obtains samples to confirm the presence of sandy hydrate reservoirs. The disadvantages of this method are: 1. Reduced drilling efficiency. Generally, hydrate drilling involves two separate trips: logging while drilling and coring. After completing the pilot hole logging operation, a coring trip is conducted based on logging characteristics to determine if the sample is from a sandy hydrate reservoir. This prevents immediate confirmation at the drilling site, hindering timely adjustments and improvements to the drilling plan based on the results. 2. The ambiguity of logging results means that the drilled sample cannot guarantee that it is a sandy, pore-filled hydrate. According to actual drilling results in the South China Sea, wells with logging characteristics such as "low natural gamma, high resistivity, and high P-wave velocity" were found to contain silty mudstone or fractured hydrates with increased carbonate mineral components after coring, or even carbonate rocks without hydrates. This situation resulted in a significant waste of time and increased operating costs. Summary of the Invention

[0004] To address at least one technical problem existing in the prior art, the present invention provides a method and apparatus for rapid on-site identification of marine sandy pore-filled hydrates, thereby reducing time and operational costs.

[0005] To achieve the above objectives, the technical solution of the present invention is as follows:

[0006] In a first aspect, the present invention provides a method for rapid on-site identification of marine sandy pore-filled hydrates, comprising:

[0007] Input the gamma value and deep resistivity value of the logging data of the drilling well location in the field, and divide several reservoir segments that are lower than the low natural gamma cutoff value and have a deep resistivity greater than the resistivity baseline value according to the determined low natural gamma cutoff value and resistivity baseline value.

[0008] Create Figure 1, which is a diagram showing the intersection of deep resistivity and longitudinal wave velocity.

[0009] Based on the divided reservoir segments, the deep resistivity value RD and the longitudinal wave velocity value Vp corresponding to each reservoir segment are obtained, and the RD value and Vp value are projected onto the first map to obtain the projection position.

[0010] The type of hydrate in the reservoir segment can be determined based on the area where the projection is located in Plate 1.

[0011] Furthermore, the first diagram includes three regions: a, b, and c. If more than 90% of the projected points are in region a, the reservoir segment is determined to be a sandy pore-filled hydrate. If more than 90% of the projected points are in region b, the reservoir segment is determined to be a fractured hydrate. If more than 90% of the projected points are in region c, the reservoir segment is determined to be a calcareous clayey silt.

[0012] Furthermore, the horizontal axis of the first chart represents deep resistivity, with a scale range of 0.1-1000, and uses a logarithmic scale; the vertical axis represents longitudinal wave velocity, with a scale range of 0-4000.

[0013] Furthermore, the regions a, b, and c are divided by two dividing lines, namely dividing line one and dividing line two.

[0014] The area enclosed by the first boundary line is region a, and the first boundary line is the line connecting the coordinates of the corner points A, B, C, D, G, and H.

[0015] The region enclosed by the second boundary line and the first boundary line is region b, and the second boundary line is the line connecting corner points D, E, and F.

[0016] Furthermore, the coordinates of the corner points A, B, C, D, E, F, G, and H are as follows: A(8.5, 4000), B(8.5, 2475), C(5.2, 2475), D(5.2, 2080), E(5.2, 1620), F(1000, 1620), G(400, 2080), H(1000, 3393).

[0017] Furthermore, the method for rapid on-site identification of marine sandy pore-filling hydrates also includes:

[0018] The density and neutron values ​​of the two hydrate types of stations in regions a and b are obtained and projected onto plate two to obtain the projection positions; plate two is a density and neutron intersection diagram.

[0019] The type of hydrate can be determined by the area in Plate 2 where the projection is located, to determine whether the hydrate is a sandy pore-filled hydrate.

[0020] Furthermore, the horizontal axis of the second plate represents density, with a scale range of 1-2.6, and the vertical axis represents neutron porosity, with a scale range of 0.2-1; the coordinates of the decomposition point I are (1.6, 0.75), and two straight lines parallel to the x-axis and y-axis are drawn with the coordinates of point I. The region enclosed by the two straight lines parallel to the x-axis and y-axis is region d.

[0021] If the projected location is within region d, then the reservoir segment is determined to be a sandy pore-filled hydrate.

[0022] Furthermore, the low natural gamma cutoff value is 65 API, and the resistivity baseline value is 1-1.5 ohm.

[0023] Secondly, the present invention provides a rapid on-site identification device for marine sandy pore-filled hydrates, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the steps of any of the methods described above.

[0024] In a second aspect, the present invention provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the steps of any of the methods described above.

[0025] Compared with the prior art, the advantages of this invention are as follows:

[0026] 1) The type of hydrate can be preliminarily determined by obtaining the deep resistivity RD and longitudinal wave velocity Vp from the logging data of the natural gas hydrate drilling site; by adding the density RHON and neutron porosity TNPH values, the sandy pore-filled hydrate can be quickly and accurately identified without other data.

[0027] 2) The rapid identification chart established under the plane rectangular coordinate system does not require professional software for data processing. It can identify sandy pore-filling hydrates on site in the first time, saving time and with extremely high efficiency.

[0028] 3) The sandy pore-filled hydrate identification chart is simple to use, requires no professional knowledge, and can be widely used by anyone, making it applicable to a wider range of people. Attached Figure Description

[0029] Figure 1 This is a flowchart of the rapid on-site identification method for marine sandy pore-filling hydrates provided in Embodiment 1 of the present invention;

[0030] Figure 2 The figure shows the intersection of deep resistance and longitudinal wave velocity.

[0031] Figure 3 Template for the intersection of resistivity and sound velocity in the field;

[0032] Figure 4 This is a diagram showing the intersection of density and neutron porosity on a plate.

[0033] Figure 5 Template for on-site density-neutron porosity cross-section diagram;

[0034] Figure 6 This is a schematic diagram of the low-gamma interval in the four wells in the experimental example;

[0035] Figure 7 The diagram shown in the experimental example is a graph depicting the intersection of deep resistance and longitudinal wave velocity.

[0036] Figure 8 The diagram shows the intersection of density and neutron porosity in the experimental example.

[0037] Figure 9 The actual hydrates corresponding to the low gamma layers in the four wells in the experimental example;

[0038] Figure 10 This is a schematic diagram of the composition of the rapid on-site identification device for marine sandy pore-filled hydrates provided in Embodiment 2 of the present invention. Detailed Implementation

[0039] The technical solution of the present invention will be further described below with reference to the accompanying drawings and embodiments.

[0040] Implementation 1:

[0041] See Figure 1 As shown in this embodiment, the method for rapid on-site identification of marine sandy pore-filling hydrates mainly includes the following steps:

[0042] 101. Input the gamma value and deep resistivity value of the well logging data of the field drilling location, and divide several reservoir segments below the low natural gamma cutoff value and with deep resistivity greater than the resistivity baseline value according to the determined low natural gamma cutoff value and resistivity baseline value.

[0043] The primary condition for sandy hydrate reservoirs is that the logging curves exhibit low gamma and high resistivity characteristics. Therefore, in this step, the cutoff values ​​for these two parameters are determined first.

[0044] In one specific embodiment, the values ​​of the two parameters, the low natural gamma cutoff and the resistivity baseline, are determined as follows:

[0045] (1) Determination of the low natural gamma cutoff value: Based on hydrate drilling in the southeastern Qiongdong sea area, sandy pore-filled hydrate samples were obtained from wells QDN-W01-2019 and QDN-W03-2019, with corresponding average natural gamma values ​​of 65 API and 62 API, respectively. A higher natural gamma value generally indicates a higher clay content and a lower sand content. Therefore, the upper limit threshold of natural gamma for sandy hydrate samples was set at 65 API.

[0046] (2) Determination of resistivity baseline: The resistivity value of hydrate-bearing sediments is greater than that of saturated water layers. Therefore, based on the hydrate drilling results, the deep resistivity baseline value of aquifers is generally around 1-1.5 ohms.

[0047] After determining the cutoff values ​​of the above two parameters, the gamma value and deep resistivity value of the logging data of the drilling site are obtained, and the reservoir segments unit1, unit2...unitN with natural gamma below 65 API and deep resistivity greater than the resistivity baseline are delineated.

[0048] 102. Create Figure 1, which is a diagram showing the intersection of deep resistivity and longitudinal wave velocity.

[0049] In one specific embodiment, the manufacturing process of Figure 1 is as follows:

[0050] like Figure 2 As shown in the figure, Plate 1 selects typical wells with low gamma-ray intervals in three sea areas: Xisha, Shenhu, and Qiongdongnan. Core samples were obtained from all these wells, and the data in the figure are reliable.

[0051] In the figure, △ represents the lower natural gamma ray due to the increased carbonate mineral composition (mainly calcite and dolomite). Gray △ represents the low gamma ray interval of well XH-W06-2016 in the Xisha Islands (average natural gamma ray: 49 API), and purple △ represents the low gamma ray interval of well SH-W19-2015 in the Shenhu Sea Area (average natural gamma ray: 53 API). XRD analysis of the cores from both wells indicates that the low gamma ray is due to the increased calcium content, not the increased sand content; the cores are all calcareous clayey silt. Green △ represents the low gamma ray interval of well QDN-W03-2019 in the Qiongdongnan Sea Area (average natural gamma ray: 51 API), and the core is confirmed to be carbonate rock.

[0052] In the figure, ○ represents the low natural gamma mainly caused by the increase in sandy mineral components. Three wells, QDN-W05-2021, QDN-W10-2021 and QDN-W01-2019, were confirmed by core samples to be sandy pore-filled hydrates.

[0053] In the figure, ж represents the lower natural gamma caused by fracture-filled hydrates (which replace part of the framework). Wells QDN-W03-2021 and QDN-W01-2019 are fracture-type hydrate zones (among which well QDN-W01-2019 contains massive hydrates, a special type of fracture-type hydrate), and core samples confirm that the lithology is mainly silty mudstone.

[0054] From Figure 1, it can be seen that the low gamma points caused by calcareous clayey silt and carbonate rocks (symbol △) fall in region c. More than 90% of the points (symbol ○) of the target reservoir sandy pore-filling hydrates are in region a inside boundary line 1. More than 90% of the points of fractured hydrates (symbol ж) fall in region b, which is enclosed by boundary line 1 and boundary line 2. Some fractured hydrate points fall in region a.

[0055] Regions a, b, and c are divided by two boundary lines, and the coordinates of their corner points A, B, C, D, E, F, G, and H are as follows:

[0056] A(8.5,4000), B(8.5,2475), C(5.2,2475), D(5.2,2080), E(5.2,1620), F(1000,1620), G(400,2080), H(1000,3393).

[0057] Quickly create a diagram at the drilling site, namely the intersection of deep resistivity and P-wave velocity. The horizontal axis represents deep resistivity, with a scale range of 0.1-1000, using a logarithmic scale; the vertical axis represents P-wave velocity, with a scale range of 0-4000. Based on the coordinates of the corner points A, B, C, D, E, F, G, and H, draw two red and blue dividing lines (colors can be freely chosen) to distinguish the three different regions a, b, and c. Figure 3 As shown.

[0058] 103. Based on the identified reservoir segments unit1, unit2...unitN, obtain the deep resistivity value RD and P-wave velocity value Vp corresponding to each reservoir segment, and project the RD value and Vp value onto the first drawing to obtain the projection position;

[0059] 104. Determine the hydrate type of the reservoir segment based on the area where the projection is located in Plate 1.

[0060] In one specific embodiment, if more than 90% of the projected locations are in region a, the reservoir segment can be determined to be a sandy pore-filled hydrate; if most of the projected locations are in region b, the reservoir segment can be determined to be a fractured hydrate; if most of the projected locations are in region c, the reservoir segment can be determined to be a calcareous clayey silt.

[0061] In addition, to further and more accurately identify sandy pore-filled hydrates, the rapid on-site identification method for marine sandy pore-filled hydrates provided in this embodiment also includes the following steps:

[0062] The density and neutron values ​​of the two hydrate types of sites in regions a and b are obtained and projected onto Plate 2.

[0063] According to Plate 1, some fractured hydrate sites were found to fall within region a. Therefore, to further and more accurately identify sandy pore-filled hydrates, the density and neutron values ​​of the two hydrate types can be projected onto the density-neutron-porosity intersection plate 2.

[0064] Specifically, the production process for Plate 2 is as follows:

[0065] like Figure 4 As shown, the data points in Figure 2 and Figure 1 are from the same well and the same low-gamma layer, and therefore both have core evidence, making the data highly reliable.

[0066] like Figure 4 As shown in Figure 2, it can be seen that there is a clear boundary between fractured hydrates (symbol ж) and sandy pore-filled hydrates (symbol ○). The data points of sandy pore-filled hydrates only fall in region d and are completely separated from the data points of fractured hydrates. The coordinates of the decomposition point I are (1.6, 0.75).

[0067] Therefore, at the drilling site, if the projection point of the RD value and Vp value, determined according to Chart 1, is located in region c, this step is unnecessary; if it is in region a or region b, this step is performed to quickly create a density-neutron porosity intersection diagram, i.e., Chart 2. The horizontal axis of Chart 2 represents density, with a scale range of 1-2.6, and the vertical axis represents neutron porosity, with a scale range of 0.2-1. Using the coordinates of the aforementioned decomposition point I, draw two straight lines parallel to the x-axis and y-axis, as shown below. Figure 5 As shown, the region enclosed by the two straight lines parallel to the x-axis and y-axis is region d.

[0068] If the density neutron value of the projection is located within region d, it can be further identified as a sandy pore-filled hydrate.

[0069] The following experimental example will further verify and illustrate this method:

[0070] Experimental example:

[0071] Four representative wells with low gamma characteristics from different sea areas in the Xisha, Shenhu, and Qiongdongnan regions of the South China Sea were selected, and core samples were obtained from all of them. Figure 6As shown in (a)-6(d), they are XH-W03-2016, SH-W18-2015, QDN-W03-2019, and QDN-W08-2018, respectively.

[0072] (1) Identify the depth intervals of these four wells that simultaneously meet the requirements of low gamma (GR≤65API) and above the resistivity baseline, and represent them with dashed rectangles;

[0073] (2) Based on the chart production method determined in step 102 above, the resistivity and P-wave velocity corresponding to the low-gamma intervals of the above four wells are projected onto chart one, and the results are as follows. Figure 7 As shown.

[0074] according to Figure 7 It can be determined that wells XH-W03-2016 and SH-W18-2015 are calcareous clayey silt (symbol △), QDN-W08-2018 is a fractured hydrate (symbol ж), and only QDN-W03-2019 is a sandstone pore-filling hydrate (symbol ○).

[0075] (3) To further verify the porous hydrate in the QDN-W03-2019 sandstone, a second plot was created according to step 105 above. The density and neutron porosity values ​​corresponding to the fractured hydrate in the QDN-W08-2018 well and the low-gamma interval in the QDN-W03-2019 well were projected onto the second plot. The projection results are as follows: Figure 8 As shown.

[0076] According to such Figure 8 The projection results shown further confirm that the low-gamma interval of QDN-W03-2019 is a sandstone pore-filling hydrate.

[0077] (4) Comparison of image recognition results with core samples

[0078] pass Figure 7 The resistivity vs. acoustic velocity plot shown in Figure 1 indicates that wells XH-W03-2016 and SH-W18-2015 have low gamma rays and consist of calcareous clayey silt, while well QDN-W08-2018 has fractured hydrates. Figure 9 (a), (b), and (d) are consistent with the actual core analysis results; combined with Figure 8 As shown in Figure 2 (neutron porosity vs. density cross-plot), well QDN-W03-2019 is identified as a sandy pore-filled hydrate, consistent with... Figure 9 (c) The core analysis results are consistent.

[0079] In summary, compared with the prior art, the present invention has the following technical advantages:

[0080] 1) The type of hydrate can be preliminarily determined by obtaining the deep resistivity RD and longitudinal wave velocity Vp from the logging data of the natural gas hydrate drilling site; by adding the density RHON and neutron porosity TNPH values, the sandy pore-filled hydrate can be quickly and accurately identified without other data.

[0081] 2) The rapid identification chart established under the plane rectangular coordinate system does not require professional software for data processing. It can identify sandy pore-filling hydrates on site in the first time, saving time and with extremely high efficiency.

[0082] 3) The sandy pore-filled hydrate identification chart is simple to use, requires no professional knowledge, and can be widely used by anyone, making it applicable to a wider range of people.

[0083] Example 2:

[0084] See Figure 10 As shown, the rapid on-site identification device for marine sandy porous hydrates provided in this embodiment includes a processor, a memory, and a computer program stored in the memory and executable on the processor, such as a rapid on-site identification program for marine sandy porous hydrates. When the processor executes the computer program, it implements the steps of Embodiment 1 described above, for example... Figure 1 The steps are shown.

[0085] For example, the computer program can be divided into one or more modules / units, which are stored in the memory and executed by the processor to complete the present invention. The one or more modules / units can be a series of computer program instruction segments capable of performing specific functions, which describe the execution process of the computer program in the rapid on-site identification device for marine sandy porous hydrates.

[0086] The rapid on-site identification device for marine sandy porous hydrates can be a desktop computer, laptop, handheld computer, or cloud server, etc. The rapid on-site identification device for marine sandy porous hydrates may include, but is not limited to, a processor and memory. Those skilled in the art will understand that... Figure 10 This is merely an example of a rapid on-site identification device for marine sandy porous hydrates and does not constitute a limitation on the device. It may include more or fewer components than shown in the figure, or combine certain components, or different components. For example, the rapid on-site identification device for marine sandy porous hydrates may also include input / output devices, network access devices, buses, etc.

[0087] The processor referred to can be a Central Processing Unit (CPU), or other general-purpose processors, digital signal processors (DSPs), application-specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. A general-purpose processor can be a microprocessor or any conventional processor.

[0088] The memory can be an internal storage unit of the rapid on-site identification device for marine sandy porous hydrates, such as the hard drive or memory of the device. The memory can also be an external storage device, such as a plug-in hard drive, SmartMediaCard (SMC), Secure Digital (SD) card, or Flash Card equipped on the device. Furthermore, the memory can include both internal storage units and external storage devices. The memory is used to store the computer program and other programs and data required by the rapid on-site identification device for marine sandy porous hydrates. The memory can also be used to temporarily store data that has been output or will be output.

[0089] Example 3:

[0090] This embodiment provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the steps of the method described in Embodiment 1.

[0091] The computer-readable medium shown can be any means that can contain, store, communicate, propagate, or transmit a program for use in or in conjunction with an instruction execution system, apparatus, or device. More specific examples (a non-exhaustive list) of computer-readable media include: an electrical connection having one or more wires (electronic device), a portable computer disk drive (magnetic device), random access memory (RAM), read-only memory (ROM), erasable and editable read-only memory (EPROM or flash memory), fiber optic devices, and portable optical disc read-only memory (CDROM). Additionally, the computer-readable medium can even be paper or other suitable media on which the program can be printed, for example, by optically scanning the paper or other medium, then editing, interpreting, or otherwise processing it as necessary to obtain the program electronically, and then storing it in computer memory.

[0092] For those skilled in the art, various other corresponding changes and modifications can be made based on the technical solutions and concepts described above, and all such changes and modifications should fall within the protection scope of the claims of this invention.

[0093] The above embodiments are merely illustrative of the technical concept and features of the present invention, and are intended to enable those skilled in the art to understand the content of the present invention and implement it accordingly. They should not be construed as limiting the scope of protection of the present invention. All equivalent changes or modifications made based on the essence of the content of the present invention should be covered within the scope of protection of the present invention.

Claims

1. A method for rapid identification of sand pore-filling hydrate in a marine area, characterized in that, The method comprises the following steps: inputting gamma value and deep resistivity value of well logging data of field drilling well site, dividing several reservoir layers which are lower than the low natural gamma cutoff value and deep resistivity value greater than the resistivity baseline value according to the determined low natural gamma cutoff value and resistivity baseline value; making a chart one, which is a deep resistivity and longitudinal wave velocity intersection chart; according to the divided several reservoir layers, obtaining deep resistivity value RD and longitudinal wave velocity value Vp corresponding to each reservoir layer, and projecting RD value and Vp value into the chart one to obtain a projection position; judging hydrate type of the reservoir layer according to the projection position in the area of the chart one; the chart one comprises a, b, c three areas; if more than 90% of the projected points are in area a, it is judged that the reservoir layer is sand pore filling type hydrate; if more than 90% of the projected points are in area b, it is judged that the reservoir layer is fissure type hydrate; if more than 90% of the projected points are in area c, it is judged that the reservoir layer is calcareous clayey silt; the method further comprises the following steps: obtaining density value and neutron value of the two hydrate type stations in area a and area b to project into chart two to obtain a projection position; the chart two is a density and neutron intersection chart; judging whether the hydrate type is sand pore filling type hydrate according to the projection position in the area of the chart two; the horizontal coordinate of the chart two is density, the scale range is 1-2.6, and the vertical coordinate is neutron porosity, the scale range is 0.2-1; the decomposition point I coordinate is (1.6, 0.75), two straight lines parallel to the x axis and the y axis are drawn according to the I point coordinate, and the area surrounded by the two straight lines is area d; if the projection position is in area d, it is determined that the reservoir layer is sand pore filling type hydrate.

2. The method for rapidly identifying a hydrate in a sand pore filling type of marine area according to claim 1, wherein the horizontal coordinate of the chart one is deep resistivity, the scale range is 0.1-1000, and logarithmic coordinate scale is adopted; the vertical coordinate is longitudinal wave velocity, the scale range is 0-4000.

3. The method for rapidly identifying a hydrate in a sand pore filling type of marine sand according to claim 2, wherein the areas a, b and c are divided by two boundary lines, which are boundary line one and boundary line two respectively; the area surrounded by the boundary line one is area a, and the boundary line one is the line connecting the corner point coordinates A, B, C, D, G and H; the area surrounded by the boundary line two and the boundary line one is area b, and the boundary line two is the line connecting the corner point coordinates D, E and F.

4. The method for rapidly identifying a hydrate of the sand pore filling type in a sea area according to claim 3, characterized by, the corner point coordinates A, B, C, D, E, F, G and H are as follows: A (8.5, 4000), B (8.5, 2475), C (5.2, 2475), D (5.2, 2080), E (5.2, 1620), F (1000, 1620), G (400, 2080), and H (1000, 3393).

5. The method for rapidly identifying a hydrate of sand pore filling type in a sea area according to claim 1, wherein the low natural gamma cutoff value is 65 API, and the resistivity baseline value is 1-1.5 ohm.

6. A device for rapid identification of sand pore-filling hydrate in a sea area, comprising a memory, a processor, and a computer program stored in the memory and capable of running on the processor, characterized in that, the processor executes the computer program to realize the steps of the method according to any one of claims 1 to 5.

7. A computer-readable storage medium storing a computer program, wherein the computer program comprises the following steps of: the computer program is executed by the processor to realize the steps of the method according to any one of claims 1 to 5.

Images