CN111077577B - 井地联合储层描述方法及装置 - Google Patents
井地联合储层描述方法及装置 Download PDFInfo
- Publication number
- CN111077577B CN111077577B CN201811228362.0A CN201811228362A CN111077577B CN 111077577 B CN111077577 B CN 111077577B CN 201811228362 A CN201811228362 A CN 201811228362A CN 111077577 B CN111077577 B CN 111077577B
- Authority
- CN
- China
- Prior art keywords
- seismic data
- model
- ground
- anisotropic
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 51
- 238000012545 processing Methods 0.000 claims abstract description 125
- 238000010521 absorption reaction Methods 0.000 claims abstract description 50
- 238000009792 diffusion process Methods 0.000 claims abstract description 48
- 238000000926 separation method Methods 0.000 claims abstract description 20
- 230000005012 migration Effects 0.000 claims abstract description 14
- 238000013508 migration Methods 0.000 claims abstract description 14
- 238000012937 correction Methods 0.000 claims description 41
- 238000007781 pre-processing Methods 0.000 claims description 9
- 238000001514 detection method Methods 0.000 claims description 6
- 230000010354 integration Effects 0.000 claims 1
- 238000012958 reprocessing Methods 0.000 abstract description 2
- 238000001914 filtration Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000005284 excitation Effects 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 4
- 230000003068 static effect Effects 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 239000010410 layer Substances 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 239000000523 sample Substances 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000002344 surface layer Substances 0.000 description 1
- 230000003313 weakening effect Effects 0.000 description 1
Images
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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
-
- 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
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (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
本发明实施例提供一种井地联合储层描述方法及装置，该方法包括：获取地震数据；根据井中地震数据进行数据处理，得到速度模型、各向异性模型；根据井中地震数据进行波场分离，得到下行纵波；根据下行纵波，得到球面扩散因子、吸收衰减因子以及反褶积算子；根据球面扩散因子、吸收衰减因子以及反褶积算子对地面地震数据进行补偿处理，得到补偿处理后的地面地震数据；根据补偿处理后的地面地震数据，更新速度模型以及各向异性模型；根据更新后的速度模型以及各向异性模型对补偿处理后的地面地震数据进行各向异性叠前深度偏移处理，再根据重新处理后的地面地震数据得到用于描述储层的描述信息，从而实现对储层的精细描述。
Description
技术领域
本发明涉及地震观测领域，尤其涉及一种井地联合储层描述方法及装置。
背景技术
地震勘探是指人工激发地震波，通过观测和分析人工地震产生的地震波在地下的传播规律，描述地下储层的性质和形态的地球物理勘探方法。
现有技术中经常使用地面地震数据实现对地下储层的描述，然而因为地面地震数据因近地表起伏、吸收衰减和干扰严重，地震波场特征受复杂构造和火成岩体等影响而失真，无法准确获取地层真实信息，从而无法实现对储层精确的描述。
发明内容
本发明实施例提供一种井地联合储层描述方法，以克服解决地面地震数据对进行储层描述时无法准确获取地层真实信息的问题。
第一方面，本发明实施例提供一种井地联合储层描述方法，包括：
获取地震数据，所述地震数据包括地面地震数据与井中地震数据；
根据所述井中地震数据进行数据处理，得到速度模型、各向异性模型，所述速度模型用于指示地震波在井中各地层中的传播速度，所述各向异性模型用于指示各地层的各向异性；
根据所述井中地震数据进行波场分离，得到下行纵波；
根据所述下行纵波，得到球面扩散因子、吸收衰减因子以及反褶积算子；
根据所述球面扩散因子、吸收衰减因子以及反褶积算子对所述地面地震数据进行补偿处理，得到补偿处理后的地面地震数据；
根据所述补偿处理后的地面地震数据，更新所述速度模型以及所述各向异性模型；
根据更新后的速度模型以及各向异性模型对所述补偿处理后的地面地震数据进行各向异性叠前深度偏移处理，得到用于描述储层的描述信息。
在一种可能的实现方式中，所述根据所述球面扩散因子、吸收衰减因子以及反褶积算子对所述地面地震数据进行补偿处理，得到补偿处理后的地面地震数据，包括：
依次根据所述球面扩散因子对所述地面地震数据进行球面扩散补偿处理、根据所述吸收衰减因子对所述地面地震数据进行吸收衰减补偿处理以及根据所述反褶积算子对所述地面地震数据进行反褶积处理，得到补偿处理后的地面地震数据。
在一种可能的实现方式中，所述根据所述补偿处理后的地面地震数据，更新所述速度模型以及所述各向异性模型，包括：
对所述补偿处理后的地面地震数据进行动校正，得到剩余动校正量；
根据所述剩余动校正量更新所述速度模型以及所述各向异性模型。
在一种可能的实现方式中，所述根据所述井中地震数据进行数据处理，得到速度模型、各向异性模型，包括：
对所述井中地震数据进行预处理，得到预处理后的井中地震数据；
根据所述井中地震数据获取直达波信息和各向异性信息，所述直达波信息包括直达波的初至时间和深度，所述各向异性信息包括各炮点至各检波点的各向异性时差；
根据所述直达波信息，得到地震波在井中各地层中的传播时间和深度，根据所述地震波在井中各地层中的传播时间和深度，得到地震波在各地层中的传播速度，得到速度模型；
根据所述各向异性信息，获取各深度点的弱各向异性参数，根据各深度点的弱各向异性参数，得到所述各向异性模型。
第二方面，本发明实施例提供一种井地联合储层描述装置，包括：
获取模块，用于获取地震数据，所述地震数据包括地面地震数据与井中地震数据；
处理模块，用于根据所述井中地震数据进行数据处理，得到速度模型、各向异性模型，所述速度模型用于指示地震波在井中各地层中的传播速度，所述各向异性模型用于指示各地层的各向异性；
分离模块，用于根据所述井中地震数据进行波场分离，得到下行纵波；
所述处理模块还用于，根据所述下行纵波，得到球面扩散因子、吸收衰减因子以及反褶积算子；
所述处理模块还用于，根据所述球面扩散因子、吸收衰减因子以及反褶积算子对所述地面地震数据进行补偿处理，得到补偿处理后的地面地震数据；
更新模块，用于根据所述补偿处理后的地面地震数据，更新所述速度模型以及所述各向异性模型；
所述处理模块还用于，根据更新后的速度模型以及各向异性模型对所述补偿处理后的地面地震数据进行各向异性叠前深度偏移处理，得到用于描述储层的描述信息。
在一种可能的实现方式中，所述处理模块具体用于：
依次根据所述球面扩散因子对所述地面地震数据进行球面扩散补偿处理、根据所述吸收衰减因子对所述地面地震数据进行吸收衰减补偿处理以及根据所述反褶积算子对所述地面地震数据进行反褶积处理，得到补偿处理后的地面地震数据。
在一种可能的实现方式中，所述更新模块具体用于：
对所述补偿处理后的地面地震数据进行动校正，得到剩余动校正量；
根据所述剩余动校正量更新所述速度模型以及所述各向异性模型。
在一种可能的实现方式中，所述处理模块还具体用于：
对所述井中地震数据进行预处理，得到预处理后的井中地震数据；
根据所述井中地震数据获取直达波信息和各向异性信息，所述直达波信息包括直达波的初至时间和深度，所述各向异性信息包括各炮点至各检波点的各向异性时差；
根据所述直达波信息，得到地震波在井中各地层中的传播时间和深度，根据所述地震波在井中各地层中的传播时间和深度，得到地震波在各地层中的传播速度，得到速度模型；
根据所述各向异性信息，获取各深度点的弱各向异性参数，根据各深度点的弱各向异性参数，得到所述各向异性模型。
第三方面，本发明实施例提供一种井地联合储层描述设备，包括：
存储器，用于存储程序；
处理器，用于执行所述存储器存储的所述程序，当所述程序被执行时，所述处理器用于执行如上第一方面以及第一方面各种可能的设计中任一所述的方法。
第四方面，本发明实施例提供一种计算机可读存储介质，包括指令，当其在计算机上运行时，使得计算机执行如上第一方面以及第一方面各种可能的设计中任一所述的方法。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的井地联合储层描述方法的流程图一；
图2为本发明实施例提供的井地联合储层描述方法的流程图二；
图3为本发明实施例提供的井地联合储层描述装置的结构示意图；
图4为为本发明实施例提供的井地联合储层描述设备的硬件结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
图1为本发明实施例提供的井地联合储层描述方法的流程图一，如图1所示，该方法包括：
S101、获取地震数据，地震数据包括地面地震数据与井中地震数据。
在本发明实施例中，获取地震数据的方法优选为垂直地震剖面方法，垂直地震剖面方法是一种地震观测方法，在地表附近的一些点上激发地震波，在井中不同深度的检波点上放置检波器，通过对检波器获取到的地震数据进行分析从而解释解决地质问题，其中激发地震波的点称为炮点。
具体的，垂直地震剖面方法对应的观测系统可以分为两种，一种是零偏移距观测系统，其得到的地震数据为零偏移距地震数据，零偏移距观测系统是指，当地震波激发点和井口的水平距离小于等于预设值时建立的地震数据观测系统，另一种是非零偏移距观测系统，其得到的数据为非零偏移距地震数据，非零偏移距观测系统是指，当地震波激发点和井口的水平距离大于预设值时建立的地震数据观测系统，本发明实施例的观测系统优选为零偏移距观测系统，对应得到的地震数据为零偏移距地震数据。
具体的，地震数据为井中放置的检波器检测得到的数据资料，例如可以包括初至直达波、一次反射波、上行波、下行波等，地震数据具体还可以包括井中地震数据和地面地震数据，井中地震数据是指在井中设置的检波器检测得到的地震数据，而地面地震数据是指在井口或地面放置的检波器检测得到的地震数据。
S102、根据井中地震数据进行数据处理，得到速度模型、各向异性模型，速度模型用于指示地震波在井中各地层中的传播速度，各向异性模型用于指示各地层的各向异性。
通过检波器获取得到地震数据之后，需要首先对井中地震数据进行数据处理，具体的，井中不同深度对应的地层的构造以及成分都不相同，因此地震波在井中不同深度对应的各地层的传播速度也不相同，对井中地震数据进行数据处理得到地震波在井中各地层的传播速度，建立起速度模型，在本发明实施例中速度模型用于指示地震波在井中各地层中的传播速度。
进一步地，还需要对井中地震数据进行处理建立各向异性模型，各向异性是指物质的全部或部分化学、物理等性质随着方向的改变而有所变化，在不同的方向上呈现出差异的性质，上述介绍了井中不同深度对应的地层的构造以及成分都不相同，因此井中不同深度的对应的各地层的各向异性也不同，通常采用各向异性参数来描述各向异性，各向异性参数例如可以为强度、渗透率、传播速度等，通过对井中地震数据进行数据处理得到井中不同深度对应的各地层的各向异性参数，从而建立起各向异性模型。
S103、根据井中地震数据进行波场分离，得到下行纵波。
波场分离是对井中地震数据进行分析中的一个关键环节，垂直地震剖面方法获取到的井中地震数据包含了大量的波场信息，并且各种波场相互交织，需要对复杂波长进行高效精确的波场分离得到相对准确的单一波场才有利于地震数据的下一步处理。
地震波可以分为上行波和下行波，上行波是指来自接收点下方的波场信息，而下行波是指来自接收点上方的波场信息，上行波又具体可以分为上行横波、上行纵波，下行波具体可以分为下行横波和下行纵波，其中横波纵波是波的两种类型，波就是振动的传播，在本发明实施例中就是地震波通过介质传播。横波是质点的振动方向与波的传播方向垂直，纵波是质点的振动方向与波的传播方向平行的波。
通过对井中地震数据进行波场分离处理，可以得到下行纵波。
S104、根据下行纵波，得到球面扩散因子、吸收衰减因子以及反褶积算子。
S105、根据球面扩散因子、吸收衰减因子以及反褶积算子对地面地震数据进行补偿处理，得到补偿处理后的地面地震数据。
其中，地震波在地下介质中传播时，由于球面扩散和地层吸收衰减，会导致地震波的振幅减弱，在地震数据处理的过程中需要对地震数据进行补偿，从而消除因球面扩散和地层吸收衰减导致的地震波的振幅减弱带来的影响。
进一步地，地层对地震波的吸收作用相当于对地震波进行了一次低通滤波，通常将该滤波称为大地滤波，此时得到的地震数据由于受复杂子波的作用和干扰的影响，分辨能力较低，地质界面上各反射波互相叠加之后相当复杂，难以解释，因此还需要对地震数据进行反褶积处理，反褶积的基本作用是压缩地震记录中的地震子波，同时可以压制鸣震和多次波，因而反褶积处理可以明显提高地震数据的垂直分辨率。
具体的，在对地震波进行补偿处理之前，需要首先对下行纵波进行处理得到球面扩散因子、吸收衰减因子以及反褶积算子，其中球面扩散因子就是用来消除球面扩散对地震波带来的影响的因子，吸收衰减因子就是用来消除地层吸收衰减对地震波带来影响的因子，反褶积算子就是用来进行反褶积处理的算子。
接着利用球面扩散因子、吸收衰减因子以及反褶积算子对地面地震数据进行补偿处理，消除球面扩散、地层吸收衰减以及褶积作用对地震波产生的影响，得到补偿处理后的地面地震数据。
S106、根据补偿处理后的地面地震数据，更新速度模型以及各向异性模型。
具体的，通过上述步骤得到了速度模型以及各向异性模型，但是上述步骤中是根据井中地震数据处理后得到的，此时得到的速度模型以及各向异性模型未经补偿处理，对于地震波在各地层中的传播速度以及在各地层中的各向异性的描述会存在一定的偏差，因此根据补偿处理后的地面地震数据，进一步更新速度模型以及各向异性模型，使得速度模型以及各向异性模型对于速度和各向异性的指示更加精确。
S107、根据更新后的速度模型以及各向异性模型对补偿处理后的地面地震数据进行各向异性叠前深度偏移处理，得到用于描述储层的描述信息。
利用更新后的速度模型以及各向异性模型对补偿处理后的地面地震数据进行各向异性叠前深度偏移，具体的，各向异性叠前深度偏移是用于确定地下复杂构造形态和位置的技术手段，得到用于描述储层的描述信息，描述信息例如可以为相干和曲率属性、蚂蚁体属性等，再进一步利用重新处理后的地面地震数据进行叠后波阻抗反演，叠后波阻抗反演是指利用地震信息得到地下地质信息的过程，从而利用得到的储层的描述信息实现对地下储层的精细描述。
本发明实施例提供的井地联合储层描述方法，包括：获取地震数据，地震数据包括地面地震数据与井中地震数据；根据井中地震数据进行数据处理，得到速度模型、各向异性模型，速度模型用于指示地震波在井中各地层中的传播速度，各向异性模型用于指示各地层的各向异性；根据井中地震数据进行波场分离，得到下行纵波；根据下行纵波，得到球面扩散因子、吸收衰减因子以及反褶积算子；根据球面扩散因子、吸收衰减因子以及反褶积算子对地面地震数据进行补偿处理，得到补偿处理后的地面地震数据；根据补偿处理后的地面地震数据，更新速度模型以及各向异性模型；根据更新后的速度模型以及各向异性模型对补偿处理后的地面地震数据进行各向异性叠前深度偏移处理，得到用于描述储层的描述信息。能够反映出更加精确的时深关系，并利用速度模型以及各向异性模型对地面地震数据重新处理以弥补地面地震的不足，提高地震数据成像的精度，从而实现对储层更加准确精细的描述。
在上述实施例的基础上，结合图2对井地联合储层描述方法进行进一步地详细介绍，图2为本发明实施例提供的井地联合储层描述方法的流程图二，如图2所示，该方法包括：
S201、获取地震数据，地震数据包括地面地震数据与井中地震数据。
此步骤的实现方式与S101类似，此处不再赘述。
S202、对井中地震数据进行预处理，得到预处理后的井中地震数据。
垂直地震剖面方法获得的原始井中地震数据存在干扰强、存在大量表层和层间多次波影响等问题，因此必须对原始的井中地震数据进行预处理，使之变成可处理的地震数据，具体的，对井中地震数据进行数据解编、成像基准面确定、观测系统建立、数据分选、野外静校正和去噪处理。
进一步地，数据解编是指将采集的井中地震数据转换为处理器可识别的格式，成像基准面确定是指在进行地震数据分析时，首先要将所有的地震数据都基于同一个基准面来进行处理，因此基准面的确定也就影响着地震数据处理的结果。其次需要进行观测系统的建立，目的是为模拟野外，定义一个相对坐标系，将野外的地震激发点以及检波器放置点放到这个相对的坐标系中，并进一步进行数据分选处理，数据分选是指将地震数据对应的地震道排列到一起，其中地震道是指在每个观测点上记录地震波，都必须经过检波器、放大系统和记录系统三个基本环节，它们连在一起总称为地震道。
其次还需要进行野外静校正处理，具体的，需要假设地面水平，并且将地震波激发点以及检波器校正至同一水平面上，该水平面即为静校正基准面，其目的是将地形、低速带以及爆炸深度等因素对地震波传播时间带来的影响加以消除，从而满足地表水平、表层介质均匀的几何地震学假设条件。再进一步进行去噪处理，从而去除噪声对地震波产生的干扰。在完成一系列数据预处理之后，得到预处理后的井中地震数据。
S203、根据井中地震数据获取直达波信息和各向异性信息，直达波信息包括直达波的初至时间和深度，各向异性信息包括各炮点至各检波点的各向异性时差。
S204、根据直达波信息，得到地震波在井中各地层中的传播时间和深度，根据地震波在井中各地层中的传播时间和深度，得到地震波在各地层中的传播速度，得到速度模型。
S205、根据各向异性信息，获取各深度点的弱各向异性参数，根据各深度点的弱各向异性参数，得到各向异性模型。
在得到预处理后的井中地震数据之后，对预处理后的井中地震数据进行进一步地数据处理得到速度模型以及各向异性模型。
具体的，对井中地震数据的直达波中拾取初至信息，其中地震波可分为直达波、一次反射波和多次反射波，直达波是指在均匀地层中由震源直接传播到检波点的地震波。直达波信息中包含直达波的初至时间和该直达波到达的对应地层的深度，进一步利用时间和深度的关系求取井中各地层对应的平均速度以及层速度，其中平均速度是指地震波垂直穿过一组水平层状介质各层的总厚度与总的传播时间之比，层速度是指两个地层之间的速度，再以地层的大致范围为约束建立起速度模型，因此速度模型用于指示地震波在井中各地层中的传播速度。
可选地，利用求取的层速度制作合成记录，并进行精细的井震标定，将地震剖面上的同相轴与地下地质界面建立准确的对应关系，为构造解释和储层精细描述奠定基础。
进一步地，需要求取各向异性模型，各向异性模型首先需要在具有垂直对称轴的横向各向同性介质弱各向异性假设条件下，具体的，上述介绍了各向异性，其中六方各向异性介质又称为横向各向同性介质，是地震勘探中使用最广的介质，由于横向各向同性介质具有一个无限次的对称轴，表现出一定程度上的各向同性性质，所以它的这种性质被称为弱各向异性。
在该假设条件下，利用井中地震数据求取初至时差和动校正速度，其中动校正是将炮点和检波点之间距离不同的各地震道上，来自同一界面同一点的反射波的到达时间校正为共中心点处的回声时间，从而正确反映地层的真深度，进行动校正所用的速度即为动校正速度。
在该基础上求其各炮点至各检波点的各向异性时差，其中各向异性时差是指在各向异性介质中地震波沿不同方向的传播速度会存在差异，求取各炮点至各检波点的各向异性时差后，在最小平方的原则下，求取井中各地层对应的各深度点的各向异性参数，并在速度模型的约束下得到各向异性模型，因此各向异性模型用于指示各地层的各向异性。
S206、根据井中地震数据进行波场分离，得到下行纵波。
具体的，地震波具有不同波场其偏振方向也不同的特点，首先通过坐标旋转和射线追踪将接收的波场初步分离为上行纵波—下行横波波场，以及下行纵波—上行横波波场，其中坐标旋转是指检波器的水平分量方位是随机的，因此接收到的记录在振幅和相位上出现差别很大的非一致性随机变化，所以需要进行坐标旋转进行方位校正处理，其中射线追踪方法作为一种快速有效的波场近似计算方法能够有效的实现波场分离。
进一步地，基于上行波场和下行波场的传播方向不同，再利用F-K滤波实现最终波场分离，其中F-K滤波是一种应用广泛的二维滤波方法。
S207、根据下行纵波，得到球面扩散因子、吸收衰减因子以及反褶积算子。
S208、依次根据球面扩散因子对地面地震数据进行球面扩散补偿处理、根据吸收衰减因子对地面地震数据进行吸收衰减补偿处理以及根据反褶积算子对地面地震数据进行反褶积处理，得到补偿处理后的地面地震数据。
对波场分离得到的下行纵波进行处理，首先基于地震波不同频率成分在传播过程中呈现球面扩散和吸收衰减差异，即低频吸收衰减小而高频吸收衰减大的特性，将下行纵波转换到频率域进行低频滤波处理。
具体的，针对某频率范围内的地震数据分析其频率随深度的变化，从而得到一系列不同深度点的球面扩散因子的值，之后统计出一个全局的球面扩散因子值作为该深度段的球面扩散因子。进一步地，并通过计算震源谱与接收谱频率差异求取吸收衰减因子，将下行纵波沿初至方向拉平并叠加，得到一个单道子波数据即为反褶积算子。
在计算得到球面扩散因子、吸收衰减因子以及反褶积算子之后，依次根据球面扩散因子对地面地震数据进行球面扩散补偿处理、根据吸收衰减因子对地面地震数据进行吸收衰减补偿处理以及根据反褶积算子对地面地震数据进行反褶积处理，得到补偿处理后的地面地震数据。
S209、对补偿处理后的地面地震数据进行动校正，得到剩余动校正量。
S210、根据剩余动校正量更新速度模型以及各向异性模型。
对补偿处理后的地面地震数据进行动校正，动校正在上述实施例中以及进行了介绍，此处不再赘述，得到剩余动校正量，其中剩余动校正是指在进行动校正的过程中，动校正速度难以精确地提取，因此动校正也不可避免的会引入一些误差，为了解决这些误差就需要进行剩余动校正，剩余动校正量即为进行剩余动校正所需要的动校正量。
根据剩余动校正量更新上述获得的速度模型以及各向异性模型，使得速度模型以及各向异性模型能够更加精确的描述井中不同深度各地层对应的速度以及各向异性参数。
S211、根据更新后的速度模型以及各向异性模型对补偿处理后的地面地震数据进行各向异性叠前深度偏移处理，得到用于描述储层的描述信息。
根据处理得到的地面地震数据和更新后的速度模型以及各向异性模型对地面地震进行各向异性叠前深度偏移，得到各向异性叠前深度偏移处理后的地面地震数据，从中提取均方根振幅属性、相干和曲率属性以及断层自动识别(Automatic Fault Extraction，AFE)和蚂蚁体属性，其中均方根振幅属性能够反映储层平面展布，相干和曲率属性能够反映储层断裂，AFE和蚂蚁体属性能够反映储层裂缝，并进一步进行叠后波阻抗反演实现对储层的精细描述。
本发明实施例提供的井地联合储层描述方法，在上述实施例的基础上，包括：获取地震数据；对井中地震数据进行预处理；根据井中地震数据获取直达波信息和各向异性信息；根据直达波信息，得到地震波在井中各地层中的传播时间和深度，根据地震波在井中各地层中的传播时间和深度，得到地震波在各地层中的传播速度，得到速度模型；根据各向异性信息，获取各深度点的弱各向异性参数，根据各深度点的弱各向异性参数，得到各向异性模型。根据井中地震数据进行波场分离，得到下行纵波。根据下行纵波，得到球面扩散因子、吸收衰减因子以及反褶积算子。依次根据球面扩散因子对地面地震数据进行球面扩散补偿处理、根据吸收衰减因子对地面地震数据进行吸收衰减补偿处理以及根据反褶积算子对地面地震数据进行反褶积处理，得到补偿处理后的地面地震数据。对补偿处理后的地面地震数据进行动校正，得到剩余动校正量。根据剩余动校正量更新速度模型以及各向异性模型。根据更新后的速度模型以及各向异性模型对补偿处理后的地面地震数据进行各向异性叠前深度偏移处理，得到用于描述储层的描述信息。能够反映出更加精确的时深关系，并利用速度模型以及各向异性模型对地面地震数据重新处理以弥补地面地震的不足，提高地震数据成像的精度，从而实现对储层更加准确精细的描述。
在实现对储层精细描述之后，可利用描述储层实现钻井轨迹优化、老井治理、井位部署等一系列应用。
图3为本发明实施例提供的井地联合储层描述装置的结构示意图。如图3所示，该井地联合储层描述装置30包括：获取模块301、处理模块302、分离模块303以及更新模块304。
获取模块301，用于获取地震数据，地震数据包括地面地震数据与井中地震数据；
处理模块302，用于根据井中地震数据进行数据处理，得到速度模型、各向异性模型，速度模型用于指示地震波在井中各地层中的传播速度，各向异性模型用于指示各地层的各向异性；
分离模块303，用于根据井中地震数据进行波场分离，得到下行纵波；
处理模块302还用于，根据下行纵波，得到球面扩散因子、吸收衰减因子以及反褶积算子；
处理模块302还用于，根据球面扩散因子、吸收衰减因子以及反褶积算子对地面地震数据进行补偿处理，得到补偿处理后的地面地震数据；
更新模块304，用于根据补偿处理后的地面地震数据，更新速度模型以及各向异性模型；
处理模块302还用于，根据更新后的速度模型以及各向异性模型对补偿处理后的地面地震数据进行各向异性叠前深度偏移处理，得到用于描述储层的描述信息。
可选地，处理模块302还具体用于：依次根据球面扩散因子对地面地震数据进行球面扩散补偿处理、根据吸收衰减因子对地面地震数据进行吸收衰减补偿处理以及根据反褶积算子对地面地震数据进行反褶积处理，得到补偿处理后的地面地震数据。
可选地，更新模块304具体用于：对补偿处理后的地面地震数据进行动校正，得到剩余动校正量；
根据剩余动校正量更新速度模型以及各向异性模型。
可选地，处理模块302还具体用于：
对井中地震数据进行预处理，得到预处理后的井中地震数据；
根据井中地震数据获取直达波信息和各向异性信息，直达波信息包括直达波的初至时间和深度，各向异性信息包括各炮点至各检波点的各向异性时差；
根据直达波信息，得到地震波在井中各地层中的传播时间和深度，根据地震波在井中各地层中的传播时间和深度，得到地震波在各地层中的传播速度，得到速度模型；
根据各向异性信息，获取各深度点的弱各向异性参数，根据各深度点的弱各向异性参数，得到各向异性模型。
本发明实施例提供的装置，可用于执行上述方法实施例的技术方案，其实现原理和技术效果类似，本实施例此处不再赘述。
图4为为本发明实施例提供的井地联合储层描述设备的硬件结构示意图，如图4所示，本发明实施例提供的井地联合储层描述设备40包括：处理器401以及存储器402；其中
存储器402，用于存储计算机执行指令；
处理器401，用于执行存储器存储的计算机执行指令，以实现上述实施例中井地联合储层描述装置所执行的各个步骤。具体可以参见前述方法实施例中的相关描述。
可选地，存储器402既可以是独立的，也可以跟处理器401集成在一起。
当存储器402独立设置时，该语音交互设备还包括总线403，用于连接存储器402和处理器401。
本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质中存储有计算机执行指令，当处理器执行所述计算机执行指令时，实现如上第一终端所执行的路线同步方法。
在本发明所提供的几个实施例中，应该理解到，所揭露的设备和方法，可以通过其它的方式实现。例如，以上所描述的设备实施例仅仅是示意性的，例如，所述模块的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个模块可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或模块的间接耦合或通信连接，可以是电性，机械或其它的形式。
上述以软件功能模块的形式实现的集成的模块，可以存储在一个计算机可读取存储介质中。上述软件功能模块存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)或处理器(英文：processor)执行本申请各个实施例所述方法的部分步骤。
应理解，上述处理器可以是中央处理单元(英文：Central Processing Unit，简称：CPU)，还可以是其他通用处理器、数字信号处理器(英文：Digital Signal Processor，简称：DSP)、专用集成电路(英文：Application Specific Integrated Circuit，简称：ASIC)等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合发明所公开的方法的步骤可以直接体现为硬件处理器执行完成，或者用处理器中的硬件及软件模块组合执行完成。
存储器可能包含高速RAM存储器，也可能还包括非易失性存储NVM，例如至少一个磁盘存储器，还可以为U盘、移动硬盘、只读存储器、磁盘或光盘等。
总线可以是工业标准体系结构(Industry Standard Architecture，ISA)总线、外部设备互连(Peripheral Component，PCI)总线或扩展工业标准体系结构(ExtendedIndustry Standard Architecture，EISA)总线等。总线可以分为地址总线、数据总线、控制总线等。为便于表示，本申请附图中的总线并不限定仅有一根总线或一种类型的总线。
上述存储介质可以是由任何类型的易失性或非易失性存储设备或者它们的组合实现，如静态随机存取存储器(SRAM)，电可擦除可编程只读存储器(EEPROM)，可擦除可编程只读存储器(EPROM)，可编程只读存储器(PROM)，只读存储器(ROM)，磁存储器，快闪存储器，磁盘或光盘。存储介质可以是通用或专用计算机能够存取的任何可用介质。
本领域普通技术人员可以理解：实现上述各方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成。前述的程序可以存储于一计算机可读取存储介质中。该程序在执行时，执行包括上述各方法实施例的步骤；而前述的存储介质包括：ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (8)
1.一种井地联合储层描述方法，其特征在于，包括：
获取地震数据，所述地震数据包括地面地震数据与井中地震数据；
根据所述井中地震数据进行数据处理，得到速度模型、各向异性模型，所述速度模型用于指示地震波在井中各地层中的传播速度，所述各向异性模型用于指示各地层的各向异性；
根据所述井中地震数据进行波场分离，得到下行纵波；
根据所述下行纵波，得到球面扩散因子、吸收衰减因子以及反褶积算子；
根据所述球面扩散因子、吸收衰减因子以及反褶积算子对所述地面地震数据进行补偿处理，得到补偿处理后的地面地震数据；
根据所述补偿处理后的地面地震数据，更新所述速度模型以及所述各向异性模型；
根据更新后的速度模型以及各向异性模型对所述补偿处理后的地面地震数据进行各向异性叠前深度偏移处理，得到用于描述储层的描述信息；
所述根据所述补偿处理后的地面地震数据，更新所述速度模型以及所述各向异性模型，包括：
对所述补偿处理后的地面地震数据进行动校正，得到剩余动校正量；
根据所述剩余动校正量更新所述速度模型以及所述各向异性模型。
2.根据权利要求1所述的方法，其特征在于，所述根据所述球面扩散因子、吸收衰减因子以及反褶积算子对所述地面地震数据进行补偿处理，得到补偿处理后的地面地震数据，包括：
依次根据所述球面扩散因子对所述地面地震数据进行球面扩散补偿处理、根据所述吸收衰减因子对所述地面地震数据进行吸收衰减补偿处理以及根据所述反褶积算子对所述地面地震数据进行反褶积处理，得到补偿处理后的地面地震数据。
3.根据权利要求1所述的方法，其特征在于，所述根据所述井中地震数据进行数据处理，得到速度模型、各向异性模型，包括：
对所述井中地震数据进行预处理，得到预处理后的井中地震数据；
根据所述井中地震数据获取直达波信息和各向异性信息，所述直达波信息包括直达波的初至时间和深度，所述各向异性信息包括各炮点至各检波点的各向异性时差；
根据所述直达波信息，得到地震波在井中各地层中的传播时间和深度，根据所述地震波在井中各地层中的传播时间和深度，得到地震波在各地层中的传播速度，得到速度模型；
根据所述各向异性信息，获取各深度点的弱各向异性参数，根据各深度点的弱各向异性参数，得到所述各向异性模型。
4.一种井地联合储层描述装置，其特征在于，包括：
获取模块，用于获取地震数据，所述地震数据包括地面地震数据与井中地震数据；
处理模块，用于根据所述井中地震数据进行数据处理，得到速度模型、各向异性模型，所述速度模型用于指示地震波在井中各地层中的传播速度，所述各向异性模型用于指示各地层的各向异性；
分离模块，用于根据所述井中地震数据进行波场分离，得到下行纵波；
所述处理模块还用于，根据所述下行纵波，得到球面扩散因子、吸收衰减因子以及反褶积算子；
所述处理模块还用于，根据所述球面扩散因子、吸收衰减因子以及反褶积算子对所述地面地震数据进行补偿处理，得到补偿处理后的地面地震数据；
更新模块，用于根据所述补偿处理后的地面地震数据，更新所述速度模型以及所述各向异性模型；
所述处理模块还用于，根据更新后的速度模型以及各向异性模型对所述补偿处理后的地面地震数据进行各向异性叠前深度偏移处理，得到用于描述储层的描述信息；
所述更新模块具体用于：
对所述补偿处理后的地面地震数据进行动校正，得到剩余动校正量；
根据所述剩余动校正量更新所述速度模型以及所述各向异性模型。
5.根据权利要求4所述的装置，其特征在于，所述处理模块具体用于：
依次根据所述球面扩散因子对所述地面地震数据进行球面扩散补偿处理、根据所述吸收衰减因子对所述地面地震数据进行吸收衰减补偿处理以及根据所述反褶积算子对所述地面地震数据进行反褶积处理，得到补偿处理后的地面地震数据。
6.根据权利要求4所述的装置，其特征在于，所述处理模块还具体用于：
对所述井中地震数据进行预处理，得到预处理后的井中地震数据；
根据所述井中地震数据获取直达波信息和各向异性信息，所述直达波信息包括直达波的初至时间和深度，所述各向异性信息包括各炮点至各检波点的各向异性时差；
根据所述直达波信息，得到地震波在井中各地层中的传播时间和深度，根据所述地震波在井中各地层中的传播时间和深度，得到地震波在各地层中的传播速度，得到速度模型；
根据所述各向异性信息，获取各深度点的弱各向异性参数，根据各深度点的弱各向异性参数，得到所述各向异性模型。
7.一种井地联合储层描述设备，其特征在于，包括：
存储器，用于存储程序；
处理器，用于执行所述存储器存储的所述程序，当所述程序被执行时，所述处理器用于执行如权利要求1至3中任一所述的方法。
8.一种计算机可读存储介质，包括指令，当其在计算机上运行时，使得计算机执行如权利要求1至3中任一所述的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811228362.0A CN111077577B (zh) | 2018-10-22 | 2018-10-22 | 井地联合储层描述方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811228362.0A CN111077577B (zh) | 2018-10-22 | 2018-10-22 | 井地联合储层描述方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111077577A CN111077577A (zh) | 2020-04-28 |
| CN111077577B true CN111077577B (zh) | 2022-02-01 |
Family
ID=70309703
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811228362.0A Active CN111077577B (zh) | 2018-10-22 | 2018-10-22 | 井地联合储层描述方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111077577B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112180465B (zh) * | 2020-11-03 | 2022-06-21 | 中国石油化工股份有限公司 | 一种井震联合的地层展布确定方法 |
| CN115567102B (zh) * | 2022-09-23 | 2024-06-11 | 中国联合网络通信集团有限公司 | 衰减确定方法、装置、设备及存储介质 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102736107A (zh) * | 2011-04-07 | 2012-10-17 | 中国石油天然气股份有限公司 | 一种能量约束非均质储层厚度识别系统 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN100383557C (zh) * | 2004-06-25 | 2008-04-23 | 大庆油田有限责任公司 | 一种提高地震分辨率的方法 |
| US9562983B2 (en) * | 2014-04-17 | 2017-02-07 | Saudi Arabian Oil Company | Generating subterranean imaging data based on vertical seismic profile data |
| WO2016105765A1 (en) * | 2014-12-23 | 2016-06-30 | Halliburton Energy Services, Inc. | Determining p-wave azimuthal anisotropy from walkaround vsp with offset dependent slowness corrections |
| CN107765306B (zh) * | 2016-08-19 | 2019-08-27 | 中国石油化工股份有限公司 | 一种vsp初始速度建模方法及装置 |
| CN107643542B (zh) * | 2017-08-28 | 2019-07-02 | 中煤科工集团西安研究院有限公司 | 基于井地联合的煤矿采空区三维弹性波层析探测方法 |
| CN108375794B (zh) * | 2018-01-22 | 2020-06-09 | 上海锦迪软件开发有限公司 | 基于对称观测的vsp缝洞绕射成像技术方法 |
-
2018
- 2018-10-22 CN CN201811228362.0A patent/CN111077577B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102736107A (zh) * | 2011-04-07 | 2012-10-17 | 中国石油天然气股份有限公司 | 一种能量约束非均质储层厚度识别系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111077577A (zh) | 2020-04-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109425896B (zh) | 白云岩油气储层分布预测方法及装置 | |
| CN109581501B (zh) | 用于沙漠区深度域速度建模的方法 | |
| US11243317B2 (en) | Optimal survey design | |
| Sairam et al. | Influence of local site effects in the Ahmedabad mega city on the damage due to past earthquakes in northwestern India | |
| CN102176054A (zh) | 近地表综合信息处理解释方法 | |
| CN102819038B (zh) | 一种碳酸盐岩内油水识别的方法及系统 | |
| Brodic et al. | Three-component seismic land streamer study of an esker architecture through S-and surface-wave imaging | |
| CN111077577B (zh) | 井地联合储层描述方法及装置 | |
| Malagnini et al. | Estimating absolute site effects | |
| CN107179553B (zh) | 基于双压实规律的沙漠区表层静校正方法 | |
| CN109471162A (zh) | 层间多次波处理方法、系统、电子设备及可读介质 | |
| US20100118655A1 (en) | Progressive 3D vertical seismic profiling method | |
| CN105425300B (zh) | 一种剩余静校正方法 | |
| EP4390466A1 (en) | Vsp-based level calibration method and apparatus for depth-domain seismic profile | |
| CN112630841B (zh) | 一种微地震事件检测分析方法 | |
| Wang | Full-azimuth, high-density, 3D single-point seismic survey for shale gas exploration in a loess plateau area, Southeast of Ordos Basin, China | |
| Maufroy et al. | Travel time inversion from ground level to gallery: protocol for the characterization of P‐wave seismic signature in a fractured‐porous Urgonian platform at hectometric scale | |
| Gosar | A microtremor HVSR study of the seismic site effects in the area of the town of Brežice (SE Slovenia) | |
| CN105989240A (zh) | 砂层预测方法及装置 | |
| CN112824939B (zh) | 地震数据的剩余静校正方法、装置、设备及可读存储介质 | |
| CN104237939B (zh) | 一种三维地震观测系统多频率偏移噪声获取方法及装置 | |
| CN108535779B (zh) | 一种用于地震勘探数据的优化方法 | |
| CN113267810B (zh) | 地震勘探全深度速度建模方法及装置 | |
| CN110716230B (zh) | 一种井地联合微地震定位方法 | |
| CN115774288A (zh) | 横波矢量反演建模方法及装置 |
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 |