CN109490965A - 一种定量评价地层非均匀性的方法及装置 - Google Patents
一种定量评价地层非均匀性的方法及装置 Download PDFInfo
- Publication number
- CN109490965A CN109490965A CN201811196745.4A CN201811196745A CN109490965A CN 109490965 A CN109490965 A CN 109490965A CN 201811196745 A CN201811196745 A CN 201811196745A CN 109490965 A CN109490965 A CN 109490965A
- Authority
- CN
- China
- Prior art keywords
- wave
- velocity
- quality factor
- random medium
- shear wave
- 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.)
- Granted
Links
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/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
- G01V2210/6222—Velocity; travel time
-
- 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/626—Physical property of subsurface with anisotropy
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)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明提出一种定量评价地层非均匀性的方法及装置,所述方法包括:获取深度区间内波形数据;进行深度校正,利用波形相干叠加法处理得到地层的纵波速度和横波速度;计算波形数据的散射衰减品质因子和中心频率;获取平均纵波速度、横波速度和横波速度的扰动量;建立随机介质模型;选取不同的相关距离,利用时域有限差分数值模拟所述随机介质模型的声场相应特征,并计算在不同相关距离下所述随机介质模型的散射衰减品质因子;建立所述随机介质模型的散射衰减品质因子图版;将实测数据散射衰减品质因子与散射衰减品质因子图版进行对比,获得所述深度区间的相关距离。本发明通过获取岩石的相关距离来定量评价岩石的非均匀性,属于油气资源勘探领域。
Description
技术领域
本发明属于油气资源勘探领域,具体涉及一种定量评价地层非均匀性的方法及装置。
背景技术
地球从地壳、地幔到地核,普遍存在非均匀性,而非均匀性往往意味着地质构造运动以及油气和矿产资源的存在,因此利用现有技术手段对地层非均匀性进行定量评价至关重要。
不同性质的非均匀性对地震波有不同的影响,速度和密度的非均匀性可引起波形变化、相位/走时和振幅的起伏以及直达波的视衰减,并导致尾波和前驱波。例如:岩石层的非均匀性会产生纵波(P)和横波(S)尾波、以及二者之间转换波的尾波。当地下介质存在裂缝,溶孔(如缝洞型碳酸盐岩储层)等非均匀地质体时,也会对地震波传播产生散射。不同尺度和不同组成的非均匀体会引起不同形式和不同强度的声波散射,因此可以从散射波场中判断介质的非均匀性并推断非均匀体的分布和性质。目前,尾波已成为定量评价岩石大尺度(地震波尺度)非均匀性的重要技术手段,利用尾波的衰减等信息可以获取岩石的相关距离,进而根据相关距离定量评价岩石的非均匀性。而小尺度(测井尺度)岩石的非均匀性的研究较少,目前仅是根据偶极声波测井数据定性识别岩石的非均匀性,无法有效的定量评价。
发明内容
针对现有井周小尺度非均匀性评价存在的不足,本发明提出一种利用偶极声波测井数据对地层非均匀性定量评价的方法及装置,为油气资源勘探领域地层评价提供一种实用可行的方法和技术手段。
本发明一种定量评价地层非均匀性的方法,包括如下技术方案:
S1、在深度区间内进行阵列声波测井,获取所述深度区间内单极、偶极阵列声波测井波形数据;
S2、对所述波形数据进行深度校正,获得处理深度位置处的单极、偶极阵列声波测井全波列数据,利用波形相干叠加法处理得到地层的纵波速度和横波速度;
S3、计算所述偶极阵列声波测井全波列数据的散射衰减品质因子和中心频率;
S4、对所述纵波速度和横波速度做统计分析,获取平均纵波速度、横波速度和横波速度的扰动量;
S5、利用所述平均纵波速度、横波速度和横波速度的扰动量建立随机介质模型;
S6、选取不同的相关距离,使所述相关距离与波数的乘积ka在预设数值区间范围内,利用时域有限差分数值模拟所述随机介质模型的声场相应特征,并计算在不同相关距离下所述随机介质模型的散射衰减品质因子;
S7、建立不同ka时所述随机介质模型的散射衰减品质因子图版;
S8、将步骤S3所述偶极阵列声波测井数据的散射衰减品质因子与步骤S7所述随机介质模型的散射衰减品质因子图版进行对比,获取偶极阵列声波测井数据的ka值,进而获得所述深度区间的相关距离。
在上述技术方案中,所述步骤S2中,所述利用波形相干叠加法处理得到地层的纵波速度和横波速度的具体过程为:
对整个声波波形或者波形中的一个时段以及给定的速度区间按公式1计算出二维相关函数Corr(v,T):
其中,Xm(t)是N个声波测井仪器接收换能器阵列中的第m个接收换能器,d是声波测井仪器接收换能器之间的间隔,T是时间窗Tw的位置,v是速度区间中的速度值;
计算当相关函数Corr(v,T)取极大值时对应的v值,求出纵波速度和横波速度。
在上述技术方案中,步骤S6所述相关距离为所述随机介质模型的自相关长度参数,横向、纵向的自相关长度相等。
在上述技术方案中,步骤S6所述利用时域有限差分数值模拟所述随机介质模型的声场相应特征的过程中,声源的激发频率为步骤S3中所述中心频率。
本发明还提供一种定量评价地层非均匀性的装置,所述装置包括:
测井模块:用于在深度区间内进行阵列声波测井,获取所述深度区间内单极、偶极阵列声波测井波形数据;
处理模块:用于对所述波形数据进行深度校正,获得处理深度位置处的单极、偶极阵列声波测井全波列数据,利用波形相干叠加法处理得到地层的纵波速度和横波速度;
计算模块:用于计算所述偶极阵列声波测井全波列数据的散射衰减品质因子和中心频率,对所述纵波速度和横波速度做统计分析,获取平均纵波速度、横波速度和横波速度的扰动量;
模拟模块:用于利用所述平均纵波速度、横波速度和横波速度的扰动量建立随机介质模型,选取不同的相关距离,使所述相关距离与波数的乘积ka在预设数值区间范围内,利用时域有限差分数值模拟所述随机介质模型的声场相应特征,并计算在不同相关距离下所述随机介质模型的散射衰减品质因子;
对比模块:用于建立不同ka时所述随机介质模型的散射衰减品质因子图版,将计算模块所述偶极阵列声波测井数据的散射衰减品质因子与所述随机介质模型的散射衰减品质因子图版进行对比,获取偶极阵列声波测井数据的ka值,进而获得所述深度区间的相关距离。
本发明采用偶极声波测井数据来对地层非均匀性进行定量评价,可产生两方面的效益:其一是与已有的非均匀性评价方法相比,该方法可以依据数值模拟图版获得地层的相关距离,可对地层非均匀性进行定量评价;其二是该方法可以对测井尺度下的地层非均匀性进行定量评价,对于地层多尺度非均匀性评价具有重要的学术和应用价值,是地震尺度非均匀性评价的有益补充。
附图说明
为了更清楚地说明本发明的技术方案,下面将对本发明技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的方法流程示意图;
图2为本发明实施例提供的装置结构示意图;
图3为本发明提供的某压裂井阵列声波测井实测数据的统计分析图;
图4位本发明实施例提供的随机散射衰减品质因子图版。
具体实施方式
下面结合附图对本发明实施例中的技术方案进行清楚、完整地描述,显然,下面所描述的实施例仅仅是本发明一部分实施例,而非全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
请参见图1,本发明一种定量评价地层非均匀性的方法,包括如下步骤:
S1、在深度区间内进行阵列声波测井,获取所述深度区间内单极、偶极阵列声波测井波形数据;
选取合适的阵列声波测井仪器及相关设备进行声波测井,获取深度区间内的声波阵列测井波形数据。
S2、进行深度校正,获得处理深度位置处的单极、偶极阵列声波测井全波列数据,利用波形相干叠加法处理得到地层的纵波速度和横波速度;
阵列声波测井获取到的声波的波形数据需要做一系列的数据预处理才能进行进一步的计算,由于仪器采样点与测量点之间有个距离偏移量,需要进行深度校正才能获得处理深度位置处的正确波形数据,波形数据预处理还可以包括增益控制、滤波、均衡化处理等。经过数据预处理后,利用波形相干叠加法快速有效第确定地层的弹性波速度剖面,得到地层纵波速度和横波速度。
可选的,步骤S2所述利用波形相干叠加法处理得到地层的纵波速度和横波速度的具体过程为:
对整个声波波形或者波形中的一个时段以及预先给定的速度区间按计算出二维相关函数Corr(v,T):
其中,Xm(t)是N个声波测井仪器接收换能器阵列中的第m个接收换能器,d是声波测井仪器接收换能器之间的间隔,T是时间窗Tw的位置,v是速度区间中的速度值;计算当相关函数Corr(v,T)取极大值时对应的v值,求出纵波速度和横波速度。
S3、计算偶极阵列声波测井数据的散射衰减品质因子和中心频率;
散射衰减品质因子的计算可使用传统方法如频谱比法、中心频率偏移法、振幅衰减法等,按需选取。
S4、对所述纵波速度和横波速度做统计分析,获取平均纵波速度、横波速度和横波速度的扰动量;
S5、利用所述平均纵波速度、横波速度和横波速度的扰动量建立随机介质模型;
假定纵波和横波速度的相对扰动是相同的,而密度的相对扰动与其有线性关系,从而可以只用一个参数比如横波速度的相对扰动量来描述随机介质在小尺度上的非均匀性。随机介质模型的自相关函数包括指数型、高斯型、混合型等,可根据实际需求选择合适的自相关函数用来建立随机介质模型,自相关函数的参数通常有两个,分别为介质在水平方向、深度方向的自相关长度a、b,它们分别描述随机介质在水平方向、深度方向上的非均匀异常的平均尺度。比如指数椭圆形自相关函数的公式为:
其中a、b为分别为x、z方向上的自相关长度,r为粗糙度。在本实施例中,取a=b,即在此随机介质模型中的相关距离即为a。
S6、选取不同的相关距离,使所述相关距离与波数的乘积ka在指定数值区间内,利用时域有限差分数值模拟所述随机介质模型的声场相应特征,并计算在不同相关距离下所述随机介质模型的散射衰减品质因子;
可选的,步骤S6所述相关距离为所述随机介质模型的自相关长度参数,横向、纵向的自相关长度相等。
相关距离即为步骤S6所述相关距离为所述随机介质模型的自相关长度参数a,且a=b,波数根据频率计算得到,设k为波数,则相关距离与波数的乘积ka=k*a,预设的数值区间取[0.01,10],即选取的相关距离与波数的乘积应落在区间[0.01,10]内。
可选的,步骤S6所述利用时域有限差分数值模拟所述随机介质模型的声场相应特征的过程中,声源的激发频率为步骤S3中所述中心频率。
S7、建立不同ka时所述随机介质模型的散射衰减品质因子图版;
S8、将步骤S3所述偶极阵列声波测井数据的散射衰减品质因子与步骤S7所述随机介质模型的散射衰减品质因子图版进行对比,获取偶极阵列声波测井数据的ka值,进而获得所述深度区间的相关距离。
参见图2,本发明还提供一种定量评价地层非均匀性的装置,所述装置包括:
测井模块210:用于在深度区间内进行阵列声波测井,获取所述深度区间内单极、偶极阵列声波测井波形数据;
处理模块220:用于进行深度校正,获得处理深度位置处的单极、偶极阵列声波测井全波列数据,利用波形相干叠加法处理得到地层的纵波速度和横波速度;
计算模块230:用于计算所述偶极阵列声波测井全波列数据的散射衰减品质因子和中心频率,对所述纵波速度和横波速度做统计分析,获取平均纵波速度、横波速度和横波速度的扰动量;
模拟模块240:用于利用所述平均纵波速度、横波速度和横波速度的扰动量建立随机介质模型,选取不同的相关距离,使所述相关距离与波数的乘积ka在预设数值区间范围内,利用时域有限差分数值模拟所述随机介质模型的声场相应特征,并计算在不同相关距离下所述随机介质模型的散射衰减品质因子;
对比模块250:用于建立不同ka时所述随机介质模型的散射衰减品质因子图版,将计算模块230所述偶极阵列声波测井数据的散射衰减品质因子与所述随机介质模型的散射衰减品质因子图版进行对比,获取偶极阵列声波测井数据的ka值,进而获得所述深度区间的相关距离。
下面结合图3给出的某压裂井阵列声波测井实测数据的统计分析图对本发明提出的方法作出充分说明。
图3给出了压裂前、后测井曲线及横波时差统计分析图。其中第一道给出了深度区间,第二道给出了井径和自然伽马曲线,第三道给出了压裂前、后的横波慢度,第四道给出了压裂前、后的各向异性大小。从图中可以看出,在压裂层段(X635m~X775m)压裂后的慢度明显增加,与压裂前几乎为零的各向异性(左曲线)相比,压裂后压裂层段井壁附近的岩石呈现较强的各向异性(右曲线),这是压裂后井壁附近的形成的裂缝引起的。
为了分析其非均匀性大小,根据本发明提出的方法,选取压裂层段时差数据统计了压裂后的横波慢度的分布区间,压后的横波慢度扰动约为6%~7%,计算得到的实测数据散射衰减品质因子QS=262。根据统计得到的平均纵波速度、均横波速度及横波扰动量建立随机介质模型,利用有限差分数值模拟了不同相关距离时的偶极子声波测井阵列波形,模拟时采用的声源中心频率为3kHz,所需地层弹性参数可以由测井资料获得,地层弹性参数扰动量为6%。根据所述随机介质模型,计算不同相关距离时的散射衰减,并建立图版。图4给出了所建立的随机介质模型的散射衰减品质因子图版,即散射衰减与波数与相关距离的乘积ka(ka=k*a,其中k为波数,a为相关距离)的对应关系。从图中可以看出当ka≈1时随机介质的散射衰减最大,这与理论预测结果一致。实测数据的散射衰减品质因子QS=262,它对应散射衰减的值约为0.004,比较接近于最大衰减值,对应的ka≈1.5,由此估计出来的压裂后地层的相关距离为0.2m。根据相关距离即可定量评价地层的非均匀性。
以上所述,以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (5)
1.一种定量评价地层非均匀性的方法,其特征在于,所述方法包括:
S1、在深度区间内进行阵列声波测井,获取所述深度区间内单极、偶极阵列声波测井波形数据;
S2、对所述波形数据进行深度校正,获得处理深度位置处的单极、偶极阵列声波测井全波列数据,利用波形相干叠加法处理得到地层的纵波速度和横波速度;
S3、计算所述偶极阵列声波测井全波列数据的散射衰减品质因子和中心频率;
S4、对所述纵波速度和横波速度做统计分析,获取平均纵波速度、横波速度和横波速度的扰动量;
S5、利用所述平均纵波速度、横波速度和横波速度的扰动量建立随机介质模型;
S6、选取不同的相关距离,使所述相关距离与波数的乘积ka在预设数值区间范围内,利用时域有限差分数值模拟所述随机介质模型的声场相应特征,并计算在不同相关距离下所述随机介质模型的散射衰减品质因子;
S7、建立不同ka时所述随机介质模型的散射衰减品质因子图版;
S8、将步骤S3所述偶极阵列声波测井数据的散射衰减品质因子与步骤S7所述随机介质模型的散射衰减品质因子图版进行对比,获取偶极阵列声波测井数据的ka值,进而获得所述深度区间的相关距离。
2.根据权利要求1所述一种定量评价地层非均匀性的方法,其特征在于,所述步骤S2中,所述利用波形相干叠加法处理得到地层的纵波速度和横波速度的具体过程为:
对整个声波波形或者波形中的一个时段以及给定的速度区间按公式1计算出二维相关函数Corr(v,T):
其中,t是时间,Xm(t)是N个声波测井仪器接收换能器阵列中的第m个接收换能器,d是声波测井仪器接收换能器之间的间隔,T是时间窗Tw的位置,v是所述速度区间中的速度值;
计算当相关函数Corr(v,T)取极大值时对应的v值,求出纵波速度和横波速度。
3.根据权利要求1所述一种定量评价地层非均匀性的方法,其特征在于,步骤S6所述相关距离即为所述随机介质模型的自相关长度参数,横向、纵向的自相关长度相等。
4.根据权利要求1所述一种定量评价地层非均匀性的方法,其特征在于,步骤S6所述利用时域有限差分数值模拟所述随机介质模型的声场相应特征的过程中,声源的激发频率为步骤S3中所述中心频率。
5.一种定量评价地层非均匀性的装置,其特征在于,所述装置包括:
测井模块:用于在深度区间内进行阵列声波测井,获取所述深度区间内单极、偶极阵列声波测井波形数据;
处理模块:用于对所述波形数据进行深度校正,获得处理深度位置处的单极、偶极阵列声波测井全波列数据,利用波形相干叠加法处理得到地层的纵波速度和横波速度;
计算模块:用于计算所述偶极阵列声波测井全波列数据的散射衰减品质因子和中心频率,对所述纵波速度和横波速度做统计分析,获取平均纵波速度、横波速度和横波速度的扰动量;
模拟模块:用于利用所述平均纵波速度、横波速度和横波速度的扰动量建立随机介质模型,选取不同的相关距离,使所述相关距离与波数的乘积ka在预设数值区间范围内,利用时域有限差分数值模拟所述随机介质模型的声场相应特征,并计算在不同相关距离下所述随机介质模型的散射衰减品质因子;
对比模块:用于建立不同ka时所述随机介质模型的散射衰减品质因子图版,将计算模块所述偶极阵列声波测井数据的散射衰减品质因子与所述随机介质模型的散射衰减品质因子图版进行对比,获取偶极阵列声波测井数据的ka值,进而获得所述深度区间的相关距离。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811196745.4A CN109490965B (zh) | 2018-10-15 | 2018-10-15 | 一种定量评价地层非均匀性的方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811196745.4A CN109490965B (zh) | 2018-10-15 | 2018-10-15 | 一种定量评价地层非均匀性的方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109490965A true CN109490965A (zh) | 2019-03-19 |
CN109490965B CN109490965B (zh) | 2020-09-01 |
Family
ID=65690240
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811196745.4A Active CN109490965B (zh) | 2018-10-15 | 2018-10-15 | 一种定量评价地层非均匀性的方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109490965B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111236933A (zh) * | 2020-01-15 | 2020-06-05 | 中国科学院声学研究所 | 利用声波测井仪器对地层径向速度进行分步反演的方法 |
CN111624661A (zh) * | 2020-04-30 | 2020-09-04 | 中国石油大学(北京) | 基于微震散射波品质因子的压裂效果评价方法及装置 |
CN112034516A (zh) * | 2020-08-17 | 2020-12-04 | 中国石油大学(华东) | 一种评价井周横波三维各向异性的方法 |
CN112444871A (zh) * | 2019-08-30 | 2021-03-05 | 中国石油化工股份有限公司 | 基于散射波地震响应特征的裂缝间距定量分析方法和设备 |
CN112444865A (zh) * | 2019-08-28 | 2021-03-05 | 中国石油化工股份有限公司 | 勘探缝洞储层的方法、***、装置及存储介质 |
CN112859163A (zh) * | 2021-01-07 | 2021-05-28 | 中国石油大学(北京) | 利用散射波确定压裂区地层品质因子变化的方法和装置 |
CN112986247A (zh) * | 2021-04-29 | 2021-06-18 | 中国科学院地质与地球物理研究所 | 基于矿物含量和分布的岩石脆性分析方法、***及设备 |
CN113187457A (zh) * | 2021-05-27 | 2021-07-30 | 中海油田服务股份有限公司 | 一种确定储层压裂效果的方法、装置和存储介质 |
CN113406697A (zh) * | 2020-03-17 | 2021-09-17 | 中国石油天然气股份有限公司 | 储层的横波品质因子获取方法和装置 |
CN113866829A (zh) * | 2021-09-28 | 2021-12-31 | 中国石油大学(北京) | 一种三维散射声波远探测扫描成像测井方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102478670A (zh) * | 2010-11-24 | 2012-05-30 | 中国石油天然气集团公司 | 一种利用地震衰减属性进行储层流体性质预测的方法 |
WO2012134657A2 (en) * | 2011-03-31 | 2012-10-04 | Chevron U.S.A. Inc. | System and method for processing seismic data |
CN103233727A (zh) * | 2013-05-13 | 2013-08-07 | 中国石油大学(华东) | 一种反演地层横波速度径向剖面的方法 |
CN106814397A (zh) * | 2016-12-21 | 2017-06-09 | 长江大学 | 一种多参数联合反演计算岩石散射衰减的方法 |
-
2018
- 2018-10-15 CN CN201811196745.4A patent/CN109490965B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102478670A (zh) * | 2010-11-24 | 2012-05-30 | 中国石油天然气集团公司 | 一种利用地震衰减属性进行储层流体性质预测的方法 |
WO2012134657A2 (en) * | 2011-03-31 | 2012-10-04 | Chevron U.S.A. Inc. | System and method for processing seismic data |
CN103233727A (zh) * | 2013-05-13 | 2013-08-07 | 中国石油大学(华东) | 一种反演地层横波速度径向剖面的方法 |
CN106814397A (zh) * | 2016-12-21 | 2017-06-09 | 长江大学 | 一种多参数联合反演计算岩石散射衰减的方法 |
Non-Patent Citations (2)
Title |
---|
苏远大 等: "正交偶极声波测井资料在评价地层各向异性中的应用", 《石油物探》 * |
黄杰 等: "二维随机介质中地震波的散射衰减分析", 《地球物理学进展》 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112444865B (zh) * | 2019-08-28 | 2024-05-17 | 中国石油化工股份有限公司 | 勘探缝洞储层的方法、***、装置及存储介质 |
CN112444865A (zh) * | 2019-08-28 | 2021-03-05 | 中国石油化工股份有限公司 | 勘探缝洞储层的方法、***、装置及存储介质 |
CN112444871A (zh) * | 2019-08-30 | 2021-03-05 | 中国石油化工股份有限公司 | 基于散射波地震响应特征的裂缝间距定量分析方法和设备 |
CN111236933A (zh) * | 2020-01-15 | 2020-06-05 | 中国科学院声学研究所 | 利用声波测井仪器对地层径向速度进行分步反演的方法 |
CN111236933B (zh) * | 2020-01-15 | 2021-11-09 | 中国科学院声学研究所 | 利用声波测井仪器对地层径向速度进行分步反演的方法 |
CN113406697A (zh) * | 2020-03-17 | 2021-09-17 | 中国石油天然气股份有限公司 | 储层的横波品质因子获取方法和装置 |
CN113406697B (zh) * | 2020-03-17 | 2023-09-26 | 中国石油天然气股份有限公司 | 储层的横波品质因子获取方法和装置 |
CN111624661A (zh) * | 2020-04-30 | 2020-09-04 | 中国石油大学(北京) | 基于微震散射波品质因子的压裂效果评价方法及装置 |
CN111624661B (zh) * | 2020-04-30 | 2021-04-20 | 中国石油大学(北京) | 基于微震散射波品质因子的压裂效果评价方法及装置 |
CN112034516B (zh) * | 2020-08-17 | 2021-11-12 | 中国石油大学(华东) | 一种评价井周横波三维各向异性的方法 |
CN112034516A (zh) * | 2020-08-17 | 2020-12-04 | 中国石油大学(华东) | 一种评价井周横波三维各向异性的方法 |
CN112859163A (zh) * | 2021-01-07 | 2021-05-28 | 中国石油大学(北京) | 利用散射波确定压裂区地层品质因子变化的方法和装置 |
CN112986247B (zh) * | 2021-04-29 | 2021-08-03 | 中国科学院地质与地球物理研究所 | 基于矿物含量和分布的岩石脆性分析方法、***及设备 |
CN112986247A (zh) * | 2021-04-29 | 2021-06-18 | 中国科学院地质与地球物理研究所 | 基于矿物含量和分布的岩石脆性分析方法、***及设备 |
US11282183B1 (en) | 2021-04-29 | 2022-03-22 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Rock brittleness analysis method and system based on mineral content and distribution and device |
CN113187457A (zh) * | 2021-05-27 | 2021-07-30 | 中海油田服务股份有限公司 | 一种确定储层压裂效果的方法、装置和存储介质 |
CN113866829A (zh) * | 2021-09-28 | 2021-12-31 | 中国石油大学(北京) | 一种三维散射声波远探测扫描成像测井方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN109490965B (zh) | 2020-09-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109490965A (zh) | 一种定量评价地层非均匀性的方法及装置 | |
CN103258091B (zh) | 非常规油气藏水平井段三维岩体力学模型建立的方法及装置 | |
CN104502997B (zh) | 一种利用裂缝密度曲线预测裂缝密度体的方法 | |
CN105182424B (zh) | 一种基于斑块饱和模型定量预测储层孔隙度的方法和装置 | |
CN106855636A (zh) | 基于碳酸盐岩储层露头的原型地质模型地震正演方法 | |
CN105277982B (zh) | 一种泥页岩总有机碳含量地震预测方法 | |
CN104345345A (zh) | 一种页岩储层总有机碳toc含量预测方法 | |
Cox et al. | Surface wave benchmarking exercise: methodologies, results, and uncertainties | |
CN102636812B (zh) | 一种获得碳酸盐岩储层储集空间体积的方法 | |
CN104677734A (zh) | 一种多尺度信息融合的致密地层脆性指数测量方法 | |
CN102213769A (zh) | 一种利用三维垂直地震剖面资料确定各向异性参数的方法 | |
CN106772673B (zh) | 一种页岩储层地应力预测建模*** | |
CN104407381B (zh) | 从井中弹性波速的径向变化获取地层岩石脆裂性质的方法 | |
CN109884710B (zh) | 针对激发井深设计的微测井层析成像方法 | |
CN103728659A (zh) | 一种提高地下岩溶探测精度的方法 | |
CN104730574B (zh) | 构建近地表结构模型的方法 | |
CN107728204A (zh) | 基于叠前纵波各向异性的裂缝预测方法及*** | |
CN106908837A (zh) | 一种压裂裂缝形态和裂缝高度确定方法 | |
CN105588883A (zh) | 三维岩石力学参数获取方法和*** | |
CN112727423A (zh) | 精细评价海上砂岩储层水力压裂裂缝高度的方法 | |
CN104166159A (zh) | 四维微地震监测的裂缝形态处理方法和*** | |
Cai et al. | Intelligent calculation method of relative sonic attenuation and its application to fracture evaluation in tight sandstone reservoir | |
CN104484573B (zh) | 一种确定地层刚性系数的方法 | |
RU2278401C1 (ru) | Способ микросейсмического мониторинга пространственного распределения источников эмиссии и рассеянного излучения и устройство для его осуществления | |
CN106908838A (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 |