静校正新方法
折射波剩余静校正方法

*北京市海淀区学院路中国地质大学地下信息探测技术与仪器教育部重点实验室,100083本文于2005年2月21日收到,修改稿于同年10月2日收到。
・处理方法・折射波剩余静校正方法段云卿*(中国地质大学地下信息探测技术与仪器教育部重点实验室)摘 要段云卿.折射波剩余静校正方法.石油地球物理勘探,2006,41(1):32~35山地、沙漠及其他复杂地表地区地震资料的线性散射噪声和随机噪声很强,有效反射信号弱,资料信噪比较低,静校正问题严重,使用常规剩余静校正方法难以见效。
本文利用折射波信噪比高的特点,将反射波剩余静校正方法应用于折射波资料处理,通过交互手段,逐段估算折射波的速度,用合适的速度对地震记录进行线性动校正,在共炮点或共中心点道集上,用相关方法计算各道与模型道时差,再用统计方法计算出炮点和检波点剩余静校正量。
将该方法应用于信噪比较低、反射波剩余静校正方法难以奏效的复杂地表区,获得良好处理效果。
关键词 剩余静校正 折射波法 共炮点道集 共中心点道集 复杂地表区 模型道1 引言静校正是地震资料处理中至关重要的一环。
我国西部地区地表条件极为复杂,静校正问题尤为严重。
如在沙漠、戈壁、黄土塬或山地等复杂地表区,地形起伏大,表层岩性变化非常剧烈,低降速带厚度变化大,激发和接收条件复杂,近地表条件纵横向千差万别。
近地表地形和低降速带的影响导致地震反射资料不能准确成像,也造成地下构造发生扭曲。
因此,研究复杂地表区静校正方法,对于提高地震勘探精度、降低勘探风险及节约勘探成本有着重要的意义。
本文基于反射波剩余静校正思路,提出一种实现折射波剩余静校正的方法,从而较好地解决了信噪比较低、反射波剩余静校正方法难以奏效地区的静校正问题。
2 方法实现思路静校正的常规步骤为:首先对地震资料进行野外静校正;随后进行折射波静校正;在动效正之后,再进行反射波剩余静校正。
通过这些处理,可初步解决长、中、短波长静校正问题。
但在山地、沙漠及其他复杂地表区,由于线性散射噪声和随机噪声强,有效反射信号弱,地震资料的信噪比往往较低,因此采用常规剩余静校正方法不能建立准确的模型道而达不到预期处理效果。
层析静校正应用第一部分

砂20
砂40
砂37 英东101
砂33
砂34
砂38
砂39
15
层析静校正应用第一部分
西部 中部 东部
英东三维
单尺度层析
微测井
速度(m/s)
深 度 (m)
商业化层析
时间(ms)
深 度 (m)
单尺度层析
层析静校正应用第一部分
16
单尺度层析
商业化软件层析
单尺度层析
层析静校正应用第一部分
17
单尺度层析
尺度1(1*1)
尺度2(2*2)
尺度3(3*3) 尺度4(5*5)
尺度5(7*7)
尺度6(11*11) 尺度7(13*13) 尺度8(19*19) 尺度9(23*23) 尺度10(25*75)
层析静校正应用第一部分
23
层析正反演初始模型
dx=30m dz=10m
24
层析静校正应用第一部分
多尺度第1次迭代分解尺度1
2600m/s
长波长 区域
表层调查模型内插方法
层析静校正应用第一部绿分山层析反演方法
12
商业化层析软件
单尺度层析
单尺度层析
层析静校正应用第一部分
13
单尺度层析
商业化层析软件
单尺度层析
层析静校正应用第一部分
14
单尺度层析
柴达木盆地英东三维是2011年重点勘探项目,该区地表起伏剧烈,相对高 差100m以上,低降速带巨厚,表层结构复杂,静校正问题比较突出。我们采用 自主研发的层析静校正方法较好解决了该区的静校正问题。
层析静校正应用第一部分
42
多尺度第1次迭代累计合成尺度18
层析静校正应用第一部分
第三节动校正与静校正

2、折射波静校正
基准面静校正需要风化层速度和厚度的信息。
但是,野外测量工作有时不能准确地提供这些信 息,由于风化层的速度低于下伏地层的速度,因 此地震记录上能够记录到来自风化层底界的折射 波。
一般情况下,折射波先于地下反射波到达地表, 我们能够比较容易地从地震记录中识别折射波, 进而拾取到折射波的初至时间。
①在O点激发,S点接收,记录下 来的反射波并不是来自O点正下 方,也不是来自S点正下方。在 界面水平的情况下,该反射来自 OS/2处M点的正下方。
②在M点自激自 收,R点的反射 时间为:
t0M
2h V
,
(t0
M
表示是M点的t0时)
若按上述条件在O点激发,S点接收,当然仍可接收到来自R点
的反射,但反射时间:tORS比t0M大,tORS
不同,动校正量不同,即动校正量随时间而 变,这就是动校正中所谓“动”的含义。
三、动校正的实现
1、计算动校正量 tij
2、从 t0 j tij 对应的存储单元搬到与 t0 j 对应的
存储单元中。这样就实现了某道对应时刻 t0 j 的动
校正。 显然,实现动校正,要进行两个循环:
先 t0 j 循环;后 xi 循环。
2
改写成各样点的动校正量的计算公式为:
tij
xi 2 2t0 Vj t0 j 2
(i 1,2, , M ; j 1,2, , N )
式中:M-道集内总道数;N-每道的总样点 数。Vt0j-t0j时刻的叠加速度。
显然 tij 既是 t0 j 的函数,又是 xi 的函数。
对任一道来说,深、浅层反射波 t0 j 和 xi
正量
j
h v0
因为井深校正总是向时间增大的方向校正(所有的校正都是减去 校正量),故此式前面取负号。
探析复杂地表区地震资料静校正方法

探析复杂地表区地震资料静校正方法
吕秋玲
【期刊名称】《石化技术》
【年(卷),期】2017(024)010
【摘要】随着我国石油勘探的发展,勘探区域向地表复杂化方向发展,对地震资料的校正变得越来越重要,静校正方法成为了地震勘探的核心技术.针对地震资料较正问题,结合我国地表复杂性的现状,对现有的高程静校正、模型法静校正、层析静校正等方法进行阐述和介绍,并对折射波静校正、反射波静校正等剩余静校正方法进行总结,最后得出最适合我国现状的静校正方法.研究表明:选择一种合适的静校正方法有利于我国地震勘探的发展,但没有任何一种静校正方法可以适用于我国所有地区,每个区域应根据实际情况选择最合适的静校正方法.
【总页数】1页(P215)
【作者】吕秋玲
【作者单位】中石化石油物探技术研究院江苏南京 211104
【正文语种】中文
【相关文献】
1.复杂地表区三维地震资料成图方法研究与应用效果分析 [J], 张永华
2.复杂地表区地震资料静校正方法探讨 [J], 赵万书;赵兴
3.塔里木盆地沙漠区低信噪比地震资料静校正和去噪方法 [J], 孔剑冰;庄道川;高雁;郑宇清
4.无微测井情况下的可控震源三维地震资料静校正解决方法——以951D-Doshan
区三维为例 [J], 蒋琳;刘宜文
5.江苏油田复杂探区地震资料处理静校正方法研究 [J], 盛道显;谢纯华;孙颖
因版权原因,仅展示原文概要,查看原文内容请购买。
两种静校正方法在准东煤田三维地震勘探的应用

两种静校正方法在准东煤田三维地震勘探的应用耿春明(山东省煤田地质局物探测量队,山东泰安271000)摘要准噶尔盆地东缘沙漠、戈壁、山地并存,地表条件极其复杂,静校正工作是地震勘探资料处理的关键.本文以准东地区石钱滩、奥塔三维地震勘探实例为依据,介绍高程静校正和初至折射静校正两种静校正方法在准东煤田三维地震勘探的应用。
关键词准东盆地静校正三维地震勘探应用中图分类号P631.4文献标识码B1高程静校正高程静校正不考虑低速带,仅对物理点与固定基准面高程差进行校正。
校正方法是直接给出基准面和地表速度,以高程差求出校正量。
这种方法只适用于没有低速带或者低速带结构横向变化不大的地区。
高程静校正的优点是计算效率高,较好的解决长波长静矫正的问题;其缺点是适应条件有限,解决短波长静校正问题的能力有限。
2初至折射静校正初至折射静校正是基于层状地表模型的反演方法。
在近地表结构存在明显的折射层,单炮记录初至中有较明显的折射波时,通过对初至的拾取,建立表层速度模型。
折射波法能够较好的解决短波长和长波长静校正量。
其缺点是初至拾取要耗费大量的时间和人力。
3应用实例3.1高程静校正在石钱滩地区的应用勘查区位于克拉麦里山南麓的戈壁荒漠平原区,*收稿日期:2011-08-10作者简介:耿春明(1981-),男,山东泰安人。
2003年本科毕业于山东理工大学计算机科学与技术专业获学士学位,2008年取得山东科技大学地质工程硕士学位。
物探工程师。
现在山东省煤田地质局物探测量队从事地震资料的采集处理工作。
地势总趋势北高南低,南部和北部稍高,地貌形态为残丘状剥蚀平原与戈壁,三维区海拔507.21 524.41m,比高17.2m,三维区内地形平坦。
勘查区内含煤地层为侏罗系,主要含煤岩组为中统西山窑组,中上统石树沟群下亚群和下统八道湾组中也含有不可采的薄煤。
勘查区内无煤层露头(如图1)。
图1石钱滩三维区地形示意图本区绝大部分无潜水,只在测区的东南角发现有潜水,潜水面约在地表以下2 3m处,区内第四系地层广泛分布,多以砂、砾为主,混合有粘土、亚粘土和砂土,地表为风蚀残余戈壁砾石,部分地段松散的第四系地层对地震波的高频成分有严重的吸收衰减作用,单炮资料初至不好。
初至与静校正

初至拾取与静校正1、初至拾取初至拾取一直是地震勘探中的一个基础而又非常重要的问题。
初至拾取是地表低速带静校正的基础,特别是再地表速度教委复杂的地区,这种静校正更为重要也更为复杂。
准确拾取地震波走时还是地震层析成像方法中计算速度的重要依据。
迄今为止,已经提出了许多初至拾取的方法,如相关法,能力比法,最大振幅法及神经网络法等。
这些方法的一个共同特点是只利用了初至的波形和振幅特征。
拾取是在共炮点道集上进行的,拾取的初至范围要依反演的地质区域内降速层的最大速度和深度而定。
而近偏移距道初至时间的准确拾取是层析静校正方法控制表层速度的关键。
2、准确拾取初至初至时间是折射波静校正的基础数据,直接影响到延迟时间、折射层速度以及低降速带速度或厚度的计算精度。
为了提高初至拾取精度,可以采用以下措施;1对单炮记录进行增益、滤波处理。
改善初至显示质量;2显示地表高程、炮检距曲线。
为初至拾取提供参考。
同时,还可以监视观测系统的正确性;3对单一类型的震源,应保证全区初至拾取位置的一致性;4对单炮记录首先应用野外静校正量,然后进行初至拾取。
同时显示线性动校正,以保证初至拾取的质量,提高初至拾取的效率。
3、静校正的方法近地表静校正方法有很多种,常用的主要有高程静校正、模型法、沙丘曲线法、折射波近地表校正法、层析反演法以及初至波剩余静校正等6种方法。
(1)高程静校正:它不考虑低速带, 仅对物理点与基准面的高程差进行校正。
严格地讲, 高程静校正只适合不存在低速带或者低速带结构横向没有变化的地区。
如果某一地区有低速带, 但低速带横向上的变化对静校正量的影响仅仅是高频的, 或者低速带的剥离和填充对静校正量的影响仅仅是高频的, 可以在精细剩余静校正的基础上使用高程静校正。
高程静校正的优点是计算效率高,较好地解决了长波长静校正问题。
其缺点在于适应条件有限, 解决短波长静校正问题能力比较差。
(2)模型法近地表校正:通过小折射、微测井等常规低速带调查方法得到离散点, 内插出空间速度结构, 进而计算出静校正量的一种方法。
层析反演约束近地表模型静校正方法
层析反演约束近地表模型静校正方法[摘要]以沙漠区野外近地表调查资料为基础,结合低、降速层的分层情况及常用的几种野外静校正方法的原理,依据各种野外静校正方法的优势和实际资料处理中所存在的问题,提出了一种新的针对沙漠区资料野外静校正的方法——利用层析反演约束近地表模型方法。
该方法充分利用了资料的初至信息与野外近地表调查数据资料,获得更为准确的野外静校正量,更好地解决了沙漠区野外静校正问题,在实际资料处理中效果明显[关键词]沙漠区资料;资料处理;野外静校正;近地表模型;层析反演约束中图分类号:td353.5 文献标识码:a 文章编号:1009-914x (2013)14-0181-01静校正技术是沙漠区地震资料处理的一项关键技术,静校正问题的解决程度,直接影响资料处理的最终效果。
沙漠区地表条件复杂,沙丘起伏剧烈,沙山沙墙起伏不断。
地震资料处理中针对沙漠区地震资料的野外静校正方法有近地表模型静校正、折射波静校正、层析反演静校正以及最新提出的近地表数据约束条件下的层析反演静校正等。
不同的沙漠区,其地震资料特点不一样,上述方法针对不同地区地震资料都有相对理想的处理效果。
根据实际。
资料处理情况,分析提出了一种新的解决沙漠区静校正的方法——层析反演约束近地表模型静校正。
1 方法介绍实现该方法首先要消除低降速带的影响,然后用替换速度充填速度校正到基准面上。
由于沙丘起伏剧烈,沙丘高差可达一百多米,野外施工困难,对表层结构调查数据带来不利影响,加上表层结构调查数据控制点少,在沙丘起伏大的地区不能真实反映近地表情况,简单的利用表层调查数据所建立近地表模型不能很好地解决沙漠区的野外静校正问题。
层析反演静校正是利用地震初至波射线的走时和路径反演介质速度结构,根据正演初至时间与实际初至时间的误差修正速度模型,经反复迭代,最终达到误差精度。
该方法求取静校正量时采用射线法计算炮点到检波点的旅行时,从而得到基准面校正量。
其优点是不受地表及近地表结构纵横向变化的影响,但在实际资料应用中发现在信噪比极低的沙漠区出现串层现象,不能得到真实的地下构造。
折射静校正步骤
划分的时候,跟2D不一样的是,这要划分为几个方形。其他与 2D一样。
5、模型的速度分析和延迟时的计算。这一步骤包括两个方面: 一是计算出折射层的速度。此软件有三种算法建立速度模型,它们是:交互速度 分析,简单速度分析和CMP速度分析。交互速度分析适应于折射层视倾角(道间时差) 较大并且校正量较大的2D和3D数据。简单速度分析方法计算速度较快,但不应用在 折射模型有奇异形状的地区,对于海洋资料,这是唯一合适的方法。CMP速度分析用 于各种情况,CMP之间的斜率就显示了在那个位置的速度,网格较大的3D资料计算 出的速度比较稳定。当然也可以给定一个常速度或者是输入一个速度文件。你做的时 候选择第一个:交互速度分析。 二是求取折射层的延迟时。这个软件给出了三种计算延迟时的方法:高斯-赛德 尔(Gauss-Seidel)交互算法、扩展广义互换法(EGRM)、混合算法。高斯-赛德 尔(Gauss-Seidel)交互算法适用于2D和3D数据,这种方法几乎应用于任何情况, 特别是当折射层速度场稳定、可靠时。扩展广义互换法(EGRM)适用于拾取的折射 层靠近炮点的情况,能较好的解决折射层的长波长分量。混合算法是集扩展广义互换 法的长波长特性和高斯-赛德尔(Gauss-Seidel)交互算法的短波长特性于一身的方 法。扩展广义互换法(EGRM)和混合算法只适用于2D数据。你做的时候选择第一个, 也只能选择第一个。
3又由于延迟时是已经计算出的所以风化层底界面平滑以后相应的风化层的速度也不是开始的初始值而是与风化层底界面相关的一系列变化的数值
绿山软件使用方法
1、预备工作 2、建库GO2 3、拾取 4、分区BRANCH 5、模型的速度分析和延迟时的计算 FATHANAL
6、模型建立并计算静校正量 FATHMODL
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
2.
3.
总之,是永远算法要求资料信噪比较高并且反射界面最好水平。
延迟时差法方法原理
转换波检波点延迟时等于转换波与纵波初至时差和纵波检 波点延迟时的和。 波点延迟时的和。
从激发点A到接收点D的纵波 与转换波初至时间差为:
∆t =[TAB +TBE +TED]−[TAB +TBC +TCD] = TCE +TED −TCD
∆ t = TED − TEF + TEF + TCE − TCD = (TED − TEF ) − (TCD − TCF )
d p = TCD − TCF , d s = TED − TEF
p
d
s
= ∆t + d
spΒιβλιοθήκη 根据延迟时定义可知:d 为纵波延迟时,类似地把 d 定义为转换波 延迟时; 为转换波与纵波的初至时差。 因此可以得出如下结论:转换波检波点延迟时等于转换波与纵波初 至时差和纵波检波点延迟时的和。
6.什么是静校正?为什么要作静校正?这种校正为什么叫“静” 校正? 7.静校正和动校正相比,有什么相同和不同之处? 8.静校正的基本假设条件是什么?如不满足会出现什么现象? *试举1~2例说明实际生产中可能发生的不满足基本假设条件的 情况。 9.什么是剩余静校正?为什么要进行剩余静校正? 10.长波长剩余静校正量和短波长剩余静校正量的分类依据是 什么?它们分别会给地震资料的处理与解释带来什么样的影响? 11.如何求取CDP道集中各道的相对剩余静校正量?它使用的是 什么记录?为什么互相关法只能求取相对剩余静校正量?
上节要点
1、折射静校正的目的是什么? 2、折射静校正和野外一次静校正的关系是什么? 3、延迟时(时间深度)的定义? 4、如何求取延迟时? 5、目前静校正方法的发展方向是什么?
前述方法的局限性
最大能量法-----模拟退火法 即使解决了低信噪比大静校正量的计算问题,但是无法 避免的一个问题是,由于常规静校正是基于地表一致性 假设条件的,随着勘探程度的深入,实际工作情况越来 越不满足该条件,从而出现“静校而不静”的问题。 波动方程静校正可以解决这个问题。
式中 s是 i 位置炮点静校正量,g j是 j 位置接收点静 i c 校正量,k 是共中心点位置 k 剩余正常时差,它随 y 炮检距 h的平方而变化, k是由于构造引起的静校 正误差量。 炮点静校正量 s 采用P波的常规方法求取。 i
根据地表一致性假设条件,在构造平缓的情况下,检波点产生的静校正 量远大于动校正的剩余时差和地层倾斜产生的静校正误差之和。
随着地震勘探的不断深入,勘探方式也由单纯的纵波勘 探向多波多分量勘探发展。
在多波勘探中,由于转换波勘探成本较低,在实际生产 中已经得到了推广。但是仍然有很多问题急需解决,转 换波静校正就是其中之一。
5)转换波静校正 )
转换波资料特点:
信噪比较同一地区纵波资料要低 转换波静校正量较同一位置纵波静校正量要大2---10倍 转换波初至不易识别和拾取
静校正技术现在已不是一项单纯的资料处 理技术,它和地震波传播理论、野外采集技 术等都有密切关系,目前应该说还没有根 本解决这一世界性难题,有些问题有待于 你们走上工作岗位后去发现和解决。
第三章思考题
1.什么是动校正?为什么要作动校正?这种校正为什么 叫“动”校正? 2.动校正量Δt和界面的埋藏深度及炮检距有什么关系? 速度取值不准会给动校正带来什么影响?当速度取值偏 小或偏大时,校正后的反射波时距曲线是什么形式的? 3.为什么说成组“搬家”法动校正是一滤波器?它是一 什么滤波器?为什么? 4.动校正引起的波形拉伸畸变与哪些量有关?能否根本 避免?为什么?生产实际中是如何处理这种畸变的? 5*.如何利用动校正量板表中的参数来判断拉伸畸变的 严重程度?
方法存在的问题
1.
方法推导的假设条件为构造平缓,仅当反射界面水平时,取得 好的效果。而地下界面倾斜角度较大或者存在一定构造时,难 以计算得到好的结果。 方法使用标准道进行互相关计算静校正量,如果初始标准道不 合适或者由于静校正量较大经过叠加后标准道变得不能代表实 际情况时,计算效果较差。 在信噪比较低的资料中,应用互相关算法计算静校正量,效果 较差。
12.统计分离法求取炮点和检波点的绝对剩余静校正量的基本依据是什 么?它需要使用哪些形式的地震记录? 13.折射静校正主要用于消除哪类静校正量?它的主要依据是什么? 14.时间深度(延迟时)的定义? 15.目前采用的非地表一致性静校正方法是什么方法?应用条件是什么? 16.转换波静校正量为什么比纵波静校正量更难求取? 另外应注意: 另外应注意: 1)关于动校正量表的制作及应用也要求掌握,对课堂上所讲例题一定要 理解,要自己会制作量板表、会根据量板表完成“搬家”; 2)静校正的实现,特别是当静校正量是负值时如何“搬家”,应搞清 楚。
宣汉地区纵波典型记录
宣汉地区转换波典型记录
宣汉地区纵波典型记录
广元地区纵波典型记录
广元地区转换波波典型记录
目前转换波静校正的研究现状
1)共接收点叠加道集相干法(利用反射波进行静校正) P.W.Cary和D.W.S.Eaton(1993)年提出,来解决短波长的较 大的静校正量
2)延迟时时差法(利用折射初至进行静校正) 李彦鹏在 2002 年 SEG 年会上提出
v( x, y, z,ω) = M (x, y, z,ω) / ρ(x, y, z )
可以解决的问题
非地表一致性静校正
补偿地震波在穿过低降速带时,其能量和频率的的衰减,提高地震 资料分辨率和信噪比。
波动方程静校正存在的问题
需要精确的近地表速度和厚度信息,才可以完成正确的 延拓。 在三维地区特别是复杂的三维地区,搞清楚近地表的速 度和厚度信息存在很大困难。 波动方程静校正还处于一个发展阶段。
4)波动方程静校正 波动方程静校正
思路
通过波动方程延拓的方法,把地表接收到的地震波场延拓到高 速层顶。(相当于把检波器埋置在高速层顶,从而完成静校正工作)
∆Z
V,Q
基于粘滞声波介质的延拓方法可归纳为两类: 单道相移延拓方法 解频率–空间域粘滞声波方程延拓方法
基于水平层状介质假设下的单道延拓
U (t + ∆t, ω) = U (t, ω)e−i ( Kω −iLω ) ∆t = U (t, ω)eiKω∆t e Lω∆t
g
j
>>
c
k
h
2 ij
+
y
k
因此经P波炮点静校正后有 :
t
ij
≈
g
i
那么每一道的静校正量可由下式计算:
g j = max[ Pj (t ) ⊗ G j (t )]
式中为 g j 第个 j 接收点的静校正量, j (t)为第 j 个 G 共接收点(CRP)的叠加道。⊗代表相关,j (t)为参考 P 道。
P j (t ) =
j k = j − n +1
∑
G k (t + g k )
其中 n 为参考道的叠加道数
正演模型(黑色线为地表,红色线为地下反射界面)
Vs =1000m/s Vp= 2000m/s 倾角 = 11.5o
进行炮点静校正量得到的共检波点转换波时间剖面
对静校正量进行构造改正后得到的共检波点叠加剖面
K (ω ) = ω v (ω ) = K ω − iL ω
复体积模量
M ( x, y , z , ω ) = M R ( x, y , z , ω ) + iM I (x, y , z , ω )
Q −1 = M I (x , y , z , ω ) M R (x , y , z , ω )
介质的吸收特性-品质因子 复速度
共检波点叠加道相关法基本原理 由于在共检波点道集中各道检波点静校正分量相同。那么对共检 波点道集应用炮点静校正,进行动校正然后叠加,形成的叠加道 就体现了对应检波点的静校正量。 根据地表一致性条件,任意一道转换波总的静校正 量都可分解为4部分:
t ij = s i + g j + c k h ij2 + y k
频率-空间域标准线性体粘滞声波方程
ω2 ∂ 2 P(ω , x, y , z ) ∂ 2 P (ω , x, y, z ) ∂ 2 P (ω , x, y , z ) + + + P (ω , x, y , z ) = 0 ∂x 2 ∂y 2 ∂z 2 M ( x, y , z , ω ) ρ
把速度表示成复数,此时可以定义复波数K(w)