A seismic profile data visualization method and system based on adaptive mosaicking

Abstract

The present invention relates to a method and system for visualizing seismic profile data based on adaptive mosaicking, comprising: coloring grid vertices and pixels between vertices corresponding to seismic samples of different amplitudes with different colors; calculating, in a pixel shader, the sample point on the next trace most relevant to the current seismic sample point; mapping the relative offset of the most relevant sample point in the time direction into a tilt pattern and rendering it to a texture; and fully considering the tilt patterns of adjacent sample points on the same trace when filling vertex indices to achieve adaptive mosaicking of vertices. The mosaicking method proposed in the present invention can be adaptively selected based on the tilt patterns of seismic sample points, effectively resolving mosaicking errors in mesh mosaicking during three-dimensional visualization of two-dimensional seismic profiles, effectively improving the visualization of seismic microstructures, and facilitating the interpretation and analysis of seismic data. This method has important implications for fields such as oil and gas exploration and geological exploration.

Description

Technical Field

[0001] The present invention relates to the field of three-dimensional visualization of seismic profiles, and in particular to a method and system for visualizing seismic profile data based on adaptive mosaicking. Background Art

[0002] Seismic data profile visualization technology plays a key role in the interpretation process of seismic exploration. Visualization of seismic data that is intuitive, accurate, and has distinct features can greatly facilitate geological researchers to make accurate analysis of relevant data and draw correct conclusions, which is of great significance to fields such as oil and gas exploration and geological exploration.

[0003] Experts and scholars at home and abroad have conducted extensive research on the visualization of seismic profile data and have made significant progress. Seismic data visualization techniques are primarily categorized into two types: two-dimensional (2D) and three-dimensional (3D). 2D visualization converts seismic profile data into a 2D image for display. This can be categorized into four types: waveform display, waveform + variable area display, grayscale display, and color-variable density display. While these traditional visualization methods demonstrate the spatial variation of seismic data to some extent, they suffer from a limited dynamic range, significant loss of detail, and inability to intuitively display microscopic geological structures, hindering seismic interpretation and geological analysis. 3D visualization of seismic data has been extensively studied and has achieved significant progress. This technology visualizes seismic data in a 3D space, addressing the shortcomings of 2D visualization methods. However, these methods typically visualize the 3D data volume without extracting specific seismic profiles for grayscale or color-variable density display. Consequently, these visualization methods do not substantially improve the ability to visualize seismic microstructures. Summary of the Invention

[0004] In view of the above problems, the purpose of the present invention is to provide a seismic profile data visualization method and system based on adaptive mosaicking, which effectively solves the mosaicking error problem existing in the grid mosaicking in the three-dimensional visualization of two-dimensional seismic profiles.

[0005] To achieve the above objectives, the present invention adopts the following technical solutions: a method for visualizing seismic profile data based on adaptive mosaicking, which includes: coloring grid vertices and pixels between vertices corresponding to seismic samples of different amplitudes with different colors; calculating the sample point on the next channel that is most relevant to the current seismic sample point in a pixel shader; mapping the relative offset of the most relevant sample point in the time direction into a tilted manner and rendering it to a texture, and when filling in the vertex index, fully considering the tilt manner of adjacent sample points on the same channel to achieve adaptive mosaicking of vertices.

[0006] Furthermore, the grid vertices and pixels between the vertices corresponding to earthquake sample points of different amplitudes are colored with different colors using a linear progressive palette coloring method.

[0007] Furthermore, the linear progressive palette coloring method includes:

[0008] Create a one-dimensional floating-point texture. Each texture element is a color vector. When the texture coordinate U changes from 0 to 1.0, the color vector first gradually changes from pure red to pure white and then to pure blue.

[0009] After creating a linear progressive palette, the two-dimensional positions of the rasterized and interpolated seismic sample points are normalized into texture coordinates in a programmable pixel shader to sample the floating-point texture storing the profile data, and the amplitude values ​​at the vertices and any positions between the vertices are obtained.

[0010] The amplitude of the obtained seismic sample points is normalized into texture coordinates, and then the texture coordinates are used to sample the palette texture to achieve accurate coloring of the grid vertices and the pixels inside the triangular units between the vertices.

[0011] Furthermore, the method of calculating the sample point on the next line that is most correlated with the current earthquake sample point in the pixel shader includes: connecting the vertices corresponding to the earthquake sample points in the three-dimensional space into triangular units to form a surface mesh; calculating the correlation coefficient between a sample point sequence of a certain length centered on the current earthquake sample point and a sequence of the same length on the next line within a certain sliding window in the pixel shader, and finding the sliding value corresponding to the maximum correlation coefficient.

[0012] Furthermore, the tilting mode used for tessellation of a quadrilateral with the current vertex as a control point is determined according to the sliding value, and the determined tilting mode is stored in a texture element by rendering to a texture, and is used to fill the mesh vertex index.

[0013] Furthermore, the tilted tessellation method adopted for the quadrilateral with the current vertex as the control point is determined according to the sliding value, including: if the sliding value is greater than or equal to 0, downward tilted tessellation is adopted, represented by -1, and the quadrilateral with the current vertex as the control point adopts downward tilted tessellation; conversely, if the sliding value is less than 0, the quadrilateral with the current vertex as the control point adopts upward tilted tessellation, represented by +1, and the tilted method digitized in this way is then stored in the texture element by rendering to the texture.

[0014] Furthermore, the pixel shader is a GPU pixel shader.

[0015] A seismic profile data visualization system based on adaptive mosaicking includes: a coloring module that colors grid vertices and pixels between vertices corresponding to seismic samples of different amplitudes with different colors; a correlation module that calculates the sample point on the next channel that is most relevant to the current seismic sample point in a pixel shader; and a mosaicking module that maps the relative offset of the most relevant sample point in the time direction into a tilted pattern and renders it to a texture. When filling in vertex indices, the tilt patterns of adjacent sample points on the same channel are fully considered to achieve adaptive mosaicking of vertices.

[0016] A computer-readable storage medium storing one or more programs, wherein the one or more programs include instructions that, when executed by a computing device, cause the computing device to perform any of the above methods.

[0017] A computing device comprises: one or more processors, a memory, and one or more programs, wherein the one or more programs are stored in the memory and configured to be executed by the one or more processors, and the one or more programs include instructions for executing any of the above methods.

[0018] The present invention has the following advantages due to the adoption of the above technical solution:

[0019] 1. This paper proposes a new low-tilt-correlation adaptive mosaicking method, which effectively solves the mosaicking error problem existing in the 3D visualization of 2D seismic profiles. It can effectively improve the display capability of microstructures and facilitate the interpretation and analysis of seismic data. It is of great significance for oil and gas exploration, geological exploration and other fields.

[0020] 2. The seismic data profile visualization technology used in the present invention plays a key role in the interpretation link in seismic exploration. The intuitive, accurate, and clearly layered visualization of seismic data can greatly facilitate geological researchers to make accurate analysis of relevant data and draw correct conclusions, which is of great significance to the fields of oil and gas exploration, geological exploration, etc. BRIEF DESCRIPTION OF THE DRAWINGS

[0021] Figure 1 is a flow chart of a method for visualizing seismic profile data based on adaptive mosaicking in one embodiment of the present invention;

[0022] Figure 2 This is a diagram showing the principle of calculating the cross-correlation function value in one embodiment of the present invention;

[0023] Figure 3a This is a schematic diagram of an upper edge point mosaic in one embodiment of the present invention;

[0024] Figure 3b This is a schematic diagram of the lower edge point mosaic in one embodiment of the present invention;

[0025] Figure 4a Schematic diagram of the triangle unit intersection phenomenon caused by downward tilted mosaic in one embodiment of the present invention;

[0026] Figure 4b This is a schematic diagram of an upward tilted mosaic in one embodiment of the present invention;

[0027] Figure 5a Schematic diagram of a hole phenomenon caused by tilting at points other than the upper and lower edges in one embodiment of the present invention;

[0028] Figure 5b This is a schematic diagram of patching an upward-inclined leak in one embodiment of the present invention;

[0029] Figure 5c Schematic diagram of oblique seam repair of a downward-sloping leak in one embodiment of the present invention;

[0030] Figure 6 It is a schematic diagram of the structure of a computing device in one embodiment of the present invention. DETAILED DESCRIPTION

[0031] To make the objectives, technical solutions, and advantages of the embodiments of the present invention more clear, the technical solutions of the embodiments of the present invention will be clearly and completely described below in conjunction with the accompanying drawings of the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the described embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field are within the scope of protection of the present invention.

[0032] It should be noted that the terms used herein are only for describing specific embodiments and are not intended to limit the exemplary embodiments according to the present application. As used herein, unless the context clearly indicates otherwise, the singular form is also intended to include the plural form. In addition, it should be understood that when the terms "comprise" and / or "include" are used in this specification, they indicate the presence of features, steps, operations, devices, components and / or combinations thereof.

[0033] The present invention provides a method and system for visualizing seismic profile data based on adaptive mosaicking. First, a pixel shader calculates the sample point on the next trace that is most relevant to the current sample point. Then, its relative offset in the time direction is mapped to a tilt pattern and rendered to a texture. When populating vertex indices, the tilt patterns of adjacent sample points on the same trace are fully considered to achieve adaptive mosaicking of vertices. The mosaicking method proposed in the present invention can be adaptively selected based on the tilt patterns of seismic sample points, effectively resolving mosaicking errors that can occur in mesh mosaicking during three-dimensional visualization of two-dimensional seismic profiles. This method can effectively improve the visualization of seismic microstructures and facilitate the interpretation and analysis of seismic data, with significant implications for oil and gas exploration, geological exploration, and other fields.

[0034] In one embodiment of the present invention, Figure 1 As shown, a method for visualizing seismic profile data based on adaptive mosaicking is provided. In this embodiment, the method includes the following steps:

[0035] 1) Color the grid vertices and pixels between vertices corresponding to earthquake samples of different amplitudes with different colors;

[0036] 2) Calculate the sample point on the next trace that is most relevant to the current earthquake sample point in the pixel shader;

[0037] 3) The relative offset of the most relevant sample points in the time direction is mapped into a tilted manner and rendered to the texture. When filling the vertex index, the tilt of adjacent sample points on the same channel is fully considered to achieve adaptive tessellation of vertices.

[0038] In the above step 1), the grid vertices and pixels between the vertices corresponding to earthquake sample points of different amplitudes are colored with different colors using a linear progressive palette coloring method.

[0039] The linear progressive palette shading method includes the following steps:

[0040] 1.1) Create a one-dimensional floating-point texture where each texture element is a color vector. When the texture coordinate U changes from 0 to 1.0, the color vector first gradually changes from pure red to pure white and then to pure blue.

[0041] 1.2) After creating a linear progressive palette, the rasterized and interpolated two-dimensional positions (x, y) of the seismic sample points are normalized into texture coordinates in a programmable pixel shader to sample the floating-point texture storing the profile data, obtaining the amplitude values ​​at the vertices and any positions between them;

[0042] 1.3) The obtained seismic sample amplitudes are normalized into texture coordinates, and then the palette texture is sampled using the texture coordinates to achieve accurate coloring of the mesh vertices and the pixels inside the triangular units between the vertices.

[0043] In step 2) above, the pixel shader calculates the sample point on the next trace that is most relevant to the current earthquake sample point, including the following steps:

[0044] 2.1) The vertices corresponding to the earthquake sample points in three-dimensional space are connected into triangular units to form a surface grid;

[0045] 2.2) In the pixel shader, calculate the correlation coefficient between a sample sequence of a certain length centered on the current earthquake sample and a sequence of the same length on the next trace within a sliding window (sliding from top to bottom), and find the sliding value corresponding to the maximum correlation coefficient; wherein the pixel shader is a GPU pixel shader.

[0046] The principle of calculating the cross-correlation function value of two sample point sequences on adjacent seismic traces in GPU pixel shader is as follows: Figure 2 shown.

[0047] The calculation formulas for the cross-correlation function and correlation coefficient are as follows:

[0048]

[0049]

[0050] Wherein, formula (1) is to calculate the cross-correlation function value R of sequence X and Y XY The formula of (τ) is as follows: the subscript XY represents two sequences, and the order of X and Y means that sequence X is not moved and sequence Y is moved by τ units; R' in formula (2) XY (τ) represents the correlation coefficient of sequences X and Y; τ is the shift parameter, T is the sequence period, N represents the maximum value of the element number in sequences X and Y, and k represents the element number, k∈[0,N].

[0051] After calculating the correlation coefficient, find the maximum value and use the corresponding τ value to determine whether the current earthquake sample is tilted upward, downward, or not tilted. Upward tilt is represented by 1, downward tilt is represented by -1, and a small tilt is considered not tilted and represented by 0. In other words, the τ value corresponding to the maximum correlation coefficient is inverted to obtain the tilt mode of the earthquake sample. Then, the tilt mode is rendered into the R component of the texture element of the tilt mode texture by rendering it to the texture, and the other components are rendered as 1.0.

[0052] In step 2.2) above, the tilting method used for tessellation of the quadrilateral with the current vertex as the control point is determined according to the sliding value, and the determined tilting method is stored in the texture element by rendering to the texture to fill the mesh vertex index.

[0053] In this embodiment, the tilted tessellation method used for the quadrilateral with the current vertex as the control point is determined according to the sliding value, specifically:

[0054] If the sliding value is greater than or equal to 0, downward tilt tessellation is adopted, represented by -1, and the quadrilateral with the current vertex as the control point (the vertex in the upper left corner of the quadrilateral) adopts downward tilt tessellation. Conversely, if the sliding value is less than 0, the quadrilateral with the current vertex as the control point adopts upward tilt tessellation, represented by +1. The tilt mode digitized in this way is then stored in the texture element by rendering to the texture.

[0055] The low-tilt correlation adaptive mosaic method mentioned above has certain shortcomings. It is necessary to make each vertex mosaic according to the tilt mode of the corresponding sample point as much as possible (for boundary vertices and areas where holes occur, in order to ensure the integrity and continuity of the grid, the vertex mosaic mode may be inconsistent with the tilt mode of the seismic sample point). To avoid the shortcomings of this low-tilt correlation adaptive mosaic method, the present invention proposes a new low-tilt correlation adaptive mosaic method, which divides the tilt mode of the seismic sample point into two cases: upward tilt and downward tilt (including a tilt value of 0).

[0056] The top edge samples are mosaicked with downward tilt, as shown below Figure 3a 、 Figure 3b As shown. Among them, Figure 3a is a schematic diagram of the upper edge point mosaic, Figure 3b The figure below shows the mosaicking diagram for the lower edge points. T and B are the seismic points to be mosaicked. They are located at the upper edge (with the smallest time-depth value) and the lower edge (with the largest time-depth value) of the seismic profile, respectively. This unified mosaicking method is used for the edge points to avoid situations where the top two or bottom two seismic points on the same seismic trace have opposite inclinations. Since there are no mosaicking triangles in the rectangular area between them, this would create a gap. Next, adaptive mosaicking can be performed on the remaining non-upper and lower edge points.

[0057] When the tilt mode of non-upper and lower edge points is upward tilt, in order to avoid the four sample points in the same rectangle being repeatedly mosaicked and causing triangle overlap, first check whether the tilt mode of the adjacent seismic sample points above it is downward tilt. If not, then perform upward tilt mosaic; if so, do not mosaic, such as Figure 4a 、 Figure 4b As shown. Among them, point B is the current sample point, Figure 4a This is the intersection phenomenon that occurs when the tilt mode of the current sample point A is downward tilt. At this time, point B is not mosaicked, that is, no vertex index is assigned to B, D, and C to form a triangle BDC; Figure 4b This is the case when point A is upwardly inclined.

[0058] When the tilt mode of non-upper and lower edge points is downward tilt, downward tilt mosaic is performed first. Then, in order to avoid loopholes, check whether the tilt mode of the adjacent sample point above it is upward tilt. If so, mosaic the rectangular area between the upper earthquake sample point and the current sample point to fill the loophole. If the tilt mode of the upper sample point is not upward tilt, no patching is performed. Figures 5a to 5c As shown. Where E is the current sample point, Figure 5a This is the leak phenomenon that occurs when point B is tilted upward. Figure 5b and Figure 5c There are two methods that can be used to patch holes.

[0059] In the above mosaic method proposed in this embodiment, the mosaic method of the four sample points in each quadrilateral unit is determined by the tilt method of the upper left corner sample point, which can maximize the true adaptive mosaicking of each seismic sample point according to its tilt method.

[0060] In one embodiment of the present invention, a seismic profile data visualization system based on adaptive mosaicking is provided, comprising:

[0061] The coloring module uses different colors to color the grid vertices and pixels between vertices corresponding to earthquake samples of different amplitudes;

[0062] The correlation module calculates the sample point on the next trace that is most correlated with the current earthquake sample point in the pixel shader;

[0063] The mosaic module maps the relative offset of the most relevant sample points in the time direction into a tilted manner and renders it to the texture. When filling the vertex index, the tilt manner of adjacent sample points on the same channel is fully considered to achieve adaptive mosaicking of vertices.

[0064] The system provided in this embodiment is used to execute the above-mentioned method embodiments. Please refer to the above-mentioned embodiments for specific processes and detailed contents, which will not be repeated here.

[0065] like Figure 6FIG. 1 is a schematic diagram of the structure of a computing device provided in one embodiment of the present invention. The computing device may be a terminal and may include: a processor, a communications interface, a memory, a display screen, and an input device. The processor, communications interface, and memory communicate with each other via a communications bus. The processor is used to provide computing and control capabilities. The memory includes a non-volatile storage medium and internal memory. The non-volatile storage medium stores an operating system and a computer program. When the computer program is executed by the processor, it implements a visualization method. The internal memory provides an environment for the operation of the operating system and computer program in the non-volatile storage medium. The communications interface is used to communicate with an external terminal via wired or wireless communication. The wireless communication may be achieved via Wi-Fi, a management network, NFC (near field communication), or other technologies. The display screen may be a liquid crystal display or an electronic ink display. The input device may be a touch screen layer covering the display screen, or may be a key, trackball, or touchpad provided on the computing device housing, or may be an external keyboard, touchpad, or mouse. The processor can call the logic instructions in the memory to execute the following method: color the grid vertices and pixels between the vertices corresponding to earthquake samples of different amplitudes with different colors; calculate the sample point on the next track that is most relevant to the current earthquake sample point in the pixel shader; map the relative offset of the most relevant sample point in the time direction into a tilt mode and render it to the texture; when filling the vertex index, fully consider the tilt mode of adjacent sample points on the same track to achieve adaptive tessellation of vertices.

[0066] In addition, the logical instructions in the above-mentioned memory can be implemented in the form of a software functional unit and can be stored in a computer-readable storage medium when sold or used as an independent product. Based on this understanding, the technical solution of the present invention, or the part that contributes to the prior art, or the part of the technical solution, can be embodied in the form of a software product. The computer software product is stored in a storage medium and includes several instructions for enabling a computer device (which can be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the method described in each embodiment of the present invention. The aforementioned storage medium includes: various media that can store program codes, such as a USB flash drive, a mobile hard disk, a read-only memory (ROM), a random access memory (RAM), a magnetic disk or an optical disk.

[0067] Those skilled in the art will understand that Figure 6The structure shown in the figure is only a block diagram of a part of the structure related to the solution of the present application, and does not constitute a limitation on the computing device to which the solution of the present application is applied. The specific computing device may include more or fewer components than shown in the figure, or combine certain components, or have a different component arrangement.

[0068] In one embodiment of the present invention, a computer program product is provided, which includes a computer program stored on a non-transitory computer-readable storage medium, and the computer program includes program instructions. When the program instructions are executed by a computer, the computer can execute the methods provided by the above-mentioned method embodiments, for example, including: coloring the grid vertices and pixels between the vertices corresponding to seismic samples of different amplitudes with different colors; calculating the sample point on the next track that is most relevant to the current seismic sample point in a pixel shader; mapping the relative offset of the most relevant sample point in the time direction into an oblique manner and rendering it to a texture, and when filling in the vertex index, fully considering the oblique manner of adjacent sample points on the same track to achieve adaptive mosaicking of vertices.

[0069] In one embodiment of the present invention, a non-transitory computer-readable storage medium is provided, which stores server instructions. The computer instructions enable a computer to execute the methods provided in the above embodiments, for example, including: coloring the grid vertices and pixels between the vertices corresponding to earthquake samples of different amplitudes with different colors; calculating the sample point on the next track that is most relevant to the current earthquake sample point in a pixel shader; mapping the relative offset of the most relevant sample point in the time direction into an oblique manner and rendering it to a texture, and when filling in the vertex index, fully considering the oblique manner of adjacent sample points on the same track to achieve adaptive mosaicking of vertices.

[0070] The above embodiment provides a computer-readable storage medium, whose implementation principle and technical effects are similar to those of the above method embodiment, and will not be repeated here.

[0071] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the steps in the process. Figure 1 a process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0072] These computer program instructions may also be stored in a computer readable memory that can direct a computer or other programmable data processing device to work in a specific manner, so that the instructions stored in the computer readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 a process or multiple processes and / or boxes Figure 1 The function specified in one or more boxes.

[0073] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operational steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing the instructions executed on the computer or other programmable device for implementing the process. Figure 1 a process or multiple processes and / or boxes Figure 1 A step that specifies a function in one or more boxes.

[0074] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the various embodiments of the present invention.

Claims

1. A seismic profile data visualization method based on adaptive mosaicking, characterized in that: include: The grid vertices and pixels between vertices corresponding to earthquake samples of different amplitudes are colored with different colors; In the pixel shader, calculate the sample point on the next trace that is most relevant to the current earthquake sample point; The relative offset of the most relevant sample points in the time direction is mapped to a tilted pattern and rendered to the texture. When filling the vertex index, the tilt patterns of adjacent sample points on the same channel are fully considered to achieve adaptive tessellation of vertices. The grid vertices and pixels between the vertices corresponding to earthquake samples of different amplitudes are colored with different colors using a linear progressive palette coloring method; The linear progressive palette shading method includes: Create a one-dimensional floating-point texture. Each texture element is a color vector. When the texture coordinate U changes from 0 to 1.0, the color vector first gradually changes from pure red to pure white and then to pure blue. After creating a linear progressive palette, the rasterized and interpolated two-dimensional positions of the seismic sample points are normalized into texture coordinates in a programmable pixel shader. The floating-point texture storing the profile data is sampled to obtain the amplitude values ​​at any vertex and any position between vertices. Normalize the amplitude of the seismic sample points obtained into texture coordinates, and then use the texture coordinates to sample the palette texture to achieve accurate coloring of the mesh vertices and the pixels inside the triangular units between the vertices; The step of calculating the sample point on the next trace that is most relevant to the current earthquake sample point in the pixel shader includes: The vertices corresponding to the seismic sample points in three-dimensional space are connected into triangular units to form a surface grid; In the pixel shader, calculate the correlation coefficient between a sample sequence of a certain length centered on the current earthquake sample point and a sequence of the same length on the next trace within a sliding window, and find the sliding value corresponding to the maximum correlation coefficient; Determine the tilting mode of the quadrilateral with the current vertex as the control point according to the sliding value, store the determined tilting mode into the texture element by rendering it into the texture element, and use it to fill the mesh vertex index; Determining, according to the sliding value, an inclined tessellation method of a quadrilateral with the current vertex as a control point, including: If the sliding value is greater than or equal to 0, the downward tilt tessellation is adopted, which is represented by -1. The quadrilateral with the current vertex as the control point adopts the downward tilt tessellation. Conversely, if the sliding value is less than 0, the quadrilateral with the current vertex as the control point adopts the upward tilt tessellation, which is represented by +1. The tilt mode digitized in this way is then rendered into the texture and stored in the texture element. The top edge sample points are mosaicked with downward tilt; when the tilt mode of non-upper and lower edge points is upward tilt, first check whether the tilt mode of the adjacent seismic sample point above it is downward tilt. If not, mosaicked with upward tilt; if so, no mosaicking is performed; when the tilt mode of non-upper and lower edge points is downward tilt, first mosaicked with downward tilt, check whether the tilt mode of the adjacent sample point above it is upward tilt. If so, mosaicked the rectangular area between the seismic sample point above and the current sample point to fill the gap. If the tilt mode of the sample point above is not upward tilt, no mosaicking is performed.

2. The seismic profile data visualization method according to claim 1, wherein: The pixel shader is a GPU pixel shader.

3. A seismic profile data visualization system based on adaptive mosaicking, used to implement the seismic profile data visualization method according to claim 1 or 2, characterized in that: include: The coloring module uses different colors to color the grid vertices and pixels between vertices corresponding to earthquake samples of different amplitudes; The correlation module calculates the sample point on the next trace that is most correlated with the current earthquake sample point in the pixel shader; The mosaic module maps the relative offset of the most relevant sample points in the time direction into a tilted manner and renders it to the texture. When filling the vertex index, it fully considers the tilt of adjacent sample points on the same channel to achieve adaptive mosaic of vertices.

4. A computer-readable storage medium storing one or more programs, characterized in that: The one or more programs include instructions that, when executed by a computing device, cause the computing device to perform any one of the methods of claims 1 to 2 .

5. A computing device, characterized in that include: One or more processors, a memory, and one or more programs, wherein the one or more programs are stored in the memory and configured to be executed by the one or more processors, and the one or more programs include instructions for executing any one of the methods according to claims 1 to 2.

Images