薄互层型油气储层同步挤压变换域分析方法
2018-08-01高静怀刘乃豪张茁生姜秀娣陈树民
高静怀,刘乃豪,吕 奇,张茁生,姜秀娣,陈树民
(1.西安交通大学电子与信息工程学院,西安交通大学数学与地球物理探测研究中心,海洋石油勘探国家工程实验室,陕西西安710049;2.南京电子技术研究所,中国电子科技集团公司智能感知技术重点实验室,江苏南京210029;3.西安交通大学数学与统计学院信息科学系,陕西西安710049;4.中国海洋石油总公司研究总院,海洋石油勘探国家工程实验室,北京100029;5.大庆油田有限责任公司勘探开发研究院,黑龙江大庆163712)
薄互层储层是陆相沉积盆地中一种重要的储层类型,在我国鄂尔多斯盆地、松辽盆地等盆地中,薄互层储层发育,因此,这类储层具有重要研究意义。薄互层由多个薄层叠合而成。所谓薄层是相对于地震波波长来说的,如果地层的厚度远小于地震波波长,就称其为薄层。松辽盆地中、浅层地层的单层厚度一般仅为几米,而地震波的波长为百米量级。鄂尔多斯盆地单层厚度仅为几米到十几米,地震波波长也为百米量级。因此,在此类盆地的地震勘探中,地震波无法分辨单层的厚度及反射界面的位置。
来自薄互层的地震信号是一种包含频率和幅度快速变化的非平稳信号,时频分析作为研究非平稳信号的一种有力工具[1],已经被广泛应用于地震信号处理和解释领域。例如,PARTYKA等[2]将短时Fourier变换(STFT)用于墨西哥湾某区块实际地震资料的河道厚度刻画以及不连续性检测。CHAKRABORTY等[3]提出了基于小波变换(CWT)的地震信号谱分析方法。高静怀等[4-5]研究了利用小波变换对地震信号进行属性分析的方法,并构造出适合于地震信号分析及瞬时属性分析的三参数小波[6]。ODEBEATU等利用S变换[7]检测与含气饱和度相关的异常现象[8]。高静怀等[9]提出一种广义S变换,并将其用于薄互层的地震响应分析。LI等[10]利用Wigner-Ville分布(WVD)刻画了碳酸盐储层。LIU等[11]利用局域化属性检测河道和低频异常。时频分析还被广泛用于地震瞬时属性的提取,TANER等[12]首先提出了复地震道分析方法,并给出了基于复地震道提取的地震瞬时属性的定义、物理意义及其在地震解释中的应用。CASTAGNA等[13]和MARFURT[14]将瞬时谱分别用于碳氢检测和地质结构刻画。HAN等[15]提出了结合经验模态分解(EMD)和瞬时频率的地震信号分析方法,该方法能够更精细地刻画地质结构。
STFT由于窗函数在时间域及频率域宽度固定,所以一旦选定窗函数,它在时间和频率两个方向上的分辨率就不变。CWT克服了STFT的缺点,在分析信号的缓变分量时利用长时窗以获得较高的频率分辨率;分析快变分量时采用短时窗以获得较高的时间分辨率。CHAKRABORTY等[3]论证了CWT相比于STFT在地震信号谱分解上的优势。S变换用频率控制窗的宽度,不仅具有小波变换的优点,而且还能够直接得到时间-频率谱[7]。WVD虽然可以得到更集中的时频表示,但是在分析多分量信号时,会产生交叉项的干扰,带来分析误差。为了提高时频分辨率,AUGER等在KODERA等[16]的研究基础上,提出了利用STFT和WVD等时频分布的相位信息进行能量重排的时频分析方法[17]。重排时频分析方法已经得到了广泛的应用,PENG等[18]将基于CWT的重排方法用于旋转机械缺陷的特征提取,WU等[19]将基于WVD的重排方法用于地震资料的谱分解。然而这种时频重排方法没有给出重构公式,限制了其应用。
综上所述,各种时频分析方法虽然在地震信号分析中得到了成功的应用,然而它们都有一定的局限性[16]。为了得到具有更高时频域分辨率的分析方法,HUANG等[20]提出了EMD的分析方法(称为Huang变换),与上述时频分析方法不同,Huang变换是将待分析信号分解为若干个本征模态函数之和,它在故障检测、地震信号分析等应用中具有显著效果[21],但该变换缺乏严格的理论基础。
DAUBECHIES等[22]在研究EMD的基础上,系统地提出了同步挤压小波变换,并给出了一套完整的数学理论,可实现待分析信号的本征模态函数分解。同步挤压变换已经被应用于语音信号、心电图、古气候学、降噪、故障检测及地震数据解释等领域[22-26]。LI等[26]提出了广义同步挤压变换,HERRERA等[27]研究了一般的同步挤压变换在地震信号解释中的应用,WANG等[28]将同步挤压变换应用于储层特征刻画。THAKUR等[29]将同步挤压变换推广到STFT以及用更一般的波形函数来替代余弦函数,HUANG等[30]提出了同步挤压S变换,并用于地震信号谱分解。
本文聚焦薄互层型储层,该类储层的地震响应信号往往具有频率和幅度变化都较快的特点,常用的母小波函数(也称基本小波函数),如Morlet小波等,分析这种信号时分辨率较低,相比之下高静怀等人提出的三参数小波[6]变换对薄层具有更强的刻画能力。然而,与所有的小波变换类似,三参数小波变换随着尺度的减小,小波函数频域展宽,导致频率分辨率降低。为了进一步提高三参数小波变换的分辨能力,本文提出了同步挤压三参数小波变换,进而提出了基于同步挤压三参数小波变换的薄互层组及其内部结构的谱分解分析方法。将该方法用于合成数据、薄互层模型合成数据及实际地震资料,并与常用方法的结果进行了对比以说明方法的有效性。
1 同步挤压三参数小波变换
1.1 连续小波变换
式中,小波变换系数Ws(a,b)≠0。
(1)
1.3.1 计算瞬时频率
德国波龙是致力于尖端测量与测试技术的高科技企业，拥有50年的行业经验。波龙的测量技术“德国制造”协助各行业提高生产率和产品品质。主要提供激光对刀仪，接触式对刀仪，触发式工件测头，模拟量工件测头，粗糙度测头，孔径测量仪，3D测量软件等。
其中：xij=0，表示Xi和Yj不匹配；xij=1，表示Xi和Yj匹配和分别为双边主体的匹配竞争度。稀缺性资源具有更强的竞争力，在双边匹配中，数量较少的一方占据优势地位。当M>N时，Y方主体占优势，且ωX<ωy；当m>
1.2 小波函数的选取
由公式(12)计算得到瞬时频率,意味着能量应该往这个频率上挤压(集中),将小波变换系数累加到该频率成分上,即(a,b)→(ω,b),从而进行能量重排,那么如何将对应在同一个频率成分上的小波系数累加到一起呢？
Morlet小波的时域表达式为:
(2)
此次大直径水力旋流器的更换，为后续粉精煤系统的改造提供了实验数据和参考依据。大直径水力旋流器结构简单，现场管理方便，避免了小直径水力旋流器组因单个旋流器堵塞而造成溢流跑粗的现象，降低了生产事故率。
(3)
(4)
1.2.2 三参数小波
肠体成品表面红色有光泽，肠体表面干燥饱满，成型性好，切片坚实，组织紧密细致有弹性，黑米分布均匀，口感细腻，有肉香味和米香味，嚼劲适中，咸味适中，特有滋味，鲜香可口，无异味。
三参数小波是针对地震信号的特点而提出的,可以很好地匹配待分析信号,Morlet小波、改进的Morlet小波和匹配地震子波小波(BMSW)均为其特例,它在地震信号分析中表现出很多优势,并被广泛应用[5-6]。
三参数小波的时域表达式为:
(5)
式中,Λ=(σ,τ,β)为参数集;σ为解析小波的调制频率;τ可以控制衰减快慢,即控制小波的频率域宽度;β为能量延迟因子;其它参数如下。
三参数小波的频域表达式为:
(9)
1.2.3Morlet小波与三参数小波的对比
AMSCs的分离、培养 采集患者和正常人脂肪组织各50 ml，提取AMSCs，培养于AMSCs专用培养基，传代培养至第3(P3)代备用。
信息化教学资源本身具有很强的丰富性和多样性，它能把原本枯燥乏味的教学内容生动形象地呈现给学生，在一定程度上促进学生感性思维的发展，有效挖掘学生的学习潜能，在无形中促使学生把所学内容和实际生活紧密结合起来，从而真正起到优化小学课堂教学结构的作用。
1) Morlet小波要求调制频率较大,一般取σ>5.33,这就导致待分析信号中低频分量的时间分辨率较低,为了提高低频分量的时间分辨率,需降低调制频率,但是此时Morlet小波在时域包络会由单峰变成双峰,对待分析信号做小波变换会使信号在多个位置局域化,小波变换得到的结果产生一些假象。
严格灵活按照《熔炼铅顶吹炉标准化操作规程》进行粗铅及炉渣的排放。做好顶吹炉与炉前的协调工作，确保所有的铅和铅冰铜从铅口排出。若渣含铅过高，可打开金属口再次排放金属，或根据实际情况延迟炉渣的排放继续烟化，确保烟化后炉渣含铅量达到排放标准。
2) Morlet小波只有调制频率σ一个参数可调,而且要求σ比较大时才是严格解析的小波函数,自由度很低,缺乏灵活性。在小波变换的实际应用中,为了刻画奇异性,需要度量待分析信号的局部正则性,此时小波的消失矩非常重要。我们就需要根据实际需要来选取合适的小波函数,不仅希望能调节调制频率,还希望能够控制小波的消失矩等特性,因此希望小波函数能有更大的自由度。
三参数小波由3个参数来控制,有着很高的自由度,不仅可以控制调制频率,还可以控制小波包络的衰减速度,此外,在调制频率很低的时候,依然保持时域的单峰特性,具有良好的时频局域化特性。
通过公式(17),在某个固定的时刻b,且小波变换系数Ws(a,b)≠0时,计算其瞬时频率ωs(a,b),将所有瞬时频率都为某一频率ω的小波变换系数通过公式(17)累加,就完成了能量的重分配,得到了挤压后的时频分布。
1.3 同步挤压三参数小波变换
在小波变换中,由于基本小波函数在时频域支集(即Heisenberg盒)为一矩形,因此待分析信号在小波变换域能量分布在某个区域内,这不利于准确刻画快速时变信号的时频特性。为减小基本小波函数的影响,通常在小波变换域进行能量重排,以减少能量扩散带来的不利影响。常用的能量重排方法都需要解决两个基本问题:能量往哪个频率成分集中？能量如何往某个频率上集中？DAUBECHIES等[22]在小波变换的基础上,通过同步挤压操作,来解决上述两个基本问题。同步挤压操作分为两个步骤:第一步,计算其瞬时频率ωs(a,b),即重排准则;第二步,通过ωs(a,b)来对小波变换系数Ws(a,b)做频率重排。
Morlet小波的Heisenberg盒的面积最小,具有最佳的时频联合域分辨率,在实际中得到了广泛的应用,但也存在着一定的局限性,主要为:
(10)
(11)
从公式(11)可以看到,假如小波ψ(t)的峰值频率为ξM,则小波变换的结果将在尺度a=ξM/ω0处取到最大值,并以这个能量最大的尺度为中心形成一个尺度带,造成能量扩散,为了得到更集中的时频分布,需要进行同步挤压操作,具体如下:
瞬时频率有很多种定义方式,DAUBECHIES等[22]提出的信号s(t)的瞬时频率定义为:
(12)
对于任意待分析信号s(t)∈L2(R),t为时间,R为实数集合,L2表示平方可积函数空间,则其小波变换定义为:
根据公式(12),对于一个单频余弦信号的小波变换会形成一个尺度带,瞬时频率可表示为:
(13)
从公式(13)可以看出,公式(12)定义的瞬时频率就是余弦信号的频率。至此,我们可归纳为:对于一个余弦信号,它的小波变换得到的时间-尺度域结果会在某个能量最大的尺度邻域形成一个尺度带,但是这些尺度对应的小波系数通过公式(12)计算出来的瞬时频率都为余弦信号的频率ω0,因此我们可以设想将这些尺度的能量都集中到ω0上,这就解决了前文提到的能量往哪个频率成分集中的问题。
还可以看到,公式(12)是定义在Ws(a,b)≠0的位置,如果某个信号有若干个频率成分,那么得到的小波变换系数会形成若干个尺度带,公式(12)在每个固定的时刻,对相应的尺度都能算出一个瞬时频率,于是可以在一个固定的时刻算出若干个瞬时频率。但这和通过Hilbert变换等方法定义的瞬时频率不同,Hilbert变换等方法是在每个固定的时刻只能算出一个瞬时频率,因此,这种方法类似于EMD方法,它适用于分析多分量信号。
1.3.2 时间-尺度域到时间-频率域的映射
1.2.1Morlet小波
在2.3单因素试验结果上，对四氢呋喃用量（X1）、KOH 甲醇溶液质量浓度（X2）、提取温度（X3）3 个因素进行响应曲面试验设计。设响应曲面因素与水平及编码值见表2，响应曲面法优化稻谷中叶黄素提取方法见表3。
对于给定的信号s(t)∈L2(R)的小波变换Ws(a,b),若选取的小波ψ(t)是解析小波,则有表达式:
在计算机网络技术专业课程体系中，广泛引入企业课程并认证，将行业企业标准引入相关课程标准。典型行业与企业认证证书的技术要求核心要素融入专业课程内容当中，学生学习完相关课程后即可考取相关证书。
(14)
(15)
因为sa(b)是s(b)的解析信号,故有:
(16)
由公式(16)可以看出,小波变换系数再乘以因子a-1得到的结果和原信号的解析信号只差一个常数因子,因此可以容易地得到其反变换,即公式(16)所示。同样,如果将Ws(a,b)a-1都分配给公式(12)中相应的瞬时频率成分上,则就得到挤压后的时频分布,且存在简单严格的反变换。由此,我们得出时间-尺度域到时间-频率域的映射如下:
(17)
采用统计学软件SPSS15.0进行数据分析处理,其中,计数资料采用X2检验,以百分比表示,计量资料采用t检验,以均数±标准差表示,P<0.05表示差异突出,具有统计学意义。
1.3.3 阈值的选取
(18)
ση=median·
(19)
式中,ση为与噪声水平相关的量,可利用小波变换的前nv个尺度来估计;median代表取中值;0.6745是针对高斯噪声的正则化因子。由于有N个时刻,所以实际应用中,可以对N个时刻求得的ση取均值,即可得到最优的阈值。
2 模型算例分析
2.1 合成信号
将本文方法用于合成信号,以验证其有效性。图1a 中合成信号[27,32]的表达式为:
(20)
该信号的前半段为一余弦信号,其频率成分不随时间变化;后半段的频率随时间而变化,为典型的时变信号,这为准确刻画其频率分布和频率变化带来了极大困难。此外,在0.20s处有一个调制频率为100Hz的Morlet小波原子,在1.00s,1.53s和1.56s处分别有一个50Hz的Ricker子波,三参数小波的参数选为Λ=(3,0.5,0)[6]。基于Morlet小波变换、Morlet小波同步挤压变换(后文简称为同步挤压变换)、本文提出的同步挤压三参数小波变换(下文称为本文方法)计算得到的时频结果分别如图1b至图1d所示。与连续小波变换方法相比,看到图1c和图1d两个同步挤压变换的结果,有着更高的时频分辨率。此外,从红色矩形处可以看到,在频率突变和频率变化较快的地方,相比于同步挤压变换结果,本文方法具有更高的时频分辨率,时频局域性更好,可以准确地刻画时变信号的频率分布和频率变化情况。从红色箭头所指处可以看出,与同步挤压变换相比,本文方法法可以更加准确地刻画单个Ricker子波;即使当两个子波到达时间相近时,本文方法法仍能清晰地分辨两个Ricker子波。
群众路线是中国特色社会主义政治发展路径的现实指向。在社会主义政治话语体系中“无产阶级政党—群众”的能动性关系不同于西方政治话语体系中的“政党（精英）—大众”的互动关系。从理论维度来看，在无产阶级政党的视野中“群众”既有主体性，也是党领导的对象。从实践维度来讲，群众路线是党的工作方法，体现为民主决策的过程和政策执行过程的动员能力，更是权力运行的监督机制。本文从中国共产党历届领导人对群众路线的阐释出发，集中探讨中国特色社会主义政治发展视域下群众路线的实践图景。
图1 合成信号及其时频结果a 合成信号; b 基于Morlet小波变换的时频结果; c 基于同步挤压变换的时频结果; d 基于本文方法的时频结果
2.2 薄互层模型
我们将不同时频分析方法用于经典的薄互层模型,以验证本文方法的有效性和准确性。薄互层模型的反射系数序列和合成地震记录如图2所示,图2a中反射系数的幅度大小都为0.5,正负相间,时间间隔从左往右依次为1,2,3,…,9,10,9,…,3,2,1ms,薄层起始位置为200ms,时间采样间隔为1ms。与主频为50Hz的Ricker子波卷积得到合成地震记录,如图2b所示。将不同时频分析方法用于该合成地震记录,得到的时频结果如图3所示。当选取Morlet小波为基本小波函数进行小波变换时,由于Morlet小波的峰值频率较高,对于低频来说,对应的尺度较大,所以其时间分辨率较低,主能量带的中心频率(给定时刻主能量带中能量最大值对应的频率)对地层厚度的变化不灵敏,导致对薄互层的时-频响应特性(调谐频率随着层的变薄而升高)刻画不够清晰准确。
图2 反射系数序列(a)和合成地震记录(b)
图3 图2b合成地震记录的时频结果a 基于Morlet小波变换; b 基于三参数小波变换; c 基于同步挤压变换; d 基于本文方法
图3c 中同步挤压变换结果也可以看到,尤其对于低频来说,时间分辨率很低,而此时的同步挤压效果也不好,不足以刻画薄层的时频特征。图3b和图3d中,三参数小波的参数选取为Λ=(3,1,0),此时图3b 三参数小波变换结果的锥形条纹较为清晰,但仍然不能清晰地刻画薄互层的局部结构特征。图3d为本文方法的时频结果,不仅可以清楚地反映主能量带,这与薄互层的时-频响应特性是一致的,而且可以清晰地揭示薄互层的局部结构,即锥形条纹。我们提取了图3b和图3d中的单频(160Hz)结果,如图4所示,其中蓝线是原始反射系数序列,可见本文方法提取的单频分量的位置与反射系数序列的位置吻合较好,从而可以准确地确定薄互层模型中反射系数的位置。但当反射系数间隔过小(小于3ms)时,本文方法提取的单频分量的结果也不能准确地刻画薄层的位置。
图4 自图3时频结果中提取的单频(160Hz)分量a 基于三参数小波变换; b 基于本文方法
3 实际资料应用
将本文方法应用于渤海地区某实际地震数据。首先,将本文方法用于某二维地震叠后剖面,共400道,时间采样点数500,时间采样间隔2ms,如图5所示,图中绿色椭圆指示了河道所在位置。将不同的时频分析方法应用于该地震数据,并提取30Hz的单频结果,如图6所示。图6a中Morlet小波变换的单频切片过于粗糙,不能有效地识别和定位绿色椭圆处的河道,而图6b和图6c中同步挤压变换和本文方法均可以较为准确地定位河道。但对比图6b和图6c 可见,本文方法对于河道的位置和河道展布的刻画,比同步挤压变换刻画得更加清晰和准确,这是因为本文方法具有更高的时频分辨率,可以得到更稀疏的时空表示。
图5 二维实际地震剖面(绿色椭圆指示了河道)
将不同的时频分析方法应用于渤海地区某三维实际地震数据,对比其对河道等地质构造的刻画能力。图7是该三维地震数据的沿层切片,该层河流相薄砂体较为发育。原始地震数据分辨率较低,不利于河道的期次划分,河道之间的叠置关系也不够清楚。我们选取不同的时频分析方法,提取单频30Hz沿层切片,如图8所示。图8b为基于三参数小波变换提取的单频切片结果,可见其对河道的刻画比图8a中Morlet小波变换结果更为清晰。图8c为同步挤压Morlet小波变换提取的单频切片结果,虽然分辨率略高于三参数小波变换的结果,但仍不能清晰地刻画河道叠置现象。图8d 为本文方法提取的单频切片结果,对河道的刻画更为清晰,也可以较为准确地刻画叠置的河道,如红色箭头和红色矩形指示位置。
图6 单频30Hz剖面a 基于Morlet小波变换; b 基于同步挤压变换; c 基于本文方法
图7 某工区三维原始地震数据体沿层切片
图8 单频30Hz沿层切片结果(红色箭头和红色矩形指示了河道分布情况)a 基于Morlet小波变换; b 基于三参数小波变换; c 基于同步挤压变换; d 基于本文方法
4 讨论与结论
本文提出了同步挤压三参数小波变换,并用于薄互层型油气储层刻画,该变换可以得到更稀疏的时频表示,有利于精确刻画薄互层等地质结构。
在我妈的审美观里，只有拥有像林黛玉那样生着樱桃小嘴的女孩，以后才嫁得出去。于是她抱着我去找医生，很认真地问能不能把我的嘴稍微缝上一点儿，结果被医生轰出了办公室。
合成地震记录算例表明,本文方法具有更高的时频分辨率,可以得到更稀疏的时频表示。薄互层模型算例表明,本文方法不仅可以准确地刻画薄互层的时-频响应特性,而且可以清楚地反映薄互层的局部结构;本文方法提取的单频结果可以准确地定位薄互层模型中的反射界面。实际地震资料算例表明,本文方法具有更稀疏的时空表示,可以更加准确地刻画河道的位置及展布情况,对于河道局部细节的刻画更为清晰。
“他们不是我的朋友，”威尔放下盘子，“同一派别的人不一定都能和睦相处。况且，爱德华和迈拉在约会，我可不想当电灯泡。”
同步挤压三参数小波变换可对薄互层组的整体特性(图3d中的主能量带的变化趋势,即先降低后升高)进行精细刻画,该变化趋势与地层的韵律及沉积相密切相关,需要进一步研究。同步挤压三参数小波变换也可对薄互层内部结构进行精细表征(图3d 中的垂向条纹),这些与薄互层整体特征结合,可表征薄互层的横向变换,为储层预测提供帮助。
下一步的工作将根据我国松辽盆地和鄂尔多斯盆地的实际情况,对薄互层类型进行分类,建立数据库,系统地研究其时频域的特征,为储层预测提供依据。
