A method and device for preprocessing and optimizing basic data of seismic interpretation system

Abstract

The present invention provides a method and device for preprocessing and optimizing basic data of a seismic interpretation system, the method comprising establishing different user environments, including establishing various parameter files; preprocessing Omega stacking velocity data; assisting Z-Map_Plus drawing software; optimizing GeoFrame basic data; providing data query, unit conversion, azimuth calculation, parameter extraction, and data management functions. The present application preprocesses and optimizes various data to achieve the purpose of achieving smooth docking of various application software, and adopts data segment overlap recognition technology to eliminate data redundancy and improve the access efficiency of GeoFrame old work areas.

Description

Technical Field

[0001] The invention belongs to the field of oil and gas field exploration, and in particular, relates to a method and a device for preprocessing and optimizing basic data of a seismic interpretation system. Background Art

[0002] LandMark Depth Team variable speed mapping software cannot directly load the stacked velocity spectrum format file output by Omega processing software. There are two main reasons: first, the coordinates of the spectrum points are not expressed in geodetic coordinates, but in pile numbers; second, there are abnormal data in the velocity spectrum TV data pair. The layer or section data output by GeoFrame interpretation software is in the time domain, while Z-Map-Plus structural mapping software requires depth domain interpretation data and survey network data. That is, the connection between Omega and LandMarkDepth Team, GeoFrame and Z-Map-Plus application software is not smooth. GeoFrame has been used in the unit for 30 years. The data sources of some old work areas are quite complex, and the data redundancy in the work areas is serious, which seriously affects the efficiency of data access in the work areas.

[0003] Manually modifying the velocity spectrum output by Omega is very cumbersome and time-consuming due to the large amount of data. The time domain layer conversion can be converted into the depth domain using variable speed mapping software. However, sometimes the lack of velocity bodies or the small mapping range makes time-depth conversion inconvenient. The additional measurement network on the structural map can be obtained by outputting the standard interpretation grid data, but it is inconvenient for users. Various data redundancies can be solved by continuous manual modification, but the time spent is unbearable for users. Summary of the invention

[0004] The embodiment of the present application provides a method and device for preprocessing and optimizing basic data of a seismic interpretation system, which preprocesses and optimizes various data to achieve smooth docking of various application software, and adopts data segment overlap recognition technology to eliminate data redundancy and improve the access efficiency of GeoFrame old work areas.

[0005] In a first aspect, an embodiment of the present application provides a method for preprocessing and optimizing basic data of a seismic interpretation system, comprising:

[0006] Establish different user environments, including establishing various parameter files;

[0007] Preprocess Omega stacking velocity data;

[0008] Auxiliary Z-Map_Plus drawing software;

[0009] Optimize GeoFrame basic data;

[0010] Provides data query, unit conversion, azimuth calculation, parameter extraction and data management functions.

[0011] Among them, preprocessing Omega stacking velocity data includes:

[0012] Filter outliers in the velocity spectrum through the constraints of the Dix formula, and output the velocity spectrum in VELF or XYTV format to provide auxiliary support for LandMark Depth Team velocity spectrum loading;

[0013] The formula for converting pile number into geodetic coordinates is:

[0014] x=x0-co2d*ae*cos(az2d)+an*sin(az2d)

[0015] y=y0+co2d*ae*sin(az2d)+an*cos(az2d)

[0016] Among them: x0, y0 are the coordinates of the origin of the two-dimensional survey network, az2d is the azimuth of the main survey line, an, ae are the stake numbers of the main survey line and the connecting survey line respectively, x, y are the geodetic coordinates, and co2d is the direction of progress of the main survey line.

[0017] Among them, the auxiliary Z-Map_Plus drawing software includes:

[0018] The time domain layer data / faults are converted into the depth domain to generate a 2D / 3D survey network with a controllable plane range.

[0019] Among them, the optimization of GeoFrame basic data includes:

[0020] The interpretation data merging technology of multiple working areas is adopted. On the same survey line interpretation section, there are section faults that completely overlap or basically overlap; fault trajectory segments overlap on the same layer; and the recognition technology of whether there is overlap is judged by calculating the average distance between two data segments, and redundant data segments are automatically deleted.

[0021] Among them, it provides data query, unit conversion, azimuth calculation, parameter extraction, and data management functions, including:

[0022] Provides seismic line coordinate query, time-depth data query, unit conversion, azimuth calculation, GeoFrame 3D survey area parameter extraction, 3D survey area coordinate system and time-depth velocity table management functions;

[0023] The basic formula for checking the 2D pile number from the coordinates (x, y) is:

[0024] ae=(x-x0)*sin(az2d)+(y-y0)*cos(az2d)

[0025] an=-(x-x0)*cos(az2d)+(y-y0)*sin(az2d)

[0026] In a 3D measurement area, the formula for checking the 3D line number from the coordinates (x, y) is:

[0027] LineM=l1+(l2-l1) / (xx4-xx1)*(xt-xx1)

[0028] Cdpm=Cdp1+(Cdp2-Cdp1) / (yy2-yy1)*(yt-yy1)

[0029] Wherein: xt is the x coordinate of the normal projection of (x,y) on the line connecting points 1 and 4 of the 3D basic survey line coordinates; yt is the y coordinate of the normal projection of (x,y) on the line connecting points 1 and 2 of the 3D basic survey line coordinates; l1 and l2 are the main survey line numbers of the border; Cdp1 and Cdp2 are the connecting survey lines of points 1 and 2; xx4 and xx1 are the x coordinates of points 4 and 1; yy2 and yy1 are the y coordinates of points 2 and 1.

[0030] In a second aspect, the present application provides a device for preprocessing and optimizing basic data of a seismic interpretation system, comprising:

[0031] The establishment unit is used to establish different user environments, including establishing various parameter files;

[0032] A preprocessing unit, used for preprocessing Omega stacking velocity data;

[0033] Auxiliary unit, used to assist Z-Map_Plus mapping software;

[0034] Optimization unit, used to optimize GeoFrame basic data;

[0035] The tool unit is used to provide data query, unit conversion, azimuth calculation, parameter extraction, and data management functions.

[0036] Wherein, the pre-processing unit is used for:

[0037] Filter outliers in the velocity spectrum through the constraints of the Dix formula, and output the velocity spectrum in VELF or XYTV format to provide auxiliary support for LandMark Depth Team velocity spectrum loading;

[0038] The formula for converting pile number into geodetic coordinates is:

[0039] x=x0-co2d*ae*cos(az2d)+an*sin(az2d)

[0040] y=y0+co2d*ae*sin(az2d)+an*cos(az2d)

[0041] Among them: x0, y0 are the coordinates of the origin of the two-dimensional survey network, az2d is the azimuth of the main survey line, an, ae are the stake numbers of the main survey line and the connecting survey line respectively, x, y are the geodetic coordinates, and co2d is the direction of progress of the main survey line.

[0042] Wherein, the optimization unit is used for:

[0043] The interpretation data merging technology of multiple working areas is adopted. On the same survey line interpretation section, there are section faults that completely overlap or basically overlap; fault trajectory segments overlap on the same layer; and the recognition technology of whether there is overlap is judged by calculating the average distance between two data segments, and redundant data segments are automatically deleted.

[0044] In a third aspect, an embodiment of the present application provides a computer-readable storage medium having a computer program stored thereon, which implements the steps of any of the above methods when executed by a processor.

[0045] In a fourth aspect, an embodiment of the present application provides a computer device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the steps of any of the above methods when executing the program.

[0046] The method and device for preprocessing and optimizing basic data of the seismic interpretation system in the embodiment of the present application have the following beneficial effects:

[0047] The method for preprocessing and optimizing the basic data of the seismic interpretation system of this application includes establishing different user environments, including establishing various parameter files; preprocessing Omega stacking velocity data; assisting Z-Map_Plus drawing software; optimizing GeoFrame basic data; providing data query, unit conversion, azimuth calculation, parameter extraction, and data management functions. This application preprocesses and optimizes various data to achieve smooth docking of various application software, and adopts data segment overlap recognition technology to eliminate data redundancy and improve the efficiency of accessing GeoFrame old work areas. BRIEF DESCRIPTION OF THE DRAWINGS

[0048] Figure 1 A schematic flow chart of a method for preprocessing and optimizing basic data of a seismic interpretation system according to an embodiment of the present application;

[0049] Figure 2 A schematic diagram of the structure of a device for preprocessing and optimizing basic data of a seismic interpretation system according to an embodiment of the present application;

[0050] Figure 3 This is a schematic diagram of the structure of a computer device according to an embodiment of the present application. DETAILED DESCRIPTION

[0051] The present application is further described below in conjunction with the accompanying drawings and embodiments.

[0052] In the following introduction, the terms "first" and "second" are used for descriptive purposes only and should not be understood as indicating or implying relative importance. The following introduction provides multiple embodiments of the present invention, and different embodiments can be replaced or combined, so this application can also be considered to include all possible combinations of the same and / or different embodiments recorded. Therefore, if one embodiment includes features A, B, and C, and another embodiment includes features B and D, then this application should also be considered to include embodiments containing one or more of all other possible combinations of A, B, C, and D, even though the embodiment may not be clearly recorded in the following text.

[0053] The following description provides examples and does not limit the scope, applicability or examples set forth in the claims. Changes may be made to the functions and arrangements of the elements described without departing from the scope of the present application. Various processes or components may be appropriately omitted, substituted or added to each example. For example, the described method may be performed in an order different from the order described, and various steps may be added, omitted or combined. In addition, features described in some examples may be combined in other examples.

[0054] like Figure 1 As shown, the present application provides a method for preprocessing and optimizing basic data of a seismic interpretation system, including: S101, establishing different user environments, including establishing various parameter files; S102, preprocessing Omega stacking velocity data; S103, assisting Z-Map_Plus drawing software; S104, optimizing GeoFrame basic data; S105, providing data query, unit conversion, azimuth calculation, parameter extraction, and data management functions.

[0055] This application uses software programming technology to pre-process and optimize various data to achieve smooth docking of various application software, and adopts data segment overlap recognition technology to eliminate data redundancy and improve the efficiency of accessing GeoFrame old work areas.

[0056] The invention discloses a data preprocessing and optimization scheme for an interpretation system, the main functions of which are preprocessing Omega stacking velocity spectrum data, assisting Z-Map-Plus drawing and optimizing GeoFrame work area basic data.

[0057] The present invention solves the problem of converting pile numbers into coordinates through software programming technology, and identifies data anomalies through the dix formula, automatically deletes the anomalies, thereby facilitating the LandMarkDepthTeam superimposed variable speed spectrum loading and reducing the workload of editing anomalies. The present invention provides a simple time-depth conversion function according to user needs, and generates a measurement network data file, which facilitates Z-Map-Plus to construct a map. The present invention solves the optimization of well deviation data, profile faults, and plane fault trajectory data through software programming technology, eliminating data redundancy.

[0058] The technical idea of ​​this application is to analyze the problems existing in the loading of Omega stacking velocity data, the requirements for interface data of Z-Map-Plus drawing software, and various potential problems existing in the low efficiency of accessing GeoFrame old work areas, and to seek technical solutions to various problems. Ultimately, various data are preprocessed and optimized through software programming technology to achieve smooth docking of various application software, and data segment overlap recognition technology is used to eliminate data redundancy and improve the access efficiency of GeoFrame old work areas.

[0059] The action process of the present invention is described below:

[0060] The first step is to establish different user environments, including the creation of various parameter files.

[0061] (1) Modify the .cshrc file in the user's system home directory

[0062] Add the command to set the locale variable in the .cshrc file in the user's system home directory. For example, for user kt1gf:

[0063] setenv PATH / patch / bin:$PATH

[0064] alias ktgy_opti' / patch / bin / ktgy_opti_x'

[0065] alias ktha_opti' / patch / bin / ktha_opti_x'

[0066] The executable file of this product is stored in the / path / bin directory.

[0067] The contents of the ktgy_opti_x file are:

[0068] #!csh

[0069] setenv DATAPAR / patch / ktgy / datapar /

[0070] setenv THTAB / patch / ktgy / thtab /

[0071] / patch / bin / data_opti_x

[0072] The user starts different commands (ktgy_opti, ktha_opti) and dynamically sets the environment variables: DATAPAR, THTAB to enter the usage environment of the relevant depression work area.

[0073] (2) Create various parameter files in the specified directory

[0074] Except for the time-depth table files which are placed in the directory specified by $THTAB, other basic data files are placed in the directory specified by ($DATAPAR).

[0075] ①2D survey network basic coordinates survey2d.par

[0076] The first line of the file header, each 2D survey area occupies one line, parameter 1 is the horizontal coordinate of the pile origin (columns 1-12), parameter 2 is the vertical coordinate of the pile origin (columns 14-24), parameter 3 is the azimuth of the main survey line (columns 28-34), parameter 4 is the direction of the main survey line (left is 1, right is -1) (columns 37-38), parameter 5 is the name of the 2D survey area (columns 47-56). For example:

[0077]

[0078] ②3D survey network basic coordinates survey3d.par

[0079] Each survey area has 5 lines. The first line is the survey area name. The second line is the line number, CDP, x and y coordinates of the first control point in the survey area. The third line is the line number, CDP, x and y coordinates of the second control point in the survey area. The fourth line is the line number, CDP, x and y coordinates of the third control point in the survey area. The fifth line is the line number, CDP, x and y coordinates of the fourth control point in the survey area. The data are separated by spaces. For example:

[0080] libao3d

[0081]

[0082] ③Time depth table index file thtab.index

[0083] Each time-depth table occupies one line, consisting of the full path and time-depth table name, organizing the time-depth file name.

[0084] ④Time depth table file data format

[0085] Each pair of round-trip time (ms) and depth (m) occupies one line, and the data pairs are separated by spaces.

[0086] The above parameter files are mainly established by the tools in step 5.

[0087] The second step is to preprocess the Omega stacking velocity data

[0088] Supports 2D / 3D Omega velocity spectrum data input, filters outliers in the velocity spectrum through the constraints of the Dix formula, outputs VELF or XYTV velocity spectrum format, and provides auxiliary support for LandMark Depth Team velocity spectrum loading. This involves the conversion of stake numbers into geodetic coordinates.

[0089] The formula for converting from pile number to coordinate is:

[0090] x=x0-co2d*ae*cos(az2d)+an*sin(az2d) (1)

[0091] y=y0+co2d*ae*sin(az2d)+an*cos(az2d) (2)

[0092] Among them: x0, y0 are the coordinates of the origin of the two-dimensional survey network, az2d is the azimuth of the main survey line, an, ae are the stake numbers of the main survey line and the connecting survey line respectively, x, y are the geodetic coordinates, and co2d is the direction of progress of the main survey line.

[0093] The third step is to use the auxiliary Z-Map_Plus drawing software

[0094] The time domain layer data / faults are converted into the depth domain to generate a 2D / 3D survey network with a controllable plane range, and support bidirectional conversion of GeoFrame and Z-Map-Plus fault trajectory data formats.

[0095] Step 4: GeoFrame basic data optimization

[0096] It supports the optimization of well deviation data, profile faults, and plane fault trajectory data, eliminates data redundancy, and can delete invalid or obsolete survey area interpretation layer data. It is particularly effective in sorting out the interpretation of old work areas and can significantly improve the user's work efficiency.

[0097] Due to the use of the well data work area copy function, there may be multiple versions of well deviation data in a well in the GeoFrame work area. The software technology automatically deletes redundant versions with different version names but the same data, and prompts multiple version information for users to choose. The interpretation data merging technology of multiple work areas is used. On the same survey line interpretation section, there are section faults that completely overlap or basically overlap; fault track segments overlap on the same layer. The recognition technology that determines whether there is overlap by calculating the average distance between two data segments automatically deletes redundant data segments.

[0098] Step 5. Tools

[0099] It provides functions such as seismic line coordinate query, time-depth data query, unit conversion, azimuth calculation, GeoFrame 3D survey area parameter extraction, 3D survey area coordinate system and time-depth velocity table management.

[0100] The basic formula for checking the 2D stake number (main survey line an, contact survey line ae) from the coordinates (x, y) is:

[0101] ae=(x-x0)*sin(az2d)+(y-y0)*cos(az2d) (3)

[0102] an=-(x-x0)*cos(az2d)+(y-y0)*sin(az2d) (4)

[0103] In a 3D measurement area, the formula for checking the 3D line number (Linem, Cdpm) from the coordinates (x, y) is:

[0104] LineM=l1+(l2-l1) / (xx4-xx1)*(xt-xx1) (5)

[0105] Cdpm=Cdp1+(Cdp2-Cdp1) / (yy2-yy1)*(yt-yy1) (6)

[0106] Where: xt is the x coordinate of the normal projection of (x,y) on the line connecting points 1 and 4 of the 3D basic survey line coordinates. yt is the y coordinate of the normal projection of (x,y) on the line connecting points 1 and 2 of the 3D basic survey line coordinates. l1, l2 are the main survey line numbers of the frame. Cdp1, Cdp2 are the connecting survey lines of points 1 and 2. xx4, xx1 are the x coordinates of points 4 and 1. yy2, yy1 are the y coordinates of points 2 and 1.

[0107] In a 3D measurement area, the formula for looking up coordinates (x, y) from the 3D line number (Linem, Cdpm) is:

[0108] px1=xx1+(xx4-xx1) / (l2-l1)*(Linem-l1) (8)

[0109] py1=yy1+(yy4-yy1) / (l2-l1)*(Linem-l1) (9)

[0110] px2=xx2+(xx3-xx2) / (l2-l1)*(Linem-l1) (10)

[0111] py2=yy2+(yy3-yy2) / (l2-l1)*(Linem-l1) (11)

[0112] x=px1+(px2-px1) / (cdp2-cdp1)*(Cdpm-Cdp1) (12)

[0113] y=py1+(py2-py1) / (cdp2-cdp1)*(Cdpm-Cdp1) (13)

[0114] Note: The coordinates of the range points of a 3D survey area are controlled by the coordinates of the 4 corner points. The 1st and 2nd points are on the 1st border line. Relative to the position of point 1, point 2 is in the direction of increasing the connecting survey line. The 3rd and 4th points of the border are determined clockwise or counterclockwise.

[0115] Except for the first step which is mandatory at the beginning, the other steps are optional and independent. Once the first step is prepared, it is not necessary to repeat it later unless new parameter changes occur.

[0116] The present invention uses software programming technology to modify the Omega stacking velocity data file, standardize the format, and remove abnormal points, opening up the channel from Omega velocity spectrum data to LandMark Depth Team. The present invention realizes the smooth docking of GeoFrame interpretation software and Z-Map-Plus drawing software. The present invention realizes the automatic optimization of digital oblique data, profile fault, and plane fault trajectory data, significantly improving the work efficiency of the old work area.

[0117] The specific implementation of this application can be written in Fortran77 language, run on Linux platform, and has the following features:

[0118] (1) Preprocess the Omega format stacking velocity spectrum data to make it seamlessly connected with the Depth Team.

[0119] (2) It provides convenient time-depth conversion for GeoFrame software interpretation data. The fault trajectory data formats of GeoFrame and Z-Map-Plus can be converted to each other. At the same time, 3D survey networks are generated according to user needs, which enhances the connection between interpretation application software and drawing software.

[0120] (3) The data overlap segment recognition technology is used to optimize the basic data of the GeoFrame work area, eliminate data redundancy to the maximum extent, and significantly improve user work efficiency.

[0121] (4) The software runs directly in the user's environment, and hundreds of basic parameters are managed in the form of files (most parameters are obtained by the system tool software), making it convenient and fast for users to use.

[0122] Beneficial effects achieved by this application:

[0123] (1) The Omega format stacked velocity spectrum data preprocessed by this system not only changes the pile number into coordinates and is directly loaded in Depth Team, but also removes abnormal velocity spectrum points, saving users a lot of time in editing abnormal velocity spectrum points.

[0124] (2) It provides an organic bridge for data transfer between the GeoFrame interpretation system and the Z-Map-Plus drawing software.

[0125] (3) By batch-exporting the basic data of GeoFrame work areas, using overlapping segment recognition technology to optimize data, and eliminating data redundancy to the maximum extent, the optimized data is then stored back in the interpretation work area, which significantly speeds up the access efficiency of GeoFrame old work area data. For example, the data access efficiency of the North Slope work area is more than 5 times higher, which means that the work efficiency of a single user is increased by 5 times. For a team of 10 people, the work pace is accelerated by more than 50 times.

[0126] Example 1: Preprocessing of Yiyang 2D Omega format velocity spectrum data

[0127] 1. Copy the 2D line Omega format velocity spectrum output to the working directory

[0128] 2. Create a velocity spectrum 2D measurement line list file using the list command

[0129] ls 2015*>lines_2d_stk_list.dat

[0130] 3. Create the survey line coordinate file yy2d_coor.dat through 2D survey line coordinate extraction technology.

[0131] 4. Start the command ktyy_Opti in the current working directory

[0132] Enter the 2D Omega format velocity spectrum data preprocessing module through the 1-2 menu steps.

[0133] ⑴ Input velocity spectrum 2D line list file: lines_2d_stk_list.dat

[0134] ⑵Set whether the velocity spectrum file contains x,y coordinates (1 / 0): 0

[0135] ⑶ Output VELF format 2D velocity profile file: lines_2d_stk_VELF.dat

[0136] Velocity spectrum file format before preprocessing

[0137]

[0138] The format of velocity spectrum file after preprocessing:

[0139]

[0140] Note: The differences before and after the velocity spectrum file preprocessing are: first, the survey line stake number (250, 1972) is replaced with the corresponding geodetic coordinates (20475161, 3118367); second, the abnormal points of the velocity spectrum are removed.

[0141] Example 2: Preprocessing of Huangma high-precision 3D Omega format velocity spectrum data

[0142] Enter the 3D Omega format velocity spectrum data preprocessing module through menu steps 1-3.

[0143] 1. Input Omega format 3D velocity spectrum file: 13HM3D_PSTM_MIG_VEL_GWH_2013.dat

[0144] 2. Output 3D velocity map file: hm_gfj3D_Vel_VELF.dat

[0145] 3. Set the output velocity spectrum format (VELF / XYTV) (1 / -1): -1

[0146] 4. Set 2D pile number measurement network code: 1

[0147] Example 3: Sha'an Dahuazhuang T25 time-depth conversion

[0148] Enter the 3D Omega format velocity spectrum data preprocessing module through the 2-1 menu steps.

[0149] 1. Enter the time domain layer file name: sn_dhz_t25_h.p701

[0150] 2. Set the file format (standard format S, non-standard format N): S

[0151] 3. Output depth domain file name (xyzt…): sn_dhz_t25_h.dat

[0152] 4. Time depth table code number: 2

[0153]

[0154] Example 4: Generation of Libao 3D measurement network data file

[0155] Enter the 3D Omega format velocity spectrum data preprocessing module through menu steps 2-5.

[0156] 1. Select 3D measurement area code: 2

[0157] 2. Set the start, end and interval of the main measurement line: 30 890 10

[0158] 3. Set the start, end and interval of the contact measurement line: 200 1150 10

[0159] 4. Set the shape range: determined by module prompts or defined by the user

[0160] 5. Set edge: determined by module prompt or user defined

[0161] 6. Set two or more points mode (1 / 2): 1

[0162] 7. Output 3 test network file name: libao3d_lines.dat

[0163]

[0164] Example 5: Optimization of Baixiep data in the northern slope area of ​​Gaoyou Sag

[0165] 1. Well deviation data optimization

[0166] Enter the well deviation data optimization submenu through the 3-1 menu step. Well deviation data optimization is divided into two steps: the first is well deviation data format conversion, the second is data optimization.

[0167] (1) Well deviation data format conversion

[0168] Enter the well deviation data format conversion module through the 3-1-1 menu steps and follow the steps below:

[0169] ① Input the format file controlled by am_multi_devsur_order4uwi.ctl: bxp_dev_2016.dat

[0170] ② Output GF-LMK format file: bxp_dev_gf-lmk.dat

[0171]

[0172] (2) Well deviation data optimization

[0173] Enter the well deviation data format conversion module through the 3-1-2 menu steps and follow the steps below:

[0174] ① Input the well deviation data file in GF-LMK format: bxp_dev_gf-lmk.dat

[0175] ②Set version control file: bxp_dev_version.dat

[0176] ③ Output target well deviation data file: bxp_dev_obje.dat

[0177] ④ Output multi-version log files: bxp_dev_log.dat

[0178] Note: This module should be run at least twice. For the first run, the initial version control file can be empty. Before running again, copy the log file bxp_dev_log.dat output for the first time to the version control file bxp_dev_version.dat, and then modify this file to ensure that the version record of each well is unique. The data format of the version control file and its log file is the same. After running again, the final target well deviation data file bxp_dev_obje.dat can be obtained.

[0179] 2. Optimization of profile and fault data

[0180] Enter the section fault data optimization module through the 3-2 menu steps and follow the steps below:

[0181] (1) Enter the fault file name in the format of faultm_gf.ifdf:

[0182] bxp_f_20160512_12may_165004.flt

[0183] (2) Input the associated survey area list file: bxp_surveys_list.dat

[0184] (3) Set the maximum error (10): 10

[0185] (4) Output faultm_gf.ifdf format target fault file name: bxp_f.flt

[0186] (5) Output redundant log file: bxp_f_error.dat

[0187] 3. Plane fault trajectory data optimization

[0188] Enter the section fault data optimization module through the 3-3 menu steps and follow the steps below:

[0189] (1) Input the fault trajectory file name in fbm_gf.ifdf format: bxp_fb.p701

[0190] (2) Set the maximum error (10): 10

[0191] (3) Output the target fault trajectory file in fbm_gf.ifdf format: bxp_fb_opti.p701

[0192] (4) Output redundant log file: bxp_fb_log.dat

[0193] 4. Layer data filtering

[0194] Enter the layer data filtering module through the 3-4 menu steps and follow the steps below:

[0195] (1) Input the layer data file name in 3d_ci7m_gf.ifdf format: bxp_h.p701

[0196] (2) Set the 3D survey area list file name: bxp_h_surveys_list.dat

[0197] (3) Output 3d_ci7m_gf.ifdf format target layer file name: bxp_h_obje.p701

[0198] (4) Output redundant log file: bxp_h_log.dat

[0199] 5. Evaluation of optimization effect

[0200] Comparison table of data file sizes before and after optimization of the North Slope work area

[0201] Serial number Data Objects Before optimization(MB) After optimization (MB) Front-to-back ratio 1 Horizon 1600.7 1291.5 1.24 2 Cross-section plane 88 76.8 1.14 3 Plane fault trace 1666 22 75.7 4 Well deviation trajectory 20.2 7.7 2.62

[0202] According to the ratio of data file sizes before and after optimization, the main redundant data in the work area are: plane fault trajectory, well deviation trajectory, layer position, and profile fault.

[0203] Comparison table of display response time before and after optimization of the interpretation system in the North Slope area

[0204] Serial number Data Objects Before optimization(s) After optimization(s) Front-to-back ratio 1 Section Display Time 60 11 5.4 2 Basemap display time 120 20 6

[0205] The response time ratio before and after optimization indicates that after data optimization, the response speed of the work area profile or base map display has increased by more than 5 times. This basically means that the user's work efficiency has increased by 5 times, and the efficiency is even greater for multi-person work teams.

[0206] like Figure 2 As shown, the present application provides a device for preprocessing and optimizing basic data of a seismic interpretation system, including: an establishment unit 201, used to establish different user usage environments, including establishing various parameter files; a preprocessing unit 202, used to preprocess Omega stacking velocity data; an auxiliary unit 203, used to assist Z-Map_Plus drawing software; an optimization unit 204, used to optimize GeoFrame basic data; a tool unit 205, used to provide data query, unit conversion, azimuth calculation, parameter extraction, and data management functions.

[0207] The pre-processing unit is used for:

[0208] Filter outliers in the velocity spectrum through the constraints of the Dix formula, and output the velocity spectrum in VELF or XYTV format to provide auxiliary support for LandMark Depth Team velocity spectrum loading;

[0209] The formula for converting pile number into geodetic coordinates is:

[0210] x=x0-co2d*ae*cos(az2d)+an*sin(az2d)

[0211] y=y0+co2d*ae*sin(az2d)+an*cos(az2d)

[0212] Among them: x0, y0 are the coordinates of the origin of the two-dimensional survey network, az2d is the azimuth of the main survey line, an, ae are the stake numbers of the main survey line and the connecting survey line respectively, x, y are the geodetic coordinates, and co2d is the direction of progress of the main survey line.

[0213] Among them, the optimization unit is used to:

[0214] The interpretation data merging technology of multiple working areas is adopted. On the same survey line interpretation section, there are section faults that completely overlap or basically overlap; fault trajectory segments overlap on the same layer; and the recognition technology of whether there is overlap is judged by calculating the average distance between two data segments, and redundant data segments are automatically deleted.

[0215] In the present application, the device embodiment for preprocessing and optimizing basic data of a seismic interpretation system is substantially similar to the method embodiment for preprocessing and optimizing basic data of a seismic interpretation system. For relevant matters, please refer to the introduction to the method embodiment for preprocessing and optimizing basic data of a seismic interpretation system.

[0216] Those skilled in the art can clearly understand that the technical solutions of the embodiments of the present invention can be implemented with the help of software and / or hardware. "Unit" and "module" in this specification refer to software and / or hardware that can independently complete or cooperate with other components to complete specific functions, where the hardware can be, for example, FPGA (Field-Programmable Gate Array), IC (Integrated Circuit), etc.

[0217] The embodiment of the present invention also provides a computer-readable storage medium, on which a computer program is stored, and when the program is executed by a processor, the method steps of preprocessing and optimizing the basic data of the seismic interpretation system are implemented. The computer-readable storage medium may include, but is not limited to, any type of disk, including a floppy disk, an optical disk, a DVD, a CD-ROM, a micro drive, and a magneto-optical disk, a ROM, a RAM, an EPROM, an EEPROM, a DRAM, a VRAM, a flash memory device, a magnetic card or an optical card, a nano system (including a molecular memory IC), or any type of medium or device suitable for storing instructions and / or data.

[0218] Figure 3 This is a schematic diagram of the structure of the computer device according to the embodiment of the present application. Figure 3 As shown, the computer device of the present application is, for example, a laptop computer, a desktop computer, a workbench, a personal digital assistant, a server, a blade server, a mainframe computer, and other suitable computers. The computer device of the present application includes a processor 401, a memory 402, an input device 403, and an output device 404. The processor 401, the memory 402, the input device 403, and the output device 404 can be connected via a bus 405 or other means. The memory 402 stores a computer program, which can be run on the processor 401, and when the processor 401 executes the program, the method steps of preprocessing and optimizing the basic data of the seismic interpretation system described above are implemented.

[0219] The input device 403 is, for example, a touch screen, a keypad, a mouse, a track pad, a touchpad, a pointer, one or more mouse buttons, a trackball, a joystick, etc. The output device 404 may include a display device, an auxiliary lighting device (e.g., an LED), and a tactile feedback device (e.g., a vibration motor), etc. The display device may include, but is not limited to, a liquid crystal display (LCD), a light emitting diode (LED) display, a plasma display, and a touch screen.

[0220] In the several embodiments provided in the present application, it should be understood that the disclosed devices and methods can be implemented in other ways. The device embodiments described above are only schematic. For example, the division of the units is only a logical function division. There may be other division methods in actual implementation, such as: multiple units or components can be combined, or can be integrated into another system, or some features can be ignored or not executed. In addition, the coupling, direct coupling, or communication connection between the components shown or discussed can be through some interfaces, indirect coupling or communication connection of devices or units, which can be electrical, mechanical or other forms.

[0221] In the embodiments of the present invention, all functional units may be integrated into one processing unit, or each unit may be separately used as a unit, or two or more units may be integrated into one unit; the above-mentioned integrated units may be implemented in the form of hardware or in the form of hardware plus software functional units.

[0222] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. For those skilled in the art, the present invention may have various modifications and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included in the protection scope of the present invention.

Claims

1. A method for preprocessing and optimizing basic data of a seismic interpretation system, characterized in that: include: Establish different user environments, including establishing various parameter files; Preprocess Omega stacking velocity data; including: Filter outliers in the velocity spectrum through the constraints of the Dix formula, and output the velocity spectrum in VELF or XYTV format to provide auxiliary support for LandMarkDepth Team velocity spectrum loading; The formula for converting pile number into geodetic coordinates is: x=x0-co2d*ae*cos(az2d)+an*sin(az2d) y=y0+co2d*ae*sin(az2d)+an*cos(az2d) Among them: x0, y0 are the coordinates of the origin of the two-dimensional survey network, az2d is the azimuth of the main survey line, an, ae are the stake numbers of the main survey line and the connecting survey line respectively, x, y are the geodetic coordinates, and co2d is the direction of the main survey line; Auxiliary Z-Map_Plus mapping software; including: Convert time domain layer data / faults into depth domain to generate 2D / 3D measurement network with controllable plane range; Optimize GeoFrame basic data; including: Adopting the interpretation data merging technology of multiple working areas, on the same survey line interpretation section, there are section faults that completely overlap or basically overlap; fault track segments overlap on the same layer; by calculating the average distance between two data segments to determine whether they overlap, the redundant data segments are automatically deleted; Provides data query, unit conversion, azimuth calculation, parameter extraction and data management functions.

2. The method for preprocessing and optimizing basic data of a seismic interpretation system according to claim 1, characterized in that: Provides data query, unit conversion, azimuth calculation, parameter extraction, and data management functions, including: Provides seismic line coordinate query, time-depth data query, unit conversion, azimuth calculation, GeoFrame 3D survey area parameter extraction, 3D survey area coordinate system and time-depth velocity table management functions; The basic formula for checking the 2D pile number from the coordinates (x, y) is: ae=(x-x0)*sin(az2d)+(y-y0)*cos(az2d) an=-(x-x0)*cos(az2d)+(y-y0)*sin(az2d) In a 3D measurement area, the formula for checking the 3D line number from the coordinates (x, y) is: LineM=l1+(l2-l1) / (xx4-xx1)*(xt-xx1) Cdpm=Cdp1+(Cdp2-Cdp1) / (yy2-yy1)*(yt-yy1) Wherein: xt is the x coordinate of the normal projection of (x,y) on the line connecting points 1 and 4 of the 3D basic survey line coordinates; yt is the y coordinate of the normal projection of (x,y) on the line connecting points 1 and 2 of the 3D basic survey line coordinates; l1 and l2 are the main survey line numbers of the border; Cdp1 and Cdp2 are the connecting survey lines of points 1 and 2; xx4 and xx1 are the x coordinates of points 4 and 1; yy2 and yy1 are the y coordinates of points 2 and 1.

3. A device for preprocessing and optimizing basic data of a seismic interpretation system, characterized in that: include: The establishment unit is used to establish different user environments, including establishing various parameter files; A preprocessing unit, used for preprocessing Omega stacking velocity data; Filter outliers in the velocity spectrum through the constraints of the Dix formula, and output the velocity spectrum in VELF or XYTV format to provide auxiliary support for LandMarkDepth Team velocity spectrum loading; The formula for converting pile number into geodetic coordinates is: x=x0-co2d*ae*cos(az2d)+an*sin(az2d) y=y0+co2d*ae*sin(az2d)+an*cos(az2d) Among them: x0, y0 are the coordinates of the origin of the two-dimensional survey network, az2d is the azimuth of the main survey line, an, ae are the stake numbers of the main survey line and the connecting survey line respectively, x, y are the geodetic coordinates, and co2d is the direction of the main survey line; Auxiliary unit, used to assist Z-Map_Plus mapping software; Convert time domain layer data / faults into depth domain to generate 2D / 3D measurement network with controllable plane range; Optimization unit, used to optimize GeoFrame basic data; By using the interpretation data merging technology of multiple work areas, on the same survey line interpretation section, there are section faults that completely overlap or basically overlap; fault track segments overlap on the same layer; By calculating the average distance between two data segments, the recognition technology determines whether they overlap and automatically deletes redundant data segments. The tool unit is used to provide data query, unit conversion, azimuth calculation, parameter extraction, and data management functions.

4. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the program is executed by a processor, the steps of the method described in any one of claims 1 to 2 are implemented.

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

Images