一种确定储层压裂效果的方法、装置和存储介质制造方法及图纸

本申请实施例公开了一种确定储层压裂效果的方法、装置和存储介质，包括：获取目标深度范围内储层压裂前和压裂后的单极阵列波形数据；分别根据储层压裂前和压裂后的单极阵列波形数据计算得到对应的径向相对速度剖面；根据压裂前和压裂后的径向相对速度剖面的差异确定储层压裂效果。本申请实施例通过对比压裂前和压裂后的径向相对速度剖面的差异准确的确定储层压裂改造效果。

【技术实现步骤摘要】
一种确定储层压裂效果的方法、装置和存储介质
本申请实施例涉及但不限于声波测井
，尤其涉及一种确定储层压裂效果的方法、装置和存储介质。
技术介绍
目前，对非常规储层的勘探程度不断提高，非常规储层具有致密、非均质性强等特征，因此，压裂成为了有效开发非常规储层的关键环节。压裂环节中对压裂效果的评价至关重要，最有效的压裂技术为水力压裂，评价压裂效果主要的测井方法包括井温测量方法、同位素测井方法和注硼中子寿命检测方法，但目前这些评价地层压裂效果的方法总体上可靠性和实用性较差。声波各向异性在套管井中的应用已证实是一种评价压裂效果和确定压裂裂缝特征的有效方法，该方法直观准确，但是该方法有其适用性，也存在很多局限性，会受压裂缝网形态、地层相对于井眼的角度等影响，且手段单一，无法精细刻画压裂对井壁附近地层的改造作用，已经不能满足非常规储层压裂效果评价的需求。亟待需要实现一种新的储层压裂效果评价方法，能够不受压裂缝网形态、斜井等影响的确定储层压裂改造效果的方法，以满足非常规储层压力效果的评价。
技术实现思路
以下是对本文详细描述的主题的概述。本概述并非是为了限制权利要求的保护范围。本公开提供了一种确定储层压裂效果的方法、装置和存储介质，应用压裂前、后的单极波形数据获取压裂前、后井壁附近地层的径向速度剖面，通过对比压裂前、后径向相对速度剖面的差异确定储层压裂改造效果。一方面，一种确定储层压裂效果的方法，所述方法包括：获取目标深度范围内储层压裂前和压裂后的单极阵列波形数据；分别根据储层压裂前和压裂后的单极阵列波形数据计算得到对应的径向相对速度剖面；根据压裂前和压裂后的径向相对速度剖面的差异确定储层压裂效果。一种示例性的实施例中，所述分别根据储层压裂前和压裂后的单极阵列波形数据计算得到对应的径向相对速度剖面，包括：对储层压裂前和压裂后的单极阵列波形数据分别执行以下操作：S1.采用慢度-时间相关函数提取地层纵波时差；S2.根据单极阵列波形数据确定标定参数；S3.根据所确定的标定参数、所提取的地层纵波时差计算径向相对速度剖面。一种示例性的实施例中，所述慢度-时间相关函数为:上述公式中，Xm(t)是N个接收换能器阵列中的第m个接收换能器的单极阵列波形数据，d是接收换能器之间的间隔距离，T'是时间窗Tw的起始时间点，s是时差区间中的时差值。一种示例性的实施例中，所述根据单极阵列波形数据确定标定参数，包括：对所获取的单极阵列波形数据，应用模拟指数函数计算单极阵列波形的走时Tn，其中，n表示接收器编号，n＝1,2,3，…；对所获取的单极阵列波形数据，应用积分法计算第1个接收器上声波的走时曲线并定义为参考R1走时曲线；当判断参考R1走时曲线与单极阵列波形T1走时曲线满足预设条件时，确定该仪器直径和流体时差作为标定参数。一种示例性的实施例中，所述对所获取的单极阵列波形数据，应用模拟指数函数计算单极阵列波形的走时Tn，包括：步骤5.1根据所获取的单极阵列波形数据计算第一波峰或波谷所对应的时间和振幅；步骤5.2根据所获取的单极阵列波形数据计算第一波峰或波谷右边所对应的第一过零点的时间；步骤5.3应用模拟指数函数计算包括第一波峰或波谷的半周期波形，其中，模拟指数函数为：f(A，t，α，T)＝Atαsin(2πt/T)上述公式中，A表示第一波峰或波谷所对应的振幅，t表示时间，α表示指数变量，T表示第一波峰或波谷对应的波形周期；步骤5.4对上述模拟指数函数求导计算，得到α与t0和T的函数关系：α＝-2πt0/(Ttan(2πt0/T))其中，t0表示第一波峰或波谷所对应的时间；步骤5.5将波形的坐标系按照坐标变换公式进行坐标变换，其中，坐标变换公式如下：f(t)＝f(t+(tzx-T/2))；其中，tzx表示第一波峰或波谷右边第一过零点的时间；步骤5.6利用最小二乘法确定符合预设条件的半周期波形，并根据确定的半周期波形确定起跳时间，其中，该半周期波形是应用模拟指数函数所模拟的波形；步骤5.7重复执行步骤5.1至5.6，计算出每个接收器上单极波形的走时Tn，其中n表示接收器的编号，n＝1,2,3，…。一种示例性的实施例中，所述积分法的计算公式为：上述公式中，v(z)是提取的声速曲线，s是声源的深度位置，r1是第一接收器的深度位置，TTf是声波在井眼流体中的传播时间，TTref是参考R1走时。一种示例性的实施例中，所述根据所确定的标定参数、所提取的地层纵波时差计算径向相对速度剖面，包括：利用走时层析技术建立井壁周围地层沿井眼轴向和径向的二维速度剖面，得到径向纵波速度剖面；根据该径向纵波速度剖面和所提取的地层纵波时差，计算径向相对速度剖面。一种示例性的实施例中，所述利用走时层析技术建立井壁周围地层沿井眼轴向和径向的二维速度剖面，得到径向纵波速度剖面包括：步骤8.1根据所确定的标定参数，应用射线追踪法和初始迭代模型，计算声波走时tn：上述公式中，Sn(vn)为依赖于地层速度模型vn(r，z)的声波传播的最短路径，vn(r，z)为根据所提取的纵波时差建立的地层速度模型，其中，该初始迭代模型是根据二维速度剖面建立的模型，该模型包括m*s个网格点；步骤8.2利用最小二乘法计算出满足实测走时||Tn-tn||min的速度模型Vn(r，z)，其中，Vn(r,z)∈vn(r,z)；步骤8.3将步骤8.2计算得到的速度模型作为初始迭代模型，重复执行步骤8.1和步骤8.2，当模型中每个网格点全部计算完成后得到径向纵波速度剖面。一种示例性的实施例中，所述根据该径向纵波速度剖面和所提取的地层纵波时差，计算径向相对速度剖面包括：根据该径向纵波速度剖面和所提取的地层纵波时差，采用径向相对速度剖面的计算公式得到径向相对速度剖面；其中，该径向相对速度剖面的计算公式为：Δvp'＝(vpi-vpw)/vpi×100％上述公式中，Δvp'为径向深度上纵波速度vpw与地层纵波速度vpi的相对差别的百分比，即径向相对速度变化量。另一方面，本公开还提供了一种确定储层压裂效果的装置，包括：存储器和处理器；所述存储器用于保存用于确定储层压裂效果的程序，所述处理器用于读取执行所述用于确定储层压裂效果的程序，执行上述实施例中任一项所述的确定储层压裂效果的方法。另一方面，本公开还提供了一种存储介质，所述存储介质中存储有确定储层压裂效果的程序，所述程序被设置为在运行时执行上述实施例中任一项所述的确定储层压裂效果的方法。本申请实施例公开一种确定储层压裂效果的方法、装置和存储介质，包括：获取目标深度范围内储层压裂前和压裂后的单极阵列波形数据；分别根据储层压裂前和压裂后的单极阵列波形数据计算得到对应的径向相对速度剖面。通过本公开的方案，可以根据压本文档来自技高网...

【技术保护点】
1.一种确定储层压裂效果的方法，其特征在于，所述方法包括：/n获取目标深度范围内储层压裂前和压裂后的单极阵列波形数据；/n分别根据储层压裂前和压裂后的单极阵列波形数据计算得到对应的径向相对速度剖面；/n根据压裂前和压裂后的径向相对速度剖面的差异确定储层压裂效果。/n

【技术特征摘要】
1.一种确定储层压裂效果的方法，其特征在于，所述方法包括：
获取目标深度范围内储层压裂前和压裂后的单极阵列波形数据；
分别根据储层压裂前和压裂后的单极阵列波形数据计算得到对应的径向相对速度剖面；
根据压裂前和压裂后的径向相对速度剖面的差异确定储层压裂效果。


2.根据权利要求1所述的确定储层压裂效果的方法，其特征在于，所述分别根据储层压裂前和压裂后的单极阵列波形数据计算得到对应的径向相对速度剖面，包括：
对储层压裂前和压裂后的单极阵列波形数据分别执行以下操作：
S1.采用慢度-时间相关函数提取地层纵波时差；
S2.根据单极阵列波形数据确定标定参数；
S3.根据所确定的标定参数、所提取的地层纵波时差计算径向相对速度剖面。


3.根据权利要求2所述的确定储层压裂效果的方法，其特征在于，所述慢度-时间相关函数为:



上述函数中，Xm(t)是N个接收换能器阵列中的第m个接收换能器的单极阵列波形数据，d是接收换能器之间的间隔距离，T'是时间窗Tw的起始时间点，s是时差区间中的时差值。


4.根据权利要求3所述的确定储层压裂效果的方法，其特征在于，所述根据单极阵列波形数据确定标定参数，包括：
对所获取的单极阵列波形数据，应用模拟指数函数计算单极阵列波形的走时Tn，其中，n表示接收器编号，n＝1,2,3，…；
对所获取的单极阵列波形数据，应用积分法计算第1个接收器上声波的走时曲线并定义为参考R1走时曲线；
当判断参考R1走时曲线与单极阵列波形T1走时曲线满足预设条件时，确定测量仪器直径和流体时差作为标定参数。


5.根据权利要求4所述的确定储层压裂效果的方法，其特征在于，所述对所获取的单极阵列波形数据，应用模拟指数函数计算单极阵列波形的走时Tn，包括：
步骤5.1根据所获取的单极阵列波形数据计算第一波峰或波谷所对应的时间和振幅；
步骤5.2根据所获取的单极阵列波形数据计算第一波峰或波谷右边所对应的第一过零点的时间；
步骤5.3应用模拟指数函数计算包括第一波峰或波谷的半周期波形，其中，模拟指数函数为：
f(A，t，α，T)＝Atαsin(2πt/T)
上述公式中，A表示第一波峰或波谷所对应的振幅，t表示时间，α表示指数变量，T表示第一波峰或波谷对应的波形周期；
步骤5.4对上述模拟指数函数求导计算，得到α与t0和T的函数关系：
α＝-2πt0/(Ttan(2πt0/T))
其中，t0表示第一波峰或波谷所对应的时间；
步骤5.5将波形的坐标系按照坐标变换公式进行坐标变换，其中，坐标变换公式如下：
f(t)＝f(t+(tzx-T/2))；
其中，tzx表示第一波峰或波谷中右边第一过零点的时间；
步骤5.6利用最小二乘法确定符合预设条件的半周期波形，并根据确定的半...

【专利技术属性】
技术研发人员：王文文，钱玉萍，杨玉卿，徐大年，孙大君，张宏伟，祁晓，侯振学，成家杰，林龙生，范川，
申请(专利权)人：中海油田服务股份有限公司，
类型：发明
国别省市：天津;12