Method and system for evaluating movable water saturation of a tight sand gas reservoir

Abstract

The application discloses a method and system for evaluating movable water saturation of a tight sandstone gas reservoir. The method comprises the following steps: using resistivity, acoustic time difference, permeability and water saturation, and combining a movable water saturation model to complete the evaluation of the movable water saturation of the tight sandstone gas reservoir. The method does not need mercury injection experiment data, and can be used by using conventional logging data, and has good universality. The method realizes movable water saturation interpretation without increasing logging series, is simple to operate, and has wide application range.

Description

TECHNICAL FIELD

[0001] The application belongs to the technical field of natural gas development, and relates to a movable water saturation evaluation method and system for a tight sandstone gas reservoir. BACKGROUND

[0002] For a tight sandstone gas reservoir, movable water saturation is a critical parameter for gas phase and gas-water two-phase percolation, and is also an important parameter for evaluating an oil and gas reservoir, predicting productivity and formulating a development plan, so it is of great significance to study the movable water saturation of a tight sandstone gas reservoir.

[0003] Chinese Patent 201910221661.X discloses a logging saturation calculation method for a sandstone oil and gas layer, and describes that the main reason for forming a low resistance of an oil (gas) layer is the equivalent additional conductivity caused by the conductive network formed by a small amount of movable water or bound water distributed in effective pores or movable pores and the water in the inner wall of the hydrophilic rock pore throat, and the conductivity of the bound water itself. After revision, the calculation model or method is suitable for saturation calculation of conventional resistivity oil (gas) layers and low-resistance oil (gas) layers, significantly improves the calculation result of the hydrocarbon saturation of a low-resistance oil (gas) layer. It is suitable for the case where the saturation calculation result of the oil (gas) layer is low, and improves the oil (gas) saturation calculation result, becoming the standard of the resistivity calculation saturation model of the sandstone oil (gas) layer.

[0004] Chinese Patent 201710890856.4 discloses an array induction and lateral logging method for quantitatively evaluating movable water saturation of a sandstone reservoir. The Arpich formula is used as a water saturation calculation model, and the water saturation of the reservoir is calculated by using the lateral deep resistivity and the induction deep resistivity respectively. The water saturation calculated by the induction deep resistivity is corrected for the bound water saturation. The difference between the induction water saturation corrected for the bound water saturation and the water saturation calculated by the lateral deep resistivity is the increment D. The ratio of the increment D to the gas saturation calculated by the lateral deep resistivity can be used to quantitatively evaluate the movable water saturation. The method realizes quantitative evaluation of the movable water saturation of the reservoir based on the array induction deep resistivity and the lateral logging deep resistivity of the conventional logging data, can quickly and accurately realize continuous quantitative evaluation of the movable water saturation of the reservoir, and can reliably identify the fluid of the reservoir.

[0005] The equivalent additional conductivity caused by the conductive network formed by the water in the inner wall of the rock pore throat in Chinese patent 201910221661.X can obviously improve the calculation result of the hydrocarbon saturation of low-resistance oil (gas) layer, but the result can only explain the oil (gas) saturation and water saturation, and cannot explain the movable water saturation and irreducible water saturation in the water saturation. The difference between the water saturation calculated by the induced water saturation and the lateral deep resistivity and the ratio of the gas saturation calculated by the lateral deep resistivity can be used to quantitatively evaluate the movable water saturation, and the movable water saturation can be explained, and the method is suitable for oil and gas wells with deep and shallow lateral resistivity logging and array induction logging. Array induction is not implemented in the general well logging series of conventional exploration and development wells in general oil and gas fields, and therefore has limitations. SUMMARY

[0006] In view of the problems in the prior art, the present application provides a tight sandstone gas reservoir movable water saturation evaluation method and system, thereby solving the technical problem that the movable water saturation of conventional exploration and development wells and logging cannot be determined in the prior art.

[0007] The present application is realized by the following technical solutions:

[0008] A tight sandstone gas reservoir movable water saturation evaluation method, comprising the following steps:

[0009] The resistivity, acoustic time difference, permeability and water saturation are utilized, and a movable water saturation model is combined to complete the evaluation of the movable water saturation of the tight sandstone gas reservoir;

[0010] The movable water saturation model is

[0011] S MW = S i -S W

[0012] S W = a*RT + b*AC + c*K + d

[0013] In the formula, S MW is the movable water saturation, S i is the water saturation, S W is an irreducible water saturation model, RT is the resistivity, AC is the acoustic time difference, K is the permeability; a, b, c and d are constants.

[0014] Preferably, the establishment of the irreducible water saturation model is specifically:

[0015] S201: Obtain a plurality of groups of logging maximum capillary radius, minimum capillary radius, critical capillary radius, capillary radius, displacement pressure difference and water phase viscosity;

[0016] S202: Obtain a plurality of sets of the tortuosity fractal dimension and the pore fractal dimension of the well logging according to the plurality of sets of the maximum capillary radius, the minimum capillary radius, the critical capillary radius and the capillary radius;

[0017] S203: Obtain a plurality of sets of the irreducible water saturation by using the plurality of sets of the maximum capillary radius, the minimum capillary radius, the critical capillary radius, the capillary radius, the displacement pressure difference, the water phase viscosity, the tortuosity fractal dimension and the pore fractal dimension, and combining the irreducible water saturation model established according to the fractal theory.

[0018] S204: Establish the irreducible water saturation model according to the plurality of sets of the irreducible water saturation and the well logging parameters.

[0019] Preferably, the irreducible water saturation model established according to the fractal theory is as follows:

[0020]

[0021] In the formula, S' is the irreducible water saturation, %; r is the maximum capillary radius, μm; r is the minimum capillary radius, μm; r is the critical capillary radius, μm; r is the capillary radius, μm; D is the tortuosity fractal dimension; D is the pore fractal dimension; Δp is the displacement pressure difference, MPa; μ is the water phase viscosity, MPa·s; and A, m and n are constants. w max min c τ p w

[0022] Preferably, the maximum capillary radius, the minimum capillary radius, the critical capillary radius and the capillary radius are obtained by the mercury injection experiment.

[0023] Preferably, the plurality of sets of the tortuosity fractal dimension and the pore fractal dimension of the well logging are obtained according to the plurality of sets of the maximum capillary radius, the minimum capillary radius, the critical capillary radius and the capillary radius according to the fractal theory.

[0024] Preferably, the well logging parameters include lithology, physical property and gas-bearing property.

[0025] The above-described tight sandstone gas reservoir movable water saturation evaluation method is applied to the establishment of a reservoir classification identification chart.

[0026] A tight sandstone gas reservoir movable water saturation evaluation system comprises:

[0027] ​​​​​​​The data processing module is used for evaluating the movable water saturation of the tight sandstone gas reservoir by using resistivity, acoustic travel time, permeability and water saturation, and combining with a movable water saturation model, wherein the movable water saturation model is

[0028] S MW = S i -S W

[0029] S W = a*RT + b*aC + c*K + d

[0030] In the formula, S MW is the movable water saturation, S i is the water saturation, S W is a bound water saturation model, RT is the resistivity, AC is the acoustic travel time, K is the permeability, and a, b, c and d are constants.

[0031] A terminal device comprises a memory, a processor, and a computer program stored in the memory and executable on the processor, and the processor implements the steps of the above method when executing the computer program.

[0032] A computer readable storage medium stores a computer program, and the computer program implements the steps of the above method when executed by a processor.

[0033] Compared with the prior art, the present application has the following beneficial technical effects:

[0034] Based on the above technical scheme, the tight sandstone gas reservoir movable water saturation evaluation method disclosed in the present application utilizes the parameters such as resistivity, acoustic travel time, permeability and water saturation of conventional logging and a movable water saturation model established in combination, and can effectively realize the evaluation of the movable water saturation of conventional wells, and the evaluation method does not need mercury injection experiment data, but utilizes conventional logging data, and has good universality. The method realizes the movable water saturation interpretation without increasing logging series, and is simple to operate and has wide application range. BRIEF DESCRIPTION OF DRAWINGS

[0035] In order to more clearly illustrate the technical solutions of the embodiments of the present application, the following will briefly introduce the drawings needed to be used in the embodiments, and it should be understood that the following drawings only show some embodiments of the present application, and should not be regarded as a limitation to the scope, and for those skilled in the art, other related drawings can also be obtained without creative labor on the basis of these drawings.

[0036] Figure 1 It is a flowchart of the tight sandstone gas reservoir movable water saturation evaluation method in embodiment 1 of the present application.

[0037] Figure 2 This is a schematic structural diagram of a system for evaluating movable water saturation of a tight sandstone gas reservoir in Example 3 of the present invention;

[0038] Figure 3 is the movable water saturation S established in Example 5 of the present invention MW Explanatory diagram;

[0039] Figure 4 This is a cross-plot of KH and movable water saturation in the application of the movable water saturation obtained in Example 5 of the present invention in establishing a reservoir classification and identification chart. DETAILED DESCRIPTION

[0040] To make the objectives, technical solutions, and advantages of the embodiments of the present invention more clear, the technical solutions of the embodiments of the present invention will be clearly and completely described below in conjunction with the accompanying drawings of the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Generally, the components of the embodiments of the present invention described and shown in the drawings herein can be arranged and designed in various different configurations.

[0041] Therefore, the following detailed description of the embodiments of the present invention provided in the accompanying drawings is not intended to limit the scope of the invention as claimed, but rather merely represents selected embodiments of the present invention. All other embodiments derived by persons of ordinary skill in the art based on the embodiments of the present invention without creative effort shall fall within the scope of protection of the present invention.

[0042] It should be noted that similar reference numerals and letters denote similar items in the following drawings, and therefore, once an item is defined in one drawing, it does not need to be further defined or explained in subsequent drawings.

[0043] In the description of the embodiments of the present invention, it should be noted that if the terms "upper," "lower," "horizontal," "inner," etc. appear, the orientation or positional relationship indicated is based on the orientation or positional relationship shown in the accompanying drawings, or the orientation or positional relationship in which the inventive product is typically placed when in use. These terms are merely for the convenience of describing the present invention and simplifying the description, and do not indicate or imply that the device or component referred to must have a specific orientation, be constructed, or operate in a specific orientation. Therefore, they should not be construed as limitations on the present invention. In addition, the terms "first," "second," etc. are used only to distinguish descriptions and should not be construed as indicating or implying relative importance.

[0044] In addition, if the term "horizontal" appears, it does not mean that the component must be absolutely horizontal, but can be slightly tilted. For example, "horizontal" only means that its direction is more horizontal than "vertical", and does not mean that the structure must be completely horizontal, but can be slightly tilted.

[0045] In the description of the embodiments of the present application, it should also be noted that unless specifically defined and limited, if the terms "arrange", "install", "connect", "connect" appear, they should be understood in a broad sense, for example, they can be fixedly connected, or can be detachably connected, or integrally connected; can be mechanically connected, or can be electrically connected; can be directly connected, or can be indirectly connected through an intermediate medium, or can be connected inside two elements. For those skilled in the art, the specific meaning of the above terms in the present application can be understood according to the specific circumstances.

[0046] The present application will be described in further detail below with reference to the accompanying drawings:

[0047] Embodiment 1

[0048] As Figure 1 shown, the present application discloses a method for evaluating movable water saturation of tight sandstone gas reservoir, comprising the following steps:

[0049] S1: using the resistivity, acoustic time difference, permeability and water saturation, and combining the movable water saturation model to complete the evaluation of the movable water saturation of the tight sandstone gas reservoir;

[0050] The movable water saturation model is

[0051] S MW =S i -S W

[0052] S W =a*RT+b*AC+c*K+d

[0053] In the formula, S MW is the movable water saturation, S i is the water saturation, S W is the irreducible water saturation model, RT is the resistivity, AC is the acoustic time difference, K is the permeability; a, b, c and d are constants. Among them, a, b, c and d are obtained by fitting specific parameters of a specific gas field or block. The parameters obtained by regression for different gas fields may be different. a, b, c and d are obtained by fitting known mercury injection test data and logging test data.

[0054] Among them, the establishment of the irreducible water saturation model is specifically:

[0055] S201: Obtain a plurality of sets of logging maximum capillary radius, minimum capillary radius, critical capillary radius, capillary radius, displacement pressure difference and water phase viscosity; wherein the maximum capillary radius, the minimum capillary radius, the critical capillary radius and the capillary radius are obtained by mercury injection experiment.

[0056] S202: Obtain a plurality of sets of the tortuosity fractal dimension and the pore fractal dimension of the well from the plurality of sets of the maximum capillary radius, the minimum capillary radius, the critical capillary radius, and the capillary radius according to a fractal theory.

[0057] S203: Obtain a plurality of sets of the irreducible water saturation by using the maximum capillary radius, the minimum capillary radius, the critical capillary radius, the capillary radius, the displacement pressure difference, the water phase viscosity, the tortuosity fractal dimension, and the pore fractal dimension, and combining the irreducible water saturation model established according to the fractal theory.

[0058] S204: Establish the irreducible water saturation model according to the plurality of sets of the irreducible water saturation and the logging parameters, wherein the logging parameters include lithology, physical property, and gas-bearing property.

[0059] The irreducible water saturation model established according to the fractal theory is as follows:

[0060]

[0061] In the formula, S' is the irreducible water saturation, %; rmax is the maximum capillary radius, μm; rmin is the minimum capillary radius, μm; rc is the critical capillary radius, μm; rcap is the capillary radius, μm; D is the tortuosity fractal dimension; Dp is the pore fractal dimension; Δp is the displacement pressure difference, MPa; μ is the water phase viscosity, MPa·s; A, m, and n are constants. w max min c τ p w A, m, and n are obtained by fitting specific parameters of a certain gas field or block, and the parameters obtained by regression for different gas fields may be different. A, m, and n are obtained by fitting known mercury injection test data and irreducible water saturation data.

[0062] Embodiment 2

[0063] In another preferred embodiment of the present application, a method for evaluating the movable water saturation of a tight sandstone gas reservoir can include the following steps:

[0064] S1: Obtain the resistivity, acoustic travel time, permeability, and water saturation of a well to be measured;

[0065] S2: Complete the evaluation of the movable water saturation of the tight sandstone gas reservoir by using the resistivity, acoustic travel time, permeability, and water saturation, and combining a movable water saturation model.​​​​​​

[0066] The movable water saturation model is

[0067] S MW = S i - S W

[0068] S W = a*RT + b*AC + c*K + d

[0069] In the formula, S MW is the movable water saturation, S i is the water saturation, S W is the irreducible water saturation model, RT is the resistivity, AC is the acoustic time difference, K is the permeability; a, b, c and d are constants. a, b, c and d are fitted by known mercury injection test data and logging test data.

[0070] The application further discloses application of the movable water saturation evaluation method for the tight sandstone gas reservoir in the application in establishing a reservoir classification identification chart.

[0071] Embodiment 3

[0072] As Figure 2 shown, the application further discloses a movable water saturation evaluation system for a tight sandstone gas reservoir, which comprises:

[0073] The data processing module is used for evaluating the movable water saturation of the tight sandstone gas reservoir by using the resistivity, the acoustic time difference, the permeability and the water saturation and combining the movable water saturation model.

[0074] The movable water saturation obtained in the application is a reservoir property independent of porosity, permeability, gas saturation and reservoir effective thickness, and can be used as one of the reservoir evaluation parameters for the tight sandstone gas reservoir, so that the tight sandstone water-bearing gas reservoir evaluation is more comprehensive and more accurate.

[0075] Embodiment 4

[0076] The movable water saturation evaluation system for the tight sandstone gas reservoir in the application can further comprise the following components:

[0077] The data acquisition module is used for acquiring the resistivity, the acoustic time difference, the permeability and the water saturation of a well to be measured.

[0078] The data processing module is used for evaluating the movable water saturation of the tight sandstone gas reservoir by using the resistivity, the acoustic time difference, the permeability and the water saturation and combining the movable water saturation model.

[0079] A result output module is configured to output the movable water saturation of the tight sand gas reservoir.

[0080] Embodiment 5

[0081] Further to explain the technical scheme of the present application, the following embodiments are used for explanation.

[0082] A movable water saturation interpretation method for a tight sand gas reservoir includes the following steps:

[0083] (1) Obtain a tight sand gas reservoir irreducible water saturation formula. As a porous medium, the tight sand has a fractal characteristic in pore structure. The fractal theory is used to obtain the irreducible water saturation formula.

[0084] The fractal theory is used to obtain the irreducible water saturation formula, and the formula is as follows:

[0085]

[0086] In the formula, S' w is the irreducible water saturation according to the fractal theory, %; r max is the maximum capillary radius, μm; r min is the minimum capillary radius, μm; r c is the critical capillary radius, μm; r is the capillary radius, μm; D τ is the tortuosity fractal dimension; D p is the pore fractal dimension; Δp is the displacement pressure difference, MPa; μ w is the water phase viscosity, MPa·s.

[0087] (2) Obtain the core irreducible water saturation. According to a tight sand reservoir irreducible water saturation prediction method based on the fractal theory, a large amount of mercury injection experiment data is used to calculate the irreducible water saturation.

[0088] (3) Establish the relationship between the calculated irreducible water saturation and the logging parameters. The calculated irreducible water saturation is fitted with different logging parameters (lithology, physical property, gas bearing property, and other logging parameters related to water saturation) to establish multiple regression formulas, and the optimal regression formula is selected as the irreducible water saturation interpretation formula. The established irreducible water saturation fitting regression formula is shown in Table 1, and the formula with the highest correlation coefficient is selected as the irreducible water saturation interpretation formula, i.e., the optimal formula 11.

[0089] Table 1 Irreducible water saturation fitting regression formula table

[0090]

[0091] In the formula, S WFor irreducible water saturation model, RT is resistivity, Ω.m, AC is acoustic travel time, mu s / m, K is permeability, mD, GR is natural gamma, API, For porosity, %, SH is shale content, %, DEN is density, g / cm 3 .

[0092] (4) obtain movable water saturation interpretation formula. Combined with the original water saturation interpretation formula, subtract the calculated irreducible water saturation, obtain the movable water saturation, and obtain the movable water saturation interpretation formula;

[0093] S MW = S i -S W

[0094] S W =-0.0793*RT+0.3159*AC-12.4143*K-35.21

[0095] The movable water saturation is used for gas well logging interpretation, and the movable water saturation S MW The interpretation is as shown in Figure 3 , so that the reservoir evaluation of tight sandstone water-bearing gas reservoirs is more comprehensive and accurate.

[0096] At the same time, in one application of the present application, the formation coefficient (KH-permeability x gas layer thickness) can be introduced, associated with the movable water saturation, to obtain a reservoir classification identification chart for classification and evaluation of the reservoir. The evaluation results are shown in Figure 4 , and Figure 4 It can be seen that according to the movable water saturation and the formation coefficient, the gas wells are classified, so that both the static geological characteristics of the gas wells and the influence of water production on the production capacity of the gas wells are considered. After the introduction of the movable water saturation, the dynamic and static coincidence rate of the gas wells is increased by 24.3%, which is of great significance for fine classification and management of the gas wells; at the same time, by using the reservoir classification identification chart, the effect of perforation optimization and old well layer checking and perforation supplementing can be effectively guided, and the effect is remarkable.

[0097] An embodiment of the terminal device provided by the present application provides a schematic diagram of the terminal device. The terminal device of the embodiment comprises a processor, a memory, and a computer program stored in the memory and executable on the processor. The processor implements the steps in each of the method embodiments when executing the computer program. Alternatively, the processor implements the functions of each module / unit in each of the device embodiments when executing the computer program.

[0098] 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 application.

[0099] The terminal device can be a desktop computer, a notebook computer, a palm computer, a cloud server, and the like. The terminal device can include, but is not limited to, a processor and a memory.

[0100] The processor can be a central processing unit (CPU), and can also be other general-purpose processors, a digital signal processor (DSP), an application-specific integrated circuit (ASIC), a field-programmable gate array (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, and the like.

[0101] The memory can be used to store the computer program and / or modules, and the processor can realize various functions of the terminal device by running or executing the computer program and / or modules stored in the memory, and calling data stored in the memory.

[0102] The modules / units integrated in the terminal device, if realized in the form of software function units and sold or used as independent products, can be stored in a computer-readable storage medium. Based on this understanding, all or part of the processes in the above-mentioned embodiments can also be completed by a computer program instructing related hardware, and the computer program can be stored in a computer-readable storage medium. When the processor executes the computer program, the steps of the above-mentioned various method embodiments can be realized. The computer program includes computer program code, which can be in the form of source code, object code, executable files, or some intermediate forms, etc. The computer-readable medium can include any entity or device capable of carrying the computer program code, recording medium, U disk, mobile hard disk, magnetic disk, optical disk, computer memory, read-only memory (ROM), random access memory (RAM), electrical carrier signals, telecommunication signals, and software distribution medium, etc. It should be noted that the contents included in the computer-readable medium can be appropriately increased or decreased according to the requirements of legislation and patent practice in the jurisdiction, for example, in some jurisdictions, according to legislation and patent practice, the computer-readable medium does not include electrical carrier signals and telecommunication signals.

[0103] The above are merely preferred embodiments of the present invention and are not intended to limit the present invention. Those skilled in the art will readily appreciate that various modifications and variations of the present invention are possible. Any modifications, equivalent substitutions, or improvements made within the spirit and principles of the present invention shall be included within the scope of protection of the present invention.

Claims

1. A method for evaluating movable water saturation in tight sandstone gas reservoirs, characterized in that: The following steps are involved: The movable water saturation of the tight sandstone gas reservoir is evaluated by using resistivity, acoustic transit time, permeability and water saturation in combination with a movable water saturation model; The movable water saturation model is Where, is the movable water saturation, is the water saturation, is the irreducible water saturation model, is the resistivity, is the time difference of sound waves, is the permeability; 、 、 as well as is a constant; The establishment of the irreducible water saturation model is specifically as follows: S201: Obtaining the maximum capillary radius, minimum capillary radius, critical capillary radius, capillary radius, displacement pressure difference, and water phase viscosity of several groups of well logging data; S202: Obtaining a plurality of sets of tortuosity fractal dimensions and pore fractal dimensions of well logging according to a plurality of sets of maximum capillary radius, minimum capillary radius, critical capillary radius, and capillary radius; S203: using the plurality of sets of the maximum capillary radius, the minimum capillary radius, the critical capillary radius, the capillary radius, the displacement pressure difference, the water phase viscosity, the tortuosity fractal dimension, and the pore fractal dimension, in combination with an irreducible water saturation model established based on fractal theory, to obtain a plurality of sets of irreducible water saturations; S204: Establishing the irreducible water saturation model according to the plurality of irreducible water saturation groups and logging parameters.

2. The method for evaluating movable water saturation of a tight sandstone gas reservoir according to claim 1, wherein: The irreducible water saturation model established based on fractal theory is: Where: is the irreducible water saturation established based on fractal theory, %; is the maximum capillary radius, ; is the minimum capillary radius, ; is the critical capillary radius, ; is the capillary radius, ; is the tortuosity fractal dimension; is the pore fractal dimension; is the displacement pressure difference, ; is the viscosity of the water phase, ; 、 as well as is a constant.

3. The method for evaluating movable water saturation of a tight sandstone gas reservoir according to claim 1, wherein: The maximum capillary radius, the minimum capillary radius, the critical capillary radius and the capillary radius are obtained through mercury injection experiments.

4. The method for evaluating movable water saturation of a tight sandstone gas reservoir according to claim 1, wherein: According to fractal theory, several groups of tortuosity fractal dimensions and pore fractal dimensions of the well logging are obtained through several groups of maximum capillary radius, minimum capillary radius, critical capillary radius and capillary radius.

5. The method for evaluating movable water saturation of a tight sandstone gas reservoir according to claim 1, wherein: The logging parameters include lithology, physical properties and gas content.

6. Application of the method for evaluating movable water saturation of tight sandstone gas reservoirs according to any one of claims 1 to 5 in establishing a reservoir classification identification chart.

7. A system for evaluating movable water saturation in tight sandstone gas reservoirs, characterized in that: The method for evaluating movable water saturation of a tight sandstone gas reservoir according to any one of claims 1 to 5 comprises: Data processing module: The data processing module is used to evaluate the movable water saturation of the tight sandstone gas reservoir using resistivity, acoustic wave transit time, permeability and water saturation in combination with the movable water saturation model. The movable water saturation model is Where, is the movable water saturation, is the water saturation, is the irreducible water saturation model, is the resistivity, is the time difference of sound waves, is the permeability; 、 、 as well as is a constant.

8. A terminal device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: When the processor executes the computer program, the steps of the method according to any one of claims 1 to 5 are implemented.

9. A computer-readable storage medium storing a computer program, characterized in that: When the computer program is executed by a processor, the steps of the method according to any one of claims 1 to 5 are implemented.

Images