基于压缩感知的非规则地震勘探观测系统设计与数据重建-周松.pdf

《基于压缩感知的非规则地震勘探观测系统设计与数据重建-周松.pdf》由会员分享，可在线阅读，更多相关《基于压缩感知的非规则地震勘探观测系统设计与数据重建-周松.pdf（9页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第56卷第5期2017年9月石油物探GEOPHYSICAL PROSPECTING F()R PETROLEUMVoL56，No5Sep，2017周松，吕尧，吕公河，等基于压缩感知的非规则地震勘探观测系统设计与数据重建J石油物探，2017，56(5)：617625ZHOU Song，LV Yao，LV Gonghe，et a1Irregular seismic geometry design and data reconstruction based on compressive sensingI-JGeophysical Prospecting for Petroleum，2017，56(5

2、)：617625基于压缩感知的非规则地震勘探观测系统设计与数据重建周 松1，吕 尧2，吕公河1，舒国旭2，石太昆2，霍守东2(1中石化石油工程地球物理有限公司，北京100020；2中国科学院地质与地球物理研究所，北京100029)摘要：提出了基于压缩感知的非规则观测系统优化设计和地震数据重建方法。基于贪心序贯策略，利用逐个选取采样点位置的方法来降低感知矩阵的最大互相关值，从而完成非规则观测系统的优化；利用L。正则化和L。正则化混合迭代的方法求解欠定采样矩阵，重建地震记录，获得良好的效果。进行了TFT工区实际观测系统的设计，鉴于接收电缆长度的限制，添加最大道间距和最小道间距作为观测系统设计过程中

3、的约束，得到了符合实际勘探要求的非规则观测系统的优化设计。该方法为三维高密度地震勘探的大规模开展提供了有效的方案。关键词：压缩感知；稀疏；非规则；观测系统设计；数据重建中图分类号：P631文章编号：10001441(2017)05061709文献标识码：ADOI：103969jissn10001441201705001Irregular seismic geometry design and data reconstruction based oncompressive sensingZHOU Son91，LV Ya02，LV Gonghel，SHU Guoxu2，SHI Taikun2。HU

4、O Shoudon92(1Sinopec Petroleum Engineering Geophysics Co，Ltd，Beijing 100020，China；2Institute of Geology and Geophysics，Chinese Academy of Sciences，Beijing 100029，China)Abstract：This paper proposes a method of irregular seismic geometry optimal design and data reconstruction based on compressive sens

5、ingTo optimize the irregular geometry design in seismic data acquisition，the method selects sampling points one by oneaccording tO a greedy sequential strategy，which could decrease the maximum cross-correlation value of the sensing matrixTo reconstruct seismic data effectively，the method solves the

6、underdetermined sampling matrix through an Lo and L1 regularizationmixed iterative algorithmWhen applied to the irregular seismic geometry design in the TFT exploration area，the method uses themaximum and the minimum group intervals as constraints tO overcome the limitation of the receiver cable len

7、gthThe results showthat the proposed method provides an effective way for the development of high-density 3D seismic explorationKeywords：compressive sensing，sparse，irregular，seismic geometry design，seismic data reconstruction随着我国油气资源勘探向中西部发展，地球物理的研究也逐渐转向地表构造复杂的区域。在这些区域进行勘探，不论是资料采集，还是后期的资料处理和解释，都对地球物

8、理人员提出了更高的要求。对于地震勘探来说，在复杂地区施工，往往会遇到峭壁、河流、沟壑、村庄和工业区等不易甚至无法布置检波器的情况，从而给数据的采集带来了一定的困难，甚至导致无法采集的情况发生。如果因为这些原因导致原始数据的大量缺失，就会严重影响整个工区的勘探质量和成像效果，从而给资料解释带来困难。收稿日期：2017-0124；改回日期：2017-0615。作者简介：周松(1963一)，男，高级工程师，从事地球物理数据采集、非常规勘探以及物探技术管理等工作。基金项目：中国石油化工股份有限公司科技攻关项目“基于压缩感知的地震勘探采集技术研究”及中石化石油工程技术服务有限公司科技项目“压缩感知技术在

9、地震勘探中的应用研究”(SGl652K)联合资助。This research is financially supported by the project of the Ministry of Science and Technology of SINOPEC and the project of Sinopec OilfieldService Corporation(Grant NoSGl652K)万方数据618 石油物探 第56卷经典信号采集理论一般都会遵循ShannonNyquist采样定理，即：若要使采集到的数据能不失真地保持原信号中的信息，采样频率必须是原信号频带宽度的两倍以上。地

10、球物理的野外资料采集也都遵循该定理。为了描述目标体的精细特征，往往需要采集更精确的信号，因此，需要通过在空间上加密采样点数和在时间上减小采样间隔来实现。这样的采集方式大大增加了采集的数据量，从而导致采集成本的急剧增高。基于信号稀疏性的采样理论1。2的压缩感知(Compressed Sensing，CS)在图像压缩、无线通信、模式识别和医疗成像等研究领域受到了格外的关注。2006年，CANDES等3提出基于压缩感知的信号恢复技术，大大推动了压缩感知在地球物理勘探方面的应用。压缩感知技术的发展和应用基于大部分信号都可以进行稀疏表示。在稀疏域对信号进行重建，可以获得原始的真实信号。压缩感知技术突破了

11、Shannon-Nyquist采样定理的局限，大大节省了采集数据所占的空间，成为提高地震勘探效率的重要方法。目前，压缩感知理论在地球物理勘探领域得到广泛应用4。1 5|。基于压缩感知理论，HERRMANN等41提出了一种多震源的采集框架；MOLDOVEANU5探索了海上数据采集的随机观测方法；MOSHER等6。7提出并完善了一种在观测系统设计中基于约束条件的更优选择激发点和接收点位置的非均匀采样方法；陈生昌等8提出了一种针对地球物理数据的高效采集方法，并通过归一化各列数据的感知矩阵的互相关系数最小作为约束条件来设计最优的采集方案。基于压缩感知的地震采集技术，获得的数据需要进行重建9。1 0|。

12、HENNENFENT等11提出了利用抖动欠采样进行波场重建的方法，并基于曲波变换进行了地震数据的重建；wU等12提出了Dreamlet域地震成像和数据重建的方法；LI等13提出了基于插值压缩感知的数据重建方法(interpolated CompressiveSensing，interpolated CS)，设计了非均匀优化采样观测系统(NonUniform Optimal Sampling，NUOS)，并通过地震数据重建，分别在海上和陆上取得了良好的试验效果；王华忠等14详细讨论了压缩感知技术在地震勘探应用中的问题和难点，从不同方面展示了压缩感知的有效应用。本文通过降低感知矩阵的最大互相关值，

13、从而进行非规则观测系统的优化设计；并利用L。正则化和L，正则化混合迭代的方法进行欠定采样矩阵的求解，完成了地震数据的重建。1 基于压缩感知的非均匀观测系统设计方法压缩感知理论的基本假设是目标信号具有稀疏性或者可压缩性，也就是说，目标信号或者其在某个变换域中只有有限个成分不等于o(对应稀疏性)，或者只有有限个成分远大于0(对应可压缩性)。若目标信号有K个成分不等于0，则称该信号为K稀疏的。假设具有稀疏性或者可压缩性的目标信号为X，稀疏变换采用傅里叶变换F，则正交稀疏基为傅里叶基函数或分窗傅里叶基函数，于是有：xFHs (1)式中：s为信号X在傅里叶域的稀疏表示；()H表示共轭转置。这样采集得到的

14、数据可以看作是采样函数或采样矩阵与目标信号相乘的结果。记采样矩阵为西，当进行满采样时，有一I，其中，J表示单位矩阵。当不能满足满采样条件(如采样数不足)时，m为J抽出的若干列组成的矩阵，只有采样位置对应的列向量被保留，则有采样数据：Y一嘶一FHs一弧 (2)式中：y一鲫H记为感知矩阵。根据压缩感知理论，成功重建K稀疏的目标信号的充分必要条件是感知矩阵甲满足有限等距准则(Restricted Isometry Property，RIP)条件，即对于任意2K稀疏向量l，均有(1一艿)iIv惦|甲v惦(1+3)lv孵成立。然而，求解满足RIP条件的矩阵是个NP-hard问题，取而代之，可采用一种可计

15、算的准则来进行采样设计。设矩阵y的列向量为lf，列向量间的最大互相关值为：lf，。Hlf，川 ，。户一警可百盯_可瓦面 JAMAIL-RAD等16指出，感知矩阵列向量间的最大互相关值即为非规则采样归一化频谱的最大非零频率振幅。在频率域，采样数据Y的频谱是感知矩阵V的频谱和目标信号的稀疏表达s的频谱之间卷积作用的结果。肚就是因为非规则采样造成傅里叶基(稀疏基)的正交性被破坏所引起的最大频谱泄漏，压制口使感知矩阵甲的频谱近似Delta函数，否则，感知矩阵y的频谱出现多个峰，与稀疏表达s卷万方数据第5期 周 松等基于压缩感知的非规则地震勘探观测系统设计与数据重建 619积作用后，将在频率域(稀疏域)

16、产生“假频”噪声。CANDES等17指出，口越小，则信号在非规则采样后能够重建的概率就越高。由于掣一面FH，其中F是固定的，因此，在采样数目不足时，可以通过改变采样矩阵o，即改变非规则采样点的分布来降低最大互相关值，从而优化采样矩阵，提高目标信号重建的概率，即：西一arg min“(垂) (4)JAMAIL-RAD等16给出了口达到下界时的一种确定性采样方式，然而这种采样方式只在一些特定的条件下才能成立。注意到该优化问题是非凸的，具体优化时只要能够找到局部最优解就足够满足需求18-19。这里我们采用贪心序贯策略来优化采样点的位置，进而构建采样矩阵。具体策略如下：1)首先确定目标信号的网格点数目

17、和间距，将所有目标信号的网格点作为候选采样点；2)遍历所有候选采样点，计算将其加入采样矩阵后的肛值，取使1值最小的候选点作为新的采样点，更新采样矩阵；3)重复步骤2)的操作，直到采样点数达到规定数目；4)在初步选定规定数目的采样点位置后，利用时间抖动(Jitter)方法D8-19，以(4)式为目标，进行采样微调，最终确定采样矩阵。以一维简单情况为例，假设在751个网格点上，选取375个点(近50)作为采样点。当规则采样时，采样点分布及其归一化后的傅里叶谱如图1所示。可以看出，欠采样的规则采样函数会产生很强的频谱泄漏，造成假频混叠。同样的网格和采样点数，当随机采样时，采样点分布及其归一化后的傅里

18、叶谱如图2所示。可以看到，随机产生的非规则采样降低了采样矩阵与傅里叶矩阵的相关性，没有假频产生，但是在整个频率范围都出现了很低的频谱泄漏，最大幅值为0794 9lo。保持网格和采样点数目不变，采用贪心序贯策略获得优化的非规则采样，采样点分布及其归一化后的傅里叶谱如图3所示。图3中采样点非均匀分布，压制了假频，降低了整个频率范围内的频谱泄漏，频谱泄漏的最大幅值0538210 ，比随机采样的频谱泄漏小。由压缩感知原理可知，更小的频谱泄漏表明由非规则采集而造成的漫布在整个频带内的杂乱噪声更弱，这样采样频谱就更加逼近规则满采样时的频谱，更有利于目标信号的重构，优化后的非规则采样比随机采样所得到的数据能

19、够完美重建的概率更高。当采样点数较少或者采样率较低时，这种差异更加明显。而当采样点较多时，采样条件往往会受其它因素如地形、实际施工条件限制等影响，此时可在采样点优化设计过程中加入这些约束，从而得到满足施工要求的采集方案，这是随机采样很难实现的。另外，对于某些二次施工工区，采样设计时还需要考虑其它指导目标如地下照明、数值模拟恢复信号的信噪比等，此时可在优化目标(4)式中加入其它约束项，然后优图2随机采样点分布(a)及其频谱(b)万方数据620 第56卷化完成满足要求的非规则采样点设计，这也是随机采样很难实现的。地震信号的稀疏性已经被广泛应用。对于三维地震勘探的数据采集情况来说，假设真实信号为D。

20、-，其时间维度为N，与非规则采集的时间维度相同，而空间维度为NN，N，与非规则采集数据重建后的规则数据所具有的空间维度对应。选择傅里叶变换作为稀疏变换，先将均匀满采样的时间方向做傅里叶变换，真实信号D。l由t-xy域变换至f-xy域。由于采样网格只与空间有关，故对所有频率切片使用相同的非规则采样设计。为了降低公式的维度，不失一般性，取其第q个频率切片，记为S。，再对|sa做空间方向上的傅里叶变换，将其变换至f-k域，记为Sa，则有：S。=F，x雪9 Fx。 (5)式中：F，和F。分别为37方向和Y方向上的傅里叶变换矩阵。引入向量化操作1 6|，将二维目标信号转化为一维信号。向量化算子：vec(

21、A)一口 (6)表示将任意矩阵A各列叠加排列成一个很长的列向量n。向量化算子具有如下性质：vec(ABc)一(cT o A)vec(B) (7)式中：A，B，c是任意3个连乘有意义的矩阵；o表示克罗内克积。通过以上计算，可以将f-k域中第q个频率切片的数据向量化，得到真实信号为D reaI向量化的稀疏表示s 9。在实际非规则采样时，时间上均匀满采样，与真实信号的时间维度保持一致，均为N，。在空间上，使用基于压缩感知的非规则优化采集方案，实际的采样点数为M，MN。也就是说，要应用一个采样算子G，对稀疏表示s 9进行采样，采样算子G实际上就是从N维单位矩阵JN中取M行而得到，其决定了实际非规则采集

22、的采样点(检波点和炮点)位置，且G的选择与这些非规则位置一一对应。于是有：Gs。一G(F?(两F?)s 9 (8)目标就是根据估计s-重建D。-。通过向量化，将二维数据转化为一维数据求解，则此时感知矩阵掣可表示为：掣一G(F?圆F?) (9)于是，可以通过最小化最大互相关值来确定采样算子G，即：Garg min1(G) (10)采用贪心序贯策略解决该问题时，需要增加候选网格的维度。考虑到非规则采样在滚动排列上的复杂性，实际施工的试验区采用了非滚动全检波点接收的方案，降低了炮点和检波点的联系。这里，将炮点和检波点的分布进行独立设计，减少了计算量。同时，假设没有地下构造的先验信息，仅凭目标信号的采

23、样网格精度和采样数目来约束设计方案。在实际观测系统设计优化时，需要根据实际情况添加一定的约束条件。特别是在设计检波点位置时，考虑到设备的实际情况，需要添加约束条件。第一是检波线距非规则，检波点非规则地分布在这些检波线上，而不是所有检波点都完全分散分布；第二是相邻检波器之间的最大距离受实际施工时电缆线长度的限制。考虑到这些因素，在设计检波点位置时，先优化检波线的位置，以(4)式为优化目标；然后固定检波线的位置，将M个采样点位置平均分布在N，N。的网格上，然后利用Jitter采样17的方式，让采样点在初始位置附近抖动取值，以(10)式为目标进行迭代优化。每次迭代只抖动其中一个采样点，而其它点保持不

24、变，取使岸最小的位置为新的采样点，更新该点位置，然后选择下一个点进行抖动，进行下一次迭代，直到抖动任何一个采样点均不能使肚变小为止1 8|。为加速迭代收敛，可以先让各采样点在平均分布位置附近随机抖动，然后以此为初始条件进行迭代优化，万方数据第5期 周 松等基于压缩感知的非规则地震勘探观测系统设计与数据重建 621同样可以达到满足要求的优化设计结果。对炮点进行设计时，考虑最大炮间距和最小炮间距，在炮间距限制范围之内进行优化设计，使得非规则炮点的位置在目标网格上完全打散分布。此时可以根据序贯思想，利用贪心算法，逐个添加炮点位置以使得当前的岸最小。每步迭代先固定之前已有的所有采样点，然后在剩余未被选

25、中的网格点中，选出使岸最小的网格点作为采样点。要加速运算，可以控制变量，每次只在一个方向取得最小值，然后换方向接着搜索，直到收敛。在所有M个采样点全部确定之后，再以(10)式为目标，进行全局微调，最终得到炮点的优化非规则采样设计。2 TFT工区的实际应用我们利用TFT工区的一块区域进行实际工区的非规则观测系统设计，所选工区范围为14km8km，检波点布设到整个工区范围，炮点区域为iokm2km。规则检波点设计为32条接收线，线距240m，每条线456道，道间距30m。非规则设计是在原始规则检波点设计的基础上，保持大致的覆盖范围，但是只取总采样点数的75，即线数32条不变，平均每条线取342道，

26、平均道间距40m，共10 944道。重建的目标网格为42条接收线，线距180m，每条线911道，道间距15m。优化后的检波点分布如图4所示，每个蓝色的十字型代表一个检波点的位置，最小点距为15ITI，最大点距为75m。同时，线距也进行了非规则优化设计，最小线距为150 m，最大线距为705m。这里由于检波线数较少，为保证重建质量，降低频谱泄漏，采样点没有落在实际的整数网格上。检波点位置归一化后的频谱如图5所示，在垂直于检波线方向上由于检波线条数较少，频谱泄漏较强，最大值为o0930，但仍然要比规则采样导致的假频小很多。实际工区的非规则炮点设计为将1800个炮点分散在10km2km的范围内，目标

27、重建网格为炮点100 908O706馨蜂0 5篁04善030201万方数据622 石油物探 第56卷距15m，炮排距90m。优化后的炮点分布如图6所示，每个红色的方框代表一个炮点位置，其归一化后的频谱如图7所示，可以看到，频谱泄漏很小，最大值为00444。最终得到的该工区优化的非规则观测系统如图8所示。妊蜊兰皿虱6优化设计的非规则炮点分布归一化波数訇7优化设计的非规则炮点频谱l 0o 90 80 7o6肇崔05篁04里0 302O l訇8优化设计的非规则观测系统(红色表示炮点位置。蓝色表示检波点位置)万方数据第5期 周松等基于压缩感知的非规则地震勘探观测系统设计与数据重建 6233 地震记录的

28、模拟与重建对于地球物理问题来说，基于压缩感知的数学重建模型可表示为：Y一叭 (11)式中：zRN表示原始地震数据；表示采样矩阵；YRM表示观测地震数据。这里，MN，即采样矩阵是一个欠定矩阵，意味着需要从欠采样的地震数据Y中重建出原始地震数据。但是，如果当原始数据X在某个域c中的系数是稀疏的，即sCx，那么方程(11)可写为：Y一弘 (12)式中：甲表示感知矩阵，掣-cIJCH。对此，可以通过求解反问题：rain 1I y一帆 (13)来获取稀疏系数s，进而通过xcHs求取原始数据00 200 300 400 500道号ax。但是，由于采样矩阵是一个欠定矩阵，可知上述问题是一个典型的NP-har

29、d问题。对此，需要对上述问题增加约束条件，本文利用L。正则化和L。正则化混合迭代的方法求解上述问题。首先通过软阈值算法求解L。范数约束问题：min|y一眺11；+圳s lI， (14)通过对(14)式的求解，可在进行下一步算法前获得一个好的迭代初值，并估算信号的稀疏水平。然后，通过硬阈值算法求解L。范数约束问题：min|Y一帆lI；+A 1I s忆 (15)反复迭代，直至获得最终结果。为了验证该算法重建地震记录的有效性，采用主频为30 Hz的Rieker子波合成500道的单炮地震记录，道间距10m，每道3001个采样点，采样间隔为2ms，如图9a所示。图9b为采样50地震道后的不完整地震记录，

30、图9c为使用本文算法恢复的地震记录，图9d为本文算法恢复的地震记录与原始模型记oo 200 300 400 500道号boo 200 300 400 500道号 道号c d图9地震记录处理结果a原始地震记录；b随机缺失的地震记录；C重建结果；d重建结果与原始记录之差05l 0234、粤蜷导1凰5OO234、擘蜡导1量5O0565O万方数据624 石油物探 第56卷录的误差。从图9c的重建地震记录与原始记录(图9a)的对比以及两者的误差(图9d)可以看出，两者同相轴清晰，连续性相同，振幅一致性准确，重建误差小。因此，该算法具有良好的数据恢复效果。图10是该地震记录的频率一波数谱。从原始地震记录和

31、重建地震记录的频率一波数谱对比可以看出，重建数据的频率一波数谱与原始记录的频率一波数谱非常接近，两者的误差很小(图lOd)，进一步表明了该算法良好的重建效果。以上实验均表明，利用基于压缩感知理论的重建方法，可以较好地恢复缺失的地震数据。)0 8 101506，： 20鞲卫睾芝25 104。 303505 413450 5(1 0归一化波数b图10地震记录处理结果频率一波数谱g原始地震记录频率一波数谱；b随机缺失的地震记录频率一波数谱；c重建结果的频率一波数谱；d重建结果与原始地震记录频率一波数谱之差4 结论与建议基于压缩感知理论，本文进行了地震勘探非规则观测系统的优化设计和地震数据的模拟重建。

32、通过降低感知矩阵的最大互相关值的方法优化设计采样点的位置，并应用于实际工区，该过程通过添加部分约束条件来使得最终的设计在实际可行的情况下达到最优解。利用L。正则化和L。正则化混合迭代的方法进行地震数据的重建，取得了良好的效果。随着勘探精度的提高，野外地震资料的采集都要求能够达到“两宽一高”的标准，这就给野外采集带来了相当大的工作量。基于压缩感知的非规则观测系统的地震数据采集和信号重建，可以达到“相同野外采集工作量的条件下获得更高分辨率的地震资料”或者“获得相同分辨率的地震资料需要较少的野外工作量”的目的。因此，利用压缩感知技术进行地震数据的采集和处理，可以大大降低野外资料采集的工作量，对目前“

33、两宽一高”要求的地震勘探具有重要意义。坚蜷芒1丑O864，0O0万方数据第5期 周松等基于压缩感知的非规则地震勘探观测系统设计与数据重建 625致谢：感谢中石化西北石油分公司在正演模型建立方面给予的大力支持，感谢中石化休斯敦研发中心的罗明秋博士在目标模型的建立和观测系统的设计方面给出的一些好的建议!参考文献1DONOHO D LCompressed sensingJIEEE Transactions on Information Theory，2006，52(4)：128913062 CANDES E J，ROMBERG J，TAO TRobust uncertainty principles

34、：exact signal reconstruction from highly incomplete frequency informationl,JIEEE Transactions on Information Theory，2006，52(2)：4895093 CANDES E J，ROMBERG JQuantitative robust uncertainty principles and optimally sparse decompositionsI-J1Foundations of Computational Mathematics，2006，6(2)：227-254r 4HE

35、RRMANN F J，ERLANGGA Y A，LIN T T YCompressive simultaneous fullwaveform simulationJGeophysics，2009，74(4)：A35一A405 MOLDOANI 7 NRandom sampling：a new strategyfor marine acquisitionJ1Expanded Abstracts of 80“Annual Internat SEG Mtg，2010：5155r 6MOSHER C C，KAPLAN S T，jANISZEWSKI F DNon-uniform optimal sam

36、pling for seismic survey designl,JExpanded Abstracts of 74山EAGE AnnualConference，2012：3333367 MOSHER C C，LI C，MORLEY L C，et a1Non-uniform optimal sampling for simultaneous source surveydesignJExpanded Abstracts of 84幽Annual InternatSEG Mtg，2014：1051098陈生昌，陈国新，王汉闯稀疏性约束的地球物理数据高效采集方法初步研究VJ石油物探，2015，54(

37、1)：2435CHEN S C。CHEN G X，WANG H CThe preliminarystudy on high efficient acquisition of geophysical datawith sparsity constraintsJGeophysical Prospectingfor Petroleum2015，54(1)：24359冯飞，王征，刘成明，等基于Shearlet变换稀疏约束地震数据重建J石油物探，2016，55(5)：682691FENG F，WANG Z，LIU C M，et a1Seismic data reconstruction based on

38、 sparse constraint in the Shearlet domainJ1Geophysical Prospecting for Petroleum，2016，55(5)：6826911,101崔永福，郭念民，吴国忱，等不规则观测系统数据规则化及在相干噪声压制中的应用l-J石油物探，2016，55(4)：524532CUI Y F，GUO N M，WU G C，et a1Regularization ofirregular geometry seismic data and its application inthe coherent noise suppressionJGeoph

39、ysical Prospecting for Petroleum，2016，55(4)：524532r11 HENNENFENT G，HERRMANN F JCurvelet reconstruction with sparsity-promoting inversionsuccessesand challengesJ1Expanded Abstracts of 69“EAGEAnnual Conference 2007：231 233121 wu R s，wu B Y，GENG YImaging in compresseddomain using dreamletsl,J1Beijing 2

40、009 InternationalGeophysical Conference and Exposition，2009：212r13LI C B，MOSHER C C，KAPLAN S TInterpolatedcompressive sensing for seismic data reconstructionJExpanded Abstracts of 82”Annual Internat SEG Mtg，2012：1-6141王华忠，冯波，王雄文，等压缩感知及其在地震勘探中的应用J石油物探，2016，55(4)：467474WANG H Z，FENG B。WANG X W，et a1Co

41、mpressedsensing and its application in seismic explorationl,JGeophysical Prospecting for Petroleum，2016，55(4)：467474E15陈祖庆，王静波基于压缩感知的稀疏脉冲反射系数反演方法研究J石油物探，2015，54(4)：459466CHEN Z Q，WANG J B-A spectral inversion method ofsparse-spike reflection coefficients based on compressed sensingJGeophysical Prosp

42、ecting for Petroleum，2015，54(4)：459466r16 JAMAIL-RAD H，KUVSHINOV B，TANG Z，et a1Deterministically subsampled acquisition geometries foroptimal reconstructionJExpanded Abstracts of 78小EAGE Annual Conference，2016：360365r171 CANDES E J，wAKIN M B An introduction to compressive sampling：a sensingsampling

43、paradigm thatgoes against the common knowledge in data acquisitionlJ1IEEE Signal Processing Magazine，2008，25(2)：2130r181 HENNENFENT G，HERRMANN F JSimply denoise：wavefield reconstruction via jittered undersamplinglJGeophysics，2008，73(3)：V19一V28r19LI C，KAPLAN S T，MOSHER C C，et a1Compressivesensing：us2014063443P20150507(编辑：朱文杰)万方数据