CN112711066B - 一种三维地震勘探炮点布设均匀性评价方法和装置 - Google Patents
一种三维地震勘探炮点布设均匀性评价方法和装置 Download PDFInfo
- Publication number
- CN112711066B CN112711066B CN201911025154.5A CN201911025154A CN112711066B CN 112711066 B CN112711066 B CN 112711066B CN 201911025154 A CN201911025154 A CN 201911025154A CN 112711066 B CN112711066 B CN 112711066B
- Authority
- CN
- China
- Prior art keywords
- grid
- shot
- statistical
- shots
- basic unit
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 57
- 238000011077 uniformity evaluation Methods 0.000 title claims description 12
- 238000011156 evaluation Methods 0.000 claims abstract description 61
- 238000001514 detection method Methods 0.000 claims description 29
- 238000010606 normalization Methods 0.000 claims description 11
- 230000008569 process Effects 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 abstract description 3
- 230000006870 function Effects 0.000 description 29
- 238000010276 construction Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 238000010304 firing Methods 0.000 description 6
- 230000006872 improvement Effects 0.000 description 5
- 230000004888 barrier function Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000005422 blasting Methods 0.000 description 2
- 238000012512 characterization method Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 239000007789 gas Substances 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 238000004445 quantitative analysis Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000005096 rolling process Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000009954 braiding Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003203 everyday effect Effects 0.000 description 1
- 238000005755 formation reaction Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Complex Calculations (AREA)
Abstract
本发明属于三维地震勘探技术领域，具体涉及一种三维地震勘探炮点布设均匀性评价方法和装置。本发明在满足地震剖面允许的最大缺口需求下确定统计网格的尺寸，同时依据确定的基础面元网格的尺寸在设定评价区域内布设统计网格，使统计网格覆盖整个设定评价区域，并统计各个统计网格内的炮点数，根据各个统计网格内的炮点数确定每个基础面元网格内的炮点数，并依据每个基础面元网格内的炮点数评价炮点均匀性。该方法及装置按照采样均匀和炮点分布均匀原则，依据实际的炮点布设情况来确定炮点布设均匀性，使得对采用人工变观设计后的系统的炮点布设均匀性评价更加准确。
Description
技术领域
本发明属于三维地震勘探技术领域，具体涉及一种三维地震勘探炮点布设均匀性评价方法和装置。
背景技术
地震勘探是利用地下介质弹性和密度的差异，通过观测和分析大地对人工激发地震波的影响，推断地下岩层的性质和形态的地球物理勘探方法。地震勘探是钻探前勘测石油与天然气资源的重要手段，在煤田和工程地质勘查、区域地质研究和地壳研究等方面，得到广泛应用。
随着国内探区油气勘探程度的不断深入，在复杂近地表和复杂构造地区寻找岩性油气田是目前的主要勘探目标。实际上在每天采区进行三维地震勘探时，由于受地形条件和建筑物的影响，特别是工区内村镇障碍更多、更大，河地表起伏大、山高坡陡地形突变点多，地表接收条件差异较大，不能按照预先设计的方案进行施工，很可能设计出的炮点、检波点位于障碍区内，或者出现不能满足安全激发的要求的情况。故在地震施工中激发点布设较为困难，较易造成点位分布不均匀甚至严重缺失；遇到较大障碍区时，因无法跨越而丢炮将造成局部地震资料空白。
针对该问题，三维地震勘探施工中通常采用变观方法加以弥补。传统变观设计通常都是人工操作，主要考虑覆盖次数、地表障碍情况，通过排列强行跨越障碍物、障碍物两端加大偏移距放炮或采用恢复性放炮的方法，来达到满足覆盖次数的目的。但是，随着地震勘探技术发展和勘探精度要求的提高，人工变观设计方法具有局限性，无法定量分析炮点变更的合理性，准确性较差，不能满足叠前偏移所要求的地震观测数据采样均匀和面元属性一致原则，难以适应实际工况要求。
授权公告号为CN102023310B的中国发明专利公开了一种三维观测系统均匀性分析方法，该方法通过计算覆盖次数分布非均匀系数、炮检距分布非均匀系数、方位角分布非均匀系数和炮检覆盖非均匀系数，从多方面分析给定的三维观测系统的覆盖次数、炮检距、方位角、炮检覆盖分布的均匀性，但该方法对三维观测系统进行评价所采用的数据为针对设计好的系统所测量、计算出的数据，而实际采集施工中变观系统的数据已变，若仍旧采用这些数据来对变观系统进行均匀性判断，便使得对变观系统均匀性(包括炮点布设均匀性)评价不准。
发明内容
本发明提供了一种三维地震勘探炮点布设均匀性评价方法和装置，用以解决现有技术方法对炮点布设均匀性评价不准确的问题。
为解决上述技术问题，本发明的技术方案和有益效果为：
本发明的一种三维地震勘探炮点布设均匀性评价方法，包括如下步骤：
1)获取设定评价区域中的检波线距、炮点距、炮线距和检波点距；并根据检波线距、炮点距、炮线距和检波点距，确定基础面元网格的尺寸，所述基础面元网格用于表征计算炮点数的基础单元，所述基础面元网格的尺寸包括横向尺寸和纵向尺寸；
2)依据地震剖面允许的最大缺口范围，确定统计网格的尺寸；
3)在设定评价区域内布设所述统计网格，使横向相邻的两个统计网格之间不重叠部分的横向距离为基础面元网格的横向尺寸，纵向相邻的两个统计网格之间不重叠部分的纵向距离为基础面元网格的纵向尺寸，直至统计网格覆盖整个设定评价区域；
4)统计各个统计网格内的炮点数，并将统计网格内的炮点数作为该统计网格内各个基础面元网格的炮点数，以得到基础面元网格在该统计网格下的炮点数；
5)将基础面元网格在各个统计网格下的炮点数相加，得到各个基础面元网格内的炮点数；
6)根据每个基础面元网格内的炮点数，对设定评价区域的炮点布设均匀性进行评价。
本发明的一种三维地震勘探炮点布设均匀性评价装置，包括存储器和处理器，所述处理器用于执行存储在存储器中的指令以实现如下方法：
1)获取设定评价区域中的检波线距、炮点距、炮线距和检波点距；并根据检波线距、炮点距、炮线距和检波点距，确定基础面元网格的尺寸，所述基础面元网格用于表征计算炮点数的基础单元，所述基础面元网格的尺寸包括横向尺寸和纵向尺寸；
2)依据地震剖面允许的最大缺口范围，确定统计网格的尺寸；
3)在设定评价区域内布设所述统计网格，使横向相邻的两个统计网格之间不重叠部分的横向距离为基础面元网格的横向尺寸，纵向相邻的两个统计网格之间不重叠部分的纵向距离为基础面元网格的纵向尺寸，直至统计网格覆盖整个设定评价区域；
4)统计各个统计网格内的炮点数，并将统计网格内的炮点数作为该统计网格内各个基础面元网格的炮点数，以得到基础面元网格在该统计网格下的炮点数；
5)将基础面元网格在各个统计网格下的炮点数相加，得到各个基础面元网格内的炮点数；
6)根据每个基础面元网格内的炮点数，对设定评价区域的炮点布设均匀性进行评价。
其有益效果：本发明在满足地震剖面允许的最大缺口需求下确定统计网格的尺寸，同时依据确定的基础面元网格的尺寸在设定评价区域内布设统计网格，使统计网格覆盖整个设定评价区域，并统计各个统计网格内的炮点数，根据各个统计网格内的炮点数确定每个基础面元网格内的炮点数，并依据每个基础面元网格内的炮点数评价炮点均匀性。该方法及装置按照采样均匀和炮点分布均匀原则，依据实际的炮点布设情况来确定炮点布设均匀性，使得对采用人工变观设计后的系统的炮点布设均匀性评价更加准确。而且，由于满足地震剖面允许的最大缺口需求，有利于减少叠前偏移噪声和采集足迹，通过统计相邻炮点间距，有效控制单条炮线上的最大炮距，达到更加均匀控制目的。
作为方法及装置的进一步改进，步骤1)中，所述基础面元网格相邻两边的尺寸分别为：检波线距和炮点距的最大公约数、炮线距和检波点距的最大公约数。
作为方法及装置的进一步改进，步骤5)和步骤6)中，将各个基础面元网格内的炮点数进行归一化处理，以得到各个基础面元网格内归一化后的炮点数；根据各个基础面元网格内归一化后的炮点数，对设定评价区域的炮点布设均匀性进行评价。
作为方法及装置的进一步改进，所述归一化处理为：将各个基础面元网格内的炮点数分别与基础面元网格面积相乘后再除以统计网格面积后得到的值作为各个基础面元网格内归一化后的炮点数。
作为方法及装置的进一步改进，为了方便直观评价设定评价区域中炮点布设均匀性，步骤6)中，得到各个基础面元网格内归一化后的炮点数后，使用曲面图将各个基础面元网格内归一化后的炮点数后进行显示，根据曲面图对设定评价区域的炮点布设均匀性进行评价。
附图说明
图1是本发明的方法实施例中的方法流程图；
图2是本发明的方法实施例中的勘探工区炮点分布图；
图3是本发明的方法实施例中的基础面元网格示意图；
图4是本发明的方法实施例中的统计网格示意图；
图5是本发明的方法实施例中的统计网格横向移动一次示意图；
图6是本发明的方法实施例中的统计网格纵向移动一次示意图；
图7是本发明的方法实施例中的网格炮数立体显示示意图。
具体实施方式
方法实施例1：
该实施例提供了一种三维地震勘探炮点布设均匀性评价方法，下面结合图1进行详细的具体说明。
步骤1，依据三维地震观测系统参数，确定用于表征计算炮点数的基础面元网格尺寸。具体的：
1-1)依据三维地震观测系统参数，获取检波线距RLI、炮点距SPI、炮线距SPLI和检波点距RI、炮线号范围、炮点号范围、每炮的炮点号和炮线号。
1-2)由式(1)和式(2)所示函数确定基础面元网格尺寸，包括基础面元网格的横向尺寸和纵向尺寸：
Gridx＝GCD(RLI，SPI) (1)
Gridy＝GCD(SPLI，RI) (2)
式中，Gridx、Gridy分别表示基础面元网格的横向尺寸、纵向尺寸，GCD()表示最大公约数函数，RLI表示检波线距，SPI表示炮点距，SPLI表示炮线距，RI表示检波点距。
步骤2，依据三维地震观测系统和工区实际炮点分布区域，确定设定评价区域；依据步骤1确定的基础面元网格尺寸，从设定评价区域的原点开始，对设定评价区域进行网格的划分。具体的：
2-1)依据三维地震观测系统和工区实际炮点分布区域，确定设定评价区域。
2-2)依据三维地震观测系统参数，设定评价区域的原点为基础面元网格起始点，并设定第一个基础面元网格的编号。
2-3)在设定评价区域范围内，依据基础面元网格的横向尺寸，沿检波线方向移动，依据基础面元网格的纵向尺寸，沿炮线方向移动，直至基础面元网格满覆盖整个设定评价区域，如图3所示，并对每个基础面元网格进行统一编号为Grid(i，j)，其中，i＝i+1，j＝j+1。
步骤3，依据地震剖面允许的最大缺口范围，确定表征炮点均匀性的统计网格尺寸。具体的：
3-1)依据地震剖面允许的最大缺口范围，确定表征评价区域内炮点均匀性的统计网格尺寸。
3-2)依据三维地震观测系统参数，确定统计网格起始点和第一个统计网格编号。
步骤4，依据步骤1-1中确定的基础面元网格尺寸和步骤3-1中确定的统计网格尺寸，获取第一个统计网格Win(1,1)范围内的炮点数。具体的：
4-1)依据基础面元网格尺寸和统计网格尺寸，确定第一个统计网格Win(1,1)范围内的基础面元网格的编号。
4-2)依据步骤1-1中获取的每炮的炮点号、炮线号和步骤3-1中确定的统计网格尺寸，确定第一个统计网格Win(1,1)范围内的炮点号范围为(SPIA～SPIB)、炮线号范围为(SPLA～SPLB)，并由式(3)所示函数确定第一个统计网格Win(1,1)范围内的炮点号在(SPIA～SPIB)并且炮线号在(SPLA～SPLB)范围内的炮点，获取第一个统计网格Win(1,1)范围内的炮点数NETWin(1,1)。
NSPI∈(SPIA～SPIB)∩NSPL∈(SPLA～SPLB) (3)
式中，NSPL表示每个炮点的炮线号，NSPI表示每个炮点的炮点号，NETWin(1,1)表示统计网格Win(1,1)范围内的炮点数。
4-3)依据步骤4-2确定的第一个统计网格Win(1,1)范围内的炮点数和步骤4-1确定的Win(1,1)范围内的基础面元网格编号，由式(4)所示函数确定Win(1,1)范围内每个基础面元网格的炮点数NETGrid(i1,j1)：
NETGrid(i1,j1)＝NETWin(1,1) (4)
式中，Grid(i1,j1)表示统计网格Win(1,1)范围内每个基础面元网格的编号，i1∈(1～10)且j1∈(1～20)；NETGrid(i1,j1)表示统计网格Win(1,1)范围内每个基础面元网格Grid(i1,j1)的炮点数。
步骤5，依据步骤1-2确定的基础面元网格尺寸和步骤3-1确定的统计网格尺寸，获取Win(1,1)横向移动一次后的统计网格Win(2,1)范围内基础面元网格的炮点数。具体的：
5-1)依据步骤1-2确定的基础面元网格尺寸和步骤3-1确定的统计网格尺寸，确定统计网格Win(2,1)范围内的基础面元网格的编号。
5-2)依据步骤1-1确定的每炮的炮点号、炮线号和步骤3-1确定的统计网格尺寸，确定统计网格Win(2,1)范围内的炮点号范围为(SPIA+1～SPIB+1)、炮线号范围为(SPLA～SPLB)，并由式(5)所示函数确定炮点号在(SPIA+1～SPIB+1)并且炮线号在(SPLA～SPLB)范围内的炮点，获取如图5所示的统计网格Win(2,1)范围内的炮点数NETWin(2,1)；
NSPI∈(SPIA+1～SPIB+1)∩NSPL∈(SPLA～SPLB) (5)
式中，NSPL表示每个炮点的炮线号，NSPI表示每个炮点的炮点号，NETWin(2,1)表示统计网格Win(2,1)范围内的炮点数。
5-3)依据步骤5-2确定的统计网格Win(2,1)范围内的炮点数和步骤5-1确定的Win(2,1)范围内的基础面元网格编号，由式(6)所示函数确定Win(2,1)范围内每个基础面元网格的炮点数NETGrid(i2,j1)：
NETGrid(i2,j1)＝NETWin(2,1) (6)
式中，Grid(i2,j1)表示统计网格Win(2,1)范围内每个基础面元网格的编号，i2∈(2～11)且j1∈(1～20)；NETGrid(i2,j1)表示统计网格Win(2,1)范围内每个基础面元网格Grid(i2,j1)的炮点数。
步骤6，依据步骤1-2确定的基础面元网格尺寸和步骤3-1确定的统计网格尺寸，获取Win(1,1)横向移动M次后的统计网格Win(M+1,1)范围内基础面元网格的炮点数。具体的：
6-1)依据步骤1-2确定的基础面元网格尺寸和步骤3-1确定的统计网格尺寸，确定统计网格Win(M+1,1)范围内的基础面元网格的编号。
6-2)依据步骤1-1确定的每炮的炮点号、炮线号和步骤3-1确定的统计网格尺寸，确定统计网格Win(M+1,1)范围内的炮点号范围为(SPIA+M～SPIB+M)、炮线号范围为(SPLA～SPLB)，并由式(7)所示函数确定炮点号在(SPIA+M～SPIB+M)并且炮线号在(SPLA～SPLB)范围内的炮点，获取如图5所示的统计网格Win(M+1,1)范围内的炮点数NETWin(M+1,1)；
NSPI∈(SPIA+M～SPIB+M)∩NSPL∈(SPLA～SPLB) (7)
式中，NSPL表示每个炮点的炮线号，NSPI表示每个炮点的炮点号，NETWin(M+1,1)表示统计网格Win(M+1,1)范围内的炮点数。
6-3)依据步骤6-2确定的统计网格Win(M+1,1)范围内的炮点数和步骤6-1确定的Win(M+1,1)范围内的基础面元网格编号，由式(8)所示函数确定Win(M+1,1)范围内每个基础面元网格的炮点数NETGrid(iM+1,j1)：
NETGrid(iM+1,j1)＝NETWin(M+1,1) (8)
式中，Grid(iM+1,j1)表示统计网格Win(M+1,1)范围内每个基础面元网格的编号，iM+1∈(M+1～M+10)且j1∈(1～20)；NETGrid(iM+1,j1)表示统计网格Win(M+1,1)范围内每个基础面元网格Grid(iM+1,j1)的炮点数。
步骤7，依据步骤1-2确定的基础面元网格尺寸和步骤3-1确定的统计网格尺寸，获取Win(1,1)纵向移动一次后的统计网格Win(1,2)范围内基础面元网格的炮点数。具体的：
7-1)依据步骤1-2确定的基础面元网格尺寸和步骤3-1确定的统计网格尺寸，确定统计网格Win(1,2)范围内的基础面元网格的编号。
7-2)依据步骤1-1确定的每炮的炮点号、炮线号和步骤3-1确定的统计网格尺寸，确定统计网格Win(1,2)范围内的炮点号范围为(SPIA～SPIB)、炮线号范围为(SPLA+1～SPLB+1)，并由式(9)所示函数确定炮点号在(SPIA～SPIB)并且炮线号在(SPLA+1～SPLB+1)范围内的炮点，获取如图5所示的统计网格Win(1,2)范围内的炮点数NETWin(1,2)；
NSPI∈(SPIA～SPIB)∩NSPL∈(SPLA+1～SPLB+1) (9)
式中，NSPL表示每个炮点的炮线号，NSPI表示每个炮点的炮点号，NETWin(1,2)表示统计网格Win(1,2)范围内的炮点数。
7-3)依据步骤7-2确定的统计网格Win(1,2)范围内的炮点数和步骤7-1确定的Win(1,2)范围内的基础面元网格编号，由式(10)所示函数确定Win(1,2)范围内每个基础面元网格的炮点数NETGrid(i1,j2)：
NETGrid(i1,j2)＝NETWin(1,2) (10)
式中，Grid(i1,j2)表示统计网格Win(1,2)范围内每个基础面元网格的编号，i1∈(1～10)且j2∈(2～21)；NETGrid(i1,j2)表示统计网格Win(1,2)范围内每个基础面元网格Grid(i1,j2)的炮点数。
步骤8，依据步骤1-2确定的基础面元网格尺寸和步骤3-1确定的统计网格尺寸，获取Win(1,2)横向移动M次后的统计网格Win(M+1,2)范围内基础面元网格的炮点数。具体的：
8-1)依据步骤1-2确定的基础面元网格尺寸和步骤3-1确定的统计网格尺寸，确定统计网格Win(M+1,2)范围内的基础面元网格的编号。
8-2)依据步骤1-1确定的每炮的炮点号、炮线号和步骤3-1确定的统计网格尺寸，确定统计网格Win(M+1,2)范围内的炮点号范围为(SPIA+M～SPIB+M)、炮线号范围为(SPLA+1～SPLB+1)，并由式(11)所示函数确定炮点号在(SPIA+M～SPIB+M)并且炮线号在(SPLA+1～SPLB+1)范围内的炮点，获取如图5所示的统计网格Win(M+1,2)范围内的炮点数NETWin(M+1,2)；
NSPI∈(SPIA+M～SPIB+M)∩NSPL∈(SPLA+1～SPLB+1) (11)
式中，NSPL表示每个炮点的炮线号，NSPI表示每个炮点的炮点号，NETWin(1+M,2)表示统计网格Win(M+1,2)范围内的炮点数。
8-3)依据步骤8-2确定的统计网格Win(M+1,2)范围内的炮点数和步骤8-1确定的Win(M+1,2)范围内的基础面元网格编号，由式(12)所示函数确定Win(M+1,2)范围内每个基础面元网格的炮点数NETGrid(iM+1,,j2)：
NETGrid(iM+1,,j2)＝NETWin(M+1,2) (12)
式中，Grid(iM+1,j2)表示统计网格Win(M+1,2)范围内每个基础面元网格的编号，i1∈(M+1～M+10)且j2∈(2～21)；NETGrid(iM+1,j2)表示统计网格Win(M+1,2)范围内每个基础面元网格Grid(iM+1,j2)的炮点数。
步骤9，依据步骤1-2确定的基础面元网格尺寸和步骤3-1确定的统计网格尺寸，获取统计网格Win(1,1)横向滚动M次，纵向滚动N次后的Win(M+1,N+1)范围内基础面元网格的炮点数，依据基础面元网格尺设制定滚动间距可以消除因为不同统计网格尺寸在滚动覆盖时引起的最高覆盖次数和最低覆盖次数的均差。具体的：
9-1)依据步骤1-2确定的基础面元网格尺寸和步骤3-1确定的统计网格尺寸，确定统计网格Win(M+1,N+1)范围内的基础面元网格的编号。
9-2)依据步骤1-1确定的每炮的炮点号、炮线号和步骤3-1确定的统计网格尺寸，确定统计网格Win(M+1,N+1)范围内的炮点号范围为(SPIA+M～SPIB+M)、炮线号范围为(SPLA+N～SPLB+N)，并由式(13)所示函数确定炮点号在(SPIA+M～SPIB+M)并且炮线号在(SPLA+N～SPLB+N)范围内的炮点，获取如图5所示的统计网格Win(M+1,N+1)范围内的炮点数NETWin(M+1,N+1)；
NSPI∈(SPIA+M～SPIB+M)∩NSPL∈(SPLA+N～SPLB+N) (13)
式中，NSPL表示每个炮点的炮线号，NSPI表示每个炮点的炮点号，NETWin(M+1,N+1)表示统计网格Win(M+1,N+1)范围内的炮点数。
9-3)依据步骤9-2确定的统计网格Win(M+1,N+1)范围内的炮点数和步骤9-1确定的Win(M+1,N+1)范围内的基础面元网格编号，由式(14)所示函数确定Win(M+1,N+1)范围内每个基础面元网格的炮点数NETGrid(iM+1,jN+1)：
NETGrid(iM+1,jN+1)＝NETWin(M+1,N+1) (14)
式中，Grid(iM+1,jN+1)表示统计网格Win(M+1,N+1)范围内每个基础面元网格的编号，i1∈(M+1～M+10)且j2∈(N+1～N+20)；NETGrid(iM+1,jN+1)表示统计网格Win(M+1,N+1)范围内每个基础面元网格Grid(iM+1,jN+1)的炮点数。
需说明的是，上述步骤4至步骤9的目的是为了在设定评价区域内布设统计网格，并使统计网格覆盖整个设定评价区域。
步骤10，按照上述步骤的原理不停重复，直至遍历完整个设定评价区域，获取得到基础面元网格在各个统计网格下的炮点数。将基础面元网格在各个统计网格下的炮点数相加，作为最终各个基础面元网格内的炮点数NETGrid(i,j)。例如，处于设定评价区域四个角的基础面元网格，即Grid(i1,j1)、Grid(iM+1,j1)、Grid(i1，jN+1)、Grid(iM+1，jN+1)，只有一次赋值，即Grid(1,1)只在统计网格Win(1,1)下对其进行了赋值，Grid(iM+1,j1)只在统计网格Win(M+1,1)下对其进行了赋值，Grid(i1，jN+1)只在统计网格Win(1,N+1)下对其进行了赋值，Grid(iM+1，jN+1)只在统计网格Win(M+1,N+1)下对其进行了赋值。而处于设定评价区域其他位置的基础面元网格，则可能在不同统计网格下对其进行了赋值。例如，Grid(i2,j1)在统计网格Win(1,1)下对其进行了赋值，且在统计网格Win(2,1)下也对其进行了赋值，最终基础面元网格Grid(i2,j1)的炮点数为两者的相加。再由式(15)所示函数确定每个基础面元网格在统计网格移动满覆盖工区范围后的表征设定评价区域内炮点分布均匀性的网格炮点归一化计数值：
A(i,j)＝NETGrid(i,j)*Gridx*Gridy/(Winx*Winy) (15)
式中，A(i,j)表示设定评价区域内的面元网格编号为(i,j)的炮点均匀性统计值；Winx表示统计网格的横向尺寸，Winy表示统计网格的纵向尺寸，Gridy表示基础面元网格的纵向尺寸；Gridx表示基础面元网格的横向尺寸。
步骤11，依据步骤10得到的每个基础面元网格的网格点炮归一化计数值，获取表征设定评价区域内炮点分布均匀性的曲面图，该曲面图的横纵坐标(X、Y坐标)分别为基础面元网格的横纵坐标，Z坐标为每个基础面元网格的网格点炮归一化计数值。通过该图可对三维地震观测系统变观的定量化分析，直观的看出空炮段范围，方便对不能满足施工需求的网格直接进行补炮或调整，以规范野外炮点布设，提高炮点调整的准确性。
下面以某地区为例，对该方法进行进一步描述。
步骤1，依据三维地震观测系统参数，确定用于表征计算炮点数的基础面元网格尺寸。具体的：
1-1)依据图2所示的三维地震观测系统，获取检波线距RLI＝200m、炮点距SPI＝50m、炮线距SPLI＝400m和检波点距RI＝25m，炮线号范围1001-4201，递增量为1、炮点号范围122-642，递增量为1、每炮的炮点号NSP和炮线号NSPL。
1-2)由式(1)和式(2)所示函数确定基础面元网格尺寸，包括基础面元网格的横向尺寸Gridx＝50m和纵向尺寸Gridy＝25m。
步骤2，依据三维地震观测系统和工区实际炮点分布区域，将图2所示的观测系统确定为设定评价区域；依据步骤1确定的基础面元网格尺寸，从设定评价区域的原点开始，对设定评价区域进行网格的划分及基础面元网格的全面覆盖。具体的：
2-1)依据三维地震观测系统和工区实际炮点分布区域，确定设定评价区域的范围为13000m×35000m＝455km2。
2-2)依据三维地震观测系统参数，设定评价区域的原点为基础面元网格起始点，并设定第一个基础面元网格的编号为Grid(1,1)。
2-3)在设定评价区域范围内，依据基础面元网格的横向尺寸50m，沿检波线方向逐次移动，依据基础面元网格的纵向尺寸25m，沿炮线方向移动，直至基础面元网格满覆盖整个设定评价区域，如图3所示，并对每个基础面元网格进行统一编号为Grid(i，j)，其中，i＝i+1，j＝j+1。
步骤3，依据如图2所示的地震剖面允许的最大缺口范围，确定表征炮点均匀性的统计网格尺寸。具体的：
3-1)依据地震剖面允许的最大缺口范围500m，确定表征评价区域内炮点均匀性的统计网格尺寸500m×500m。
3-2)依据三维地震观测系统参数，确定统计网格起始点和第一个统计网格编号。
步骤4，依据步骤1-1中确定的基础面元网格尺寸和步骤3-1中确定的统计网格尺寸，获取统计网格Win(1,1)范围内的炮点数。具体的：
4-1)依据基础面元网格尺寸25m×50m和统计网格尺寸500m×500m，确定统计网格Win(1,1)范围内的基础面元网格的编号，如图4所示，包括Grid(1,1)、Grid(2,1)、Grid(3,1)、Grid(4,1)、……、Grid(10,1)；Grid(1,2)、Grid(2,2)、Grid(3,2)、Grid(4,2)、……、Grid(10,2)；……；Grid(1,20)、Grid(2,20)、Grid(3,20)、Grid(4,20)、……、Grid(10,20)，共计200个。
4-2)依据步骤1-1中获取的每炮的炮点号、炮线号和步骤3-1中确定的统计网格尺寸500m×500m，确定统计网格Win(1,1)范围内的炮点号范围为(122～131)、炮线号范围为(1001～1020)，并由式(3)所示函数确定统计网格Win(1,1)范围内的炮点号在(122～131)并且炮线号在(1001～1020)范围内的炮点，获取如图4所示统计网格Win(1,1)范围内的炮点数NETWin(1,1)＝15。
4-3)依据步骤4-2确定的统计网格Win(1,1)范围内的炮点数和步骤4-1确定的Win(1,1)范围内的基础面元网格编号，由式(4)所示函数确定Win(1,1)范围内每个基础面元网格的炮点数NETGrid(i1,j1)＝15。
步骤5，依据步骤1-2确定的基础面元网格尺寸和步骤3-1确定的统计网格尺寸，获取Win(1,1)横向移动一次后的统计网格Win(2,1)范围内基础面元网格的炮点数。具体的：
5-1)依据步骤1-2确定的基础面元网格尺寸25m×50m和步骤3-1确定的统计网格尺寸500m×500m，确定统计网格Win(2,1)范围内的基础面元网格的编号，如图5所示，包括Grid(2,1)、Grid(3,1)、Grid(4,1)、Grid(5,1)、……、Grid(11,1)；Grid(2,2)、Grid(3,2)、Grid(4,2)、Grid(5,2)、……、Grid(11,2)；……；Grid(2,20)、Grid(3,20)、Grid(4,20)、Grid(5,20)、……、Grid(11,20)，共计200个。
5-2)依据步骤1-1确定的每炮的炮点号、炮线号和步骤3-1确定的统计网格尺寸，确定统计网格Win(2,1)范围内的炮点号范围为(123～132)、炮线号范围为(1001～1020)，并由式(5)所示函数确定炮点号在(123～132)并且炮线号在(1001～1020)范围内的炮点，获取如图5所示的统计网格Win(2,1)范围内的炮点数NETWin(2,1)＝17。
5-3)依据步骤5-2确定的统计网格Win(2,1)范围内的炮点数和步骤5-1确定的Win(2,1)范围内的基础面元网格编号，由式(6)所示函数确定Win(2,1)范围内每个基础面元网格的炮点数NETGrid(i2,j1)＝17。
步骤6，依据步骤1-2确定的基础面元网格尺寸和步骤3-1确定的统计网格尺寸，获取Win(1,1)横向移动M次后的统计网格Win(M+1,1)范围内基础面元网格的炮点数。具体的：
6-1)依据步骤1-2确定的基础面元网格尺寸25m×50m和步骤3-1确定的统计网格尺寸500m×500m，确定统计网格Win(M+1,1)范围内的基础面元网格的编号，包括Grid(M+1,1)、Grid(M+2,1)、Grid(M+3,1)、Grid(M+4,1)、……、Grid(M+10,1)；Grid(M+1,2)、Grid(M+2,2)、Grid(M+3,2)、Grid(M+4,2)、……、Grid(M+10,2)；……；Grid(M+1,20)、Grid(M+2,20)、Grid(M+3,20)、Grid(M+4,20)、……、Grid(M+10,20)共计200个。
6-2)依据步骤1-1确定的每炮的炮点号、炮线号和步骤3-1确定的统计网格尺寸500m×500m，确定统计网格Win(M+1,1)范围内的炮点号范围为(122+M～131+M)、炮线号范围为(1001～1020)，并由式(7)所示函数确定炮点号在(122+M～M+131)并且炮线号在(1001～1020)范围内的炮点，获取如图5所示的统计网格Win(M+1,1)范围内的炮点数NETWin(M+1,1)。
6-3)依据步骤6-2确定的统计网格Win(M+1,1)范围内的炮点数和步骤6-1确定的Win(M+1,1)范围内的基础面元网格编号，由式(8)所示函数确定Win(M+1,1)范围内每个基础面元网格的炮点数NETGrid(iM+1,j1)。
步骤7，依据步骤1-2确定的基础面元网格尺寸和步骤3-1确定的统计网格尺寸，获取Win(1,1)纵向移动一次后的统计网格Win(1,2)范围内基础面元网格的炮点数。具体的：
7-1)依据步骤1-2确定的基础面元网格尺寸25m×50m和步骤3-1确定的统计网格尺寸500m×500m，确定统计网格Win(1,2)范围内的基础面元网格的编号，包括Grid(1,2)、Grid(2,2)、Grid(3,2)、Grid(4,2)、……、Grid(10,2)；Grid(1,3)、Grid(2,3)、Grid(3,3)、Grid(4,3)、……、Grid(10,3)；……；Grid(1,21)、Grid(2,21)、Grid(3,21)、Grid(4,21)、……、Grid(10,21)，共计200个。
7-2)依据步骤1-1确定的每炮的炮点号、炮线号和步骤3-1确定的统计网格尺寸500m×500m，确定统计网格Win(1,2)范围内的炮点号范围为(122～131)、炮线号范围为(1002～1021)，并由式(9)所示函数确定炮点号在(122～131)并且炮线号在(1002～1021)范围内的炮点，获取如图6所示的统计网格Win(1,2)范围内的炮点数NETWin(1,2)＝15。
7-3)依据步骤7-2确定的统计网格Win(1,2)范围内的炮点数和步骤7-1确定的Win(1,2)范围内的基础面元网格编号，由式(10)所示函数确定Win(1,2)范围内每个基础面元网格的炮点数NETGrid(i1,j2)＝15。
步骤8，依据步骤1-2确定的基础面元网格尺寸和步骤3-1确定的统计网格尺寸，获取Win(1,2)横向移动M次后的统计网格Win(M+1,2)范围内基础面元网格的炮点数。具体的：
8-1)依据步骤1-2确定的基础面元网格尺寸25m×50m和步骤3-1确定的统计网格尺寸500m×500m，确定统计网格Win(M+1,2)范围内的基础面元网格的编号，包括Grid(M+1,2)、Grid(M+2,2)、Grid(M+3,2)、Grid(M+4,2)、……、Grid(M+10,2)；Grid(M+1,3)、Grid(M+2,3)、Grid(M+3,3)、Grid(M+4,3)、……、Grid(M+10,3)；……；Grid(M+1,21)、Grid(M+2,21)、Grid(M+3,21)、Grid(M+4,21)、……、Grid(M+10,21)，共计200个。
8-2)依据步骤1-1确定的每炮的炮点号、炮线号和步骤3-1确定的统计网格尺寸，确定统计网格Win(M+1,2)范围内的炮点号范围为(122+M～131+M)、炮线号范围为(1002～1021)，并由式(11)所示函数确定炮点号在(122+M～131+M)并且炮线号在(1002～1021)范围内的炮点，获取如图5所示的统计网格Win(M+1,2)范围内的炮点数NETWin(M+1,2)。
8-3)依据步骤8-2确定的统计网格Win(M+1,2)范围内的炮点数和步骤8-1确定的Win(M+1,2)范围内的基础面元网格编号，由式(12)所示函数确定Win(M+1,2)范围内每个基础面元网格的炮点数NETGrid(iM+1,j2)。
步骤9，依据步骤1-2确定的基础面元网格尺寸和步骤3-1确定的统计网格尺寸，获取统计网格Win(1,1)横向滚动M次，纵向滚动N次后的Win(M+1,N+1)范围内基础面元网格的炮点数，依据基础面元网格尺设制定滚动间距可以消除因为不同统计网格尺寸在滚动覆盖时引起的最高覆盖次数和最低覆盖次数的均差。具体的：
9-1)依据步骤1-2确定的基础面元网格尺寸25m×50m和步骤3-1确定的统计网格尺寸500m×500m，确定统计网格Win(M+1,N+1)范围内的基础面元网格的编号，包括Grid(M+1,N+1)、Grid(M+2,N+1)、Grid(M+3,N+1)、Grid(M+4,N+1)、……、Grid(M+10,N+2)；Grid(M+1,N+2)、Grid(M+2,N+2)、Grid(M+3,N+2)、Grid(M+4,N+2)、……、Grid(M+10,N+2)；……；Grid(M+1,N+19)、Grid(M+2,N+20)、Grid(M+3,N+20)、Grid(M+4,N+20)、……、Grid(M+10,N+20)，共计200个。
9-2)依据步骤1-1确定的每炮的炮点号、炮线号和步骤3-1确定的统计网格尺寸，确定统计网格Win(M+1,N+1)范围内的炮点号范围为(122+M～131+M)、炮线号范围为(1001+N～1020+N)，并由式(13)所示函数确定炮点号在(122+M～131+M)并且炮线号在(1001+N～1020+N)范围内的炮点，获取如图5所示的统计网格Win(M+1,N+1)范围内的炮点数NETWin(M+1,N+1)。
9-3)依据步骤9-2确定的统计网格Win(M+1,N+1)范围内的炮点数和步骤9-1确定的Win(M+1,N+1)范围内的基础面元网格编号，由式(14)所示函数确定Win(M+1,N+1)范围内每个基础面元网格的炮点数NETGrid(iM+1,jN+1)。
步骤10，按照上述步骤的原理不停重复，直至遍历完整个设定评价区域，获取得到基础面元网格在各个统计网格下的炮点数。将基础面元网格在各个统计网格下的炮点数相加，作为最终各个每个基础面元网格的炮点数，再由式(15)所示函数确定每个基础面元网格在统计网格移动满覆盖工区范围后的表征设定评价区域内炮点分布均匀性的网格炮点归一化计数值A(i,j)。
步骤11，依据步骤10得到的每个基础面元网格的网格炮点归一化计数值，获取表征设定评价区域内炮点分布均匀性的曲面图，如图7所示，该曲面图的横纵坐标(X、Y坐标)分别为基础面元网格的横纵坐标，Z坐标为每个基础面元网格的网格点炮归一化计数值。通过该图可对三维地震观测系统变观的定量化分析，直观的看出空炮段范围，方便对不能满足施工需求的网格直接进行补炮或调整，以规范野外炮点布设，提高炮点调整的准确性。
方法实施例2：
在方法实施例1中，在移动统计网格时，是按照一行结束后再移动另一行的顺序进行的，以使统计网格覆盖整个设定评价区域。在该实施例中，与方法实施例1的区别仅在于，不再按照该顺序移动统计网格，而是按照依次计算Win(1,1)、Win(1,2)、……、Win(1,N)、Win(2,1)、Win(2,2)、……、Win(2,N)、……、Win(M,1)、Win(M,2)、……、Win(M,N)的顺序，来使统计网格覆盖整个设定评价区域。
方法实施例3：
在方法实施例1的步骤10中，将每个基础面元网格内的炮点数进行归一化处理，并根据归一化后的基础面元网格内的炮点数对设定评价区域的炮点布设均匀性进行评价。在该实施例中，与方法实施例1的区别仅在于，不再进行归一化处理，根据实施例1步骤9中得到的每个基础面元的炮点数直接进行绘制曲面图，并根据曲面图来观察设定评价区域的炮点布设均匀性。
方法实施例4：
在方法实施例1中，基础面元网格相邻两边的尺寸分别为：检波线距和炮点距的最大公约数、炮线距和检波点距的最大公约数。在该实施例中，与方法实施例1的区别仅在于，基础面元网格相邻两边的尺寸分别为：检波线距和炮点距的最大公约数的2倍(甚至N倍，N＞2)、炮线距和检波点距的最大公约数的2倍(甚至N倍，N＞2)。但该实施例得到曲面图与方法实施例1中得到的曲面图相比，不够密集，评价不够精细。
方法实施例5：
在方法实施例1中，步骤11通过曲面图将步骤10得到的归一化数值进行显示，以更加直观地观测设定评价区域炮点布设均匀性。在该实施例中，与实施例1的区别仅在不再执行步骤11，根据步骤10得到的归一化数值直接对设定评价区域炮点布设均匀性进行评价，或者将步骤10得到的归一化数值存储在表格中，依据表格对设定评价区域炮点布设均匀性进行评价。
装置实施例1：
该实施例提供了一种三维地震勘探炮点布设均匀性评价装置，该装置包括存储器和处理器，该装置包括存储器和处理器，存储器和处理器之间直接或间接地电性连接以实现数据的传输或交互。这里的处理器可以是通用处理器，例如中央处理器CPU，也可以是其他可编程逻辑器件，例如数字信号处理器DSP，处理器用于执行存储在存储器中的指令以实现方法实施例1～5中介绍的一种三维地震勘探炮点布设均匀性评价方法，由于方法实施例1～5已对该方法做了详细说明，这里不再赘述。
尽管本发明的内容已经通过上述优选实施例作了详细介绍，但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后，对于本发明的多种修改和替代都将是显而易见的。因此，本发明的保护范围应由所附的权利要求来限定。
Claims (10)
1.一种三维地震勘探炮点布设均匀性评价方法，其特征在于，包括如下步骤：
1）获取设定评价区域中的检波线距、炮点距、炮线距和检波点距；并根据检波线距、炮点距、炮线距和检波点距，确定基础面元网格的尺寸，所述基础面元网格用于表征计算炮点数的基础单元，所述基础面元网格的尺寸包括横向尺寸和纵向尺寸；
2）依据地震剖面允许的最大缺口范围，确定统计网格的尺寸；
3）在设定评价区域内布设所述统计网格，使横向相邻的两个统计网格之间不重叠部分的横向距离为基础面元网格的横向尺寸，纵向相邻的两个统计网格之间不重叠部分的纵向距离为基础面元网格的纵向尺寸，直至统计网格覆盖整个设定评价区域；
4）统计各个统计网格内的炮点数，并将统计网格内的炮点数作为该统计网格内各个基础面元网格的炮点数，以得到基础面元网格在该统计网格下的炮点数；
5）将基础面元网格在各个统计网格下的炮点数相加，得到各个基础面元网格内的炮点数；
6）根据每个基础面元网格内的炮点数，对设定评价区域的炮点布设均匀性进行评价。
2.根据权利要求1所述的三维地震勘探炮点布设均匀性评价方法，其特征在于，
步骤1）中，所述基础面元网格相邻两边的尺寸分别为：检波线距和炮点距的最大公约数、炮线距和检波点距的最大公约数。
3.根据权利要求1所述的三维地震勘探炮点布设均匀性评价方法，其特征在于，步骤5）中，将各个基础面元网格内的炮点数进行归一化处理，以得到各个基础面元网格内归一化后的炮点数；步骤6）中，根据各个基础面元网格内归一化后的炮点数，对设定评价区域的炮点布设均匀性进行评价。
4.根据权利要求3所述的三维地震勘探炮点布设均匀性评价方法，其特征在于，所述归一化处理为：将各个基础面元网格内的炮点数分别与基础面元网格面积相乘后再除以统计网格面积后得到的值作为各个基础面元网格内归一化后的炮点数。
5.根据权利要求3所述的三维地震勘探炮点布设均匀性评价方法，其特征在于，步骤6）中，得到各个基础面元网格内归一化后的炮点数后，使用曲面图将各个基础面元网格内归一化后的炮点数进行显示，根据曲面图对设定评价区域的炮点布设均匀性进行评价。
6.一种三维地震勘探炮点布设均匀性评价装置，其特征在于，包括存储器和处理器，所述处理器用于执行存储在存储器中的指令以实现如下方法：
1）获取设定评价区域中的检波线距、炮点距、炮线距和检波点距；并根据检波线距、炮点距、炮线距和检波点距，确定基础面元网格的尺寸，所述基础面元网格用于表征计算炮点数的基础单元，所述基础面元网格的尺寸包括横向尺寸和纵向尺寸；
2）依据地震剖面允许的最大缺口范围，确定统计网格的尺寸；
3）在设定评价区域内布设所述统计网格，使横向相邻的两个统计网格之间不重叠部分的横向距离为基础面元网格的横向尺寸，纵向相邻的两个统计网格之间不重叠部分的纵向距离为基础面元网格的纵向尺寸，直至统计网格覆盖整个设定评价区域；
4）统计各个统计网格内的炮点数，并将统计网格内的炮点数作为该统计网格内各个基础面元网格的炮点数，以得到基础面元网格在该统计网格下的炮点数；
5）将基础面元网格在各个统计网格下的炮点数相加，得到各个基础面元网格内的炮点数；
6）根据每个基础面元网格内的炮点数，对设定评价区域的炮点布设均匀性进行评价。
7.根据权利要求6所述的三维地震勘探炮点布设均匀性评价装置，其特征在于，步骤1）中，所述基础面元网格相邻两边的尺寸分别为：检波线距和炮点距的最大公约数、炮线距和检波点距的最大公约数。
8.根据权利要求6所述的三维地震勘探炮点布设均匀性评价装置，其特征在于，步骤5）中，将各个基础面元网格内的炮点数进行归一化处理，以得到各个基础面元网格内归一化后的炮点数；步骤6）中，根据各个基础面元网格内归一化后的炮点数，对设定评价区域的炮点布设均匀性进行评价。
9.根据权利要求8所述的三维地震勘探炮点布设均匀性评价装置，其特征在于，所述归一化处理为：将各个基础面元网格内的炮点数分别与基础面元网格面积相乘后再除以统计网格面积后得到的值作为各个基础面元网格内归一化后的炮点数。
10.根据权利要求8所述的三维地震勘探炮点布设均匀性评价装置，其特征在于，步骤6）中，得到各个基础面元网格内归一化后的炮点数后，使用曲面图将各个基础面元网格内归一化后的炮点数进行显示，根据曲面图对设定评价区域的炮点布设均匀性进行评价。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911025154.5A CN112711066B (zh) | 2019-10-25 | 2019-10-25 | 一种三维地震勘探炮点布设均匀性评价方法和装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911025154.5A CN112711066B (zh) | 2019-10-25 | 2019-10-25 | 一种三维地震勘探炮点布设均匀性评价方法和装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112711066A CN112711066A (zh) | 2021-04-27 |
| CN112711066B true CN112711066B (zh) | 2024-02-20 |
Family
ID=75540854
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911025154.5A Active CN112711066B (zh) | 2019-10-25 | 2019-10-25 | 一种三维地震勘探炮点布设均匀性评价方法和装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112711066B (zh) |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2308993A1 (en) * | 1999-05-25 | 2000-11-25 | Baker Hughes Incorporated | Method of determining merit of offset and azimuth distributions in seismic data acquisition |
| CN102023310A (zh) * | 2009-09-19 | 2011-04-20 | 中国石油化工集团公司 | 三维观测系统均匀性分析方法 |
| CN103576213A (zh) * | 2012-07-27 | 2014-02-12 | 中国石油天然气集团公司 | 一种炮检距属性均匀性的确定方法 |
| CN103675893A (zh) * | 2012-09-04 | 2014-03-26 | 中国石油天然气集团公司 | 海上拖缆观测系统复杂模板满覆盖布设方法 |
| CN104142514A (zh) * | 2013-10-29 | 2014-11-12 | 中国石油化工股份有限公司 | 一种三维地震观测系统定量设计方法 |
| CN104142518A (zh) * | 2013-10-29 | 2014-11-12 | 中国石油化工股份有限公司 | 一种地震观测系统叠前时间偏移响应的分析方法 |
| CN106199706A (zh) * | 2015-06-01 | 2016-12-07 | 中国石油化工股份有限公司 | 三维观测系统面元属性统计方法及装置 |
| CN107561574A (zh) * | 2017-07-20 | 2018-01-09 | 中国石油天然气集团公司 | 炮点位置、接收点位置的检测方法和装置 |
| CN109298450A (zh) * | 2018-11-07 | 2019-02-01 | 中国石油天然气股份有限公司 | 地震采集观测系统属性分析方法及装置 |
| CN109655879A (zh) * | 2017-10-11 | 2019-04-19 | 中国石油化工股份有限公司 | 目的层信噪比一致性能量照明的观测系统优化方法及装置 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10386511B2 (en) * | 2014-10-03 | 2019-08-20 | Exxonmobil Upstream Research Company | Seismic survey design using full wavefield inversion |
| WO2016108896A1 (en) * | 2014-12-31 | 2016-07-07 | Landmark Graphics Corporation | Seismic elastic wave simulation for tilted transversely isotropic media using adaptive lebedev staggered grid |
-
2019
- 2019-10-25 CN CN201911025154.5A patent/CN112711066B/zh active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2308993A1 (en) * | 1999-05-25 | 2000-11-25 | Baker Hughes Incorporated | Method of determining merit of offset and azimuth distributions in seismic data acquisition |
| CN102023310A (zh) * | 2009-09-19 | 2011-04-20 | 中国石油化工集团公司 | 三维观测系统均匀性分析方法 |
| CN103576213A (zh) * | 2012-07-27 | 2014-02-12 | 中国石油天然气集团公司 | 一种炮检距属性均匀性的确定方法 |
| CN103675893A (zh) * | 2012-09-04 | 2014-03-26 | 中国石油天然气集团公司 | 海上拖缆观测系统复杂模板满覆盖布设方法 |
| CN104142514A (zh) * | 2013-10-29 | 2014-11-12 | 中国石油化工股份有限公司 | 一种三维地震观测系统定量设计方法 |
| CN104142518A (zh) * | 2013-10-29 | 2014-11-12 | 中国石油化工股份有限公司 | 一种地震观测系统叠前时间偏移响应的分析方法 |
| CN106199706A (zh) * | 2015-06-01 | 2016-12-07 | 中国石油化工股份有限公司 | 三维观测系统面元属性统计方法及装置 |
| CN107561574A (zh) * | 2017-07-20 | 2018-01-09 | 中国石油天然气集团公司 | 炮点位置、接收点位置的检测方法和装置 |
| CN109655879A (zh) * | 2017-10-11 | 2019-04-19 | 中国石油化工股份有限公司 | 目的层信噪比一致性能量照明的观测系统优化方法及装置 |
| CN109298450A (zh) * | 2018-11-07 | 2019-02-01 | 中国石油天然气股份有限公司 | 地震采集观测系统属性分析方法及装置 |
Non-Patent Citations (8)
| Title |
|---|
| 《CMP面元分析技术在弯曲测线施工中的应用》;江锋;内蒙古石油化工(第17期);第87-91页 * |
| 《Research on seismic survey design for doubly complex areas》;Zhao Hu et al;APPLIED GEOPHYSICS;第9卷(第3期);第279-285页 * |
| 《一种可行的智能三维特观设计系统的设计方法》;张立柱;辽宁工程技术大学学报;第22卷(第1期);第41-43页 * |
| 《三维地震勘探应慎用Flexi面元法》;陈学强;勘探地球物理进展;第28卷(第5期);第335-340页 * |
| 《三维观测系统均匀性技术及应用效果分析》;马峰等;石油天然气学报;第34卷(第5期);第58-64页 * |
| 《东濮凹陷三维地震重新采集技术》;蒋先勇;断块油气田;第10卷(第1期);第18-21页 * |
| 《提高地震波照明均匀性的加密炮设计新方法》;秦龙等;石油地球物理勘探;第51卷(第4期);第639-646页 * |
| 《胜利探区高密度地震采集技术应用研究》;邱燕;科技传播(第10期);第118、138页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112711066A (zh) | 2021-04-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104133245B (zh) | 一种地震资料的静校正方法及系统 | |
| Zhang et al. | Determination of statistical discontinuity persistence for a rock mass characterized by non-persistent fractures | |
| CN102494667B (zh) | 一种表征地面沉降的方法 | |
| CN103543465A (zh) | 基于目的层照明能量确定有效激发点的方法 | |
| CN107678057B (zh) | 三维地震中确定或优化检波器或设备投入量及排布的方法 | |
| CN103605157B (zh) | 衰减近地表散射波的方法 | |
| CN103454680A (zh) | Walk-away VSP观测系统垂向覆盖次数的计算方法 | |
| CN103576213B (zh) | 一种炮检距属性均匀性的确定方法 | |
| CN116703244B (zh) | 采矿沉陷区治理效果与综合评价方法 | |
| CN102169187B (zh) | 三维地震观测系统整体均衡性定量分析方法 | |
| CN102636813A (zh) | 一种用于处理物探地震数据的动校正拉伸切除方法 | |
| CN112711066B (zh) | 一种三维地震勘探炮点布设均匀性评价方法和装置 | |
| CN104866664B (zh) | 一种防渗止水帷幕体施工难度预测方法 | |
| CN100349014C (zh) | 重力勘探数据处理变密度地形校正方法 | |
| CN109542995B (zh) | 一种库区滑坡分析方法 | |
| CN114236624B (zh) | 基于电磁法估算压裂改造空间体积的方法和系统 | |
| CN105223630B (zh) | 基于地质模型的全方位观测系统参数论证方法 | |
| CN104570062B (zh) | 一种以激发为中心的vsp观测系统设计方法 | |
| CN100562764C (zh) | 自然电位连续剖面勘探方法 | |
| CN105259577B (zh) | 一种确定地层界面的角度信息的方法及装置 | |
| CN107942373A (zh) | 基于裂缝性油气储层断裂系统检测的相干算法 | |
| Gosar | Review of seismological investigations related to 1998 M w 5.6 and 2004 M w 5.2 earthquakes in Krn Mountains. | |
| CN107085236A (zh) | 最大炮检距的确定方法和装置 | |
| Haerani et al. | On the Performance of Terrestrial Laser Scanner for Volcanic and Landslide Hazard Assessment in Indonesia | |
| Rasskazov et al. | Electrical tomography-based imaging of mineral deposits with complex geology |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |