Method for identifying covering type asteroid impact crater
By identifying negative structures in three-dimensional seismic data bodies and screening features on multiple dimensions, combined with paleomorphic recovery, the problem of identifying underground cover asteroid impact craters is solved, and the accuracy and success rate of recognition are improved.
Patent Information
- Application Number
- CN202510122573.X
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2025-01-26
- Publication Date
- 2025-05-06
AI Technical Summary
The prior art is difficult to identify the impact craters of asteroids buried deep within a few thousand meters underground, and due to geological effects, the shape of the impact crater is difficult to preserve and is difficult to identify.
By obtaining three-dimensional seismic data bodies containing stratigraphic hierarchical information and fault information, negative structures are identified, and featured negative structures are screened in three dimensions: seismic profile, horizontal slices and top structural maps, restoring the paleo-terrain for further screening.
It improves the recognition accuracy of underground covered asteroid impact craters, reduces omissions, and enhances the recognition success rate.
Smart Images
Figure CN119937009A_ABST
Abstract
Description
Technical Field
[0001] The present invention belongs to the technical field of seismic exploration. Specifically, the present invention relates to a method for identifying a covered asteroid impact crater. Background Art
[0002] Meteorite impact craters are pits or ring-shaped geological structures formed by asteroids or comets hitting the surface of planets and their satellites at high speed. They are of great significance to the exploration of the formation and evolution of the Earth, paleontological changes, rock formation and mineralization, and the state of matter deep inside the Earth. Asteroid impact craters are mainly circular or elliptical depressions developed on the surface of the Earth. Their shapes are relatively clear and can be determined through field surveys and measurements.
[0003] Current methods for identifying asteroid impact craters include surface survey, iridium element determination, high-pressure quartz or high-pressure planar deformation foliations (PDFs) identification, etc.
[0004] Among them, the surface survey method is to use topographic and geomorphic data or maps to identify circular or nearly circular depressions, and use the topographic and stratigraphic rock characteristics at the edge of the depression to identify the annular raised belt and crater lip structure formed by the meteorite impact, and identify the radial fractures and shatter cone structures produced by the impact inside and around the depression as the basis for identifying the impact crater.
[0005] The iridium element determination method is to collect clay samples from the surface of circular or nearly circular depressions or shallow drill holes, and use rare element analysis instruments to determine the platinum group iridium content of the samples. When the iridium content of the sample is dozens or hundreds of times higher than the average iridium content in the area, it can be used as a basis for identifying impact craters.
[0006] The high-pressure quartz or high-pressure planar deformation foliation identification method uses the planar deformation texture identification of special minerals and special mineral crystals formed by impact high pressure. In circular or nearly circular depressions, shallow wells are drilled to obtain granite cores in the suspected impact layer. Special impact high-pressure quartz minerals, such as coesite and stishovite, are identified under a polarizing microscope. The quartz grain crystals on these rock thin sections have planar deformation foliations (quartz planar deformation foliations (PDFs)) caused by high-pressure impact, which can be used as a basis for identifying impact craters.
[0007] The existing methods for identifying impact craters on the surface rely on obtaining rock samples or stratum core samples on the surface or by drilling shallow wells for elemental content or mineral type structure analysis. However, when the impact crater is located in an underground covering area buried several thousand meters deep, the existing identification methods are no longer applicable due to the lack of intuitive data, clay or rock samples. In addition, after a long period of erosion and sedimentation, the morphology of the impact crater is difficult to preserve intact and difficult to identify.
[0008] In view of this, there is an urgent need to provide a technical solution for identifying covered asteroid impact craters, so as to realize the identification of deep underground buried impact craters. Summary of the invention
[0009] In order to at least solve one or more of the technical problems mentioned above, the present invention provides a method for identifying underground-covered asteroid impact craters, comprising: a first step, obtaining a three-dimensional seismic data body containing stratigraphic information and fault information; a second step, randomly selecting a seismic profile from the three-dimensional seismic data body to identify negative structures; a third step, respectively intercepting seismic profiles, seismic horizontal slices and top surface structural maps of the negative structures, and screening negative structures with at least one of the following characteristics: a circular or nearly circular outline, a crater lip on the edge, a broken zone inside, and radial faults in the periphery, to obtain a structure set; a fourth step, restoring the paleo-geomorphology with the nearest marker layer above the top unconformity surface of the negative structure in the structure set as a reference, screening negative structures with circular or nearly circular outlines on the paleo-geomorphology, and obtaining underground-covered asteroid impact craters.
[0010] According to one embodiment of the present invention, in the first step, a three-dimensional seismic data volume is obtained by performing horizon and fault interpretation on the three-dimensional seismic data.
[0011] According to one embodiment of the present invention, in the second step, seismic sections are intercepted from the three-dimensional seismic data volume at intervals of integer multiples of the seismic line spacing to identify negative structures.
[0012] According to one embodiment of the present invention, a structure in which the seismic in-phase axis bends downward on the seismic section is identified as a negative structure.
[0013] According to one embodiment of the present invention, in the third step, intercepting the seismic profile of the negative structure includes: taking the center line of the negative structure as the axis, rotating and selecting the seismic profiles in different directions of the negative structure.
[0014] According to one embodiment of the present invention, the angle between adjacent seismic sections is 2-45 degrees.
[0015] According to one embodiment of the present invention, in the third step, intercepting the seismic section of the negative structure includes: intercepting at least two seismic sections that are orthogonal to each other.
[0016] According to one embodiment of the present invention, in the third step, intercepting the seismic horizontal slices of the negative structure includes: intercepting at least three of the slices located at the upper part, the middle part and the lower part of the negative structure.
[0017] According to one embodiment of the present invention, in the fourth step, the ancient landform is restored by layer leveling using the marker layer as a reference.
[0018] According to another aspect of the present invention, there is provided an identification device for underground-covered asteroid impact craters, comprising: a processor for executing program instructions; and a memory storing program instructions, wherein when the program instructions are loaded and executed by the processor, the processor executes the above-mentioned identification method.
[0019] In the present invention, negative structures are screened by randomly selecting seismic sections in the three-dimensional seismic data volume, and the negative structures are respectively screened in the seismic sections, horizontal slices and top surface structures. Figure 3 Identifying the characteristics of impact craters in three dimensions can overcome the problem that the characteristics of impact craters are difficult to identify due to geological processes such as erosion and sedimentation, reduce omissions, and improve the success rate of identification. The accuracy of impact crater identification can be improved by further screening after paleo-geomorphological restoration based on the structure set. By intercepting seismic profiles at intervals of integer multiples of the spacing between seismic survey lines, the omission of negative structures can be reduced while reducing the amount of data processing. By rotating and selecting seismic profiles in different directions of the negative structure with the center line of the negative structure as the axis, the characteristics of the impact craters can be identified in multiple directions. By intercepting multiple seismic horizontal slices on the negative structure, feature identification can be performed at different depths, further reducing omissions. BRIEF DESCRIPTION OF THE DRAWINGS
[0020] By reading the following detailed description with reference to the accompanying drawings, the above and other objects, features and advantages of the exemplary embodiments of the present invention will become readily understood. In the accompanying drawings, several embodiments of the present invention are shown in an exemplary and non-limiting manner, and the same or corresponding reference numerals represent the same or corresponding parts, wherein:
[0021] Figure 1 A schematic diagram of a three-dimensional seismic survey system is shown;
[0022] Figure 2 A schematic diagram of a three-dimensional seismic data volume is shown;
[0023] Figure 3 A schematic diagram showing the steps of a method for identifying underground covered asteroid impact craters;
[0024] Figure 4 A schematic diagram of a seismic section is shown;
[0025] Figure 5 A schematic diagram of another seismic section is shown;
[0026] Figure 6 A schematic diagram of a seismic horizontal slice is shown;
[0027] Figure 7 A schematic diagram of a top surface structural diagram is shown;
[0028] Figure 8A hardware schematic diagram of a method for identifying underground covered asteroid impact craters is shown. DETAILED DESCRIPTION
[0029] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative work are within the scope of protection of the present invention.
[0030] It should be understood that the terms "include" and "comprising" used in the description and claims of the present invention indicate the presence of described features, integers, steps, operations, elements and / or components, but do not exclude the presence or addition of one or more other features, integers, steps, operations, elements, components and / or combinations thereof.
[0031] It should also be understood that the terms used in this specification of the present invention are only for the purpose of describing specific embodiments and are not intended to limit the present invention. As used in the specification of the present invention and the claims, the singular forms of "a", "an" and "the" are intended to include the plural forms unless the context clearly indicates otherwise. It should also be further understood that the term "and / or" used in the specification of the present invention and the claims refers to any combination of one or more of the associated listed items and all possible combinations, and includes these combinations.
[0032] As used in this specification and claims, the term "if" may be interpreted as "when" or "upon" or "in response to determining" or "in response to detecting," depending on the context. Similarly, the phrase "if it is determined" or "if [described condition or event] is detected" may be interpreted as meaning "upon determination" or "in response to determining" or "upon detection of [described condition or event]" or "in response to detecting [described condition or event]," depending on the context.
[0033] The specific embodiments of the present invention are described in detail below with reference to the accompanying drawings.
[0034] Figure 1 A schematic diagram of a three-dimensional seismic survey system is shown.
[0035] like Figure 1As shown, in the system 100, a plurality of mutually spaced detectors 110 for detecting seismic waves are arranged on the surface 101 of the exploration target area, forming a detector array covering the target area on a plane, and these detectors 110 are connected to the seismic information processing device by wired or wireless connection, and a plurality of seismic sources 120 are also arranged. The seismic information processing device can perform preliminary processing on the seismic data. The working process of the three-dimensional seismic exploration system is: artificially excite the seismic sources 120 located at multiple positions to generate seismic waves, which are reflected from the boundary of the stratum 102 and received by the detector array to form seismic information collected on a plane and changing with time. The seismic information received by the detector array represents certain measurements of seismic wave energy as a function of time, such as displacement, velocity, wave impedance, pressure, etc. This information can be grouped in different ways, such as traces, sets, etc., and then processed or format converted according to the corresponding relationship between time and space to form a three-dimensional seismic data body in the form of a three-dimensional array, that is, high-quality three-dimensional seismic fine imaging data is obtained, and the three-dimensional seismic data body is formed by spatial stacking of interface points. The interpretation of three-dimensional seismic data can observe the morphology of geological interfaces from different directions and study the changes of geological bodies in three-dimensional space by cutting cross sections, longitudinal sections and horizontal slices.
[0036] Figure 2 A schematic diagram of a three-dimensional seismic data volume is shown.
[0037] Figure 3 A schematic diagram showing the steps of a method for identifying underground covered asteroid impact craters.
[0038] like Figure 3 As shown, a method 300 for identifying underground-covered asteroid impact craters comprises: a first step S301, obtaining a three-dimensional seismic data body containing stratigraphic information and fault information; a second step S302, randomly selecting a seismic profile from the three-dimensional seismic data body to identify negative structures; a third step S303, respectively intercepting seismic profiles, seismic horizontal slices and top surface structural maps of negative structures, screening negative structures having at least one of the following characteristics: a circular or nearly circular outline, a crater lip at the edge, a fracture zone inside, and radial faults around, to obtain a structure set; a fourth step S304, restoring paleo-geomorphology with the nearest marker layer above the unconformity surface at the top of the negative structure in the structure set as a reference, screening negative structures with circular or nearly circular outlines on the paleo-geomorphology, and obtaining underground-covered asteroid impact craters.
[0039] Figure 2 The three-dimensional seismic data volume 200 in Figure 1The three-dimensional seismic exploration system shown in the figure is obtained, and by interpreting the three-dimensional seismic data by layer and fault, a three-dimensional seismic data body is obtained. It contains stratigraphic layer information and fault information, and is marked with a marker layer, which can be used as a basis for layer flattening and fault restoration. The interpretation process of stratigraphic layer information and fault information can be achieved by existing technology, and the present invention will not be repeated.
[0040] In the second step, taking any seismic section refers to taking a seismic section at any angle or position. By selecting seismic sections at different angles and scanning the three-dimensional seismic data volume, the structure with the seismic in-phase axis bending downward on the seismic section can be identified as a negative structure. Figure 2 In the process of identifying negative structures, seismic sections are preferably intercepted from the three-dimensional seismic data volume at intervals of integer multiples of the seismic line spacing.
[0041] The identification process of negative structures is the preliminary operation of identifying impact craters. The most obvious feature of impact craters is negative structures. By scanning and screening negative structures in 3D seismic data, the sample size for identifying other features of impact craters can be reduced, thus improving the recognition efficiency.
[0042] In the third step, at least one feature of the impact crater is identified in different dimensions, including: a circular or nearly circular outline, a crater lip at the edge, a fracture zone inside, and radial faults around the periphery. The different dimensions include seismic profiles, seismic horizontal slices, and top surface structural maps.
[0043] Since the impact craters identified by the present invention are of underground cover type, due to factors such as geological erosion and geological movement in history, the characteristics of one aspect of the impact crater are not obvious or only some characteristics are retained in a certain dimension. In the present invention, a multi-dimensional recognition method is adopted to avoid omissions as much as possible. Although the method of only recognizing some features may lead to the misidentification of other structures as impact crater features, it is possible to avoid missing real impact craters in the structure set as much as possible, and the misidentified impact craters can be left for further processing in the subsequent paleo-structural restoration steps.
[0044] The method of intercepting the seismic profile of the negative structure includes: taking the center line of the negative structure as the axis, rotating and selecting seismic profiles in different directions of the negative structure. Preferably, the angle between adjacent seismic profiles is 2-45 degrees. Preferably, at least two mutually orthogonal seismic profiles are intercepted.
[0045] Figure 4 A schematic diagram of a seismic section is shown.
[0046] When identifying impact crater features on the seismic section of the negative structure, the seismic section is rotated to scan the negative structure with the centerline of the negative structure as the axis. It is preferred that the coverage of the seismic section slightly exceeds the negative structure.
[0047] like Figure 4 As shown, in a certain seismic profile 400, the negative structure 401 is in the shape of a downward depression, and radial faults 402 exist around it.
[0048] Figure 5 A schematic diagram of another seismic section is shown.
[0049] like Figure 5 As shown, in the seismic profile 500, the negative structure 501 is in the shape of a deeper depression with a crater lip 502 at the edge.
[0050] By combining the negative structural contour lines on multiple seismic sections of the same axis, the three-dimensional shape of the negative structure 401 can be obtained, and its contour can be identified as circular or near-circular. Near-circular structures include ellipses, semicircles, semi-ellipses and other structures with circular trends.
[0051] Figure 6 A schematic diagram of a seismic horizontal slice is shown.
[0052] like Figure 6 As shown, in the seismic horizontal slice 600, the negative structure 601 is nearly circular in shape.
[0053] Intercepting the seismic horizontal slices of the negative structure includes: intercepting at least three located at the upper, middle and lower parts of the negative structure. Intercepting multiple seismic horizontal slices of the same negative structure in the depth direction can determine the morphology of the negative structure at different depths. On the seismic horizontal slices near the bottom of the negative structure, the fracture zone inside it can also be identified, and on the seismic horizontal slices near the top of the negative structure, the crater lip on the edge can also be identified.
[0054] Figure 7 A schematic diagram of the top surface structure diagram is shown.
[0055] like Figure 7 As shown, the top surface structure diagram refers to a planar structure diagram of the top or bottom surface undulating shape changes represented by contour lines, which presents the contour shape formed by the contour lines. In the top surface structure diagram 700, the contour 701 of the negative structure is nearly circular.
[0056] In the present invention, even if some features of the underground covered impact crater are lost due to geological effects such as deposition and erosion, the success rate of identification can be increased by identifying the features of the impact crater in the negative structure through the above three dimensions. The structure set is a set of negative structures with at least one feature of the impact crater identified from seismic profiles, seismic horizontal slices, and top surface structural maps. That is, in the obtained structure set, the negative structure contains at least one feature of the impact crater.
[0057] In the fourth step, each negative structure in the structure set is restored to paleo-geomorphology. The specific steps include: Figure 4 As shown, the unconformity surface 403 in the negative structure is determined, and marker layers that have been marked in the three-dimensional seismic data body are selected above the unconformity surface 403, such as marker layer 404 and marker layer 405. The lowest marker layer 404 is selected as the reference layer, and the ancient landform is restored by layer leveling.
[0058] An unconformity surface refers to a surface formed by interrupting the seismic reflection axis in a negative structure, and the surface has the characteristic of discontinuity in the age of the upper and lower strata. The top of a negative structure refers to the top of a negative structural depression. The lowest interface in the depression is the interface of the asteroid impact layer. By using the nearest marker layer above the interface as a reference layer for paleo-geomorphology restoration, the original form of the asteroid impact crater formation period can be restored to a large extent. A marker layer refers to a layer or a group of rock layers with obvious characteristics that can be used as stratigraphic comparison markers. The stratigraphic information in the three-dimensional seismic data body of the present invention includes the position of each marker layer. Since the characteristics of the marker layer are obvious, the use of the marker layer to restore the paleo-structure can ensure the authenticity of the paleo-structure.
[0059] In the process of restoring the ancient structure, the layer flattening operation and fault restoration operation are performed based on the stratigraphic information and fault information in the three-dimensional seismic data body. On the basis of restoring the ancient landform of the negative structure, the circular or nearly circular outline of the negative structure is easier to identify than other features. Based on this, the present invention adopts the method of screening the negative structure with a circular or nearly circular outline on the ancient landform, so that the screening data processing amount of the negative structure in the structural concentration is smaller, the processing speed is faster, and the processing result is more accurate.
[0060] According to another aspect of the present invention, there is provided an identification device for underground-covered asteroid impact craters, comprising: a processor for executing program instructions; and a memory storing program instructions, wherein when the program instructions are loaded and executed by the processor, the processor executes the above-mentioned identification method.
[0061] In the present invention, negative structures are screened by randomly selecting seismic sections in the three-dimensional seismic data volume, and the negative structures are respectively screened in the seismic sections, horizontal slices and top surface structures. Figure 3Identifying the characteristics of impact craters in multiple dimensions can overcome the problem of difficult-to-identify characteristics caused by geological processes such as erosion and sedimentation, reduce omissions, and improve the success rate of identification. The accuracy of impact crater identification can be improved by further screening after paleo-geomorphological restoration based on the structure set. By intercepting seismic profiles at intervals of integer multiples of the spacing between seismic survey lines, the omission of negative structures can be reduced while reducing the amount of data processing. By rotating and selecting seismic profiles in different directions of the negative structure with the center line of the negative structure as the axis, the characteristics of the impact crater can be identified in multiple directions. By intercepting multiple seismic horizontal slices on the negative structure, feature identification can be performed at different depths, further reducing omissions.
[0062] Figure 8 A hardware schematic diagram of a method for identifying underground covered asteroid impact craters is shown.
[0063] The system 800 may include a device 801 according to an embodiment of the present invention, as well as its peripheral devices and an external network, wherein the device 801 is used to execute the operation of the identification method to implement the technical solution of the aforementioned embodiment of the present invention.
[0064] like Figure 8 As shown in , the device 801 may include a CPU 8011, which may be a general-purpose CPU, a dedicated CPU, or other information processing and program execution units. Further, the device 801 may also include a large-capacity memory 8012 and a read-only memory ROM 8013, wherein the large-capacity memory 8012 may be configured to store various types of data, seismic data, model data, descending unit division parameters, etc., and various programs required for performing various operations, and the ROM 8013 may be configured to store the power-on self-test of the device 801, the initialization of various functional modules in the system, the basic input / output driver of the system, and the data required for booting the operating system.
[0065] Further, the device 801 also includes other hardware platforms or components, such as the TPU 8014, GPU 8015, FPGA 8016 and MLU 8017 shown. It is understood that although a variety of hardware platforms or components are shown in the device 800, this is only exemplary and not restrictive, and those skilled in the art can add or remove corresponding hardware according to actual needs. For example, the device 801 can include only a CPU as a well-known hardware platform and another hardware platform as a test hardware platform of the present invention.
[0066] The device 801 of the present invention also includes a communication interface 8018, so that it can be connected to a local area network / wireless local area network (LAN / WLAN) 805 through the communication interface 8018, and then connected to a local server 808 or to the Internet ("Internet") 807 through the LAN / WLAN. Alternatively or additionally, the device 801 of the present invention can also be directly connected to the Internet or a cellular network through the communication interface 8018 based on wireless communication technology, such as wireless communication technology based on the third generation ("3G"), the fourth generation ("4G") or the fifth generation ("5G") . In some application scenarios, the device 801 of the present invention can also access a server 808 of an external network and a possible database 809 as needed to obtain various known algorithms, data and modules, and can remotely store various measured data.
[0067] The peripheral devices of the device 801 may include a display device 802, an input device 803, and a data transmission interface 804. In one embodiment, the display device 802 may include, for example, one or more speakers and / or one or more visual displays, which are configured to perform voice prompts and / or image video displays on the computing process or detection results of the device of the present invention. The input device 803 may include, for example, a keyboard, a mouse, a microphone, a gesture capture camera, or other input buttons or controls, which are configured to receive input or user instructions for training data. The data transmission interface 804 may include, for example, a serial interface, a parallel interface or a universal serial bus interface ("USB"), a small computer system interface ("SCSI"), a serial ATA, a FireWire ("FireWire"), a PCI Express, and a high-definition multimedia interface ("HDMI"), etc., which are configured for data transmission and interaction with other devices or systems.
[0068] The CPU 8011, the mass storage 8012, the read-only memory ROM 8013, the TPU 8014, the GPU 8015, the FPGA 8016, the MLU 8017 and the communication interface 8018 of the device 801 of the present invention can be connected to each other through a bus 8019, and data interaction with peripheral devices can be achieved through the bus. In one embodiment, through the bus 8019, the CPU 8011 can control other hardware components in the device 801 and its peripheral devices.
[0069] In operation, the processor CPU 8011 of the device 801 of the present invention can receive the three-dimensional seismic data body through the input device 803 or the data transmission interface 804, and retrieve the computer program instructions or codes (such as various programs for constructing the stratigraphic model) stored in the memory 8012 to process the received three-dimensional seismic data body to obtain the three-dimensional geometric model of the stratigraphic layer. After the CPU 8011 determines the oil and gas reservoir by executing the program instructions, it can be displayed on the display device 802 or output by voice prompts. In addition, the device 801 can also upload the oil and gas reservoir composition to the network, such as the remote database 809, through the communication interface 8018.
[0070] It should also be understood that any module, unit, component, server, computer, terminal or device that executes instructions of the present invention examples may include or otherwise access computer-readable media, such as storage media, computer storage media or data storage devices (removable) and / or non-removable) such as disks, optical disks or tapes. Computer storage media may include volatile and non-volatile, removable and non-removable media implemented in any method or technology for storing information, such as computer-readable instructions, data structures, program modules or other data.
[0071] The computer-readable storage medium may be any suitable magnetic storage medium or magneto-optical storage medium, such as Resistive Random Access Memory (RRAM), Dynamic Random Access Memory (DRAM), Static Random Access Memory (SRAM), Enhanced Dynamic Random Access Memory (EDRAM), High-Bandwidth Memory (HBM), Hybrid Memory Cube (HMC), etc., or any other medium that can be used to store the required information and can be accessed by an application, a module, or both. Any such computer storage medium may be part of a device or accessible or connectable to a device. Any application or module described in the present invention may be implemented using computer-readable / executable instructions that may be stored or otherwise maintained by such a computer-readable medium.
[0072] Although multiple embodiments of the present invention have been shown and described herein, it will be apparent to those skilled in the art that such embodiments are provided by way of example only. Those skilled in the art may conceive of many changes, modifications, and alternatives without departing from the thought and spirit of the present invention. It should be understood that in the process of practicing the present invention, various alternatives to the embodiments of the present invention described herein may be adopted. The appended claims are intended to define the scope of protection of the present invention, and therefore cover equivalents or alternatives within the scope of these claims.
Claims
1. A method for identifying underground covered asteroid impact craters, characterized in that: include: The first step is to obtain a three-dimensional seismic data volume containing stratigraphic layer information and fault information; The second step is to randomly select a seismic section from the three-dimensional seismic data volume to identify negative structures; The third step is to intercept the seismic profile, seismic horizontal slice and top surface structural map of the negative structure respectively, and select the negative structure with at least one of the following characteristics: circular or nearly circular outline, crater lip at the edge, fracture zone inside, radial fault at the periphery, and obtain the structural set; The fourth step is to restore the paleo-geomorphology with reference to the nearest marker layer above the top unconformity surface of the negative structure in the structural concentration, and to select negative structures with circular or nearly circular outlines on the paleo-geomorphology to obtain underground covered asteroid impact craters.
2. The identification method according to claim 1, characterized in that: In the first step, the three-dimensional seismic data volume is obtained by performing horizon and fault interpretation on the three-dimensional seismic data.
3. The identification method according to claim 1, characterized in that: In the second step, seismic sections are intercepted from the three-dimensional seismic data volume at intervals of integer multiples of the seismic line spacing to identify negative structures.
4. The identification method according to claim 3, characterized in that: On the seismic section, structures with downward bending of the seismic in-phase axis are identified as negative structures.
5. The identification method according to claim 1, characterized in that: In the third step, intercepting the seismic profile of the negative structure comprises: Taking the center line of the negative structure as the axis, seismic sections in different directions are rotated and selected.
6. The identification method according to claim 5, characterized in that: The angle between adjacent seismic sections is 2-45 degrees.
7. The identification method according to claim 1, characterized in that: In the third step, intercepting the seismic section of the negative structure includes: intercepting at least two mutually orthogonal seismic sections.
8. The identification method according to claim 1, characterized in that: In the third step, intercepting the seismic horizontal slices of the negative structure includes: intercepting at least three of the slices located at the upper part, the middle part and the lower part of the negative structure.
9. The identification method according to claim 1, characterized in that: In the fourth step, the ancient landform is restored by layer leveling using the marker layer as a reference.
10. An identification device for underground covered asteroid impact craters, characterized in that: include: a processor for executing program instructions; and a memory storing program instructions, which, when loaded and executed by a processor, enables the processor to execute the identification method according to any one of claims 1 to 9.