[go: up one dir, main page]

CN111399031B - A method of collecting and processing mountain seismic data - Google Patents

A method of collecting and processing mountain seismic data Download PDF

Info

Publication number
CN111399031B
CN111399031B CN202010355435.3A CN202010355435A CN111399031B CN 111399031 B CN111399031 B CN 111399031B CN 202010355435 A CN202010355435 A CN 202010355435A CN 111399031 B CN111399031 B CN 111399031B
Authority
CN
China
Prior art keywords
virtual
seismic
seismic data
point
processing
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
Application number
CN202010355435.3A
Other languages
Chinese (zh)
Other versions
CN111399031A (en
Inventor
蒋裕强
蒋婵
付永红
周亚东
蒋增政
谷一凡
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sichuan Jierui Taike Technology Co ltd
Southwest Petroleum University
Original Assignee
Sichuan Jierui Taike Technology Co ltd
Southwest Petroleum University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Sichuan Jierui Taike Technology Co ltd, Southwest Petroleum University filed Critical Sichuan Jierui Taike Technology Co ltd
Priority to CN202010355435.3A priority Critical patent/CN111399031B/en
Publication of CN111399031A publication Critical patent/CN111399031A/en
Application granted granted Critical
Publication of CN111399031B publication Critical patent/CN111399031B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/003Seismic data acquisition in general, e.g. survey design
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/301Analysis for determining seismic cross-sections or geostructures
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/303Analysis for determining velocity profiles or travel times
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/10Aspects of acoustic signal generation or detection
    • G01V2210/16Survey configurations
    • G01V2210/161Vertical seismic profiling [VSP]
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/62Physical property of subsurface
    • G01V2210/622Velocity, density or impedance
    • G01V2210/6222Velocity; travel time

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)

Abstract

本发明公开了一种山地地震资料的采集及处理方法,涉及油气勘探领域,包括(1)根据山前带及其地形设计虚拟VSP观测系统,并进行地震数据采集;(2)对采集的地震数据进行处理,得到地震剖面及深度速度模型,通过设计虚拟VSP采集系统实施对山前带的地震资料采集,并进行相关的炮点及检波点的校正工作;利用非零偏VSP的地震资料处理方法,对相关的地震数据进行处理;再对伪VSP地震剖面数据中的CDP点坐标进行坐标转换,从而调整为在实测坐标系统中显示的地震剖面,根据本发明所采集得到的地震数据及相关处理结果,可用于复杂构造带的构造解释及储层预测、裂缝预测及提供相关深度速度模型等多项工作。

Figure 202010355435

The invention discloses a method for collecting and processing mountain seismic data, which relates to the field of oil and gas exploration. After processing the data, the seismic profile and depth velocity model are obtained, and the seismic data acquisition of the piedmont zone is carried out by designing a virtual VSP acquisition system, and the correction of the relevant shot points and detection points is carried out; seismic data processing using non-zero offset VSP method, the relevant seismic data is processed; and then the coordinates of the CDP point in the pseudo-VSP seismic profile data are converted into coordinates, so as to be adjusted to the seismic profile displayed in the measured coordinate system, and the seismic data collected according to the present invention and related The processing results can be used for structural interpretation of complex structural belts, reservoir prediction, fracture prediction, and provision of relevant depth velocity models.

Figure 202010355435

Description

一种山地地震资料的采集及处理方法A method of collecting and processing mountain seismic data

技术领域technical field

本发明涉及油气勘探领域,尤其涉及一种山地地震资料的采集及处理方法。The invention relates to the field of oil and gas exploration, in particular to a method for collecting and processing mountain seismic data.

背景技术Background technique

在地球物理勘探领域中,通过常规地震采集方法获取地震数据后,需要对地震数据进行处理、解释及反演,利用所得到的成果进行勘探井位论证。所以,地震采集所得到的地震数据及相关处理技术方法对后续的地震解释工作影响较大,决定了后续相关工作的成败。在山地地震资料处理中,地层的深度及速度的求取相当关键,它涉及到叠前深度偏移成像的成败。In the field of geophysical exploration, after acquiring seismic data through conventional seismic acquisition methods, it is necessary to process, interpret and invert the seismic data, and use the obtained results to demonstrate exploration well locations. Therefore, the seismic data obtained by seismic acquisition and related processing techniques have a great influence on the subsequent seismic interpretation work, which determines the success or failure of the subsequent related work. In the processing of mountain seismic data, it is very important to obtain the depth and velocity of the formation, which is related to the success or failure of the pre-stack depth migration imaging.

大量的勘探实践表明,在构造高陡区域,常规的地震反射成像不理想。如川东北盆缘带的某些勘探区,由于激发及接收条件恶劣,造成所得的地震资料不能准确成像。并且对地震资料进行多家及多次成像处理,无论叠前深度偏移或叠后偏移等技术,都不能使用相关处理成果进行解释工作。导致该勘探区的油气勘探工作进入停滞状态。在现今的油气勘探中,山前带得到了相关油气勘探公司的重视,并在山前带中钻获一批高产油气流。因此,山前带现在成了油气勘探中的重点研究区域。A large number of exploration practices have shown that conventional seismic reflection imaging is not ideal in high and steep structural areas. For example, in some exploration areas in the basin margin belt of northeastern Sichuan, the obtained seismic data cannot be accurately imaged due to poor excitation and reception conditions. In addition, multiple and multiple imaging processes are performed on the seismic data. Regardless of techniques such as pre-stack depth migration or post-stack migration, the relevant processing results cannot be used for interpretation work. As a result, the oil and gas exploration work in this exploration area has come to a standstill. In today's oil and gas exploration, the piedmont belt has attracted the attention of relevant oil and gas exploration companies, and a batch of high-yield oil and gas streams have been drilled in the piedmont belt. Therefore, the piedmont is now a key research area in oil and gas exploration.

总的来说,制约山前带的地震勘探,除了高昂的采集费用外,还有地震资料成像的问题。主要由于山前带的构造复杂、地层发生倒转现象。这样的地质情况导致地层的层速度求取不准确,从而影响到相关常规地震资料的处理及成像效果,如常规的叠前深度偏移处理等。因此,求取山区相对准确的深度速度模型也是相当关键的。一些发明专利如《一种基于地形因子的炮点偏移方法》提出了根据工区数字高程模型计算工区地形因子坡度值,得到工区坡度数字模型文件,根据工区坡度数字模型文件,使用双线性内插法得到每一炮点的坡度值,根据每一炮点的坡度值和设定的坡度限值,确定需要偏移的炮点,在需要偏移的炮点的纵向和横向偏移限制值范围内,选取满足所述设定的坡度限值的位置,对所述需要偏移的炮点进行偏移;发明专利如《基于地质地球物理模型的三维地震观测系统优化设计方法》提出建立先验的地质地球物理模型,设计两种或多种三维地震数据采集的观测系统,采用虚谱法三维波动方程数值模拟与三维地震物理模型模拟两种方法联合正演来完成数据采集工作,以共反射点为基础,对采集的数据进行叠前深度偏移处理,得到地震处理几种成像结果,对所述地震处理几种成像结果进行成像质量的综合对比评价,经过对比,确定其中一个最佳地震成像结果。从近年来的山地地震资料采集及处理的相关成果来看,主要存在以下一些问题:In general, in addition to the high cost of acquisition, the seismic exploration of the piedmont zone is restricted by the problem of seismic data imaging. Mainly due to the complex structure of the piedmont zone and the phenomenon of stratum inversion. Such geological conditions lead to inaccurate calculation of the layer velocity of the formation, thus affecting the processing and imaging effects of related conventional seismic data, such as conventional pre-stack depth migration processing. Therefore, it is also crucial to obtain a relatively accurate depth velocity model in mountainous areas. Some invention patents, such as "A method of shot point migration based on terrain factor", propose to calculate the slope value of the terrain factor of the working area according to the digital elevation model of the working area, and obtain the digital model file of the slope of the working area. Interpolate to get the slope value of each shot point. According to the slope value of each shot point and the set slope limit, determine the shot point that needs to be offset, and the longitudinal and lateral offset limit values of the shot point that needs to be offset. Within the range, select the position that satisfies the set slope limit, and offset the shot point that needs to be offset; the invention patent such as "Optimization Design Method of 3D Seismic Observation System Based on Geological and Geophysical Model" proposes to establish a first The three-dimensional wave equation numerical simulation of the virtual spectrum method and the three-dimensional seismic physical model simulation are combined forward modeling to complete the data acquisition work. Based on the reflection points, pre-stack depth migration is performed on the collected data to obtain several imaging results of seismic processing. Seismic imaging results. Judging from the related achievements of mountain seismic data acquisition and processing in recent years, there are mainly the following problems:

(1)常规平原上建立的采集系统不适合于山地地震资料采集,亟需要建立一套适合于山地地震资料采集的系统。(1) The acquisition system established on the conventional plain is not suitable for the acquisition of mountain seismic data, and it is urgent to establish a system suitable for the acquisition of mountain seismic data.

(2)常规的地震资料处理技术不适合于山地地震资料处理,也需要建立一套有别于常规地震资料处理的技术方法,使其适合于山地地震资料的处理工作。(2) Conventional seismic data processing technology is not suitable for mountain seismic data processing. It is also necessary to establish a set of technical methods different from conventional seismic data processing to make it suitable for mountain seismic data processing.

发明内容SUMMARY OF THE INVENTION

本发明所要解决的技术问题是:提供一种山地地震资料采集及处理的方法,以解决常规地震资料处理中的山区深度速度模型建立、地震资料成像不准确的问题,并且具有经济可行、技术可操作性强的特点。The technical problem to be solved by the present invention is to provide a method for collecting and processing mountain seismic data, so as to solve the problems of establishing a depth velocity model in a mountain area and inaccurate seismic data imaging in conventional seismic data processing, and is economically feasible and technically feasible. Features of strong operability.

本发明提供的一种山地地震资料的采集及处理方法,包括,A method for collecting and processing mountain seismic data provided by the present invention includes:

(1)根据山前带及其地形设计虚拟VSP观测系统,并进行地震数据采集,且得到相关的地震数据;(1) Design a virtual VSP observation system according to the piedmont and its topography, collect seismic data, and obtain relevant seismic data;

(2)根据虚拟VSP观测系统及相关处理技术,对采集的地震数据进行处理,得到地震剖面及深度速度模型。(2) According to the virtual VSP observation system and related processing technology, the acquired seismic data are processed to obtain the seismic profile and depth velocity model.

所述山前带是指相关地表地貌发生相对剧烈变化的区域,山前带前面一般是平缓地区,地表起伏不大,其中的山体主要是高差不大的丘陵。而山前及山体则与其面前的平原呈现高差较大的状态,高大山体区域的岩石相对出露,风化层覆盖少,植物一般不高大,山顶顶部相对较为平缓。而山前的平原或丘陵区则风化层较厚,植被相对丰富其长得高大,地形相对较平。山体的边缘一般呈高陡状态。从相关的卫星图上可以了解山前带与平原的区别及分界线。The piedmont zone refers to an area where the relevant surface landforms undergo relatively drastic changes. The front of the piedmont zone is generally a flat area with little surface fluctuation, and the mountains in it are mainly hills with small height differences. On the other hand, the front of the mountain and the mountain have a large height difference with the plain in front of it. The rocks in the tall mountain area are relatively exposed, the weathering layer is less covered, the plants are generally not tall, and the top of the mountain is relatively flat. On the other hand, the plain or hilly area in front of the mountain has thicker regolith, relatively rich vegetation, tall and relatively flat terrain. The edge of the mountain is generally high and steep. From the relevant satellite images, we can understand the difference and dividing line between the piedmont and the plain.

更进一步的,所述步骤(1)具体包括,Further, the step (1) specifically includes,

(1-1)根据山前带的地质情况及地形条件,设计地震测线;(1-1) Design seismic lines according to the geological conditions and topographic conditions of the piedmont zone;

(1-2)对目的层的射线进行正演分析,以确定野外观测系统参数;(1-2) Perform forward analysis on the rays of the target layer to determine the parameters of the field observation system;

(1-3)通过地震数据采集系统实施地震数据采集。(1-3) Implement seismic data acquisition through a seismic data acquisition system.

更进一步的,所述步骤(1-1)具体包括,Further, the step (1-1) specifically includes,

在山前带的平原地区低部位及山体的斜坡部位设置检波点(接收);Set up the detection point (reception) in the low part of the plain area of the piedmont and the slope part of the mountain;

在山体顶部设置炮点(激发)。Set a shot (excitation) on top of the hill.

或将两者的激发及接收位置进行互换。Or exchange the excitation and reception positions of the two.

该虚拟VSP观测系统主要特征在于不是在井中接收的,而是将相关的激发点设置在山前带的山体顶部表面,而接收点则设置在山下的平原地带上及山体斜坡上。这样所设计的炮点及检波点的位置关系就形成一个相对于井与地表关系的倒“L”型的非零偏VSP观测系统。The main feature of the virtual VSP observation system is that it is not received in the well, but the relevant excitation points are set on the top surface of the mountain in the piedmont, and the receiving points are set on the plains and slopes below the mountain. In this way, the positional relationship of shot points and detection points designed to form an inverted "L" type non-zero offset VSP observation system relative to the relationship between well and surface.

更进一步的,所述野外观测系统参数包括炮点及检波点的位置、移动方向和距离,根据实际的VSP观测系统设计出炮点及检波点的位置及相关的移动方向、距离等参数;Further, the parameters of the field observation system include the position, moving direction and distance of the shot point and the detection point, and the position of the shot point and the detection point and the relevant parameters such as the moving direction and distance are designed according to the actual VSP observation system;

炮点及检波点的距离为等距或不等距;The distance between the shot point and the receiver point is equidistant or unequal;

还可以使炮点及检波点的移动方向一致或非一致或设置检波点不动,移动炮点的激发接收方式。这样可以形成一个二维的地震数据采集系统。这个地震数据采集系统与常规的二维数据采集系统是类似的,经过对设计好的地震数据数据采集系统实施地震数据采集,从而得到相关的地震数据体。It is also possible to make the moving direction of the shot point and the receiver point consistent or non-uniform, or set the receiver point to remain stationary, and move the excitation and reception mode of the shot point. In this way, a two-dimensional seismic data acquisition system can be formed. This seismic data acquisition system is similar to the conventional two-dimensional data acquisition system. After the seismic data acquisition is carried out on the designed seismic data acquisition system, the relevant seismic data volume is obtained.

更进一步的,所述步骤(2)具体包括,Further, the step (2) specifically includes,

(2-1)对野外采集得到的地震资料中的炮点及检波点的位置进行校正,将其实测的坐标及高程校正到所述虚拟VSP观测系统上,从而形成相关虚拟井深接收点及地表激发点;(2-1) Correct the positions of shot points and detection points in the seismic data collected in the field, and correct the measured coordinates and elevations to the virtual VSP observation system, so as to form relevant virtual well depth receiving points and surface excitation point;

(2-2)将炮点及检波点的校正量加入到地震道头数据中,并进行处理,得到初始VSP地震剖面;(2-2) Add the corrections of shot points and receiver points to the seismic tracehead data, and process them to obtain the initial VSP seismic profile;

(2-3)对初始VSP地震剖面进行CDP坐标变换处理,得到与实测地震剖面位置相吻合的伪VSP地震剖面。(2-3) The CDP coordinate transformation is performed on the initial VSP seismic profile to obtain a pseudo-VSP seismic profile that matches the position of the measured seismic profile.

更进一步的,所述步骤(2-1)具体包括,Further, the step (2-1) specifically includes,

确定实测虚拟井口点(xo,yo,ho)及虚拟井底点(xo,yo,hi)Determine the measured virtual wellhead point (x o , yo , h o ) and virtual bottom hole point (x o , yo , hi )

将炮点及检波点分别校正到设定的虚拟地表线、虚拟井轨迹上;Correct the shot point and detection point respectively to the set virtual surface line and virtual well trajectory;

所述虚拟井轨迹为连接虚拟井口点及虚拟井底点的直线且垂直于基准面;The virtual well trajectory is a straight line connecting the virtual wellhead point and the virtual bottom hole point and is perpendicular to the reference plane;

所述虚拟地表线为井口点引出与虚拟井轨迹成垂线的直线且平行于基准面;The virtual surface line is a straight line drawn from the wellhead point that is perpendicular to the virtual well trajectory and is parallel to the reference plane;

所述虚拟井底点为高程相对较小的位置,且与山前带的位置相对较近,不影响相关的地震激发、接收的信噪比及分辨率。The virtual bottom hole is a position with a relatively small elevation, and is relatively close to the position of the piedmont, and does not affect the related seismic excitation, signal-to-noise ratio and resolution of reception.

更进一步的,所述步骤(2-2)具体包括,Further, the step (2-2) specifically includes,

根据炮点及检波点实测的坐标及高程数据将两者分别利用替换速度及距离数据,与虚拟地表线或虚拟井轨迹的垂直距离进行计算,以得到相应炮点及检波点的静校正量;According to the measured coordinates and elevation data of the shot point and the receiver point, respectively use the replacement velocity and distance data to calculate the vertical distance from the virtual surface line or the virtual well trajectory to obtain the static correction amount of the corresponding shot point and the receiver point;

利用静校正量对相应的炮点及检波点的地震数据进行校正,并将相关的炮点及检波点分别校正到虚拟地表线及虚拟井轨迹上;Correct the seismic data of the corresponding shot point and receiver point by using the static correction amount, and correct the related shot point and receiver point to the virtual surface line and virtual well trajectory respectively;

结合山区的地表地质情况,利用常规的非零偏VSP处理流程实施地震资料处理,得到初始VSP地震剖面。Combined with the surface geological conditions of the mountainous area, the conventional non-zero offset VSP processing flow is used to process the seismic data to obtain the initial VSP seismic profile.

更进一步的,所述步骤(2-3)具体包括,Further, the step (2-3) specifically includes,

建立虚拟VSP观测系统的虚拟坐标系统与实测坐标系统之间转换的数学关系式;Establish the mathematical relationship between the virtual coordinate system of the virtual VSP observation system and the measured coordinate system;

建立炮点及检波点的二维坐标关系;Establish the two-dimensional coordinate relationship between the shot point and the receiver point;

通过虚拟坐标系统与实测坐标系统之间的转换关系公式,将虚拟VSP观测系统的CDP点坐标转换到实测坐标系统中的CDP点坐标并进行相关层位的埋深及速度求取、地震剖面显示。Through the conversion formula between the virtual coordinate system and the measured coordinate system, the CDP point coordinates of the virtual VSP observation system are converted to the CDP point coordinates in the measured coordinate system, and the buried depth and velocity of the relevant horizons are obtained, and the seismic profile is displayed. .

更进一步的,所述虚拟坐标系统具体以虚拟地表线为X轴、虚拟井轨迹为H轴建立以虚拟井口点为原点的虚拟坐标系统。Further, the virtual coordinate system specifically takes the virtual surface line as the X-axis and the virtual well trajectory as the H-axis to establish a virtual coordinate system with the virtual wellhead point as the origin.

通过采用以上的技术方案,本发明的有益效果是:通过设计虚拟VSP采集系统实施对山前带的地震资料采集,并进行相关的炮点及检波点的校正工作;利用非零偏VSP的地震资料处理方法,对相关的地震数据进行处理;再对伪VSP地震剖面数据中的CDP点坐标进行坐标转换,从而调整为在实测坐标系统中显示的地震剖面。根据本发明所采集得到的地震数据及相关处理结果,可用于复杂构造带的构造解释及储层预测、裂缝预测及提供相关深度速度模型等多项工作,具有比常规地震采集所得到的地震数据所包含信息量更多及准确、经济可行、技术可靠的特点。By adopting the above technical scheme, the beneficial effects of the present invention are: by designing a virtual VSP acquisition system, the seismic data acquisition of the piedmont zone is carried out, and the correction of the relevant shot points and detection points is carried out; The data processing method is to process the relevant seismic data; then, the coordinates of the CDP points in the pseudo-VSP seismic profile data are converted into coordinates to adjust to the seismic profile displayed in the measured coordinate system. The seismic data and related processing results collected according to the present invention can be used for structural interpretation of complex structural belts, reservoir prediction, fracture prediction and provision of relevant depth velocity models, etc. The information contained is more and accurate, economically feasible, and technically reliable.

附图说明Description of drawings

本发明将通过例子并参照附图的方式说明,其中:The invention will be described by way of example and with reference to the accompanying drawings, in which:

图1为根据本发明实施例的地震采集参数设计方法的流程示意图;1 is a schematic flowchart of a method for designing seismic acquisition parameters according to an embodiment of the present invention;

图2为根据本发明实施例的相关地震采集参数及虚拟井轨迹及地表线示意图。FIG. 2 is a schematic diagram of relevant seismic acquisition parameters, virtual well trajectory and surface line according to an embodiment of the present invention.

具体实施方式Detailed ways

本说明书中公开的所有特征,或公开的所有方法或过程中的步骤,除了互相排斥的特征和/或步骤以外,均可以以任何方式组合。All features disclosed in this specification, or all disclosed steps in a method or process, may be combined in any way except mutually exclusive features and/or steps.

下面将结合附图对本发明作进一步说明。The present invention will be further described below with reference to the accompanying drawings.

本实施例的执行主体为地震采集参数设计,并建立针对该设计地震采集参数的观测系统及地震资料处理方法。具体为设计一种针对山地地震资料采集及处理的伪VSP系统,将得到的地震数据进行相关的校正后实施非零偏VSP资料处理,从而可以提供相对准确的深度及速度模型、地震剖面。具有相对经济、可行性高等特点,值得在相关勘探区的山前带地震勘探中进行推广。The execution subject of this embodiment is seismic acquisition parameter design, and an observation system and seismic data processing method for the designed seismic acquisition parameters are established. Specifically, a pseudo-VSP system is designed for the acquisition and processing of mountain seismic data, and the obtained seismic data is subjected to non-zero bias VSP data processing after relevant correction, so as to provide relatively accurate depth and velocity models and seismic profiles. It has the characteristics of relative economy and high feasibility, and it is worth popularizing in the seismic exploration of the piedmont zone in the relevant exploration area.

图1为根据本发明实施例的针对某山前带的地震采集参数设计及处理流程示意图,如图1所示,本发明提供的一种山地地震资料采集的方法,包括:1 is a schematic diagram of the design and processing flow of seismic acquisition parameters for a certain piedmont zone according to an embodiment of the present invention. As shown in FIG. 1 , a method for acquiring mountain seismic data provided by the present invention includes:

①根据山前带及其地形设计相关的虚拟VSP观测系统,并进行相关的地震数据采集且得到相关的地震数据;① Design the relevant virtual VSP observation system according to the piedmont and its topography, and collect the relevant seismic data and obtain the relevant seismic data;

②根据虚拟VSP观测系统及相关处理技术,对采集的地震数据实施地震资料处理并进行坐标变换,得到相关的地震剖面及深度速度模型;②According to the virtual VSP observation system and related processing technology, implement seismic data processing and coordinate transformation on the collected seismic data to obtain the relevant seismic profile and depth velocity model;

根据山前带及其地形设计相关的虚拟VSP观测系统,并进行相关的地震数据采集及处理后得到相关的地震数据,包括以下步骤:Design a virtual VSP observation system according to the piedmont and its topography, and collect and process the relevant seismic data to obtain the relevant seismic data, including the following steps:

在步骤1中,根据山前带及其地形设计相关的虚拟VSP观测系统,并进行相关的地震数据采集且得到相关的地震数据。主要操作为建立适合山前带的虚拟VSP观测系统。该系统的获得主要是根据对山前带区的地质情况及地形条件,设计出相关的地震测线。该测线主要是常规二维地震勘探形式的,在山前带的平原地区低部位及山体的斜坡部位设计相关的检波点位置,形成二维测线的一部分,检波点间隔在垂直方向上的投影(虚拟井轨迹上)可以是等距或不等距设置;在山体顶部设置炮点,炮点在基准面上的投影可以是等距或不等距设置。并利用相关的地震采集正演软件实施针对目的层的射线正演分析,从而确定出对应的野外观测系统参数。In step 1, a related virtual VSP observation system is designed according to the piedmont belt and its topography, and the related seismic data is collected and obtained. The main operation is to establish a virtual VSP observation system suitable for the piedmont. The acquisition of the system is mainly based on the design of the relevant seismic lines according to the geological conditions and topographic conditions of the piedmont area. The survey line is mainly in the form of conventional two-dimensional seismic exploration. The relevant detection points are designed in the low part of the plain area of the piedmont and the slope of the mountain to form a part of the two-dimensional survey line. The detection point interval is in the vertical direction. Projection (on the virtual well trajectory) can be set in equidistant or non-equidistant settings; set shot points on the top of the mountain, and the projection of shot points on the datum plane can be set in equidistant or non-equidistant settings. And use the relevant seismic acquisition forward modeling software to implement the ray forward modeling analysis for the target layer, so as to determine the corresponding field observation system parameters.

其中,山前带是指相关地表地貌发生相对剧烈变化的区域,山前带前面一般是平缓地区,地表起伏不大,其中的山体主要是高差不大的丘陵。而山前及山体则与其面前的平原呈现高差较大的状态,高大山体区域的岩石相对出露,风化层覆盖少,植物一般不高大,山顶顶部相对较为平缓。而山前的平原或丘陵区则风化层较厚,植被相对丰富其长得高大,地形相对较平。山体的边缘一般呈高陡状态。从相关的卫星图上可以了解山前带与平原的区别及分界线。Among them, the piedmont zone refers to the area where the relevant surface landforms have undergone relatively drastic changes. The front of the piedmont zone is generally a gentle area with little fluctuation of the surface, and the mountains are mainly hills with small height differences. On the other hand, the front of the mountain and the mountain have a large height difference with the plain in front of it. The rocks in the tall mountain area are relatively exposed, the weathering layer is less covered, the plants are generally not tall, and the top of the mountain is relatively flat. On the other hand, the plain or hilly area in front of the mountain has thicker regolith, relatively rich vegetation, tall and relatively flat terrain. The edge of the mountain is generally high and steep. From the relevant satellite images, we can understand the difference and dividing line between the piedmont and the plain.

其中,虚拟VSP观测系统主要特征在于不是在井中接收的,而是将相关的接收检波器设置在山前带的斜坡部位、平原低部位内,而激发点则设置在山顶区域上。这样所设计的炮点及检波点的位置关系就形成一个相对于井与地表关系的倒“L”型的非零偏VSP观测系统,并根据实际的非零VSP观测系统及正演结果设计出炮点及检波点的位置及相关的移动方向、距离等参数。本发明的所述炮点及检波点距离可以为等距或不等距离,炮点及检波点的移动方向可以一致的或非一致,也可以设置为检波点不移动,而炮点移动的激发接收方式,这样可以形成一个二维的地震数据采集系统。这个地震数据采集系统与常规的二维数据采集系统是类似的,经过对设计好的地震数据数据采集系统实施地震数据采集,从而得到相关的地震数据体。Among them, the main feature of the virtual VSP observation system is that it is not received in the well, but the relevant receiving geophones are set in the slope part of the piedmont belt and the low part of the plain, and the excitation point is set on the mountain top area. In this way, the designed positional relationship between the shot point and the receiver point forms an inverted "L" type non-zero offset VSP observation system relative to the relationship between the well and the surface. According to the actual non-zero VSP observation system and forward modeling results, a The position of the shot point and the detection point and the related moving direction, distance and other parameters. The distance between the shot point and the detection point of the present invention can be equidistant or unequal distance, the moving direction of the shot point and the detection point can be consistent or non-uniform, or it can be set so that the detection point does not move, and the excitation of the shot point moving In this way, a two-dimensional seismic data acquisition system can be formed. This seismic data acquisition system is similar to the conventional two-dimensional data acquisition system. After the seismic data acquisition is carried out on the designed seismic data acquisition system, the relevant seismic data volume is obtained.

在步骤2中,根据虚拟VSP观测系统及相关处理技术,对采集的地震数据实施地震资料处理并进行坐标变换,得到相关的地震剖面及深度速度模型,包括以下步骤:In step 2, according to the virtual VSP observation system and related processing technology, seismic data processing and coordinate transformation are performed on the acquired seismic data to obtain the relevant seismic profile and depth velocity model, including the following steps:

①对炮点及检波点位置进行校正,使其校正到相关的倒“L”型的非零偏VSP观测系统上。具体操作为对虚拟井口点(xo,yo,ho)及虚拟井底点(xi,yi,hi)、虚拟地表线数据,将连接虚拟井口点及虚拟井底点的直线设定为虚拟井轨迹,井口点引出与虚拟井轨迹成垂直的直线作为虚拟地表线。将炮点及检波点分别校正到设定的虚拟地表线、虚拟的井轨迹上。虚拟的地表线的为一条平行于基准面(平面)的直线,也就是虚拟井口(ho)的高程是一样的,该虚拟井轨迹对虚拟井口对基准面作垂线与地表面相交的位置就是所设计的检波点位置。虚拟井轨迹也为直线,该直线为虚拟井口向基准面作垂线与地表某点相交所生成。该地表点原则上是高程相对较小的位置,并且与山前带的位置相对较近,不影响相关的地震激发、接收的信噪比及分辨率。在对采集的地震数据进行处理中,也可以对相关的井口点(xo,yo,ho)三参数进行测试,优选合适的虚拟井口点(xo、yo、ho)进入相关的处理流程中。在实际操作中,根据炮点及检波点实测的坐标(xj,yj)及高程数据(hj)将两者分别利用替换速度及相关的距离数据—与相关虚拟地表线或虚拟井轨迹的垂直距离进行计算,从而得到相关炮点及检波点的静校正量。在本发明技术中,所设定的实测三维坐标系是实际地震勘探中所使用的坐标系是一致的,也就是平面上是x、y坐标轴加上一个垂直该平面的h轴,深度h定义为高于海平面为正值,低于海平面为负值。① Correct the position of shot point and detection point to make it correct to the relevant inverted "L" type non-zero offset VSP observation system. The specific operation is to connect the virtual wellhead point (x o , y o , h o ), the virtual bottom hole point ( xi , yi , hi ), and the virtual surface line data, and connect the virtual wellhead point and the virtual bottom hole point. It is set as a virtual well trajectory, and a straight line perpendicular to the virtual well trajectory is drawn from the wellhead point as a virtual surface line. Correct the shot point and detection point to the set virtual surface line and virtual well trajectory respectively. The virtual surface line is a straight line parallel to the reference plane (plane), that is, the elevation of the virtual wellhead (h o ) is the same. is the location of the designed detection point. The virtual well trajectory is also a straight line, which is generated by intersecting a vertical line from the virtual wellhead to the reference plane and a certain point on the surface. In principle, the surface point is a position with relatively small elevation and is relatively close to the position of the piedmont, which does not affect the related seismic excitation, received signal-to-noise ratio and resolution. In the processing of the acquired seismic data, the three parameters of the relevant wellhead points (x o , yo , ho ) can also be tested, and the appropriate virtual wellhead points (x o , yo , ho ) are preferably entered into the relevant wellhead points. in the processing flow. In actual operation, according to the measured coordinates (x j , y j ) and elevation data (h j ) of the shot point and the detection point, the two are respectively used to replace the speed and the related distance data—with the related virtual surface line or virtual well trajectory. The vertical distance is calculated to obtain the static correction amount of the relevant shot point and detection point. In the technology of the present invention, the set measured three-dimensional coordinate system is consistent with the coordinate system used in actual seismic exploration, that is, the x and y coordinate axes on the plane plus an h axis perpendicular to the plane, the depth h Defined as positive values above sea level and negative values below sea level.

②使用针对虚拟VSP井系统校正后的炮点及检波点的校正量加入到地震道头数据中,利用相关的非零偏VSP处理流程及软、硬件实施地震资料处理,得到地震资料处理剖面,并将伪VSP地震剖面的CDP上的虚拟坐标数据进行相关的转换计算,从而得到与实测坐标系统相吻合的伪VSP地震剖面。其中,相关操作步骤具体为利用相关的静校正量分别对炮点及检波点的地震数据进行校正,并将相关的炮点及检波点分别校正到虚拟地表线及虚拟井轨迹上。然后结合山区的地表地质情况,利用常规的非零偏VSP处理流程实施地震资料处理,从而得到一个VSP地震剖面。(2) Use the correction amounts of shot points and detection points corrected for the virtual VSP well system to be added to the seismic trace head data, and use the relevant non-zero offset VSP processing flow and software and hardware to process the seismic data to obtain the seismic data processing profile. The virtual coordinate data on the CDP of the pseudo-VSP seismic section is calculated by the relevant transformation, so as to obtain the pseudo-VSP seismic section that is consistent with the measured coordinate system. Specifically, the relevant operation steps are to use the relevant static correction amount to respectively correct the seismic data of the shot point and the detection point, and respectively correct the relevant shot point and the detection point to the virtual surface line and the virtual well trajectory. Then, combined with the surface geological conditions of the mountainous area, the conventional non-zero offset VSP processing flow is used to process the seismic data, so as to obtain a VSP seismic profile.

其中,利用相关的静校正量分别对炮点及检波点的地震数据进行校正,并将相关的炮点及检波点分别校正到虚拟地表线及虚拟的井轨迹,建立虚拟坐标系统并建立虚拟坐标系统与实测坐标系统之间转换的数学关系式。实际操作中,以虚拟地表线为X轴、虚拟井轨迹为H轴建立以虚拟井口为原点的虚拟坐标系统,并得到虚拟坐标系统与实测坐标系统的数学转换关系。通过建立虚拟坐标系统,建立了炮点及检波点的二维坐标关系,从而能够生成一个二维地震剖面。一般情况下,将虚拟地表线的坐标轴设定为X轴后,该轴也就是伪VSP地震剖面中CDP点所在的坐标轴。通过虚拟坐标系统与实测坐标系统之间的转换关系公式,可以将伪VSP地震剖面的CDP点坐标转换到实测坐标系统中的CDP点坐标并进行相关层位的埋深及速度求取、地震剖面显示等。Among them, the seismic data of shot points and receiver points are respectively corrected by the relevant static correction amount, and the relevant shot points and receiver points are respectively corrected to the virtual surface line and the virtual well trajectory, and the virtual coordinate system is established and the virtual coordinate system is established. Mathematical relationship for conversion between the system and the measured coordinate system. In actual operation, a virtual coordinate system with the virtual wellhead as the origin is established with the virtual surface line as the X axis and the virtual well trajectory as the H axis, and the mathematical conversion relationship between the virtual coordinate system and the measured coordinate system is obtained. By establishing a virtual coordinate system, the two-dimensional coordinate relationship between the shot point and the receiver point is established, so that a two-dimensional seismic section can be generated. In general, after setting the coordinate axis of the virtual surface line as the X axis, this axis is also the coordinate axis where the CDP point in the pseudo VSP seismic section is located. Through the conversion relationship formula between the virtual coordinate system and the measured coordinate system, the CDP point coordinates of the pseudo VSP seismic profile can be converted to the CDP point coordinates in the measured coordinate system, and the buried depth and velocity of the relevant horizon can be obtained, and the seismic profile can be obtained. display etc.

在本步骤中,主要是结合山区的地表地质情况,利用常规的非零偏VSP处理流程实施本发明技术所得到的地震资料资料处理,从而得到一个过虚拟井的VSP地震剖面,并将该伪VSP地震剖面的CDP点坐标进行转换,得到实测坐标系统中的坐标值,从而能在实测坐标系统中按相关CDP点进行显示。另外,可以从山区的地表地质得到相关的虚拟井轨迹上的地层、岩性、速度等信息,从而为相关VSP资料处理提供参数及标定。对处理后的伪VSP地震剖面的CDP点进行坐标变换,使其经坐标变换后地震剖面相关数据能与实测常规地震剖面一样进行显示,从而完成相关山地地震资料的采集及处理,真实的得到山前带地下的地质情况,为相关的油气勘探服务。另外,利用本发明技术得到的虚拟VSP地震资料的层位解释数据及层速度数据,可以为相关的叠前深度偏移处理提供准确的速度深度模型。In this step, the seismic data obtained by the technology of the present invention is processed by using the conventional non-zero offset VSP processing flow mainly in combination with the surface geological conditions of the mountainous area, so as to obtain a VSP seismic section passing through a virtual well, and the false The CDP point coordinates of the VSP seismic profile are converted to obtain the coordinate values in the measured coordinate system, which can be displayed according to the relevant CDP points in the measured coordinate system. In addition, the strata, lithology, velocity and other information on the relevant virtual well trajectory can be obtained from the surface geology of the mountainous area, so as to provide parameters and calibration for the relevant VSP data processing. The coordinates of the CDP points of the processed pseudo-VSP seismic profile are transformed, so that the related data of the seismic profile after the coordinate transformation can be displayed in the same way as the measured conventional seismic profile, so as to complete the acquisition and processing of the related mountain seismic data, and get the real mountain The geological conditions of the underground forezone serve related oil and gas exploration. In addition, the horizon interpretation data and horizon velocity data of the virtual VSP seismic data obtained by the technology of the present invention can provide an accurate velocity depth model for the related prestack depth migration processing.

其中,非零偏VSP处理流程主要是预处理、同深度叠加、初至拾取、频谱分析及带通滤波、震源子波整形、静态时移和排齐、波场分离、反褶积及VSPCDP叠加成像、偏移成像等。现阶段,非零偏VSP处理技术非常成熟,各种技术方法相对较多。Among them, the non-zero offset VSP processing flow mainly includes preprocessing, same-depth stacking, first-arrival picking, spectrum analysis and bandpass filtering, source wavelet shaping, static time shifting and alignment, wave field separation, deconvolution and VSPCDP stacking imaging, offset imaging, etc. At this stage, the non-zero bias VSP processing technology is very mature, and there are relatively many technical methods.

此外,根据本次发明例子,相关技术人员可以对相关的山地采集方式进行扩展,如可以扩展到山体顶部面上激发方式及平面上立体接收等方式,并设计相关地震处理系统实施地震资料处理及对相关的数据体进行坐标转换及地震资料处理。In addition, according to the example of the present invention, the relevant technical personnel can expand the related mountain acquisition methods, such as the excitation method on the top of the mountain body and the three-dimensional reception on the plane, etc., and design a related seismic processing system to implement seismic data processing and processing. Coordinate transformation and seismic data processing are performed on relevant data volumes.

下面列举具体实施例对上述步骤进行说明。The above steps are described below with reference to specific embodiments.

本实例是对某三维工区的高陡部位的海相页岩层探测实施针对山前带的地震资料采集及处理,主要是依据本发明技术流程来进行(图1)。该勘探区的二维地震资料由于各种干扰及处理方法所影响,目的层成像不清晰,不能符合相关地震资料解释的要求。因此,为说明本次发明技术的采集及处理效果,在原二维地震测线上布设本次发明的采集系统。以此来进行相关地震资料的对比,来判断常规及本次发明技术所得到的地震剖面的优劣。In this example, the acquisition and processing of seismic data for the piedmont zone is carried out for the detection of marine shale layers in a high and steep part of a three-dimensional work area, which is mainly performed according to the technical process of the present invention (Fig. 1). Due to the influence of various interference and processing methods, the 2D seismic data in this exploration area has unclear imaging of the target layer, which cannot meet the requirements for interpretation of relevant seismic data. Therefore, in order to illustrate the acquisition and processing effect of the technology of the present invention, the acquisition system of the present invention is arranged on the original two-dimensional seismic survey line. In this way, the relative seismic data are compared to judge the pros and cons of the seismic profiles obtained by the conventional and the invented technology.

在步骤1中,某勘探区针对页岩气的实际地震数据采集过程中,建立相关的二维地震观测系统实施地震数据的采集,从而形成倒“L”型的虚拟VSP采集系统(图2)。分别在山区高部位上布设炮点(激发点),山区斜坡及平原区域设置检波点,检波点的间隔为4m至10m不等,检波点线总长度为4.1km左右,垂向上为2421m。检波点布设则是顺着山前平地及缓坡向底部位布设,炮点间隔有大有小,总体炮距为10-40m距离不等,在相关范围内设计有近1292个炮点,炮点设计可以偏移虚拟地表线左右一定范围内,符合相关采集规范。根据相关的地震采集设计,完成该线上的地震数据采集工作,得到相关的地震数据。In step 1, during the actual seismic data acquisition process for shale gas in a certain exploration area, a related two-dimensional seismic observation system is established to implement seismic data acquisition, thereby forming an inverted "L"-shaped virtual VSP acquisition system (Fig. 2) . The shot points (excitation points) were arranged on the high parts of the mountainous area, and the detection points were set up on the slopes and plains of the mountainous area. The interval of the detection points was 4m to 10m. The total length of the detection point line was about 4.1km, and the vertical direction was 2421m. The layout of the detection points is along the piedmont flat and the gentle slope to the bottom. The shot point interval is large or small. The overall shot distance varies from 10 to 40m. There are nearly 1292 shot points in the relevant range. The design can be offset within a certain range to the left and right of the virtual surface line, in line with relevant acquisition specifications. According to the relevant seismic acquisition design, complete the seismic data acquisition on the line, and obtain the relevant seismic data.

在步骤2中,对上述采集的地震资料进行相关的地震资料处理工作,得到相关的地震数据用于地震资料解释。根据野外采集的地震数据,针对虚拟地表线及虚拟井轨迹实施炮点及检波点的静校正工作,得到相关的校正量。并且针对非零偏VSP地震资料处理,建立相关的技术流程。利用相关的处理软件及硬件设备,实施地震资料处理工作,得到伪VSP地震剖面数据。并对该地震剖面进行CDP坐标变换处理,得到一个与实测地震剖面位置相吻合的伪VSP地震剖面。In step 2, relevant seismic data processing is performed on the acquired seismic data to obtain relevant seismic data for interpretation of the seismic data. According to the seismic data collected in the field, static correction of shot points and detection points is carried out for the virtual surface line and virtual well trajectory, and the relevant correction amount is obtained. And for the non-zero offset VSP seismic data processing, the related technical process is established. Using relevant processing software and hardware equipment, implement seismic data processing to obtain pseudo-VSP seismic profile data. And the CDP coordinate transformation is carried out on the seismic profile, and a pseudo-VSP seismic profile is obtained which is consistent with the position of the measured seismic profile.

实际操作表明,利用本次针对山前带的野外地震资料采集及处理工作,与原二维地震剖面的对比,具有目的层反射清晰及构造成像准确的特点,优于常规二维地震采集的成果。另外,与布设在该测线上的某口页岩气井的井-震合成记录标定吻合性较好,达到了地质目的。The actual operation shows that the field seismic data acquisition and processing work for the piedmont zone, compared with the original two-dimensional seismic profile, has the characteristics of clear reflection of the target layer and accurate structural imaging, which is superior to the results of conventional two-dimensional seismic acquisition. . In addition, it is in good agreement with the calibration of the well-seismic synthetic record of a shale gas well laid on the survey line, which achieves the geological purpose.

上述说明示出并描述了本发明的一个优选实施例,但如前所述,应当理解本发明并非局限于本文所披露的形式,不应看作是对其他实施例的排除,而可用于各种其他组合、修改和环境,并能够在本文所述发明构想范围内,通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本发明的精神和范围,则都应在本发明所附权利要求的保护范围内。The foregoing description shows and describes a preferred embodiment of the present invention, but as previously mentioned, it should be understood that the present invention is not limited to the form disclosed herein, and should not be construed as an exclusion of other embodiments, but may be used in various and other combinations, modifications and environments, and can be modified within the scope of the inventive concepts described herein, from the above teachings or from skill or knowledge in the relevant art. However, modifications and changes made by those skilled in the art do not depart from the spirit and scope of the present invention, and should all fall within the protection scope of the appended claims of the present invention.

Claims (7)

1. A method for collecting and processing mountain land seismic data is characterized by comprising the following steps: comprises that
(1) Designing a virtual VSP observation system according to the mountain front zone and the terrain thereof, and acquiring seismic data; the step (1) specifically comprises the following steps,
(1-1) designing a seismic survey line according to the geological condition and the topographic condition of the mountain front zone;
the step (1-1) specifically comprises the steps of arranging wave detection points at the low part of a plain area of a mountain front zone and the slope part of a mountain body; setting a shot point at the top of the mountain; or interchanging the two;
(1-2) performing forward analysis on rays of a target layer to determine parameters of a field observation system;
(1-3) performing seismic data acquisition by a seismic data acquisition system;
(2) and processing the acquired seismic data to obtain a seismic section and a depth velocity model.
2. The method for collecting and processing mountain seismic data as claimed in claim 1, wherein: the field observation system parameters comprise positions, moving directions and distances of a shot point and a wave detection point;
the distance between the shot point and the wave detection point is equal or unequal;
the moving directions of the shot point and the wave detection point are consistent or inconsistent or the wave detection point is set to be fixed, and the shot point is moved.
3. The method for collecting and processing mountain seismic data as claimed in claim 1, wherein: the step (2) specifically comprises the following steps,
(2-1) correcting the positions of shot points and wave detection points to the virtual VSP observation system;
(2-2) adding the correction values of the shot point and the demodulator probe into the seismic channel head data, and processing to obtain an initial VSP seismic section;
and (2-3) carrying out CDP coordinate transformation processing on the initial VSP seismic section to obtain a pseudo VSP seismic section matched with the position of the actually measured seismic section.
4. The method for collecting and processing mountain seismic data as claimed in claim 3, wherein: the step (2-1) specifically comprises the steps of,
respectively correcting the shot point and the wave detection point to a set virtual earth surface line and a set virtual well track;
the virtual well track is a straight line connecting the virtual wellhead point and the virtual well bottom point and is vertical to the reference surface;
the virtual earth surface line is a straight line which is led out from a wellhead point and forms a perpendicular line with the virtual well track and is parallel to the reference surface;
the virtual well bottom point is a position with a relatively small elevation and is relatively close to the position of the mountain front zone, and the signal-to-noise ratio and the resolution of the related earthquake excitation and receiving are not influenced.
5. The method for collecting and processing mountain seismic data as claimed in claim 3, wherein: the step (2-2) specifically comprises the steps of,
calculating the vertical distance between the shot point and the demodulator probe and the virtual earth surface line or the virtual well track by respectively using the replacement speed and the distance data according to the actually measured coordinates and elevation data of the shot point and the demodulator probe so as to obtain the static correction value of the corresponding shot point and demodulator probe;
correcting the seismic data of the corresponding shot points and the corresponding demodulator probes by using the static correction values, and respectively correcting the relevant shot points and demodulator probes to the virtual earth surface line and the virtual well track;
and (3) combining the surface geological condition of the mountain area, and performing seismic data processing by using a conventional non-zero-bias VSP processing flow to obtain an initial VSP seismic profile.
6. The method for collecting and processing mountain seismic data as claimed in claim 3, wherein: the step (2-3) specifically comprises the steps of,
establishing a mathematical relation formula for conversion between a virtual coordinate system and an actually measured coordinate system of a virtual VSP observation system;
establishing a two-dimensional coordinate relationship between a shot point and a wave detection point;
and converting the CDP point coordinate of the virtual VSP observation system into the CDP point coordinate in the actual measurement coordinate system through a conversion relation formula between the virtual coordinate system and the actual measurement coordinate system, and performing buried depth and speed calculation and seismic section display of the related layer.
7. The method for collecting and processing mountain seismic data as claimed in claim 6, wherein: the virtual coordinate system specifically uses a virtual earth surface line as an X axis and a virtual well track as an H axis to establish a virtual coordinate system using a virtual well head point as an origin.
CN202010355435.3A 2020-04-29 2020-04-29 A method of collecting and processing mountain seismic data Active CN111399031B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010355435.3A CN111399031B (en) 2020-04-29 2020-04-29 A method of collecting and processing mountain seismic data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010355435.3A CN111399031B (en) 2020-04-29 2020-04-29 A method of collecting and processing mountain seismic data

Publications (2)

Publication Number Publication Date
CN111399031A CN111399031A (en) 2020-07-10
CN111399031B true CN111399031B (en) 2020-12-08

Family

ID=71431832

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010355435.3A Active CN111399031B (en) 2020-04-29 2020-04-29 A method of collecting and processing mountain seismic data

Country Status (1)

Country Link
CN (1) CN111399031B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112099086B (en) * 2020-09-16 2022-03-29 中油奥博(成都)科技有限公司 High-resolution optical fiber borehole seismic data deep frequency analysis method
CN112379424A (en) * 2020-11-02 2021-02-19 中国石油天然气集团有限公司 Graphic representation method and device of mountain two-dimensional seismic survey line observation system

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102213770A (en) * 2010-04-01 2011-10-12 罗仁泽 Combined area excitation receiving method for seismic data in low signal-to-noise ratio area
CN102313900A (en) * 2011-09-01 2012-01-11 中国石油集团川庆钻探工程有限公司 Excitation position determining method of three-dimensional seismic acquisition observation system
CN105487120A (en) * 2014-09-18 2016-04-13 中国石油化工股份有限公司 Method for determining excitation well location of mountain front zone
CN105866833A (en) * 2016-06-16 2016-08-17 中国石油集团川庆钻探工程有限公司地球物理勘探公司 VSP-CDP stacking method and three-dimensional VSP coverage number calculation method
CN107121700A (en) * 2017-02-28 2017-09-01 中国石油化工股份有限公司 The preparation method of structural map under a kind of thrust-nappe belt in front of the mountains

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7551516B2 (en) * 2005-03-09 2009-06-23 Aram Systems, Ltd. Vertical seismic profiling method utilizing seismic communication and synchronization
EP2715403A4 (en) * 2011-05-27 2015-01-28 Conocophillips Co Two-way wave equation targeted data selection for improved imaging of prospects among complex geologic structures
NO335878B1 (en) * 2013-06-10 2015-03-16 Read As Fiber optic and electrical seismic sensor cable for acquisition and transmission of information on seismic events recorded by multiple multicomponent geophones in a subsurface reservoir
CN106556860B (en) * 2015-09-29 2018-11-13 中国石油化工股份有限公司 The method and apparatus for laying VSP observation systems
CN105242315B (en) * 2015-10-20 2017-09-01 中国石油天然气集团公司 A kind of method and device of geological structure imaging
CN107290776A (en) * 2016-03-30 2017-10-24 中国石油化工股份有限公司 A kind of relief surface curved line design method centered on exciting
CN107340534A (en) * 2016-05-03 2017-11-10 中国石油化工股份有限公司 A kind of observation system shot point Optimization Design based on dimensional topography gradient slope aspect
CN109581495A (en) * 2018-10-24 2019-04-05 中国石油天然气集团有限公司 Mountainous region surface velocity model construction method and system

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102213770A (en) * 2010-04-01 2011-10-12 罗仁泽 Combined area excitation receiving method for seismic data in low signal-to-noise ratio area
CN102313900A (en) * 2011-09-01 2012-01-11 中国石油集团川庆钻探工程有限公司 Excitation position determining method of three-dimensional seismic acquisition observation system
CN105487120A (en) * 2014-09-18 2016-04-13 中国石油化工股份有限公司 Method for determining excitation well location of mountain front zone
CN105866833A (en) * 2016-06-16 2016-08-17 中国石油集团川庆钻探工程有限公司地球物理勘探公司 VSP-CDP stacking method and three-dimensional VSP coverage number calculation method
CN107121700A (en) * 2017-02-28 2017-09-01 中国石油化工股份有限公司 The preparation method of structural map under a kind of thrust-nappe belt in front of the mountains

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
《复杂山前带有偏VSP数据处理方法研究》;张山 等;《石油物探》;20120331;第51卷(第2期);3、164-171 *

Also Published As

Publication number Publication date
CN111399031A (en) 2020-07-10

Similar Documents

Publication Publication Date Title
CN111927552B (en) Sky-earth reconnaissance technology-based complex mountain area long and deep-buried tunnel reconnaissance method
CN106597533B (en) A kind of Depth Domain velocity modeling method for mountain front seism processing
CN102053270B (en) Sedimentary formation unit-based seismic facies analysis method
CN103809216B (en) A kind of resistivity data is combined speed with geological data and is built field method
CN102967882B (en) The modeling method of the interval velocity model on stratum
CN101315427A (en) A method and system for processing seismic exploration data in complex areas
CN102901985B (en) A kind of Depth Domain interval velocity modification method being applicable to relief surface
CN106094032B (en) A kind of method for building formation velocity model
CN105549087B (en) A kind of underground coal mine seam seismic exploration when walking and amplitude joint inversion method
CN109884710B (en) Micro-logging tomography method aiming at excitation well depth design
CN111399031B (en) A method of collecting and processing mountain seismic data
CN113376695A (en) Full waveform inversion method suitable for complex collapse column of coal seam floor
Bocin et al. Gravity and magnetic modelling in the Vrancea Zone, south-eastern Carpathians: Redefinition of the edge of the East European Craton beneath the south-eastern Carpathians
CN102798888B (en) Method for calculating velocity ratio of longitudinal wave to transverse wave by using non-zero wellhead distance data
CN101609163B (en) Multi-scale seismic data joint imaging method based on fluctuation theory
CN103076628A (en) Method for processing pre-stack time migration with optimized aperture
CN114398696B (en) Three-dimensional fine modeling method for alkaline granite uranium deposit
CN109085644A (en) True earth's surface imaging method when being walked based on dual-beam
CN108375794A (en) Based on the VSP fracture hole Diffraction Imaging technical methods symmetrically observed
Greenwood et al. Active seismic surveys for drilling target characterization in Ossola Valley: International Continental Scientific Drilling Program (ICDP) project Drilling the Ivrea–Verbano zonE (DIVE) phase I
CN107765341A (en) A kind of method for determining stratum residual density
Catchings et al. Faulting apparently related to the 1994 Northridge, California, earthquake and possible co-seismic origin of surface cracks in Potrero Canyon, Los Angeles County, California
CN109143398A (en) A kind of modeling method of automatic mesh chromatography Depth Domain speed
CN113589367A (en) Method for correcting structure trend near fracture based on area conservation
Xue et al. Seismic fine imaging and its significance for natural gas hydrate exploration in the Shenhu Test Area, South China Sea

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