Method, device, equipment and medium for establishing absorption attenuation field of desert surface layer

Abstract

The present disclosure provides a method, device, equipment and medium for establishing an absorption attenuation field of a desert surface layer, belonging to the technical field of seismic exploration. The method comprises: dividing the structure of the desert surface layer into a first medium structure layer and a second medium structure layer, establishing a characteristic function of the first medium structure layer to represent the variation law of the absorption attenuation parameters of the first medium structure layer, and establishing a characteristic function of the second medium structure layer to represent the variation law of the absorption attenuation parameters of the second medium structure layer. Since the structure of the desert surface layer is complex, the desert surface layer is divided into layers, and a corresponding characteristic function is established for each medium structure layer, which can improve the calculation accuracy of the Q value of the physical point in each medium structure layer, and further improve the accuracy of the established absorption attenuation field.

Description

TECHNICAL FIELD

[0001] The present disclosure relates to the technical field of seismic exploration, and particularly relates to a method and device for establishing an absorption attenuation field of a desert area surface layer, equipment and a medium. BACKGROUND

[0002] The principle of seismic exploration is to study the geological structure by studying the propagation characteristics of seismic waves in geology, so as to explore mineral resources such as oil and natural gas in geology. However, the desert area surface layer has large structure fluctuations and thick quicksand layer, resulting in serious absorption and attenuation of seismic wave energy, low signal-to-noise ratio and resolution of seismic records. By studying and characterizing the change rule of the Q (quality) value of the desert area surface layer structure, and then establishing the absorption attenuation field of the desert area surface layer, the absorption and attenuation effect of the desert area surface layer on the seismic wave energy can be weakened or even eliminated, so as to improve the seismic resolution and improve the accuracy of mineral resource exploration in the desert area.

[0003] In the related art, a single characteristic function is used to represent the change rule of the Q value of the desert area surface layer structure, and the absorption attenuation field is established according to the Q value calculated by the characteristic function.

[0004] In the process of implementing the present disclosure, the inventors found that the prior art at least has the following problems:

[0005] The structure of the desert area surface layer is complex, and the Q value of the desert area surface layer calculated by using a single characteristic function has low precision, which further reduces the precision of the established absorption attenuation field. SUMMARY

[0006] The embodiments of the present disclosure provide a method and device for establishing an absorption attenuation field of a desert area surface layer, equipment and a medium, which can improve the calculation precision of the Q value of the desert area surface layer, and further improve the precision of the established absorption attenuation field. The technical solution is as follows:

[0007] In a first aspect, a method for establishing an absorption attenuation field of a desert area surface layer is provided, the desert area surface layer includes a first medium structure layer and a second medium structure layer, the first medium structure layer is a structure layer above the water table, and the second medium structure layer is a structure layer below the water table; the desert area surface layer is provided with a plurality of physical points; and the method includes:

[0008] determining a characteristic function of the first medium structure layer and a characteristic function of the second medium structure layer, the characteristic function of the first medium structure layer being used to represent a variation law of an absorption attenuation parameter of the first medium structure layer, the characteristic function of the second medium structure layer being used to represent a variation law of an absorption attenuation parameter of the second medium structure layer, the characteristic function of the first medium structure layer being different from the characteristic function of the second medium structure layer; calculating a first absorption attenuation parameter value of each of the plurality of physical points according to the characteristic function of the first medium structure layer; calculating a second absorption attenuation parameter value of each of the plurality of physical points according to the characteristic function of the second medium structure layer; and establishing an absorption attenuation field of the desert area surface layer according to the first absorption attenuation parameter values and the second absorption attenuation parameter values of the plurality of physical points.

[0009] Optionally, the first medium structure layer has a plurality of different characteristic functions; the calculating of the first absorption attenuation parameter value of each of the plurality of physical points according to the characteristic function of the first medium structure layer comprises: obtaining a dune parameter of a position where a target physical point is located, the target physical point being one of the plurality of physical points; determining a target characteristic function of the target physical point according to the dune parameter, the target characteristic function being at least one of the plurality of characteristic functions of the first medium structure layer; and calculating the first absorption attenuation parameter value of the target physical point according to the target characteristic function.

[0010] Optionally, the dune parameter comprises a dune thickness, a slope and an elevation; the determining of the characteristic function of the first medium structure layer corresponding to the target physical point according to the dune parameter comprises: determining the target characteristic function according to at least one of a ratio of the dune thickness to the slope and an included angle between a wind direction and a tendency, the tendency being calculated according to the elevation.

[0011] Optionally, the first medium structure layer comprises a first sub-structure layer and a second sub-structure layer in a longitudinal direction, the first sub-structure layer being located above the second sub-structure layer; the determining of the target characteristic function according to at least one of the ratio of the dune thickness to the slope and the included angle between the main wind direction and the tendency comprises at least one of: in response to the ratio of the dune thickness to the slope being less than a first ratio threshold, determining the target characteristic function as a first characteristic function, the first characteristic function being used to calculate a first absorption attenuation parameter value of the first sub-structure layer; the first characteristic function being: Q(H) = A1·H α1wherein A1=0.0014S+0.4849, a1=0.0021S+0.4141, S represents the dune reference parameter of the target physical point, the dune reference parameter of the target physical point being calculated by the dune parameter of the location where the target physical point is located; in response to the ratio of the dune thickness to the slope being greater than a second ratio threshold value, and the included angle between the main wind direction and the tendency being less than a first included angle threshold value, determining the target characteristic function as a second characteristic function and a third characteristic function, the second characteristic function and the third characteristic function being respectively used to calculate the first absorption attenuation parameter value of the first substructure layer and the first absorption attenuation parameter value of the second substructure layer; the first ratio threshold value being less than the second ratio threshold value; the second characteristic function being: Q(H)=A2·H α2 wherein A2=-0.0017S+1.8517, a2=-0.00003S+0.9431, S represents the dune reference parameter of the target physical point, the dune reference parameter of the target physical point being calculated by the dune parameter of the location where the target physical point is located; the third characteristic function being: Q(H)=A3·H α3 wherein A3=0.0008S+1.7647, a3=0.0002S+0.7284, S represents the dune reference parameter of the target physical point, the dune reference parameter of the target physical point being calculated by the dune parameter of the location where the target physical point is located; in response to the ratio of the dune thickness to the slope being greater than a second ratio threshold value, and the included angle between the main wind direction and the tendency being greater than a second included angle threshold value, determining the target characteristic function as a fourth characteristic function and a fifth characteristic function, the fourth characteristic function and the fifth characteristic function being respectively used to calculate the first absorption attenuation parameter value of the first substructure layer and the first absorption attenuation parameter value of the second substructure layer; the first included angle threshold value being less than the second included angle threshold value; the fourth characteristic function being: Q(H)=A4·H α4 wherein A4=-0.0003S+1.1162, a4=0.0001S+0.9437, S represents the dune reference parameter of the target physical point, the dune reference parameter of the target physical point being calculated by the dune parameter of the location where the target physical point is located; the fifth characteristic function being: Q(H)=A5·H α5wherein A5=0.0039S+0.0479, a5=-0.0006S+1.0386, S represents the dune reference parameter of the target physical point, which is calculated by the dune parameter of the location where the target physical point is located; in response to the ratio of the dune thickness to the slope being greater than a second ratio threshold value, the included angle between the main wind direction and the tendency being greater than a first included angle threshold value and less than a second included angle threshold value, determining the target feature function as a sixth feature function, the sixth feature function being used to calculate a first absorption attenuation parameter value of the first sub-structure layer and a first absorption attenuation parameter value of the second sub-structure layer; the sixth feature function is: Q=[Q B *|Ψ-Ψ B0 |+Q Y *|Ψ-Ψ Y0 |] / |Ψ Y0 -Ψ B0 |, wherein Ψ B0 represents the first included angle threshold value, Ψ Y0 represents the second included angle threshold value, and Ψ represents the included angle between the main wind direction and the inclination; when Q B is the first absorption attenuation parameter value of the first sub-structure layer calculated by using the second feature function and Q Y represents the first absorption attenuation parameter value of the first sub-structure layer calculated by using the fourth feature function, Q represents the calculated first absorption attenuation parameter value of the first sub-structure layer; when Q B is the first absorption attenuation parameter value of the first sub-structure layer calculated by using the third feature function and Q Y represents the first absorption attenuation parameter value of the first sub-structure layer calculated by using the fifth feature function, Q represents the calculated first absorption attenuation parameter value of the second sub-structure layer; in response to the ratio of the dune thickness to the slope being greater than a first ratio threshold value and less than a second ratio threshold value, determining the target feature function as a seventh feature function, the seventh feature function being used to calculate the first absorption attenuation parameter value of the first sub-structure layer; the seventh feature function is: Q=[Q H *|P S -P S0 |+Q P *|P S -P S1 |] / |P S0 -P S1 |, wherein P S0 represents the first ratio threshold value, P S1 represents the second ratio threshold value, P S represents the ratio of the dune thickness to the slope, and Q Hrepresenting a first absorption attenuation parameter value of the first sub-structure layer calculated by using the sixth characteristic function, Q P representing a first absorption attenuation parameter value of the first sub-structure layer calculated by using the first characteristic function.

[0012] Optionally, the calculating the second absorption attenuation parameter values of the plurality of physical points according to the characteristic function of the second medium structure layer comprises: calculating the second absorption attenuation parameter value of the target physical point according to an eighth characteristic function: Q = -0.0004S + 0.4484S - 74.86, wherein S represents a dune reference parameter of the target physical point, and the dune reference parameter of the target physical point is calculated by a dune parameter of a position where the target physical point is located. 2 +0.4484S - 74.86, wherein S represents a dune reference parameter of the target physical point, and the dune reference parameter of the target physical point is calculated by a dune parameter of a position where the target physical point is located.

[0013] Optionally, the desert surface layer is provided with a plurality of measuring points, and the method further comprises: correcting the first absorption attenuation parameter value of the target physical point according to a first absorption attenuation parameter difference value of the plurality of measuring points to obtain a corrected first absorption attenuation parameter value of the target physical point, the first absorption attenuation parameter difference value being a difference value between a third absorption attenuation parameter value of the measuring point measured and a fourth absorption attenuation parameter value of the measuring point calculated; correcting the second absorption attenuation parameter value of the target physical point according to a second absorption attenuation parameter difference value of the plurality of measuring points to obtain a corrected second absorption attenuation parameter value of the target physical point, the second absorption attenuation parameter difference value being a difference value between a fifth absorption attenuation parameter value of the measuring point measured and a sixth absorption attenuation parameter value of the measuring point calculated; and establishing an absorption attenuation field of the desert surface layer according to the corrected first absorption attenuation parameter value and the corrected second absorption attenuation parameter value of the target physical point.

[0014] Optionally, the establishing the absorption attenuation field of the desert surface layer according to the corrected first absorption attenuation parameter value and the corrected second absorption attenuation parameter value of the target physical point comprises: calculating a first relative absorption attenuation total amount of the target physical point according to the corrected first absorption attenuation parameter value and the corrected second absorption attenuation parameter value of the target physical point; calculating second relative absorption attenuation total amounts of the plurality of measuring points according to the measured third absorption attenuation parameter value and the fifth absorption attenuation parameter value of the plurality of measuring points; calculating third relative absorption attenuation total amounts of the plurality of measuring points according to the calculated fourth absorption attenuation parameter value and the sixth absorption attenuation parameter value of the plurality of measuring points; obtaining relative absorption attenuation total amount difference values of the plurality of measuring points according to the second relative absorption attenuation total amounts and the third relative absorption attenuation total amounts; correcting the first relative absorption attenuation total amount of the target physical point according to the relative absorption attenuation total amount difference values of the plurality of measuring points to obtain a corrected first relative absorption attenuation total amount of the target physical point; and establishing the absorption attenuation field of the desert surface layer according to the corrected first relative absorption attenuation total amount of the target physical point.

[0015] In a second aspect, provided is a device for establishing an absorption attenuation field of a desert surface layer, the desert surface layer comprising a first medium structure layer and a second medium structure layer, the first medium structure layer being a structure layer above a water table, and the second medium structure layer being a structure layer below the water table; the desert surface layer being provided with a plurality of physical points; the device comprising: a determination module configured to determine a characteristic function of the first medium structure layer and a characteristic function of the second medium structure layer, the characteristic function of the first medium structure layer being used to represent a variation law of absorption attenuation parameters of the first medium structure layer, the characteristic function of the second medium structure layer being used to represent a variation law of absorption attenuation parameters of the second medium structure layer, and the characteristic function of the first medium structure layer being different from the characteristic function of the second medium structure layer; a first calculation module configured to calculate first absorption attenuation parameter values of the plurality of physical points according to the characteristic function of the first medium structure layer; a second calculation module configured to calculate second absorption attenuation parameter values of the plurality of physical points according to the characteristic function of the second medium structure layer; and an establishment module configured to establish the absorption attenuation field of the desert surface layer according to the first absorption attenuation parameter values and the second absorption attenuation parameter values of the plurality of physical points.

[0016] Optionally, the first medium structure layer has a plurality of different characteristic functions; the first calculation module is configured to obtain a dune parameter of a position of a target physical point, the target physical point being one of the plurality of physical points; determine a target characteristic function of the target physical point according to the dune parameter, the target characteristic function being at least one of the plurality of characteristic functions of the first medium structure layer; and calculate a first absorption attenuation parameter value of the target physical point according to the target characteristic function.

[0017] Optionally, the dune parameter includes a dune thickness, a slope and an elevation; the first calculation module is configured to determine the target characteristic function according to at least one of a ratio of the dune thickness to the slope and an included angle between a wind direction and a tendency, the tendency being calculated according to the elevation.

[0018] Optionally, the first medium structure layer includes a first sub-structure layer and a second sub-structure layer in a longitudinal direction, the first sub-structure layer being located above the second sub-structure layer; the first calculation module is configured to determine the target characteristic function as a first characteristic function in response to the ratio of the dune thickness to the slope being less than a first ratio threshold value, the first characteristic function being used to calculate a first absorption attenuation parameter value of the first sub-structure layer; the first characteristic function being: Q(H)=A1·H α1 wherein A1=0.0014S+0.4849, α1=0.0021S+0.4141, S representing a dune reference parameter of the target physical point, the dune reference parameter of the target physical point being calculated by the dune parameter of the position of the target physical point; determine the target characteristic function as a second characteristic function and a third characteristic function in response to the ratio of the dune thickness to the slope being greater than a second ratio threshold value and the included angle between the main wind direction and the tendency being less than a first included angle threshold value, the second characteristic function and the third characteristic function being used to calculate a first absorption attenuation parameter value of the first sub-structure layer and a first absorption attenuation parameter value of the second sub-structure layer respectively; the first ratio threshold value being less than the second ratio threshold value; the second characteristic function being: Q(H)=A2·H α2 wherein A2=-0.0017S+1.8517, α2=-0.00003S+0.9431, S representing a dune reference parameter of the target physical point, the dune reference parameter of the target physical point being calculated by the dune parameter of the position of the target physical point; the third characteristic function being: Q(H)=A3·H α3wherein A3=0.0008S+1.7647, a3=0.0002S+0.7284, S represents the dune reference parameter of the target physical point, the dune reference parameter of the target physical point being calculated by the dune parameter of the location where the target physical point is located; in response to the ratio of the dune thickness to the slope being greater than a second ratio threshold value, and the included angle between the main wind direction and the tendency being greater than a second included angle threshold value, determining the target feature function as a fourth feature function and a fifth feature function, the fourth feature function and the fifth feature function being respectively used to calculate the first absorption attenuation parameter value of the first sub-structure layer and the first absorption attenuation parameter value of the second sub-structure layer; the first included angle threshold value is less than the second included angle threshold value; the fourth feature function is: Q(H)=A4·H α4 wherein A4=-0.0003S+1.1162, a4=0.0001S+0.9437, S represents the dune reference parameter of the target physical point, the dune reference parameter of the target physical point being calculated by the dune parameter of the location where the target physical point is located; the fifth feature function is: Q(H)=A5·H α5 wherein A5=0.0039S+0.0479, a5=-0.0006S+1.0386, S represents the dune reference parameter of the target physical point, the dune reference parameter of the target physical point being calculated by the dune parameter of the location where the target physical point is located; in response to the ratio of the dune thickness to the slope being greater than a second ratio threshold value, the included angle between the main wind direction and the tendency being greater than a first included angle threshold value and less than a second included angle threshold value, determining the target feature function as a sixth feature function, the sixth feature function being used to calculate the first absorption attenuation parameter value of the first sub-structure layer and the first absorption attenuation parameter value of the second sub-structure layer; the sixth feature function is: Q=[Q B *|Ψ-Ψ B0 |+Q Y *|Ψ-Ψ Y0 |] / |Ψ Y0 -Ψ B0 |, wherein Ψ B0 represents the first included angle threshold value, Ψ Y0 represents the second included angle threshold value, and Ψ represents the included angle between the main wind direction and the inclination, when Q B is the first absorption attenuation parameter value of the first sub-structure layer calculated by using the second feature function, and Q Y represents the first absorption attenuation parameter value of the first sub-structure layer calculated by using the fourth feature function, Q represents the calculated first absorption attenuation parameter value of the first sub-structure layer; when Q B is the first absorption attenuation parameter value of the first sub-structure layer calculated by using the third feature function, and QY Q represents the calculated first absorption attenuation parameter value of the second substructure layer; in response to the ratio of the dune thickness and the slope being greater than a first ratio threshold and less than a second ratio threshold, determining the target characteristic function as a seventh characteristic function, the seventh characteristic function being used to calculate the first absorption attenuation parameter value of the first substructure layer; the seventh characteristic function is: Q = [Q H *|P S -P S0 |+Q P *|P S -P S1 |] / |P S0 -P S1 |, wherein P S0 represents the first ratio threshold, P S1 represents the second ratio threshold, P S represents the ratio of the dune thickness and the slope, Q H represents the first absorption attenuation parameter value of the first substructure layer calculated by the sixth characteristic function, Q P represents the first absorption attenuation parameter value of the first substructure layer calculated by the first characteristic function.

[0019] Optionally, the second calculation module is configured to calculate the second absorption attenuation parameter value of the target physical point according to an eighth characteristic function: Q = -0.0004S 2 +0.4484S-74.86, wherein S represents a dune reference parameter of the target physical point, the dune reference parameter of the target physical point being calculated by a dune parameter of a position where the target physical point is located.

[0020] Optionally, the desert area surface layer is provided with a plurality of measuring points, the establishing module is configured to correct the first absorption attenuation parameter value of the target physical point according to a first absorption attenuation parameter difference value of the plurality of measuring points, to obtain a corrected first absorption attenuation parameter value of the target physical point, the first absorption attenuation parameter difference value being a difference value between a third absorption attenuation parameter value of the measuring point measured and a fourth absorption attenuation parameter value of the measuring point calculated; correct the second absorption attenuation parameter value of the target physical point according to a second absorption attenuation parameter difference value of the plurality of measuring points, to obtain a corrected second absorption attenuation parameter value of the target physical point; the second absorption attenuation parameter difference value being a difference value between a fifth absorption attenuation parameter value of the measuring point measured and a sixth absorption attenuation parameter value of the measuring point calculated; and establish the absorption attenuation field of the desert area surface layer according to the corrected first absorption attenuation parameter value and the corrected second absorption attenuation parameter value of the target physical point.

[0021] Optionally, the establishing module is configured to calculate a first relative absorption attenuation total amount of the target physical point according to the corrected first absorption attenuation parameter value and the corrected second absorption attenuation parameter value of the target physical point; calculate a second relative absorption attenuation total amount of the plurality of measuring points according to the third absorption attenuation parameter value and the fifth absorption attenuation parameter value of the plurality of measuring points measured; calculate a third relative absorption attenuation total amount of the plurality of measuring points according to the fourth absorption attenuation parameter value and the sixth absorption attenuation parameter value of the plurality of measuring points calculated; obtain a relative absorption attenuation total amount difference value of the plurality of measuring points according to the second relative absorption attenuation total amount and the third relative absorption attenuation total amount; correct the first relative absorption attenuation total amount of the target physical point according to the relative absorption attenuation total amount difference value of the plurality of measuring points, to obtain a corrected first relative absorption attenuation total amount of the target physical point; and establish the absorption attenuation field of the desert area according to the corrected first relative absorption attenuation total amount of the target physical point.

[0022] In a third aspect, a computer device is provided, including: a processor; a memory for storing processor-executable instructions; and wherein the processor is configured to execute the method of the first aspect.

[0023] In a fourth aspect, a computer-readable storage medium is provided, which, when instructions in the computer-readable storage medium are executed by a processor of a computer device, enables the computer device to execute the method of the first aspect.

[0024] The technical scheme provided by the embodiments of the present disclosure has the following beneficial effects:

[0025] In the embodiments of the present disclosure, the desert area surface layer includes a first medium structure layer and a second medium structure layer. Since the structure of the desert area surface layer is complex, the desert area surface layer is layered, and the characteristic function of the first medium structure layer and the characteristic function of the second medium structure layer are determined respectively, so that the variation law of the absorption attenuation parameter of the first medium structure layer and the variation law of the absorption attenuation parameter of the second medium structure layer can be more accurately represented, thereby improving the calculation accuracy of the Q value of the desert area surface layer, and further improving the accuracy of the established absorption attenuation field. BRIEF DESCRIPTION OF DRAWINGS

[0026] In order to more clearly illustrate the technical solutions in the embodiments of the present disclosure, the drawings needed to be used in the embodiments will be briefly introduced. Obviously, the drawings in the following description are only some embodiments of the present disclosure, and other drawings can be obtained by those skilled in the art without creative effort.

[0027] Figure 1 is a flowchart of a method for establishing an absorption attenuation field of a desert area surface layer provided by an embodiment of the present disclosure;

[0028] Figure 2 is a flowchart of another method for establishing an absorption attenuation field of a desert area surface layer provided by an embodiment of the present disclosure;

[0029] Figure 3 is a schematic diagram of a region type of a desert area surface layer provided by an embodiment of the present disclosure;

[0030] Figure 4 is a fitting result diagram of a characteristic function of a windward area provided by an embodiment of the present disclosure;

[0031] Figure 5 is a fitting result diagram of a characteristic function of a second medium structure layer provided by an embodiment of the present disclosure;

[0032] Figure 6 is a compensation effect diagram of an absorption attenuation field provided by an embodiment of the present disclosure;

[0033] Figure 7 is a flowchart of another method for establishing an absorption attenuation field of a desert area surface layer provided by an embodiment of the present disclosure;

[0034] Figure 8 is a structural schematic diagram of an establishment device for an absorption attenuation field of a desert area surface layer provided by an embodiment of the present disclosure;

[0035] Figure 9 is a structural block diagram of a computer device provided by an embodiment of the present disclosure. DETAILED DESCRIPTION

[0036] In order to make the purpose, technical solutions and advantages of the present disclosure clearer, the present disclosure embodiments will be further described in detail below with reference to the drawings.

[0037] The absorption attenuation parameter (Q) value represents the strength of the absorption characteristics of the medium, and is an inherent characteristic of the medium. To accurately establish the absorption attenuation field of the desert area surface layer, the variation law of the Q value of the desert area surface layer in the lateral direction and in the longitudinal direction needs to be reasonably understood and represented. However, the spatial structure of the desert area surface layer changes dramatically in the lateral direction and in the longitudinal direction, and it is very difficult to establish the absorption attenuation field of the desert area surface layer.

[0038] In the present disclosure embodiments, the desert area surface layer structure includes a first medium structure layer and a second medium structure layer. The first medium structure layer is a structure layer above the phreatic surface, and the second medium structure layer is a structure layer below the phreatic surface. The phreatic water refers to the underground water within the saturated layer, buried above the first stable water-resisting layer, and having a free water surface. The free surface of the phreatic water is called the phreatic surface.

[0039] The first medium structure layer is a low velocity drop layer, which is a sand layer without water or containing water but not rich in water. Influenced by the compaction degree and water content of the sand body, the medium characteristics of the first medium structure layer are that the propagation velocity of the seismic wave changes with the depth, and the change is nonlinear, so a nonlinear characteristic function needs to be used to represent the variation relationship between the Q value and the depth in the low velocity drop layer.

[0040] The second medium structure layer is a high velocity layer, which is a water-rich sand layer or a lithologic layer with a propagation velocity of the seismic wave of 1500 m / s or more. The medium structure characteristics of the second medium structure layer are that the propagation velocity of the seismic wave basically does not change with the depth or the propagation velocity of the seismic wave changes little with the depth, and the change of the propagation time of the seismic wave with the depth is linear, and the Q value does not change with the depth.

[0041] Because the medium structure characteristics of the first medium structure layer and the second medium structure layer are different, the variation law of the Q value is different, and a single characteristic function is used to represent the variation law of the Q value of the first medium structure layer and the second medium structure layer, there is a large error, and the absorption attenuation field of the desert area surface layer is not accurate enough, which is not conducive to eliminating the absorption and attenuation effect of the desert area surface layer on the seismic wave.

[0042] Therefore, the present disclosure embodiments provide a method for establishing the absorption attenuation field of the desert area surface layer, which can establish a more accurate absorption attenuation field.

[0043] Figure 1 is a flowchart of a method for establishing the absorption attenuation field of the desert area surface layer provided by the present disclosure embodiments. The method can be executed by a computer device, such as a computer, etc. Referring to Figure 1 , the method comprises the following steps.

[0044] In step 101, a characteristic function of the first medium structure layer and a characteristic function of the second medium structure layer are determined.

[0045] The characteristic function of the first medium structure layer is used to represent the variation law of the absorption attenuation parameter of the first medium structure layer, and the characteristic function of the second medium structure layer is used to represent the variation law of the absorption attenuation parameter of the second medium structure layer. The characteristic function of the first medium structure is different from the characteristic function of the second medium structure layer.

[0046] In step 102, a first absorption attenuation parameter value of each physical point is calculated according to the characteristic function of the first medium structure layer.

[0047] The desert surface layer is provided with a plurality of physical points at intervals, and each physical point corresponds to a shot point or a receiver point. Exemplarily, each physical point is arranged at an interval of 5 meters.

[0048] In step 103, a second absorption attenuation parameter value of each physical point is calculated according to the characteristic function of the second medium structure layer.

[0049] In step 104, an absorption attenuation field of the desert surface layer is established according to the first absorption attenuation parameter value and the second absorption attenuation parameter value of each physical point.

[0050] In the embodiment of the present disclosure, the desert surface layer includes the first medium structure layer and the second medium structure layer. Since the structure of the desert surface layer is complex, the desert surface layer is layered, and the characteristic function of the first medium structure layer and the characteristic function of the second medium structure layer are determined respectively, which can more accurately represent the variation law of the absorption attenuation parameter of the first medium structure layer and the variation law of the absorption attenuation parameter of the second medium structure layer, thereby improving the calculation accuracy of the Q value of the desert surface layer, and further improving the accuracy of the established absorption attenuation field.

[0051] Figure 2 is a flow chart of another method for establishing an absorption attenuation field of a desert surface layer provided by the present disclosure. The method can be executed by a computer device, such as a computer, etc. Referring to Figure 2 , the method comprises:

[0052] In step 201, a characteristic function of the first medium structure layer is determined.

[0053] The characteristic function of the first medium structure layer is used to represent the variation law of the Q value of the first medium structure layer.

[0054] The first medium structure layer mainly includes dune-shaped sand dunes and honeycomb-shaped sand dunes and other sand dune structures.

[0055] Take the Tarim Basin as an example, more than 80% of the surface of the Taklimakan Desert in the Tarim Basin is covered by mobile aeolian dunes. The Taklimakan Desert can be divided into large sand areas, small and medium-sized sand areas, piedmont zones, and river gobi areas. For example, the Gucheng area in the middle of the Tarim Basin, under the influence of regional monsoons, presents a long and continuous dune with a height of 100-200 meters. The dune is often interspersed with 10-30 meter or even 50 meter high honeycomb dunes. On the one hand, the shapes and structures of different dunes are different, and the Q value variation law differs greatly. On the other hand, the Q value variation law also differs at different spatial positions of each dune, such as the windward and leeward areas of the dune, and the main part and wing of the dune. The shapes and structures of the surface dunes in the desert area are diverse, making it difficult to meet the demand for fine characterization of the Q value of every point of each dune under the premise of economic benefits.

[0056] It is found through research that wind is an important factor in shaping the surface morphology of the desert area. Under the influence of the main wind direction of the desert area, the dune presents a distinction between flat areas, windward areas, and leeward areas. From a horizontal perspective, compared with the windward area, the slope of the leeward area is relatively steep, and the slope of the flat area is relatively gentle. From the perspective of compaction degree, the windward area is greater than the leeward area. Whether it is a ridge-shaped dune or a honeycomb-shaped dune, it can be divided into windward areas, leeward areas, and flat areas in the horizontal direction.

[0057] Figure 3 FIG. 1 is a schematic diagram of a regional type of a surface layer of a desert area provided by an embodiment of the present disclosure. Figure 3 In FIG. 1, the region corresponding to reference numeral 1 is a flat area, the region corresponding to reference numeral 2 is a windward slope of a windward area, the region corresponding to reference numeral 3 is a leeward slope of a leeward area, the region corresponding to reference numeral 4 is an intersection area of the windward area and the leeward area, the region corresponding to reference numeral 5 is an intersection area of the leeward area and the flat area, and the region corresponding to reference numeral 6 is an intersection area of the windward area and the flat area.

[0058] The first medium structure layer can be divided into multiple sub-structure layers in the longitudinal direction, and the Q value variation law of each sub-structure layer is different. For example, the first medium structure layer includes a first sub-structure layer and a second sub-structure layer in the longitudinal direction, and the first sub-structure layer is located above the second sub-structure layer. The first sub-structure layer is an aeolian dune, which has flowability in the horizontal direction along the regional monsoon direction, and the variation characteristics of the compaction degree are primary and the variation of the water content is secondary in the horizontal direction, the structure is relatively loose, and the absorption and attenuation are the most serious. The second sub-structure layer is a stable dune, which refers to the structure from the bottom of the aeolian dune to the phreatic surface. The stable dune mainly includes the variation of the water content and the variation of the compaction degree in the longitudinal direction, and mainly includes the variation of the thickness of the dune in the horizontal direction, and the absorption and attenuation are smaller than those of the aeolian dune.

[0059] Therefore, in the embodiments of the present disclosure, the first medium structure layer has a plurality of different characteristic functions. Different regions in the lateral direction and different substructure layers in the longitudinal direction of the first medium structure layer correspond to different characteristic functions.

[0060] In some embodiments, a plurality of measuring points are arranged on the desert surface layer. The measuring point refers to a selected site on the desert surface layer. The desert surface layer of different region types includes a plurality of measuring points. Each measuring point is spaced apart by a long distance. For example, 50 measuring points are arranged, and each measuring point is spaced apart by 1 km or 2 km.

[0061] Step 201 includes determining the dune parameters of the first medium structure layer where the plurality of measuring points are located and the Q value fitting function of the plurality of measuring points, respectively; and determining a plurality of different characteristic functions of the first medium structure layer according to the dune parameters of the first medium structure layer where the plurality of measuring points are located and the Q value fitting function of the plurality of measuring points.

[0062] The dune parameters of the plurality of measuring points include the dune thickness, slope, and elevation of the measuring point. The dune thickness refers to the thickness of the first medium structure layer at the location of the measuring point; the slope refers to the inclination of the dune at the location of the measuring point; and the elevation refers to the distance between the top layer of the first medium structure layer and the bottom layer of the second medium structure layer at the location of the measuring point. In the embodiments of the present disclosure, the factors affecting the Q value of the first medium structure layer include the dune thickness, slope, and elevation. Therefore, a characteristic function with the dune thickness, slope, and elevation as independent variables and the Q value as the dependent variable can be established to represent the variation of the Q value in the first medium structure layer. In the embodiments of the present disclosure, the dune reference parameter is used as an intermediate variable to represent the three dune parameters of the dune thickness, slope, and elevation, and a characteristic function with the dune reference parameter as the independent variable and the Q value as the dependent variable is established.

[0063] The dune reference parameter at the location of each measuring point can be calculated by formula (1), which is as follows:

[0064] S=S H +(H G -1000)*3*P (1)

[0065] In formula (1), S represents the dune reference parameter, S H represents the dune thickness, H G represents the elevation, and P represents the slope. The dune thickness S H and the coordinates of each measuring point can be measured by conventional surface investigation techniques, such as single-well micro-logging technology and double-well micro-logging technology. The elevation of the grid where the measuring point is located can be obtained from the planar grid elevation map of the desert area according to the coordinates of each measuring point. The elevation is the H G of the measuring point. The slope P is calculated by formula (2).

[0066] P = (1 - SIN(2 * 0 * 3.1415927 / 180)) 2 (2)

[0067] In formula (2), P is the slope of the measuring point, and 0 is the inclination of the measuring point. The inclination of each measuring point can be calculated according to formula (3) as follows:

[0068] 0 = 180 / π * atan[(Dzx 2 + Dzy 2 )] 1 / 2 (3)

[0069] In formula (3), 0 is the inclination, Dzx = [(H 22 + 2 * H 12 + H 02 ) - (H 20 + 2 * H 10 + H 00 )] / (8 * ALx), Dzy = [(H 00 + 2 * H 01 + H 02 ) - (H 20 + 2 * H 21 + H 22 )] / (8 * ALy), H 00 is the elevation of the grid where the measuring point is located, H 01 and H 02 are the elevations of the first grid and the second grid in the x direction of the grid where the measuring point is located, H 10 and H 20 are the elevations of the first grid and the second grid in the y direction of the grid where the measuring point is located. H 11 and H 12 are the elevations of the first grid and the second grid in the y direction of the grid corresponding to H 01 , H 21 and H 22 are the elevations of the first grid and the second grid in the y direction of the grid corresponding to H 02 , ALx and ALy are the grid lengths in the x direction and the y direction respectively, and atan is the inverse tangent operation.

[0070] The Q-value fitting function of the measuring point can be obtained in the following manner: a plurality of data pairs of each measuring point can be obtained by using conventional surface investigation techniques, such as single-well micrologging technology, double-well micrologging technology, and the like, each data pair including a seismic wave propagation time and a corresponding depth; a corresponding Q-value can be calculated according to each data pair; and the Q-value fitting function of the measuring point can be fitted according to the depth and the corresponding Q-value in each data pair. The Q-value corresponding to each data pair can be calculated by using a spectral ratio method. The Q-value fitting functions of the plurality of measuring points can all be expressed as: Q(H)=A·H α wherein H is the depth of the desert surface layer, A represents the coefficient of the fitting function, and a represents the index of the fitting function. The values of A and / or a corresponding to the measuring points in different regions are different.

[0071] Exemplarily, according to the dune parameters of the first medium structure layer in which the plurality of measuring points are located and the Q-value fitting functions of the plurality of measuring points, a plurality of different characteristic functions of the first medium structure layer are determined, including: determining the dune parameters and the Q-value fitting functions corresponding to the measuring points in the target region; and linearly regressing the characteristic function of the target region according to the determined dune parameters and the coefficients and indexes of the Q-value fitting functions.

[0072] The target region is Figure 3 any one of the region types shown in the figure. The target region selects a plurality of measuring points, and exemplarily, the number of measuring points selected by the target region is 14, and each measuring point has a corresponding dune parameter and Q-value fitting function. In the embodiment of the present disclosure, the determination process of the characteristic function of the windward region of the desert surface layer is illustrated by taking the windward region as an example. Figure 4 is a fitting result figure of the characteristic function of the windward region provided by the present disclosure. Figure 4 is obtained by correspondingly displaying the coefficients and indexes of the Q-value fitting functions of the 14 measuring points in the windward region of the desert surface layer on a figure corresponding to the dune parameters of the 14 measuring points.

[0073] Figure 4 In the figure, the abscissa can be represented by the elevation, the dune thickness, and the dune reference parameter, and the ordinate can be represented by the coefficient and the index. The small square points represent the coefficients of the 14 measuring points and the elevations, the large square points represent the coefficients of the 14 measuring points and the dune thicknesses, the triangular points represent the indexes of the 14 measuring points and the elevations, the cross-shaped points represent the indexes of the 14 measuring points and the dune thicknesses, the circular points represent the indexes of the 14 measuring points and the dune reference parameters, and the cross-shaped points represent the coefficients of the 14 measuring points and the dune reference parameters. From Figure 4It can be seen that the coefficients of the 14 measuring points and the height, the coefficients and the dune thickness, the exponents of the 14 measuring points and the height, and the exponents and the dune thickness have no obvious mathematical law. However, the coefficients of the 14 measuring points and the dune reference parameters, and the exponents and the dune reference parameters both present an approximate linear relationship. The linear regression of the Q value fitting function of multiple measuring points and the corresponding dune reference parameters of the measuring points is performed, and the relationship between the coefficients and the dune reference parameters is y = -0.00025x + 1.02177, and the linear relationship between the exponents and the dune reference parameters is y = 0.00013x + 0.69177. According to the distance of each measuring point to the relationship between the coefficients and the dune reference parameters, the coefficient deviation of each measuring point can be obtained. The average value of the coefficient deviation of each measuring point can be obtained as the average coefficient deviation. According to the distance of each measuring point to the relationship between the exponents and the dune reference parameters, the exponent deviation of each measuring point can be obtained. The average value of the exponent deviation of each measuring point can be obtained as the average exponent deviation of each measuring point. In the embodiment of the present disclosure, the coefficient deviation and the exponent deviation of each measuring point are relatively small, the coefficient deviation is basically controlled within 10%, the average coefficient deviation is only 5%, the exponent deviation is basically controlled within 7%, and the average exponent deviation is only 3%. It can be known that the relationship between the coefficients and the dune reference parameters of the target region obtained by fitting can more accurately represent the change law of the coefficients and the dune reference parameters of the target region; and the relationship between the exponents and the dune reference parameters of the target region obtained by fitting can more accurately represent the change law of the exponents and the dune reference parameters of the target region.

[0074] The plurality of characteristic functions of the first dielectric structure layer are obtained by the same method as described above. In some embodiments, the first dielectric structure layer includes a first characteristic function, a second characteristic function, a third characteristic function, a fourth characteristic function, a fifth characteristic function, a sixth characteristic function, and a seventh characteristic function.

[0075] The first characteristic function is used to represent the change law of the Q value of the first sub-structure layer of the flat area. Since the second sub-structure layer of the flat area is thinner than the first sub-structure layer in the longitudinal direction, it can be ignored, and thus the first characteristic function can represent the change law of the Q value of the flat area of the first dielectric structure layer. Exemplarily, the first characteristic function is:

[0076] Q(H) = A1 · H α1 (4)

[0077] In formula (4), A1 = 0.0014S + 0.4849, and a1 = 0.0021S + 0.4141. S represents the dune reference parameter, and the dune reference parameter S is calculated by formula (1).

[0078] The second characteristic function is used to represent the variation law of the Q value of the first sub-structure layer of the leeward slope in the leeward area; and the third characteristic function is used to represent the variation law of the Q value of the second sub-structure layer of the leeward slope in the leeward area.

[0079] Exemplarily, the second characteristic function is:

[0080] Q(H)=A2·H α2 (5)

[0081] In the formula (5), A2=-0.0017S+1.8517, and α2=-0.00003S+0.9431. S represents a dune reference parameter, and the dune reference parameter S is calculated by the formula (1).

[0082] The third characteristic function is:

[0083] Q(H)=A3·H α3 (6)

[0084] In the formula (6), A3=0.0008S+1.7647, and α3=0.0002S+0.7284. S represents a dune reference parameter, and the dune reference parameter S is calculated by the formula (1).

[0085] The fourth characteristic function is used to represent the variation law of the Q value of the first sub-structure layer of the windward slope in the windward area; and the fifth characteristic function is used to represent the variation law of the Q value of the second sub-structure layer of the windward slope in the windward area.

[0086] Exemplarily, the fourth characteristic function is:

[0087] Q(H)=A4·H α4 (7)

[0088] In the formula (7), A4=-0.0003S+1.1162, and α4=0.0001S+0.9437. S represents a dune reference parameter, and the dune reference parameter S is calculated by the formula (1).

[0089] The fifth characteristic function is:

[0090] Q(H)=A5·H α5 (8)

[0091] In the formula (8), A5=0.0039S+0.0479, and α5=-0.0006S+1.0386. S represents a dune reference parameter, and the dune reference parameter S is calculated by the formula (1).

[0092] The sixth characteristic function is used to represent the variation law of the Q value of the first sub-structure layer and the second sub-structure layer of the intersection area of the leeward area and the windward area. Exemplarily, the sixth characteristic function is:

[0093] Q = [Q B *|Ψ-Ψ B0 |+Q Y *|Ψ-Ψ Y0 |] / |Ψ Y0 -Ψ B0 | (9)

[0094] In the interpolation formula (9), Ψ represents the included angle between the main wind direction and the inclination angle, Ψ B0 represents the first included angle threshold, Ψ Y0 represents the second included angle threshold, Ψ B0 and Ψ Y0 are empirical values. When Q B is the Q value calculated by using the second characteristic function, Q Y represents the Q value calculated by using the fourth characteristic function, Q represents the Q value of the first sub-structure layer of the intersection region of the leeside area and the windward area. When Q B is the Q value calculated by using the third characteristic function, Q Y represents the Q value calculated by using the fifth characteristic function, Q represents the Q value of the second sub-structure layer of the intersection region of the leeside area and the windward area.

[0095] The seventh characteristic function is used to represent the Q value variation law of the first sub-structure layer and the second sub-structure layer of the intersection region of the flat area and the leeside area or the intersection region of the flat area and the windward area. Exemplarily, the seventh characteristic function is:

[0096] Q = [Q H *|P S -P S0 |+Q P *|P S -P S1 |] / |P S0 -P S1 | (10)

[0097] In the interpolation formula (10), P S represents the ratio of the dune thickness to the slope, P S0 represents the first ratio threshold, P S1 represents the second ratio threshold, P S0 and P S1 are empirical values. Q H represents the Q value calculated by using the sixth characteristic function, Q P represents the Q value calculated by using the first characteristic function, Q represents the Q value of the first sub-structure layer of the intersection region of the flat area and the leeside area or the intersection region of the flat area and the windward area.

[0098] Since the Q value varies in different regions in the lateral direction of the first medium structure layer and in different substructure layers in the longitudinal direction, the characteristic functions of different substructure layers corresponding to different regions in the lateral direction of the first medium structure layer are respectively constructed in the embodiment of the present disclosure, so that the calculation result of the single-point Q value is more accurate, and the establishment accuracy of the absorption attenuation field is further improved.

[0099] In step 202, the characteristic function of the second medium structure layer is determined.

[0100] The characteristic function of the second medium structure layer is used to represent the variation rule of the Q value of the second medium structure layer. In the embodiment of the present disclosure, the characteristic function of the second medium structure layer includes: determining the Q value of the second medium structure layer where the plurality of measuring points are located; and fitting the characteristic function of the second medium structure layer according to the Q value of the second medium structure layer where the plurality of measuring points are located and the dune reference parameters of the plurality of measuring points. The Q value of the second medium structure layer corresponding to the plurality of measuring points can be obtained by a conventional surface investigation method.

[0101] Figure 5 is a fitting result graph of a characteristic function of a second medium structure layer provided by the embodiment of the present disclosure. As shown in the left graph, Figure 5 Overall, the Q value in the second medium structure layer corresponding to the plurality of measuring points increases with the increase of the propagation velocity of the seismic wave, but the data points are relatively scattered, the error of the fitted characteristic function is large, and therefore the accuracy of the calculated Q value is low. As shown in the right graph, Figure 5 The Q value of the plurality of measuring points has a good regular characteristic with the corresponding dune reference parameters. By fitting the Q value of the plurality of measuring points of the second medium structure layer and the corresponding dune reference parameters, the eighth characteristic function is obtained as:

[0102] Q=-0.0004S 2 +0.4484S-74.86 (11)

[0103] In formula (11), S represents the dune reference parameter, and the dune reference parameter S is calculated by formula (1).

[0104] In step 203, the first absorption attenuation parameter value of the target physical point is calculated according to the characteristic function of the first medium structure layer. Optionally, before step 203 is performed, a smoothing algorithm is further included to remove small layer dunes in the desert area surface layer.

[0105] The inclination of the dune has different levels, and the big level dune has a small dune structure. In the formation process of the flowing dune, the shape of the big level dune is stable and long-term, and the shape of the small level dune is unstable and short-term. Therefore, we should grasp the inclination characteristics and slope characteristics of the big dune. Through the smoothing algorithm, the elevation data of all physical points in the surface layer of the desert area are grid smoothed, and the small level dune is removed. The size of the grid depends on the average dune height of the target desert area surface layer. In this way, the small level dune can be removed, and only the big level dune is retained, which is convenient for identifying the region type where the target physical point is located and the characteristic function corresponding to the region type according to the dune parameters of the target physical point.

[0106] In the embodiment of the present disclosure, the surface layer of the desert area is provided with a plurality of physical points. The related content of the physical point is described in the foregoing step 102, and the detailed description is omitted here. The target physical point is one of the plurality of physical points.

[0107] In the embodiment of the present disclosure, the step 203 comprises:

[0108] Firstly, the dune parameters of the position where the target physical point is located are obtained.

[0109] The dune parameters include dune thickness, slope and elevation. The related content of the dune parameters is described in the foregoing step 201, and the detailed description is omitted here.

[0110] Secondly, the target characteristic function of the target physical point is determined according to the dune parameters.

[0111] The target characteristic function is at least one of the plurality of characteristic functions of the first medium structure layer. The plurality of characteristic functions of the first medium structure layer are described in the foregoing step 201, and the detailed description is omitted here.

[0112] In the embodiment of the present disclosure, the target characteristic function is determined according to at least one of the ratio of the dune thickness to the slope and the angle between the main wind direction and the inclination. In some examples, the target characteristic function can be determined according to the ratio of the dune thickness to the slope; in other examples, the target characteristic function can be determined according to the angle between the main wind direction and the inclination; in still other examples, the target characteristic function can be determined according to the ratio of the dune thickness to the slope and the angle between the main wind direction and the inclination. The main wind direction refers to the dominant wind direction of the desert area, which can be determined from the surface elevation map or the surface thickness map of the desert area in the past. The inclination refers to the inclination direction of the dune where the target physical point is located. The inclination of the target physical point is calculated by formula (12), and formula (12) is as follows:

[0113] Ψ = 180 / π * atan (-Dzx / Dzy) (12)

[0114] In formula (13), Ψ is the tendency, Dzx = [(H 22 + 2*H 12 + H 02 ) - (H 20 + 2*H 10 + H 00 )] / (8*ΔLx), Dzy = [(H 00 + 2*H 01 + H 02 ) - (H 20 + 2*H 21 + H 22 )] / (8*ΔLy), H 00 is the height of the grid where the measuring point is located, H 01 and H 02 are the heights of the first grid and the second grid in the x direction of the grid where the measuring point is located, H 10 and H 20 are the heights of the first grid and the second grid in the y direction of the grid where the measuring point is located. H 11 and H 12 are the heights of the first grid and the second grid in the y direction of the corresponding grid of H 01 . H 21 and H 22 are the heights of the first grid and the second grid in the y direction of the corresponding grid of H 02 , ΔLx and ΔLy are the grid lengths in the x direction and the y direction respectively, and atan is the inverse tangent operation.

[0115] In some embodiments, the target characteristic function is determined according to at least one of the ratio of the dune thickness to the slope and the angle between the main wind direction and the tendency, including at least one of:

[0116] (1) in response to the ratio of the dune thickness to the slope being less than a first ratio threshold, determining the target characteristic function as a first characteristic function.

[0117] The first ratio threshold can be represented as P S0 , when the ratio of the dune thickness to the slope is less than the first ratio threshold, it indicates that the target physical point is located in a flat area, at this time, the target characteristic function corresponding to the target physical point can be determined as the first characteristic function. The first characteristic function is used to calculate the first absorption attenuation parameter value of the first sub-structure layer.

[0118] (2) in response to the ratio of the dune thickness to the slope being greater than a second ratio threshold, and the angle between the main wind direction and the tendency being less than a first angle threshold, determining the target characteristic function as a second characteristic function and a third characteristic function.

[0119] When the ratio of the dune thickness to the slope is greater than the second ratio threshold value and the angle between the main wind direction and the tendency is less than the first angle threshold value, it is indicated that the target physical point is located on the leeward slope of the leeward area. At this time, the target characteristic function corresponding to the target physical point can be determined as the second characteristic function and the third characteristic function. The second characteristic function and the third characteristic function are respectively used to calculate the first absorption attenuation parameter value of the first sub-structure layer and the first absorption attenuation parameter value of the second sub-structure layer. The first angle threshold value is less than the second angle threshold value.

[0120] (3) In response to the ratio of the dune thickness to the slope being greater than the second ratio threshold value and the angle between the main wind direction and the tendency being greater than the second angle threshold value, the target characteristic function is determined as the fourth characteristic function and the fifth characteristic function. The first angle threshold value is less than the second angle threshold value.

[0121] When the ratio of the dune thickness to the slope is greater than the second ratio threshold value and the angle between the main wind direction and the tendency is greater than the second angle threshold value, it is indicated that the target physical point is located on the windward slope of the windward area. At this time, the target characteristic function corresponding to the target physical point can be determined as the fourth characteristic function and the fifth characteristic function. The fourth characteristic function and the fifth characteristic function are respectively used to calculate the first absorption attenuation parameter value of the first sub-structure layer and the first absorption attenuation parameter value of the second sub-structure layer.

[0122] (4) In response to the ratio of the dune thickness to the slope being greater than the second ratio threshold value and the angle between the main wind direction and the tendency being greater than the first angle threshold value and less than the second angle threshold value, the target characteristic function is determined as the sixth characteristic function.

[0123] When the ratio of the dune thickness to the slope is greater than the second ratio threshold value and the angle between the main wind direction and the tendency is greater than the first angle threshold value and less than the second angle threshold value, it is indicated that the target physical point is located in the region where the leeward area and the windward area intersect. At this time, the target characteristic function corresponding to the target physical point can be determined as the sixth characteristic function. The sixth characteristic function is used to calculate the first absorption attenuation parameter value of the first sub-structure layer and the first absorption attenuation parameter value of the second sub-structure layer.

[0124] (5) In response to the ratio of the dune thickness to the slope being greater than the first ratio threshold value and less than the second ratio threshold value, the target characteristic function is determined as the seventh characteristic function.

[0125] When the ratio of the dune thickness to the slope is greater than the first ratio threshold value and less than the second ratio threshold value, it is indicated that the target physical point is located in the region where the flat area and the leeward area intersect, or in the region where the flat area and the windward area intersect. At this time, the target characteristic function can be determined as the seventh characteristic function. The seventh characteristic function is used to calculate the first absorption attenuation parameter value of the first sub-structure layer.

[0126] Thirdly, the first absorption attenuation parameter value of the target physical point is calculated according to the target characteristic function.

[0127] The third step comprises calculating a dune reference parameter of the target physical point according to the dune parameter of the location where the target physical point is located; and calculating a first absorption attenuation parameter value of the target physical point according to the dune reference parameter of the target physical point and the target characteristic function of the target physical point determined in the second step. The dune reference parameter of the target physical point is calculated by using formula (1). In the embodiments of the present disclosure, for each target physical point, the first absorption attenuation parameter value is calculated according to the first step to the third step.

[0128] In step 204, a second absorption attenuation parameter value of the target physical point is calculated according to the characteristic function of the second medium structure layer.

[0129] The second absorption attenuation parameter value of the target physical point is calculated according to the eighth characteristic function.

[0130] In the embodiments of the present disclosure, the target physical point corresponds to at least one first absorption attenuation parameter value in the first medium structure layer and corresponds to one second absorption attenuation parameter value in the second medium structure layer. The dune reference parameter used for calculating the first absorption attenuation parameter value and the second absorption attenuation parameter value of the target physical point is the same.

[0131] In step 205, the first absorption attenuation parameter value of the target physical point is corrected according to the first absorption attenuation parameter difference values of the plurality of measuring points, to obtain a corrected first absorption attenuation parameter value of the target physical point.

[0132] For the related content of the measuring point arrangement, refer to the aforementioned step 201, and the detailed description is omitted here. The accurate third absorption attenuation parameter value of each measuring point can be measured by using the conventional surface investigation technology. The fourth absorption attenuation parameter value of each measuring point can also be calculated by using the method of the aforementioned step 203. The first absorption attenuation parameter difference value is the difference between the third absorption attenuation parameter value of the measuring point measured and the fourth absorption attenuation parameter value of the measuring point calculated.

[0133] In some embodiments, correcting the first absorption attenuation parameter value of the target physical point comprises: calculating the average value of the first absorption attenuation parameter difference values of the plurality of measuring points; and correcting the first absorption attenuation parameter value of the target physical point according to the average value of the first absorption attenuation parameter difference values of the plurality of measuring points.

[0134] Exemplarily, the sum of the first absorption attenuation parameter value of the target physical point and the average value of the third absorption attenuation parameter difference values of the plurality of measuring points is taken as the corrected first absorption attenuation parameter value of the target physical point.

[0135] In step 206, the second absorption attenuation parameter value of the target physical point is corrected according to the second absorption attenuation parameter difference values of the plurality of measuring points, to obtain a corrected second absorption attenuation parameter value of the target physical point.

[0136] Each measurement point can be measured by a conventional surface investigation technique to obtain an accurate fifth absorption attenuation parameter value. Each measurement point can also be calculated by the method of step 204 to obtain a sixth absorption attenuation parameter value. The second absorption attenuation parameter difference value is the difference between the fifth absorption attenuation parameter value of the measurement point and the sixth absorption attenuation parameter value of the measurement point.

[0137] In some embodiments, the correction of the second absorption attenuation parameter value of the target physical point comprises: calculating an average value of the second absorption attenuation parameter difference values of the plurality of measurement points, and correcting the second absorption attenuation parameter value of the target physical point according to the average value of the second absorption attenuation parameter difference values of the plurality of measurement points.

[0138] Exemplarily, the sum of the second absorption attenuation parameter value of the target physical point and the average value of the second absorption attenuation parameter difference values of the plurality of measurement points is taken as the corrected second absorption attenuation parameter value of the target physical point.

[0139] Since the characteristic functions of the first medium structure layer and the characteristic functions of the second medium structure layer are both fitted by the dune parameters and the like of the measurement points, there will be certain errors in calculating the first absorption attenuation parameter value and the second absorption attenuation parameter value of the target physical point. By correcting the first absorption attenuation parameter value and the second absorption attenuation parameter value of the target physical point, the accuracy of the calculated first absorption attenuation parameter value and the second absorption attenuation parameter value of the target physical point relative to the total absorption attenuation amount can be improved.

[0140] In step 207, an absorption attenuation field of the desert area surface layer is established according to the corrected first absorption attenuation parameter value and the corrected second absorption attenuation parameter value of the target physical point.

[0141] In some embodiments, the step 207 comprises: a first step of calculating a first relative absorption attenuation total amount of the target physical point according to the corrected first absorption attenuation parameter value and the corrected second absorption attenuation parameter value of the target physical point.

[0142] The first relative absorption attenuation total amount of the target physical point refers to the total absorption attenuation amount from the ground surface to the reference surface corresponding to the target physical point. The desert area surface layer is above the reference surface, and exemplarily, the reference surface is 800 meters away from the ground surface.

[0143] The first relative absorption attenuation total amount of the target physical point is calculated by the following formula:

[0144]

[0145] Q represents the first relative total absorption attenuation of the target physical point in formula (13). i represents the number of layers divided in the longitudinal direction of the desert area surface layer, for example, the first medium structure layer includes two substructure layers, and the second medium structure layer is one layer, so i is 3. Q i represents the Q value of the i-th layer corresponding to the target physical point, for example, at least one first Q value in the first medium structure layer and a second Q value in the second medium structure layer corresponding to the target physical point, t i represents the propagation time of the seismic wave in the i-th layer, t i It can be obtained by a conventional surface investigation method.

[0146] Secondly, the second relative total absorption attenuation of the plurality of measuring points is calculated according to the measured third absorption attenuation parameter value and the fifth absorption attenuation parameter value of the plurality of measuring points.

[0147] The formula for calculating the second relative total absorption attenuation of the plurality of measuring points is the same as formula (13).

[0148] Thirdly, the third relative total absorption attenuation of the plurality of measuring points is calculated according to the calculated fourth absorption attenuation parameter value and the sixth absorption attenuation parameter value of the plurality of measuring points.

[0149] The formula for calculating the third relative total absorption attenuation of the plurality of measuring points is the same as formula (13).

[0150] The second relative total absorption attenuation and the third relative total absorption attenuation of the measuring point both refer to the total absorption attenuation from the surface to the reference surface corresponding to the measuring point.

[0151] Fourthly, the relative total absorption attenuation difference of the plurality of measuring points is obtained according to the second relative total absorption attenuation and the third relative total absorption attenuation.

[0152] The relative total absorption attenuation difference of the plurality of measuring points is obtained by subtracting the second relative total absorption attenuation from the third relative total absorption attenuation.

[0153] Fifthly, the first relative total absorption attenuation of the target physical point is corrected according to the relative total absorption attenuation difference of the plurality of measuring points to obtain the corrected first relative total absorption attenuation of the target physical point. In some embodiments, the fifth step includes: calculating the error correction value corresponding to the target physical point according to the relative total absorption attenuation difference of the plurality of measuring points; and correcting the first relative total absorption attenuation of the target physical point according to the error correction value corresponding to the target physical point to obtain the second relative total absorption attenuation of the target physical point.

[0154] The error correction value corresponding to the target physical point is calculated by the following formula:

[0155]

[0156] In formula (14), y represents the error correction value, M represents the total number of measuring points, and Z represents the error correction value. i d represents the total difference in relative absorption attenuation at the i-th measurement point. i This represents the horizontal distance between the target physical point and the i-th measuring point.

[0157] The first relative absorption attenuation at the target physical point is corrected using the following formula:

[0158] Q = Q' + y (15)

[0159] In formula (15), Q represents the total second relative absorption attenuation of the target physical point, Q' represents the total first relative absorption attenuation of the target physical point, and y represents the error correction value of the target physical point. Since the calculated total first relative absorption attenuation of the target physical point may have errors, the error correction value of the target physical point can be obtained based on the difference in the total relative absorption attenuation of multiple measurement points and the inverse distance interpolation formula (14). This error correction value can correct the total first relative absorption attenuation of the target physical point, making the corrected total second relative absorption attenuation of the target physical point more accurate, thereby improving the accuracy of the field establishment.

[0160] The sixth step is to establish the absorption attenuation field of the desert surface based on the total relative absorption attenuation after correction of the target physical point.

[0161] In this embodiment of the disclosure, the total amount of second relative absorption attenuation of multiple physical points is input into the field construction software to obtain the absorption attenuation field of the desert surface.

[0162] To verify the accuracy of the absorption attenuation field established using the field establishment method provided in this disclosure, 42 measuring points were arranged on the desert surface. Only 13 of these points were used for calibration, while the remaining 29 points were used as a posteriori points to verify the accuracy of the field establishment. A comparison of the field establishment accuracy using the method of this disclosure and the method using related technologies is shown in Table 1.

[0163] Table 1:

[0164]

[0165]

[0166] In Table 1, the location of each measuring point corresponds to its name, used to distinguish different measuring points. The actual value represents the accurate total relative absorption attenuation obtained for each measuring point using conventional surface survey techniques. The old method referred to the total relative absorption attenuation based on the coordinates of each measuring point from... Figure 7 The total relative absorption attenuation is read from the coordinates of each measurement point. The new method refers to the total relative absorption attenuation based on the coordinates of each measurement point. Figure 6The relative absorption attenuation total amount read.

[0167] From the data in Table 1, it can be seen that in the absorption attenuation field of the present application, the average deviation of the relative absorption attenuation total amount of 42 measuring points is 3.5%, and the maximum deviation is 14.5%. In the absorption attenuation field of the related art, the average deviation of the relative absorption attenuation total amount of 42 measuring points is 17.0%, and the maximum deviation is 45.9%, which is 13.5% and 31.4% higher than that of the present application, respectively. It can be seen that the precision of the absorption attenuation field established by the field building method of the present application is higher than that of the absorption attenuation field established by the field building method of the related art.

[0168] If the 13 measuring points used for calibration do not participate in the average deviation statistics, the average deviation of the relative absorption attenuation total amount of 29 measuring points of the present application is 5%, and the average deviation of the relative absorption attenuation total amount of 29 measuring points of the related art becomes 25%. It shows that the fluctuation of the absorption attenuation field established by the method of the related art is more intense, which is more than 20% larger than the fluctuation of the absorption attenuation field established by the method of the present application.

[0169] Alternatively, in the embodiments of the present disclosure, the effect diagram of the absorption attenuation field of the desert surface layer can also be used to verify the precision of the absorption attenuation field established by the field building method provided by the present disclosure. In the effect diagram of the absorption attenuation field, as long as the coordinates of a certain point are known, the relative absorption attenuation total amount value corresponding to the point can be obtained from the diagram.

[0170] Because the sand dune region in the desert surface has higher absorption and attenuation of seismic waves, the corresponding relative absorption and attenuation value is large; while the flat region has relatively weaker absorption and attenuation of seismic waves, the corresponding relative absorption and attenuation value is small. Different colors are used to represent different relative absorption and attenuation values ​​on the desert surface. The higher the accuracy of the calculated relative absorption and attenuation values ​​at different locations on the desert surface, the more clearly the distribution of sand dunes and flat regions can be seen in the absorption and attenuation field diagram. Conversely, if the accuracy of the calculated relative absorption and attenuation values ​​at different locations on the desert surface is low, it is difficult to discern the approximate distribution of sand dunes and flat regions from the absorption and attenuation field diagram. The data in Table 1 above shows that the method of this embodiment obtains high accuracy and small error in the relative absorption and attenuation values ​​at different locations on the desert surface, while the method using related technologies obtains low accuracy and large error. Therefore, from the effect diagram of the absorption attenuation field obtained by the field construction method provided in this embodiment, it can be seen that the distribution of sand dunes and flat areas on the surface of the desert area is consistent with the distribution of sand dunes and flat areas on the elevation map or thickness map of the desert area, and wind gap traces consistent with those on the elevation map and thickness map can also be seen. From the effect diagram of the absorption attenuation field obtained by the field construction method provided in related technologies, it can be seen that the distribution of sand dunes and flat areas on the surface of the desert area is inconsistent with the distribution of sand dunes and flat areas on the elevation map or thickness map of the desert area, and wind gap traces are difficult to see. Therefore, the effect diagram of the absorption attenuation field on the surface of the desert area can also verify the accuracy of the absorption attenuation field established by the field construction method provided in this embodiment.

[0171] Figure 6 This is a diagram illustrating the effect of compensating for the absorption attenuation field established using the embodiments of this disclosure in the ancient city test area. (See attached diagram.) Figure 6 In the diagram, the leftmost curve represents the Earth's surface elevation line, the horizontal axis represents travel time, and the vertical axis represents the distance between the shot point and the receiver. Each horizontal line in the diagram represents the waveform record of the seismic waves received by each receiver. Each line consists of multiple different waveforms, and the first waveform appearing on each horizontal line is the first arrival wave. The vertical width of the first arrival wave represents its amplitude, and the horizontal width represents its dominant frequency. Under normal circumstances, as the shot-receiver distance increases, the amplitude and dominant frequency of the first arrival wave decrease systematically.

[0172] from Figure 6It can be seen that before compensation, the first arrival wave of the adjacent geophone points is arranged closely when the offset is between 200m and 300m, and the first arrival wave of the adjacent geophone points changes from close arrangement to sparse arrangement when the offset is between 300m and 400m. It is shown that, between 200m and 400m, the amplitude of the first arrival wave decreases with the increase of the offset. However, in the part marked in the figure, the first arrival wave of each geophone point changes from sparse arrangement to close arrangement when the offset is between 400m and 700m, which shows that, before compensation, the amplitude of the first arrival wave changes from small to large with the increase of the offset, and does not change according to the rule from large to small. Moreover, in the part marked in the figure, the lateral width of the first arrival wave between 400m and 500m is smaller than that between 700m and 900m. It is shown that, before compensation, the main frequency of the first arrival wave does not change according to the rule from large to small with the increase of the offset.

[0173] After the absorption and attenuation field established by the embodiment of the present disclosure is used for compensation, the first arrival wave of each geophone point changes from close arrangement to sparse arrangement when the offset is between 200m and 400m, and the first arrival wave changes from sparse arrangement to close arrangement when the offset is between 400m and 700m in the area marked in the figure, but the close degree is weaker than that before compensation. It is shown that, after compensation, the amplitude of the first arrival wave changes from small to large with the increase of the offset, but the degree of amplitude increase is weaker than that before compensation. Moreover, the lateral width of the first arrival wave changes from large to small between 200m and 1800m in the figure.

[0174] It can be seen that, before compensation, the first arrival wave is still affected by the surface layer of the desert area (the time difference effect has been eliminated), and the reliability is low if the high-frequency component is directly obtained from the record. After the absorption and attenuation field established by the embodiment of the present disclosure is used for compensation, the change characteristics of the first arrival wave are reasonable, and the influence of the surface layer of the desert area is small. Overall, the reliability of the single shot record after compensation is higher than that before compensation. Therefore, the absorption and attenuation field established by the embodiment of the present disclosure can reduce the influence of the surface layer of the desert area on the absorption and attenuation of the seismic wave to some extent, and the result is more reliable.

[0175] It should be noted that, in the embodiment of the present disclosure, the first medium structure layer is divided into two sub-structure layers for description, and in other embodiments, in order to further improve the calculation accuracy of the Q value of the first medium structure layer, the first medium structure layer can also be divided into more than two sub-structure layers, and the Q value of each sub-structure layer corresponding to the target physical point in the more than two sub-structure layers is calculated respectively.

[0176] In this embodiment, since the Q-value of the second medium structure layer varies little vertically, it is not divided into layers. However, in other embodiments, to further improve the calculation accuracy of the Q-value of the second medium structure layer, it can be divided into two or more sub-structure layers, and the Q-value of each sub-structure layer corresponding to the target physical point can be calculated separately. For example, based on the lithological characteristics of the second medium structure layer, it can be divided into two sub-structure layers: an unstable high-velocity layer and a stable high-velocity layer. The unstable high-velocity layer refers to the portion from below the water table to the top of the stable high-velocity layer, including water-bearing fine sand, weathered layers, or clay layers, etc., with very little lateral variation in the Q-value within the layer. The stable high-velocity layer mainly includes hard clay or diagenetic strata, with the Q-value within the layer remaining essentially unchanged laterally.

[0177] In this embodiment, the surface structure of the desert region is divided into a first dielectric structure layer and a second dielectric structure layer. A characteristic function is established for the first dielectric structure layer to calibrate the variation law of its absorption attenuation parameters, and a characteristic function is established for the second dielectric structure layer to calibrate the variation law of its absorption attenuation parameters. Due to the complex structure of the desert surface, dividing the desert surface into layers and establishing a corresponding characteristic function for each dielectric structure layer can improve the calculation accuracy of the Q-value of physical points in each dielectric structure layer, thereby improving the accuracy of the established absorption attenuation field.

[0178] Figure 7 This is a flowchart illustrating another method for establishing an absorption attenuation field on the surface of a desert region, provided in an embodiment of this disclosure. The method includes:

[0179] In step 301, a dune survey is conducted on the surface of the target desert area to determine the lateral regional types of the surface of the target desert area.

[0180] Sand dune surveys can identify the types of surface dunes in a target desert area. These dune types can include ridge dunes, honeycomb dunes, etc. Based on the dune types, the regional type of the surface layer in the target desert area can be determined. Figure 3 One or more of the area types shown. In some examples, the area type of the target desert area surface includes only flat areas. In other examples, the area type of the target desert area surface includes flat areas, windward areas, and leeward areas.

[0181] In this embodiment of the disclosure, step 301 further includes determining the prevailing wind direction on the surface of the target desert area. For example, the prevailing wind direction can be determined from previous surface elevation maps or surface thickness maps of the target desert area.

[0182] In step 302, the first and second media structure layers of the target desert area surface are determined in the longitudinal direction.

[0183] In some embodiments, the phreatic surface of the surface layer of the target desert area can be obtained by a conventional surface investigation method, such as single-well micro logging technology, double-well micro logging technology, etc. The part between the phreatic surface and the ground surface of the target desert area is divided into the first medium structure layer, and the part between the phreatic surface and the datum surface is divided into the second medium structure layer. The related content of the datum surface is described in the foregoing step 207, and the detailed description is omitted here.

[0184] Since the structure of the first medium structure layer is relatively complex, the first medium structure layer can be further divided into multiple sub-structure layers according to the structural characteristics of the first medium structure layer, for example, the first medium structure layer is divided into a first sub-structure layer and a second sub-structure layer. The related content of the first sub-structure layer and the second sub-structure layer is described in the foregoing step 201, and the detailed description is omitted here.

[0185] In step 303, a Q-value fitting function of the first medium structure layer in which the multiple measuring points are located is determined.

[0186] The related content of the sand dune parameters and the Q-value fitting function is described in the foregoing step 201, and the detailed description is omitted here.

[0187] In step 304, the Q-value of the second medium structure layer in which the multiple measuring points are located is determined.

[0188] The related content is described in the foregoing step 202, and the detailed description is omitted here.

[0189] In step 305, a characteristic function of the first medium structure layer is determined according to the sand dune parameters of the multiple measuring points and the Q-value fitting function of the first medium structure layer in which the multiple measuring points are located, and a characteristic function of the second medium structure layer is determined according to the sand dune parameters of the multiple measuring points and the Q-value of the second medium structure layer in which the multiple measuring points are located.

[0190] The characteristic function of the first medium structure layer has multiple, and the characteristic function of the second medium structure layer is one. The related content of determining the characteristic function of the first medium structure layer of the surface layer of the target desert area is described in the foregoing step 201, and the related content of determining the characteristic function of the second medium structure layer of the surface layer of the target desert area is described in the foregoing step 202, and the detailed description is omitted here.

[0191] In step 306, the first absorption attenuation parameter values of the multiple physical points are calculated according to the characteristic function of the first medium structure layer, and the second absorption attenuation parameter values of the multiple physical points are calculated according to the characteristic function of the second medium structure layer.

[0192] The first absorption attenuation parameter values of the plurality of physical points are calculated according to the characteristic function of the first medium structure layer, as described in the foregoing step 203; and the second absorption attenuation values of the plurality of physical points are calculated according to the characteristic function of the second medium structure layer, as described in the foregoing step 204, and the detailed description is omitted here.

[0193] In step 307, the absorption attenuation field of the desert surface layer is established according to the first absorption attenuation parameter values and the second absorption attenuation parameter values of the plurality of physical points.

[0194] The detailed description is omitted here, and refer to the foregoing steps 205 to 207.

[0195] In the embodiments of the present disclosure, the structure of the desert surface layer is divided into a first medium structure layer and a second medium structure layer, and a characteristic function of the first medium structure layer is established to calibrate the variation of the absorption attenuation parameter of the first medium structure layer, and a characteristic function of the second medium structure layer is established to calibrate the variation of the absorption attenuation parameter of the second medium structure layer. Since the structure of the desert surface layer is complex, the desert surface layer is divided into layers, and a corresponding characteristic function is established for each medium structure layer, which can improve the calculation accuracy of the Q value of the physical point in each medium structure layer, and further improve the accuracy of the established absorption attenuation field.

[0196] Figure 8 Fig. 8 is a structural block diagram of an apparatus 800 for establishing an absorption attenuation field of a desert surface layer provided by the embodiments of the present disclosure. As shown in the figure, the apparatus includes a determination module 801, a first calculation module 802, a second calculation module 803, and an establishment module 804. Figure 8

[0197] The determination module 801 is configured to determine a characteristic function of the first medium structure layer and a characteristic function of the second medium structure layer, the characteristic function of the first medium structure layer being used to represent the variation of the absorption attenuation parameter of the first medium structure layer, the characteristic function of the second medium structure layer being used to represent the variation of the absorption attenuation parameter of the second medium structure layer, and the characteristic function of the first medium structure layer being different from the characteristic function of the second medium structure layer. The first calculation module 802 is configured to calculate the first absorption attenuation parameter values of the plurality of physical points according to the characteristic function of the first medium structure layer. The second calculation module 803 is configured to calculate the second absorption attenuation parameter values of the plurality of physical points according to the characteristic function of the second medium structure layer. The establishment module 804 is configured to establish the absorption attenuation field of the desert surface layer according to the first absorption attenuation parameter values and the second absorption attenuation parameter values of the plurality of physical points.

[0198] ​Optionally, the first medium structure layer has a plurality of different characteristic functions. The first calculation module 802 is configured to obtain a dune parameter of a position of a target physical point, the target physical point being one of the plurality of physical points; determine a target characteristic function of the target physical point according to the dune parameter, the target characteristic function being at least one of the plurality of characteristic functions of the first medium structure layer; and calculate a first absorption attenuation parameter value of the target physical point according to the target characteristic function.

[0199] Optionally, the dune parameter includes a dune thickness, a slope and an elevation; and the first calculation module 802 is configured to determine the target characteristic function according to at least one of a ratio of the dune thickness to the slope and an included angle between a wind direction and a tendency, the tendency being calculated according to the elevation.

[0200] Optionally, the first medium structure layer includes a first sub-structure layer and a second sub-structure layer in the longitudinal direction, the first sub-structure layer is located above the second sub-structure layer, the first calculation module 802 is configured to determine the target characteristic function as a first characteristic function in response to that the ratio of the dune thickness to the slope is less than a first ratio threshold value, the first characteristic function is used to calculate a first absorption attenuation parameter value of the first sub-structure layer; determine the target characteristic function as a second characteristic function and a third characteristic function in response to that the ratio of the dune thickness to the slope is greater than a second ratio threshold value, and the included angle between the main wind direction and the tendency is less than a first included angle threshold value, the second characteristic function and the third characteristic function are respectively used to calculate a first absorption attenuation parameter value of the first sub-structure layer and a first absorption attenuation parameter value of the second sub-structure layer; the first ratio threshold value is less than the second ratio threshold value; determine the target characteristic function as a fourth characteristic function and a fifth characteristic function in response to that the ratio of the dune thickness to the slope is greater than the second ratio threshold value, and the included angle between the main wind direction and the tendency is greater than a second included angle threshold value, the fourth characteristic function and the fifth characteristic function are respectively used to calculate a first absorption attenuation parameter value of the first sub-structure layer and a first absorption attenuation parameter value of the second sub-structure layer; the first included angle threshold value is less than the second included angle threshold value; determine the target characteristic function as a sixth characteristic function in response to that the ratio of the dune thickness to the slope is greater than the second ratio threshold value, the included angle between the main wind direction and the tendency is greater than the first included angle threshold value and less than the second included angle threshold value, the sixth characteristic function is used to calculate a first absorption attenuation parameter value of the first sub-structure layer and a first absorption attenuation parameter value of the second sub-structure layer; determine the target characteristic function as a seventh characteristic function in response to that the ratio of the dune thickness to the slope is greater than the first ratio threshold value and less than the second ratio threshold value, the seventh characteristic function is used to calculate a first absorption attenuation parameter value of the first sub-structure layer.

[0201] Optionally, the second calculation module 803 is configured to calculate the second absorption attenuation parameter value of the target physical point according to an eighth characteristic function.

[0202] Here, each characteristic function is described in the foregoing method embodiments, and the detailed description is omitted here.

[0203] Optionally, the desert area surface layer is provided with a plurality of measuring points, and the establishing module 804 is configured to correct the first absorption attenuation parameter value of the target physical point according to a first absorption attenuation parameter difference value of the plurality of measuring points, to obtain a corrected first absorption attenuation parameter value of the target physical point, the first absorption attenuation parameter difference value being a difference value between a third absorption attenuation parameter value of the measuring point measured and a fourth absorption attenuation parameter value of the measuring point calculated; correct the second absorption attenuation parameter value of the target physical point according to a second absorption attenuation parameter difference value of the plurality of measuring points, to obtain a corrected second absorption attenuation parameter value of the target physical point; the second absorption attenuation parameter difference value being a difference value between a fifth absorption attenuation parameter value of the measuring point measured and a sixth absorption attenuation parameter value of the measuring point calculated; and establish the absorption attenuation field of the desert area surface layer according to the corrected first absorption attenuation parameter value and the corrected second absorption attenuation parameter value of the target physical point.

[0204] Optionally, the establishing module 804 is configured to calculate a first relative absorption attenuation total amount of the target physical point according to the corrected first absorption attenuation parameter value and the corrected second absorption attenuation parameter value of the target physical point; calculate a second relative absorption attenuation total amount of the plurality of measuring points according to the third absorption attenuation parameter value and the fifth absorption attenuation parameter value of the plurality of measuring points measured; calculate a third relative absorption attenuation total amount of the plurality of measuring points according to the fourth absorption attenuation parameter value and the sixth absorption attenuation parameter value of the plurality of measuring points calculated; obtain a relative absorption attenuation total amount difference value of the plurality of measuring points according to the second relative absorption attenuation total amount and the third relative absorption attenuation total amount; correct the first relative absorption attenuation total amount of the target physical point according to the relative absorption attenuation total amount difference value of the plurality of measuring points, to obtain a corrected first relative absorption attenuation total amount of the target physical point; and establish the absorption attenuation field of the desert area according to the corrected first relative absorption attenuation total amount of the target physical point.

[0205] It should be noted that the establishing device 800 of the absorption attenuation field of the desert area surface layer provided in the above embodiments is only used as an example to illustrate the division of the above functional modules in establishing the absorption attenuation field of the desert area surface layer. In actual applications, the above functions can be completed by different functional modules according to needs, that is, the internal structure of the device is divided into different functional modules to complete all or part of the functions described above. In addition, the establishing device of the absorption attenuation field of the desert area surface layer and the method for establishing the absorption attenuation field of the desert area surface layer provided in the above embodiments belong to the same concept, and the specific implementation process is described in the method embodiments, which will not be described here.

[0206] Figure 9 is a structural block diagram of a computer device provided by the embodiments of the present disclosure. As shown inFigure 9 As shown in FIG. 9, the computer device can be a computer or the like. The computer device includes a processor 901 and a memory 902.

[0207] The processor 901 can include one or more processing cores, such as a 4-core processor, an 8-core processor, or the like. The processor 901 can be implemented in the form of at least one of a DSP (Digital Signal Processing), an FPGA (Field-Programmable Gate Array), a PLA (Programmable Logic Array). The processor 901 can also include a main processor and a coprocessor. The main processor is a processor for processing data in an awake state, also referred to as a CPU (Central Processing Unit). The coprocessor is a low-power processor for processing data in a standby state. In some embodiments, the processor 901 can be integrated with a GPU (Graphics Processing Unit) that is responsible for rendering and drawing the content to be displayed on the display screen. In some embodiments, the processor 901 can further include an AI (Artificial Intelligence) processor for processing machine learning-related computing operations.

[0208] The memory 902 can include one or more computer-readable storage media, which can be non-transitory. The memory 902 can also include a high-speed random access memory, and a non-volatile memory such as one or more disk storage devices, flash storage devices. In some embodiments, the non-transitory computer-readable storage medium in the memory 902 is used to store at least one instruction for being executed by the processor 901 to implement the method for establishing the absorption and attenuation field of the surface layer of the desert area provided in the embodiments of the present application.

[0209] Those skilled in the art can understand that, Figure 9 The structure shown in FIG. 9 does not constitute a limitation on the computer device, and can include more or fewer components than those shown, or combine certain components, or adopt different component arrangements.

[0210] The embodiments of the present application also provide a non-transitory computer-readable storage medium, when the instructions in the storage medium are executed by the processor of the computer device, the computer device can execute the method for establishing the absorption and attenuation field of the surface layer of the desert area provided by the embodiments of the present application.

[0211] A computer program product containing instructions which, when run on a computer, cause the computer to perform the method for establishing an absorption attenuation field of a desert area surface layer provided by the embodiments of the present disclosure.

[0212] The above merely describes optional embodiments of the present disclosure, and is not intended to limit the present disclosure. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present disclosure shall be included in the protection scope of the present disclosure.

Claims

1. A method of establishing an absorption attenuation field of a desert surface layer, characterized by, The desert area surface layer comprises a first medium structure layer and a second medium structure layer, the first medium structure layer is a structure layer above the water table, and the second medium structure layer is a structure layer below the water table; The desert area surface layer is provided with a plurality of physical points; the method comprises: determining a characteristic function of the first medium structure layer and a characteristic function of the second medium structure layer, the characteristic function of the first medium structure layer being used to represent the variation law of the absorption attenuation parameter of the first medium structure layer, the characteristic function of the second medium structure layer being used to represent the variation law of the absorption attenuation parameter of the second medium structure layer, the characteristic function of the first medium structure layer being different from the characteristic function of the second medium structure layer; the first medium structure layer has a plurality of different characteristic functions; different regions in the lateral direction and different substructure layers in the longitudinal direction of the first medium structure layer correspond to different characteristic functions; obtaining a dune parameter of the position of a target physical point, the target physical point being one of the plurality of physical points, the dune parameter comprising a dune thickness, a slope and an elevation; determining a target characteristic function according to at least one of the ratio of the dune thickness to the slope and the included angle between the main wind direction and the tendency, the included angle between the main wind direction and the tendency being calculated according to the elevation, the target characteristic function being at least one of the plurality of characteristic functions of the first medium structure layer; and calculating a first absorption attenuation parameter value of the target physical point according to the target characteristic function; calculating second absorption attenuation parameter values of the plurality of physical points respectively according to the characteristic function of the second medium structure layer; and establishing an absorption attenuation field of the desert area surface layer according to the first absorption attenuation parameter values and the second absorption attenuation parameter values of the plurality of physical points.

2. The method of claim 1, wherein, The first medium structure layer comprises a first substructure layer and a second substructure layer in the longitudinal direction, the first substructure layer being located above the second substructure layer, and the determination of the target characteristic function according to at least one of the ratio of the dune thickness to the slope and the included angle between the main wind direction and the tendency comprises at least one of the following: In response to the ratio of the dune thickness to the slope being less than a first ratio threshold, determining the target feature function as a first feature function, the first feature function being used to calculate a first absorption attenuation parameter value of the first sub-structure layer; the first feature function is: Q(H)=A1•H α1 wherein A1=0.0014S+0.4849, α1=0.0021S+0.4141, S represents a dune reference parameter of the target physical point, the dune reference parameter of the target physical point being calculated by a dune parameter of a position where the target physical point is located; H is a depth of a desert area surface layer, A1 represents a coefficient of a fitting function, and α1 represents an index of the fitting function. in response to the ratio of the dune thickness to the slope being greater than a second ratio threshold value and the included angle between the main wind direction and the tendency being less than a first included angle threshold value, determining the target characteristic function as a second characteristic function and a third characteristic function, the second characteristic function and the third characteristic function being respectively used to calculate a first absorption attenuation parameter value of the first substructure layer and a first absorption attenuation parameter value of the second substructure layer; the first ratio threshold value being less than the second ratio threshold value; The second characteristic function is: Q(H)=A2•H α2 wherein A2=-0.0017S+1.8517, a2=-0.00003S+0.9431; A2 represents the coefficient of the fitting function, and a2 represents the index of the fitting function. The third characteristic function is: Q(H)=A3•H α3 wherein A3=0.0008S+1.7647, a3=0.0002S+0.7284; A3 represents the coefficient of the fitting function, and a3 represents the index of the fitting function. in response to the ratio of the dune thickness to the slope being greater than a second ratio threshold value and the included angle between the main wind direction and the tendency being greater than a second included angle threshold value, determining the target characteristic function as a fourth characteristic function and a fifth characteristic function, the fourth characteristic function and the fifth characteristic function being respectively used to calculate a first absorption attenuation parameter value of the first substructure layer and a first absorption attenuation parameter value of the second substructure layer; the first included angle threshold value being less than the second included angle threshold value; The fourth characteristic function is: Q(H)=A4•H α4 wherein A4=-0.0003S+1.1162, a4=0.0001S+0.9437; A4 represents the coefficient of the fitting function, and a4 represents the index of the fitting function. The fifth characteristic function is: Q(H) = A5•H α5 wherein A5 = 0.0039S + 0.0479, a5 = -0.0006S + 1.0386; A5 represents the coefficient of the fitting function, and a5 represents the index of the fitting function. determining the target feature function as a sixth feature function in response to the ratio of the dune thickness to the slope being greater than a second ratio threshold value, the angle between the main wind direction and the tendency being greater than a first angle threshold value and less than a second angle threshold value, the sixth feature function being used to calculate a first absorption attenuation parameter value of the first sub-structure layer and a first absorption attenuation parameter value of the second sub-structure layer; The sixth characteristic function is: Ψ B0 represents the first included angle threshold, Ψ Y0 represents the second included angle threshold, Ψ represents an included angle between the main wind direction and the tendency, when Q B is a first absorption attenuation parameter value of the first substructure layer calculated by using the second characteristic function and Q Y represents a first absorption attenuation parameter value of the first substructure layer calculated by using the fourth characteristic function, Q represents the calculated first absorption attenuation parameter value of the first substructure layer; when Q B is a first absorption attenuation parameter value of the first substructure layer calculated by using the third characteristic function and Q Y represents a first absorption attenuation parameter value of the first substructure layer calculated by using the fifth characteristic function, Q represents the calculated first absorption attenuation parameter value of the second substructure layer; determining the target feature function as a seventh feature function in response to the ratio of the dune thickness to the slope being greater than a first ratio threshold value and less than a second ratio threshold value, the seventh feature function being used to calculate a first absorption attenuation parameter value of the first sub-structure layer; The seventh characteristic function is: wherein P S0 represents the first ratio threshold value, P S1 represents the second ratio threshold value, P S represents the ratio of the dune thickness to the slope, Q H represents the first absorption attenuation parameter value of the first sub-structure layer calculated using the sixth characteristic function, Q P represents the first absorption attenuation parameter value of the first sub-structure layer calculated using the first characteristic function.

3. The method according to any one of claims 1 to 2, characterized in that, the calculating the second absorption attenuation parameter values of the plurality of physical points according to the feature function of the second medium structure layer comprises: calculating the second absorption attenuation parameter value of the target physical point according to an eighth feature function: Q = -0.0004S 2 +0.4484S - 74.86, wherein S represents a dune reference parameter of the target physical point, which is calculated by a dune parameter of a location where the target physical point is located.

4. The method of claim 3, wherein, the desert area surface layer is provided with a plurality of measuring points, and the method further comprises: correcting the first absorption attenuation parameter value of the target physical point according to a first absorption attenuation parameter difference value of the plurality of measuring points to obtain a corrected first absorption attenuation parameter value of the target physical point, the first absorption attenuation parameter difference value being a difference value between a third absorption attenuation parameter value of the measuring point measured and a fourth absorption attenuation parameter value of the measuring point calculated; correcting the second absorption attenuation parameter value of the target physical point according to a second absorption attenuation parameter difference value of the plurality of measuring points to obtain a corrected second absorption attenuation parameter value of the target physical point, the second absorption attenuation parameter difference value being a difference value between a fifth absorption attenuation parameter value of the measuring point measured and a sixth absorption attenuation parameter value of the measuring point calculated; and establishing an absorption attenuation field of the desert area surface layer according to the corrected first absorption attenuation parameter value and the corrected second absorption attenuation parameter value of the target physical point.

5. The method of claim 4, wherein, the establishing the absorption attenuation field of the desert area surface layer according to the corrected first absorption attenuation parameter value and the corrected second absorption attenuation parameter value of the target physical point comprises: calculating a first relative absorption attenuation total amount of the target physical point according to the corrected first absorption attenuation parameter value and the corrected second absorption attenuation parameter value of the target physical point; calculating a second relative absorption attenuation total amount of the plurality of measuring points according to the third absorption attenuation parameter value and the fifth absorption attenuation parameter value of the plurality of measuring points measured; calculating a third relative absorption attenuation total amount of the plurality of measuring points according to the fourth absorption attenuation parameter value and the sixth absorption attenuation parameter value of the plurality of measuring points calculated; obtaining a relative absorption attenuation total amount difference value of the plurality of measuring points according to the second relative absorption attenuation total amount and the third relative absorption attenuation total amount; correcting the first relative absorption attenuation total amount of the target physical point according to the relative absorption attenuation total amount difference value of the plurality of measuring points to obtain a corrected first relative absorption attenuation total amount of the target physical point; and establishing the absorption attenuation field of the desert area surface layer according to the corrected first relative absorption attenuation total amount of the target physical point.

6. A device for establishing an absorption attenuation field of a desert surface layer, characterized in that The desert area surface layer comprises a first medium structure layer and a second medium structure layer, the first medium structure layer is a structure layer above the water table, and the second medium structure layer is a structure layer below the water table. The desert area surface layer is provided with a plurality of physical points; the device comprises: A determination module is configured to determine a characteristic function of the first medium structure layer and a characteristic function of the second medium structure layer, the characteristic function of the first medium structure layer is used to represent a variation law of an absorption attenuation parameter of the first medium structure layer, the characteristic function of the second medium structure layer is used to represent a variation law of an absorption attenuation parameter of the second medium structure layer, the characteristic function of the first medium structure layer is different from the characteristic function of the second medium structure layer; the first medium structure layer has a plurality of different characteristic functions; different characteristic functions correspond to different regions in the lateral direction and different substructure layers in the longitudinal direction of the first medium structure layer; A first calculation module is configured to obtain a dune parameter of a position where a target physical point is located, the target physical point is one of the plurality of physical points, the dune parameter comprises a dune thickness, a slope and an elevation; determine a target characteristic function according to at least one of a ratio of the dune thickness to the slope and an included angle between a main wind direction and a tendency, the included angle between the main wind direction and the tendency is calculated according to the elevation, the target characteristic function is at least one of a plurality of characteristic functions of the first medium structure layer; and calculate a first absorption attenuation parameter value of the target physical point according to the target characteristic function; A second calculation module is configured to calculate second absorption attenuation parameter values of the plurality of physical points respectively according to the characteristic function of the second medium structure layer; A building module is configured to build an absorption attenuation field of the desert area surface layer according to the first absorption attenuation parameter values and the second absorption attenuation parameter values of the plurality of physical points.

7. A computer device, characterized by Comprise: A processor; A memory for storing processor-executable instructions; The processor is configured to execute the method of any one of claims 1 to 5.

8. A computer-readable storage medium, characterized in that, When the instructions in the computer-readable storage medium are executed by the processor of the computer device, the computer device can execute the method of any one of claims 1 to 5. When the instructions in the computer-readable storage medium are executed by the processor of the computer device, the computer device can execute the method of any one of claims 1 to 5.

Images