CN112379443A - 纵横波可控震源微测井表层分析***及方法 - Google Patents
纵横波可控震源微测井表层分析***及方法 Download PDFInfo
- Publication number
- CN112379443A CN112379443A CN202011210111.7A CN202011210111A CN112379443A CN 112379443 A CN112379443 A CN 112379443A CN 202011210111 A CN202011210111 A CN 202011210111A CN 112379443 A CN112379443 A CN 112379443A
- Authority
- CN
- China
- Prior art keywords
- wave
- longitudinal
- excitation device
- transverse wave
- signals
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 48
- 239000002344 surface layer Substances 0.000 title claims abstract description 24
- 238000004458 analytical method Methods 0.000 title claims abstract description 17
- 230000005284 excitation Effects 0.000 claims abstract description 129
- 238000005070 sampling Methods 0.000 claims description 49
- 238000004590 computer program Methods 0.000 claims description 17
- 238000012937 correction Methods 0.000 claims description 11
- 238000005211 surface analysis Methods 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 claims description 9
- 241001669679 Eleotris Species 0.000 claims description 8
- 230000009471 action Effects 0.000 claims description 6
- 238000010586 diagram Methods 0.000 description 13
- 238000003860 storage Methods 0.000 description 12
- 238000011835 investigation Methods 0.000 description 10
- 239000010410 layer Substances 0.000 description 8
- 230000008569 process Effects 0.000 description 8
- 230000006870 function Effects 0.000 description 7
- 238000012545 processing Methods 0.000 description 7
- 238000010276 construction Methods 0.000 description 6
- 238000012360 testing method Methods 0.000 description 6
- 238000004891 communication Methods 0.000 description 5
- 238000004519 manufacturing process Methods 0.000 description 5
- 229910000831 Steel Inorganic materials 0.000 description 4
- 239000010959 steel Substances 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 3
- 239000002390 adhesive tape Substances 0.000 description 2
- 238000010009 beating Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000010892 electric spark Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000009191 jumping Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000010079 rubber tapping Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- VVNRQZDDMYBBJY-UHFFFAOYSA-M sodium 1-[(1-sulfonaphthalen-2-yl)diazenyl]naphthalen-2-olate Chemical compound [Na+].C1=CC=CC2=C(S([O-])(=O)=O)C(N=NC3=C4C=CC=CC4=CC=C3O)=CC=C21 VVNRQZDDMYBBJY-UHFFFAOYSA-M 0.000 description 1
- 230000001360 synchronised 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
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6169—Data from specific type of measurement using well-logging
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种纵横波可控震源微测井表层分析***及方法，所述***包括横波激发装置、纵波激发装置、地震记录仪器和井下三分量检波器；其中，所述横波激发装置用于发出横波扫描信号，所述纵波激发装置用于发出纵波扫描信号；所述地震记录仪器用于在所述横波激发装置和所述纵波激发装置发出横波扫描信号和纵波扫描信号的同时控制所述井下三分量检波器采集横波扫描信号和纵波扫描信号经过地层的采样信号，本发明可提高横波激发能量以提高横波表层调查深度。
Description
技术领域
本发明涉及地球物理勘探技术领域，尤其涉及一种纵横波可控震源微测井表层分析***及方法。
背景技术
在地球物理勘探领域，随着油气勘探深入，单一纵波勘探难以解决复杂构造及薄储层分辨率问题，需要开展纵横波联合勘探。横波表层调查是横波地震采集项目一项重要工作，也是解决近地表建模和静校正计算的前提和基础。柴达木盆地三湖地区从2001年先后开展转换波二维、转换波三维、纯横波二维、纵横波二维和纵横波三维地震勘探，其中微测井是提高横波表层调查精度的有效方法，尝试了枕木左右敲击微测井、电火花激发微测井和侧向钢板敲击小折射等多种横波微测井施工方式进行横波表层调查。但是，现有的横波表层调查方法存在激发能量弱的问题，限制了横波调查深度。
发明内容
本发明的一个目的在于提供一种纵横波可控震源微测井表层分析***，提高横波激发能量以提高横波表层调查深度。本发明的另一个目的在于提供一种纵横波可控震源微测井表层分析方法。本发明的再一个目的在于提供一种计算机设备。本发明的还一个目的在于提供一种可读介质。
为了达到以上目的，本发明一方面公开了一种纵横波可控震源微测井表层分析***，包括横波激发装置、纵波激发装置、地震记录仪器和井下三分量检波器；
其中，所述横波激发装置用于发出横波扫描信号，所述纵波激发装置用于发出纵波扫描信号；
所述地震记录仪器用于在所述横波激发装置和所述纵波激发装置发出横波扫描信号和纵波扫描信号的同时控制所述井下三分量检波器采集横波扫描信号和纵波扫描信号经过地层的采样信号。
优选的，所述横波激发装置包括横波震板及设于所述横波震板底部的横波可控震源；
所述横波可控震源用于根据预设横波扫描信号形成对应的振动力，所述横波震板在所述振动力的作用下形成所述横波扫描信号。
优选的，所述纵波激发装置包括纵波震板及设于所述纵波震板底部的纵波可控震源；
所述纵波可控震源用于根据预设纵波扫描信号形成对应的振动力，所述纵波震板在所述振动力的作用下形成所述纵波扫描信号。
优选的，所述地震记录仪器通过横波触发线与所述横波激发装置电连接，所述地震记录仪器通过纵波触发线与所述纵波激发装置电连接；
所述横波激发装置发出所述横波扫描信号后，所述横波激发装置形成横波触发信号并传输至所述地震记录仪器以使所述地震记录仪器同步采集对应的采样信号；
所述纵波激发装置发出所述纵波扫描信号后，所述纵波激发装置形成纵波触发信号并传输至所述地震记录仪器以使所述地震记录仪器同步采集对应的采样信号。
优选的，所述地震记录仪器进一步用于对所述采集信号进行时移校正，并对校正后的采集信号进行旋转计算得到R和T分量。
优选的，所述地震记录仪器进一步用于通过枕木敲击形成纵波和横波，并记录对应的初至拾取位置，根据所述初至拾取位置对所述采集信号进行时移校正。
优选的，所述地震记录仪器进一步用于对采集信号的横波信号的X和Y分量进行RT旋转，得到均方根振幅最大对应角度，根据所述角度对所述X和Y分量进行旋转计算得到R和T分量。
本发明还公开了一种纵横波可控震源微测井表层分析方法，包括：
通过横波激发装置发出横波扫描信号，通过纵波激发装置发出纵波扫描信号；
在所述横波激发装置和所述纵波激发装置发出横波扫描信号和纵波扫描信号的同时控制井下三分量检波器采集横波扫描信号和纵波扫描信号经过地层的采样信号。
优选的，进一步包括预先确定所述横波扫描信号和纵波扫描信号的步骤：
预设多个不同驱动幅度、扫描长度、扫描频率和斜坡长度的预设扫描信号，所述预设扫描信号包括横波扫描信号和纵波扫描信号；
采用多个预设扫描信号分别进行激发并采集所述采样信号；
抽取采集的多个采样信号的相同方向信号分量进行分析，综合考虑旁瓣、噪音及初至清晰度从多个预设扫描信号中确定所述预设横波扫描信号和预设纵波扫描信号。
优选的，进一步包括预先确定所述横波激发装置和所述纵波激发装置的位置的步骤：
针对多个不同井口距离分别进行采样信号的采集，其中，所述井口距离为所述横波激发装置和所述纵波激发装置的中心与微测井最近边缘的距离；
抽取不同采样信号相同方向的分量信号进行井口距离校正；
综合考虑井口距离校正时间和噪音确定最终井口距离以得到所述横波激发装置和所述纵波激发装置的位置。
优选的，进一步包括预先确定采样信号的采样率和记录长度的步骤：
根据微测井井深和纵横波速度估算满足纵横波微测井采集记录最小长度；
根据用户指令设定地震记录仪器的采样率；
根据所述采样率和所述最小长度确定采样信号的记录长度。
本发明还公开了一种计算机设备，包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，
所述处理器执行所述程序时实现如上所述方法。
本发明还公开了一种计算机可读介质，其上存储有计算机程序，
该程序被处理器执行时实现如上所述方法。
本发明提出了一种纵横波可控震源进行微测井激发***，通过横波激发装置发出横波扫描信号，所述纵波激发装置发出纵波扫描信号。在所述横波激发装置和所述纵波激发装置发出横波扫描信号和纵波扫描信号的同时控制所述井下三分量检波器采集横波扫描信号和纵波扫描信号经过地层的采样信号，解决表层调查激发能量和调查深度，取得高品质纵横波记录。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1示出本发明纵横波可控震源微测井表层分析***一个具体实施例的结构图；
图2示出现有技术中采集信号X分量的示意图；
图3示出本发明横波撞钟敲击微测井表层分析***一个具体实施例R分量的示意图；
图4示出本发明纵横波可控震源微测井表层分析方法一个具体实施例的流程图；
图5示出本发明纵横波可控震源微测井表层分析方法一个具体实施例S300的流程图；
图6示出本发明纵横波可控震源微测井表层分析方法一个具体实施例S400的流程图；
图7示出本发明纵横波可控震源微测井表层分析方法一个具体实施例S500的流程图；
图8示出适于用来实现本发明实施例的计算机设备的结构示意图。
附图说明：
1井下三分量检波器、2井下三分量检波器电缆、3井壁、4井口距离、5地面、6浅层地震记录仪器、7电缆接头、8触发接头、9钢丝绳、10横波触发线、11纵波触发线、12横波可控震源、13横波震板、14纵波可控震源、15纵波震板。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
根据本发明的一个方面，本实施例公开了一种纵横波可控震源12微测井表层分析***。如图1所示，本实施例中，所述***包括横波激发装置、纵波激发装置、地震记录仪器6和井下三分量检波器1。需要说明的是，地震记录仪器6和井下三分量检波器1可采用现有技术中的设备，也根据实际需求设置具有相同功能的设备，本发明对此并不作限定。
其中，所述横波激发装置用于发出横波扫描信号，所述纵波激发装置用于发出纵波扫描信号。
所述地震记录仪器6用于在所述横波激发装置和所述纵波激发装置发出横波扫描信号和纵波扫描信号的同时控制所述井下三分量检波器1采集横波扫描信号和纵波扫描信号经过地层的采样信号。
本发明提出了一种纵横波可控震源12进行微测井激发***，通过横波激发装置发出横波扫描信号，通过所述纵波激发装置发出纵波扫描信号。在所述横波激发装置和所述纵波激发装置发出横波扫描信号和纵波扫描信号的同时控制所述井下三分量检波器1采集横波扫描信号和纵波扫描信号经过地层形成的采样信号，解决表层调查激发能量和调查深度，取得高品质纵横波记录。
在优选的实施方式中，所述横波激发装置包括横波震板13及设于所述横波震板13底部的横波可控震源12。所述横波可控震源12用于根据预设横波扫描信号形成对应的振动力，所述横波震板13在所述振动力的作用下形成所述横波扫描信号。
在优选的实施方式中，所述纵波激发装置包括纵波震板15及设于所述纵波震板15底部的纵波可控震源14。所述纵波可控震源14用于根据预设纵波扫描信号形成对应的振动力，所述纵波震板15在所述振动力的作用下形成所述纵波扫描信号。
在优选的实施方式中，所述地震记录仪器6通过横波触发线10与所述横波激发装置电连接，所述地震记录仪器6通过纵波触发线11与所述纵波激发装置电连接。
其中，所述横波激发装置发出所述横波扫描信号后，所述横波激发装置形成横波触发信号并传输至所述地震记录仪器6以使所述地震记录仪器6同步采集对应的采样信号。
进一步的，在一个具体例子中，浅层地震记录仪器6的电缆接头7通过井下三分量检波器电缆2与设于井下的井下三分量检波器1电连接。三分量检波器1上连接有钢丝绳9，可通过钢丝绳9将三分量检波器1下放或上提。
所述纵波激发装置发出所述纵波扫描信号后，所述纵波激发装置形成纵波触发信号并传输至所述地震记录仪器6以使所述地震记录仪器6同步采集对应的采样信号。
可以理解的是，可在浅层地震记录仪器6与横波激发装置和纵波激发装置间连接触发线实现电连接，利用短路触发原理，实现了浅层地震记录仪器6与横波激发装置和纵波激发装置信号的同步采集。其中，在一个具体例子中，地震记录仪器6采用GDZ48A浅层地震仪。横波可控震源和纵波可控震源用于产生横波扫描信号和纵波扫描信号的激发信号为-3V脉冲电压，GDZ48A浅层地震仪触发电压>2.5V电压可以触发。因此，可通过触发线的形式将激发信号传输至横波可控震源电连接，具体的，可裁剪2根30m炮线作为浅层地震仪与可控震源的联机触发线。在实施时，分别将纵波可控震源14和横波可控震源12箱体TB信号与炮线进行连接，连接接头可用3M防电宽胶带进行包扎固定。然后将浅层地震仪触发接头8红色(正激性)线与纵横波可控震源12联机TB信号黄色(负激性)相接，将浅层地震仪触发接头8黄色(负激性)线与纵横波可控震源12联机TB信号红色(正激性)相接，连接接头用3M防电宽胶带进行包扎，这种连接方式与实际电路连接方式相反，主要可控震源TB信号是负相位信号，因此需要反极性相接。
在优选的实施方式中，还需进一步预先确定所述横波扫描信号和纵波扫描信号。具体的，可预设多个不同驱动幅度、扫描长度、扫描频率和斜坡长度的预设扫描信号，所述预设扫描信号包括横波扫描信号和纵波扫描信号。采用多个预设扫描信号分别进行激发并采集所述采样信号。抽取采集的多个采样信号的相同方向信号分量进行分析，综合考虑旁瓣、噪音及初至清晰度从多个预设扫描信号中确定所述预设横波扫描信号和预设纵波扫描信号，从而确定最优的扫描信号，提高表层分析深度。
其中，在一个具体例子中，设计可控震源扫描信号并将它加载到可控震源(横波可控震源和纵波可控震源)的震源箱体软件(VTTI-VE464)中，该软件可控制纵、横波可控震源发出纵、横波扫描信号，采用纵、横波可控震源自震模式单台单次激发和三分量井下检波器接收微测井采集，即可采用纵、横波可控震源进行单台单次激发和三分量井下检波器接收采样信号。
进一步的，在正式生产之前，优选不同地表钻5-10m井，将三分量检波器1下到井底，固定检波器，确保检波器与井壁3耦合，将检波器接到GDZ48A浅层地震仪。每一种扫描信号进行3次激发，进行纵横波微测井数据采集，将采集得到的采集信号设置在文件中，形成的多个文件的文件号按照顺序进行记录存储。抽取相同分量(包括Z、X、Y三个方向的分量)不同扫描信号激发采集得到的采集信号形成的单道记录，按照相同波形显示参数进行记录显示，并分析初至波波峰、波谷、旁瓣、主频、振幅、噪音一致性，优选旁瓣最少、噪音最弱、初至波清晰度最好的激发扫描信号经过地层形成的采集信号进行生产采集。通过参数试验分析，影响采集信号资料的最关键参数是扫描频和驱动幅度，扫描频率3-72Hz最佳，纵波驱动幅度20-40％，横波驱动幅度30-50％，整体来说纵横波微测井驱动幅度比地面5地震数据采集小，地面5地震数据采集驱动幅度为65-70％，主要较小驱动幅度有利于减少可控震源激发产生噪音。
在优选的实施方式中，还需进一步预先确定所述横波激发装置和所述纵波激发装置的位置。具体的，针对多个不同井口距离4分别进行采样信号的采集，其中，所述井口距离4为所述横波激发装置和所述纵波激发装置的中心与微测井最近边缘的距离。抽取不同采样信号相同方向的分量信号进行井口距离4校正。综合考虑井口距离4校正时间和噪音确定最终井口距离4以得到所述横波激发装置和所述纵波激发装置的位置。需要说明的是，对信号进行井口距离校正为本领域的常规技术手段，当然也可以采用其他技术手段对信号进行井口距离校正，本发明对此并不作限定。
可以理解的是，通过不同井口距离4数据采集测试，优选最佳井口距离4。其中，在一个具体例子中，井口距离4(可控震源激发震板中心到微测井井口距离4)是纵横波微测井关键参数，常规微测井井口距离4为3m,可控震源施工微测井井口距离4受到可控震源机械震动噪音影响，必须进行试验采集分析，优选适合可控震源激发井口距离4。一般可控震源激发装置的长度为10-11m，震板中心到激发装置靠近井口边缘的距离为5m，因此井口距离4必须大于5m，主要两方面考虑：1)从安全施工要求可控震源激发装置不能覆盖井口，容易出现安全事故；二是方便野外站在井口提取三分量检波器1操作人员在井口半径1-2m有足够空间进行野外操作。分别开展井口距离4m、6m、8m、10m、20m、30m、40m、50m纵横波微测井数据采集测试，每一种对比因素激发3次，抽取相同分量(Z、X、Y)不同因素激发单道记录，按照相同波形显示参数进行记录显示，并进行井口距离4校正，选择井口校正时间最小和噪音最弱的井口距离4作为可控震源井口距离4最佳参数。通过参数试验分析，确定可控震源微测井施工井口距离4为10m,这参数既影响资料最小又兼顾安全施工要求。
在优选的实施方式中，还需进一步预先确定采样信号的采样率和记录长度。其中，根据微测井井深和纵、横波速度估算满足纵、横波微测井采集信号最小长度。根据用户指令设定地震记录仪器6的采样率。根据所述采样率和所述最小长度确定采样信号的记录长度。
其中，在一个具体例子中，可采用GDZ48A浅层地震仪微测井采集采样率和记录长度。纵横波可控震源12激发信号为未相关信号，因此其记录长度＝扫描信号长度+实际记录信号长度。
Tr＝Ts+Tm (1)
Tr＝NLr×DT (2)(2)
Ts＝NLs×DT (3)(3)
Tm＝NLm×DT (4)(4)
其中，Tr和NLr为纵、横波微测井未相关记录长度和样点长度，Ts和NLs为纵、横波可控震源12扫描信号长度和样点长度；Tm为满足目的层纵、横波采集数据最小长度和样点长度；DT浅层地震仪微测井采集采样率。
根据设计井深和纵横波速度，初步估算满足纵横波微测井采集记录最小长度。
其中，H0、H1........Hm-1、Hm分布代表不同目的层厚度；V0、V1........Vm-1、Vm分布代表不同目的层速度。根据公式(5)并给定设计最大井深、不同目的层厚度和速度初步估算纵横波微测井采集记录最小长度。
根据GDZ48A浅层地震仪记录性能，合理设计微测井采集采样率和记录长度。GDZ48A浅层地震仪一般提供采样率设计窗口，自动计算出最大采样样点个数，主要有以下三种采样率参数菜单选择：[email protected]、[email protected]、[email protected]，以柴达木盆地三湖地区近地表结构，初步估算初Tm为1000ms，根据步骤二试验扫描长度至少为8000ms，同时考虑野外施工效率，一般扫描长度越长，每一个深度点采集越长，影响野外作业效率，因此综合选择扫描长度Ts为8000，Tr记录长度至少大于9000ms，因此我们选择[email protected]进行数据采集。
然后，根据前期准备和参数设计，正式对生产井进行纵横波微测井数据采集，按照深度20米以下深度间隔3米、10-20米深度间隔2米、5-10米深度间隔1米、0-5米以下深度间隔0.5米进行数据采集，三分量检波器1下到井底，采集从井底到井口顺序逐个深度点进行数据采集，每一深度纵波和横波各采集2次，文件号按照顺序1.SG2、2.SG2、3.SG2........进行记录。每激发1次产生未相关记录进行实时监控，重点监控记录波形面貌、频率、能量及噪音，及时发现有异常记录并重新采集，确保每激发1次记录正确性。利用可控震源扫描信号与纵横波微测井记录进行相关，得到最终纵横波微测井记录。
在优选的实施方式中，所述地震记录仪器6进一步用于对所述采集信号进行时移校正，并对校正后的采集信号进行旋转计算得到R和T分量。
在优选的实施方式中，所述地震记录仪器6进一步用于通过枕木敲击形成纵波和横波，并记录对应的初至拾取位置，根据所述初至拾取位置对所述采集信号进行时移校正。具体的，常规枕木敲击记录是脉冲信号，其初至拾取位置为初至起跳点，可控震源激发产生记录为零相位信号，正常情况拾取波峰或波谷，但是受到震源机械干扰及环境噪音等影响，初至位置不容易确定。在试验阶段进行相同位置激发点和相同深度检波器接收，分布采用枕木敲击激发和可控震源激发得到记录，根据枕木敲击激发记录拾取初至位置标定到可控震源激发记录，可以直接确定可控震源激发记录拾取初至位置，根据试验分析结果表明可控震源激发记录拾取初至位置为第一个强能量的波谷。
在优选的实施方式中，所述地震记录仪器6进一步用于对采集信号的横波信号的X和Y分量进行RT旋转，得到均方根振幅最大对应角度，根据所述角度对所述X和Y分量进行旋转计算得到R和T分量。
可以理解的是，可对纵横波可控震源12记录X和Y分量进行RT旋转，直接将两个横波分量旋转校正为快慢分量。井下三分量检波在提升过程中自由旋转并无法固定检波器方向，导致记录波一致性比较差，需要对井下三分量检波器1进行RT旋转。具体的，在一个具体例子中，RT旋转技术实现步骤如下：
提取记录X分量和Y分量，初步对X分量和Y分量初至拾取，并选择初至位置50ms时窗的X分量和Y分量，图2示出采集信号的X分量。
根据RT旋转公式(6)和(7)进行旋转角度扫描，从0-180度范围和角度增量为1度，逐度计算出旋转后X分量和Y分量均方根振幅ER和ET，输出均方根振幅最大对应角度Ang；
Rk(ti)＝X(ti)×cos(Angk)+Y(ti)×sin(Angk) (6)
Tk(ti)＝X(ti)×sin(Angk)+Y(ti)×cos(Angk) (7)
其中，X(ti)和Y(ti)分布代表X分量和Y分量，ti代表选取时间窗样点对应时间，M代表时窗样点个数，Angk代表旋转角度，其中范围0-180度，k代表角度增量序号；ERk代表旋转后R分量均方位根振幅，ETk代表旋转后T分量均方位根振幅。
根据计算角度对记录X分量和Y分量进行旋转计算得到R和T分量，即得到初至一致性比较快慢波记录，如图3所示，为后续纵横波微测井记录解释提供高品质地震记录。
本发明涉及纵横波表层调查方式，是一种有效解决深井纵横波微测井激发初至波能量的激发方法。利用纵横波可控震源12激发和井下三分量接收，并直接利用浅层地震仪与可控震源联机快速采集。主要采取以下针对性技术措施：1)通过制作浅层地震仪与可控震源触发联线，解决浅层地震仪与可控震源同步采集；2)采用纵横波可控震源12自震模式单台单次激发和三分量井下检波器接收微测井采集并通过试验科学优选最佳震源激发参数；3)科学设计可控震源激发井口距离4-10m；4)利用枕木敲击记录初至拾取位置，标定纵横波可震源记录并确定最终初至拾取位置；5)对可控震源激发微测井记录进行RT数据旋转，最终得到高品质微测井记录。
本发明适应纵横波地震勘探项目表层调查，解决了长期困扰井中微测井激发能量弱问题，该方法操作性强和适用性强，纵横波微测井记录信噪比高和初至波清晰度高，拾取初至位置精准可靠，为后续纵横波表层调查提供有效方法，具有较好的应用前景。
基于相同原理，本实施例还公开了一种纵横波可控震源12微测井表层分析方法。如图4所示，本实施例中，所述方法包括：
S100：通过横波激发装置发出横波扫描信号，通过纵波激发装置发出纵波扫描信号。
S200：在所述横波激发装置和所述纵波激发装置发出横波扫描信号和纵波扫描信号的同时控制所述井下三分量检波器1采集横波扫描信号和纵波扫描信号经过地层的采样信号。
在优选的实施方式中，如图5所示，所述方法进一步包括预先确定所述横波扫描信号和纵波扫描信号的步骤S300：
S310：预设多个不同驱动幅度、扫描长度、扫描频率和斜坡长度的预设扫描信号，所述预设扫描信号包括横波扫描信号和纵波扫描信号。
S320：采用多个预设扫描信号分别进行激发并采集所述采样信号。
S330：抽取采集的多个采样信号的相同方向信号分量进行分析，综合考虑旁瓣、噪音及初至清晰度从多个预设扫描信号中确定所述预设横波扫描信号和预设纵波扫描信号。
在优选的实施方式中，如图6所示，所述方法进一步包括预先确定所述横波激发装置和所述纵波激发装置的位置的步骤S400：
S410：针对多个不同井口距离4分别进行采样信号的采集，其中，所述井口距离4为所述横波激发装置和所述纵波激发装置的中心与微测井最近边缘的距离。
S420：抽取不同采样信号相同方向的分量信号进行井口距离4校正。
S430：综合考虑井口距离4校正时间和噪音确定最终井口距离4以得到所述横波激发装置和所述纵波激发装置的位置。
在优选的实施方式中，如图7所示，所述进一步包括预先确定采样信号的采样率和记录长度的步骤S500：
S510：根据微测井井深和纵横波速度估算满足纵横波微测井采集记录最小长度。
S520：根据用户指令设定地震记录仪器6的采样率。
S530：根据所述采样率和所述最小长度确定采样信号的记录长度。
由于该方法解决问题的原理与以上***类似，因此本方法的实施可以参见***的实施，在此不再赘述。
上述实施例阐明的***、装置、模块或单元，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。一种典型的实现设备为计算机设备，具体的，计算机设备例如可以为个人计算机、膝上型计算机、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。
在一个典型的实例中计算机设备具体包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述方法。
下面参考图8，其示出了适于用来实现本申请实施例的计算机设备600的结构示意图。
如图8所示，计算机设备600包括中央处理单元(CPU)601，其可以根据存储在只读存储器(ROM)602中的程序或者从存储部分608加载到随机访问存储器(RAM))603中的程序而执行各种适当的工作和处理。在RAM603中，还存储有***600操作所需的各种程序和数据。CPU601、ROM602、以及RAM603通过总线604彼此相连。输入/输出(I/O)接口605也连接至总线604。
以下部件连接至I/O接口605：包括键盘、鼠标等的输入部分606；包括诸如阴极射线管(CRT)、液晶反馈器(LCD)等以及扬声器等的输出部分607；包括硬盘等的存储部分608；以及包括诸如LAN卡，调制解调器等的网络接口卡的通信部分609。通信部分609经由诸如因特网的网络执行通信处理。驱动器610也根据需要连接至I/O接口605。可拆卸介质611，诸如磁盘、光盘、磁光盘、半导体存储器等等，根据需要安装在驱动器610上，以便于从其上读出的计算机程序根据需要被安装如存储部分608。
特别地，根据本发明的实施例，上文参考流程图描述的过程可以被实现为计算机软件程序。例如，本发明的实施例包括一种计算机程序产品，其包括有形地包含在机器可读介质上的计算机程序，所述计算机程序包括用于执行流程图所示的方法的程序代码。在这样的实施例中，该计算机程序可以通过通信部分609从网络上被下载和安装，和/或从可拆卸介质611被安装。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
为了描述的方便，描述以上装置时以功能分为各种单元分别描述。当然，在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。
本发明是参照根据本发明实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
本领域技术人员应明白，本申请的实施例可提供为方法、***或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于***实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
以上所述仅为本申请的实施例而已，并不用于限制本申请。对于本领域技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本申请的权利要求范围之内。
Claims (13)
1.一种纵横波可控震源微测井表层分析***，其特征在于，包括横波激发装置、纵波激发装置、地震记录仪器和井下三分量检波器；
其中，所述横波激发装置用于发出横波扫描信号，所述纵波激发装置用于发出纵波扫描信号；
所述地震记录仪器用于在所述横波激发装置和所述纵波激发装置发出横波扫描信号和纵波扫描信号的同时控制所述井下三分量检波器采集横波扫描信号和纵波扫描信号经过地层的采样信号。
2.根据权利要求1所述的纵横波可控震源微测井表层分析***，其特征在于，所述横波激发装置包括横波震板及设于所述横波震板底部的横波可控震源；
所述横波可控震源用于根据预设横波扫描信号形成对应的振动力，所述横波震板在所述振动力的作用下形成所述横波扫描信号。
3.根据权利要求1所述的纵横波可控震源微测井表层分析***，其特征在于，所述纵波激发装置包括纵波震板及设于所述纵波震板底部的纵波可控震源；
所述纵波可控震源用于根据预设纵波扫描信号形成对应的振动力，所述纵波震板在所述振动力的作用下形成所述纵波扫描信号。
4.根据权利要求1所述的纵横波可控震源微测井表层分析***，其特征在于，所述地震记录仪器通过横波触发线与所述横波激发装置电连接，所述地震记录仪器通过纵波触发线与所述纵波激发装置电连接；
所述横波激发装置发出所述横波扫描信号后，所述横波激发装置形成横波触发信号并传输至所述地震记录仪器以使所述地震记录仪器同步采集对应的采样信号；
所述纵波激发装置发出所述纵波扫描信号后，所述纵波激发装置形成纵波触发信号并传输至所述地震记录仪器以使所述地震记录仪器同步采集对应的采样信号。
5.根据权利要求1所述的纵横波可控震源微测井表层分析***，其特征在于，所述地震记录仪器进一步用于对所述采集信号进行时移校正，并对校正后的采集信号进行旋转计算得到R和T分量。
6.根据权利要求5所述的纵横波可控震源微测井表层分析***，其特征在于，所述地震记录仪器进一步用于通过枕木敲击形成纵波和横波，并记录对应的初至拾取位置，根据所述初至拾取位置对所述采集信号进行时移校正。
7.根据权利要求5所述的纵横波可控震源微测井表层分析***，其特征在于，所述地震记录仪器进一步用于对采集信号的横波信号的X和Y分量进行RT旋转，得到均方根振幅最大对应角度，根据所述角度对所述X和Y分量进行旋转计算得到R和T分量。
8.一种纵横波可控震源微测井表层分析方法，其特征在于，包括：
通过横波激发装置发出横波扫描信号，通过纵波激发装置发出纵波扫描信号；
在所述横波激发装置和所述纵波激发装置发出横波扫描信号和纵波扫描信号的同时控制井下三分量检波器采集横波扫描信号和纵波扫描信号经过地层的采样信号。
9.根据权利要求8所述的纵横波可控震源微测井表层分析方法，其特征在于，进一步包括预先确定所述横波扫描信号和纵波扫描信号的步骤：
预设多个不同驱动幅度、扫描长度、扫描频率和斜坡长度的预设扫描信号，所述预设扫描信号包括横波扫描信号和纵波扫描信号；
采用多个预设扫描信号分别进行激发并采集所述采样信号；
抽取采集的多个采样信号的相同方向信号分量进行分析，综合考虑旁瓣、噪音及初至清晰度从多个预设扫描信号中确定所述预设横波扫描信号和预设纵波扫描信号。
10.根据权利要求8所述的纵横波可控震源微测井表层分析方法，其特征在于，进一步包括预先确定所述横波激发装置和所述纵波激发装置的位置的步骤：
针对多个不同井口距离分别进行采样信号的采集，其中，所述井口距离为所述横波激发装置和所述纵波激发装置的中心与微测井最近边缘的距离；
抽取不同采样信号相同方向的分量信号进行井口距离校正；
综合考虑井口距离校正时间和噪音确定最终井口距离以得到所述横波激发装置和所述纵波激发装置的位置。
11.根据权利要求8所述的纵横波可控震源微测井表层分析方法，其特征在于，进一步包括预先确定采样信号的采样率和记录长度的步骤：
根据微测井井深和纵横波速度估算满足纵横波微测井采集记录最小长度；
根据用户指令设定地震记录仪器的采样率；
根据所述采样率和所述最小长度确定采样信号的记录长度。
12.一种计算机设备，包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，
所述处理器执行所述程序时实现如权利要求8-11任一项所述方法。
13.一种计算机可读介质，其上存储有计算机程序，其特征在于，
该程序被处理器执行时实现如权利要求8-11任一项所述方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011210111.7A CN112379443A (zh) | 2020-11-03 | 2020-11-03 | 纵横波可控震源微测井表层分析***及方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011210111.7A CN112379443A (zh) | 2020-11-03 | 2020-11-03 | 纵横波可控震源微测井表层分析***及方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112379443A true CN112379443A (zh) | 2021-02-19 |
Family
ID=74577908
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011210111.7A Pending CN112379443A (zh) | 2020-11-03 | 2020-11-03 | 纵横波可控震源微测井表层分析***及方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112379443A (zh) |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4809239A (en) * | 1987-07-14 | 1989-02-28 | Schlumberger Technology Corporation | Method for evaluating parameters related to the elastic properties of subsurface earth formations |
| US4907670A (en) * | 1981-01-20 | 1990-03-13 | Anstey Nigel A | Seismic exploration using compressional and shear waves simultaneously |
| US5901113A (en) * | 1996-03-12 | 1999-05-04 | Schlumberger Technology Corporation | Inverse vertical seismic profiling using a measurement while drilling tool as a seismic source |
| CN1888933A (zh) * | 2005-06-28 | 2007-01-03 | 中国石油集团东方地球物理勘探有限责任公司 | 横波或转换横波勘探近地表表层结构调查方法 |
| CN101419292A (zh) * | 2007-10-25 | 2009-04-29 | 中国石油天然气集团公司 | 采用纵波源多分量地震数据生成横波地震剖面的方法 |
| US20090185449A1 (en) * | 2008-01-18 | 2009-07-23 | Baker Hughes Incorporated | Method for 3-C 3-D Wavefield Decomposition |
| US20140126328A1 (en) * | 2012-11-06 | 2014-05-08 | Schlumberger Technology Corporation | Methods and Systems for Improving Microseismic Event Detection and Location |
| CN108267778A (zh) * | 2018-01-04 | 2018-07-10 | 中国铁路设计集团有限公司 | 一种地层速度测试***及测试方法 |
| CN110749927A (zh) * | 2019-11-19 | 2020-02-04 | 中油奥博(成都)科技有限公司 | 光纤声波传感正交偶极声波测井***及其测量方法 |
| CN111596347A (zh) * | 2019-02-21 | 2020-08-28 | 中国石油天然气集团有限公司 | 一种快速求取表层纵横波速度比的方法及装置 |
| CN111665568A (zh) * | 2020-07-09 | 2020-09-15 | 中油奥博(成都)科技有限公司 | 基于分布式光纤声波传感技术的微测井装置及测量方法 |
-
2020
- 2020-11-03 CN CN202011210111.7A patent/CN112379443A/zh active Pending
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4907670A (en) * | 1981-01-20 | 1990-03-13 | Anstey Nigel A | Seismic exploration using compressional and shear waves simultaneously |
| US4809239A (en) * | 1987-07-14 | 1989-02-28 | Schlumberger Technology Corporation | Method for evaluating parameters related to the elastic properties of subsurface earth formations |
| US5901113A (en) * | 1996-03-12 | 1999-05-04 | Schlumberger Technology Corporation | Inverse vertical seismic profiling using a measurement while drilling tool as a seismic source |
| CN1888933A (zh) * | 2005-06-28 | 2007-01-03 | 中国石油集团东方地球物理勘探有限责任公司 | 横波或转换横波勘探近地表表层结构调查方法 |
| CN101419292A (zh) * | 2007-10-25 | 2009-04-29 | 中国石油天然气集团公司 | 采用纵波源多分量地震数据生成横波地震剖面的方法 |
| US20090185449A1 (en) * | 2008-01-18 | 2009-07-23 | Baker Hughes Incorporated | Method for 3-C 3-D Wavefield Decomposition |
| US20140126328A1 (en) * | 2012-11-06 | 2014-05-08 | Schlumberger Technology Corporation | Methods and Systems for Improving Microseismic Event Detection and Location |
| CN108267778A (zh) * | 2018-01-04 | 2018-07-10 | 中国铁路设计集团有限公司 | 一种地层速度测试***及测试方法 |
| CN111596347A (zh) * | 2019-02-21 | 2020-08-28 | 中国石油天然气集团有限公司 | 一种快速求取表层纵横波速度比的方法及装置 |
| CN110749927A (zh) * | 2019-11-19 | 2020-02-04 | 中油奥博(成都)科技有限公司 | 光纤声波传感正交偶极声波测井***及其测量方法 |
| CN111665568A (zh) * | 2020-07-09 | 2020-09-15 | 中油奥博(成都)科技有限公司 | 基于分布式光纤声波传感技术的微测井装置及测量方法 |
Non-Patent Citations (6)
| Title |
|---|
| 于宝华;王海立;吴蔚;马立新;戴勇: "三湖横波静校正技术", 中国石油学会2019年物探技术研讨会论文集, 30 September 2019 (2019-09-30), pages 835 - 838 * |
| 唐建明: "转换波三维三分量地震勘探方法技术研究", 中国博士学位论文全文数据库 基础科学辑, no. 2011, 15 January 2011 (2011-01-15) * |
| 徐峰;尹成;李忠;潘树林;熊军;吕文彪;: "井中激发地表接收的近地表结构调查方法", 石油物探, vol. 48, no. 03, 31 May 2009 (2009-05-31), pages 294 - 298 * |
| 李桂花;朱光明;: "利用三分量微测井技术调查表层纵、横波速度", 石油地球物理勘探, no. 02, 30 April 2006 (2006-04-30), pages 160 - 165 * |
| 王海立;李光耀;张立军;赵荣艳;满红霞;王婷婷;翟俊伟;: "横波表层调查与静校正技术在三湖地区的应用", 青海石油, vol. 30, no. 03, 30 September 2012 (2012-09-30), pages 7 - 15 * |
| 赵邦六;: "多分量地震勘探在岩性气藏勘探开发中的应用", 石油勘探与开发, vol. 35, no. 04, 23 August 2008 (2008-08-23), pages 397 - 409 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7590491B2 (en) | Signal integration measure for seismic data | |
| CN104280775B (zh) | 一种基于全波形矢量偏移叠加的微地震监测定位方法 | |
| KR20200014387A (ko) | 지하 구조물의 검출 | |
| US11474271B2 (en) | Methods and systems for automated sonic imaging | |
| AU3302499A (en) | Single well system for mapping sources of acoustic energy | |
| EA026344B1 (ru) | Система и способ получения и обработки сейсмических данных о полях упругих волн | |
| US11604301B2 (en) | Methods and systems for automated sonic imaging | |
| Hobiger et al. | Site characterization of Swiss strong‐motion stations: The benefit of advanced processing algorithms | |
| CN116378648B (zh) | 一种基于随钻声波前视的近钻头地层探测方法及装置 | |
| CN111257941B (zh) | 一种组合式海底地震仪方位角自动识别装置及方法 | |
| US8639442B2 (en) | Identifying invalid seismic data | |
| US20170285195A1 (en) | Integrating vertical seismic profile data for microseismic anisotropy velocity analysis | |
| CN105735971A (zh) | 一种基于弹性波的钻孔深度检测***及其检测方法 | |
| CN105259581A (zh) | 一种地震资料时深转换方法 | |
| Cheng et al. | Testing of a permanent orbital surface source and distributed acoustic sensing for monitoring of unconventional reservoirs: Preliminary results from the Eagle Ford Shale | |
| AU2011279350B2 (en) | Method for accentuating specular and non-specular seismic events from within shallow subsurface rock formations | |
| Anthony et al. | Installation and performance of the Albuquerque Seismological Laboratory small‐aperture posthole array | |
| CN112379443A (zh) | 纵横波可控震源微测井表层分析***及方法 | |
| Abbas et al. | An Open‐Access Data Set of Active‐Source and Passive‐Wavefield DAS and Nodal Seismometer Measurements at the Newberry Florida Site | |
| Hobiger et al. | Site characterization in the framework of the renewal of the Swiss strong motion network (SSMNet) | |
| CN105929445B (zh) | 微破裂向量扫描方法 | |
| Phillips et al. | Explosion source models and the scattering origin of regional phases from SPE phase 1 coda spectral ratios | |
| CN113109870B (zh) | 一种高精度浅地层速度获取方法 | |
| Abukrat et al. | Directional sources and distributed acoustic sensing deployments for near-surface characterization | |
| WO2024133179A1 (en) | Method for ambient vibration analysis |
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 |