井间及井地联合地震波CT成像方法及装置制造方法及图纸

本发明专利技术公开了一种井间及井地联合地震波CT成像方法及装置。该方法包括以下步骤：步骤1：建立井间及井地联合观测系统，其包括设置于井间的井间检波器和炮点，以及设置于地面的地面检波器；步骤2：通过井间检波器和地面检波器拾取实际初至波旅行时；步骤3：根据测井资料建立初始速度模型，根据初始速度模型计算理论旅行时和射线路径；步骤4：以实际初至波旅行时与理论旅行时的残差作为目标函数，基于SIRT算法迭代反演速度修正值，对初始速度模型进行修正。该方法的优点在于能够兼具井间及井地CT成像的优点，以提高勘探分辨率，实现空间三维探测。

【技术实现步骤摘要】
井间及井地联合地震波CT成像方法及装置
本专利技术涉及地震勘探及开发领域，更具体地，涉及一种井间及井地联合地震波CT成像方法。
技术介绍
地震层析成像是一种基于模型的数据拟合方法。与医学中的计算机辅助层析成像(ComputerTomography，CT)类似并受其启发，地震层析成像通过分析台站记录到的由天然地震或者人工地震激发的地震波信号来获取地球内部结构的信息。1976年，Aki和Lee从60个台站处提取了32个地震的走时信息，并利用这些走时信息成功地获得了美国加里福尼亚地区的地壳结构(AkiandLee,1976)。他们开创性的工作开启了利用地震层析成像方法来研究地球内部结构的新时代。根据研究范围的大小，地震层析成像可以分为局部区域层析成像、地区层析成像和全球层析成像；根据所用地震波的类型，又可分为体波层析成像和面波层析成像。面波层析成像利用的数据是具有长周期的面波，因而分辨率较低，适用于全球或者地区的成像研究；而体波层析成像利用的数据是具有更短波长的体波，分辨率相对较高，适用于局部区域、地区以及全球等多尺度结构的研究。在实际工程建设中隐伏的不良地质问题，如岩溶、陷落柱、裂隙、裂缝、断裂破碎带、软弱夹层、地下空洞和不明埋设物等，会造成工程质量隐患或施工困难、浪费，故必须进行精细的工程勘察。由于工程检测要求探测规模很小的异常体，并要求较高的分辨率，因此需要高精度的勘察手段。而在油气田资源勘探过程中，为了获取精确的速度模型，往往开展大量的钻井及测井工作。钻孔勘探的精度较高，但成本高、耗工期，不可能在工区内大量采用，层析成像方法可以弥补这一勘察手段的缺陷。目前，地震CT从观测方式上来分析可以分为井地CT、井间CT、垂直地震剖面CT，其中井地CT中，由于大多采用井中激发地面接收，射线方向主要是垂直的，故垂向分辨率不如横向分辨率，存在垂向模糊性。而相反，井间CT地震记录上的初至波主要包括透射波、饶射波，由于波传播路径段，具有较高的能量及较宽的频带。波的传播方向基本上是水平的，横向分辨率不如垂向分辨率，存在横向模糊性。因此，有必要开发一种地震波CT成像技术，能够兼具井间及井地CT的优点，以提高勘探分辨率，实现空间三维探测。本专利技术
技术介绍
部分的信息仅仅旨在加深对本专利技术的一般
技术介绍
的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
技术实现思路
本专利技术提出了一种井间及井地联合地震波CT成像方法及装置，其能够通过野外单炮记录初至波拾取获得地震波旅行时间，并依据测井等资料建立初始速度模型，求解正演射线路径和理论旅行时，以理论和实际旅行时残差为目标函数建立反演方程，逐次反演迭代计算，最终获得高精度的介质速度场，解决了速度异常体勘探以及为地震资料处理提供精确的速度模型。根据本专利技术的一方面，提出了一种井间及井地联合地震波CT成像方法，包括以下步骤：步骤1：建立井间及井地联合观测系统，所述井间及井地联合观测系统包括设置于井间的井间检波器和炮点，以及设置于地面的地面检波器；步骤2：通过所述井间检波器和地面检波器拾取实际初至波旅行时；步骤3：根据测井资料建立初始速度模型，根据所述初始速度模型计算理论旅行时和射线路径；步骤4：以所述实际初至波旅行时与所述理论旅行时的残差作为目标函数，基于SIRT算法迭代反演速度修正值，对所述初始速度模型进行修正。优选地，通过滑动时窗能量比法拾取所述实际初至波旅行时。优选地，基于SIRT算法迭代反演速度修正值，对所述初始速度模型进行修正包括：步骤401：以Sj表示初始模型中的第j个网格的慢度，第k个炮点对应Nk个接收点，其中，j＝1,2,…,M，M为总的网格数量，k＝1,2,…,NUM，NUM为总的炮点数量；步骤402：根据射线追踪方法计算第k个炮点对应的接收点的理论旅行时：其中，表示在第i次迭代中第k个炮点对应的接收点的理论旅行时，i表示迭代次数，i的初始取值为1，anj表示第n条射线在第j个网格内的路径长度；步骤403：计算第k个炮点对应的接收点的实际初至波旅行时Tn和理论旅行时之间的残差：其中，表示在第i次迭代中第k个炮点对应的接收点的实际初至波旅行时Tn和理论旅行时之间的残差；步骤404：计算每个网格的慢度修正值：其中，表示针对第k个炮点在第i次迭代中第j个网格的慢度修正值，Mj表示射线经过第j个网格的次数，即射线密度；步骤405：循环执行步骤402至404，针对所有炮点求出第j个网格的慢度修正值，并计算所有炮点的第j个网格的慢度修正值的和步骤406：对初始速度模型的每个网格的慢度值进行修正：其中，表示在第i次迭代中第j个网格的慢度值；步骤407：重复步骤402至406，对i进行迭代，直到残差小于预定值。根据本专利技术的另一方面，提出了一种井间及井地联合地震波CT成像装置，包括：井间及井地联合观测系统，所述井间及井地联合观测系统包括设置于井间的井间检波器和炮点，以及设置于地面的地面检波器；实际初至波旅行时拾取单元，用于基于所述井间检波器和地面检波器拾取实际初至波旅行时；初始速度模型建立单元，用于根据测井资料建立初始速度模型，根据所述初始速度模型计算理论旅行时和射线路径；反演单元，用于以所述实际初至波旅行时与所述理论旅行时的残差作为目标函数，基于SIRT算法迭代反演速度修正值，对所述初始速度模型进行修正。优选地，通过滑动时窗能量比法拾取所述实际初至波旅行时单元。优选地，基于SIRT算法迭代反演速度修正值，对所述初始速度模型进行修正包括以下步骤：步骤401：以Sj表示初始模型中的第j个网格的慢度，第k个炮点对应Nk个接收点，其中，j＝1,2,…,M，M为总的网格数量，k＝1,2,…,NUM，NUM为总的炮点数量；步骤402：根据射线追踪方法计算第k个炮点对应的接收点的理论旅行时：其中，表示在第i次迭代中第k个炮点对应的接收点的理论旅行时，i表示迭代次数，i的初始取值为1，anj表示第n条射线在第j个网格内的路径长度；步骤403：计算第k个炮点对应的接收点的实际初至波旅行时Tn和理论旅行时之间的残差：其中，表示在第i次迭代中第k个炮点对应的接收点的实际初至波旅行时Tn和理论旅行时之间的残差；步骤404：计算每个网格的慢度修正值：其中，表示针对第k个炮点在第i次迭代中第j个网格的慢度修正值，Mj表示射线经过第j个网格的次数，即射线密度；步骤405：循环执行步骤402至404，针对所有炮点求出第j个网格的慢度修正值，并计算所有炮点的第j个网格的慢度修正值的和步骤406：对初始速度模型的每个网格的慢度值进行修正：其中，表示在第i次迭代中第j个网格的慢度值；步骤407：重复步骤402至406，对i进行迭代，直到残差小于预定值。根据本专利技术的井间及井地联合地震波CT成像方法及装置，其优点在于：同时设置井间检波器和地面检波器，能够提高数据的观测角度，增加反演区域内射线覆盖密度和穿过异常体垂向及横向的射线，提高反演精度和改善层析成像的分辨率，实现空间三维探测；井间及井地联合地震波CT成像方法，激发点都位于井中，地震波能量传播距离短，与探测目标体接近，同时降低了低速带对波能量的吸收，地震波具有较高的频率，初至波较清晰；利用滑动时窗能量比法进行人机交互拾取野外实际初至波旅行时本文档来自技高网...

【技术保护点】
一种井间及井地联合地震波CT成像方法，包括以下步骤：步骤1：建立井间及井地联合观测系统，所述井间及井地联合观测系统包括设置于井间的井间检波器和炮点，以及设置于地面的地面检波器；步骤2：通过所述井间检波器和地面检波器拾取实际初至波旅行时；步骤3：根据测井资料建立初始速度模型，根据所述初始速度模型计算理论旅行时和射线路径；步骤4：以所述实际初至波旅行时与所述理论旅行时的残差作为目标函数，基于SIRT算法迭代反演速度修正值，对所述初始速度模型进行修正。

【技术特征摘要】
1.一种井间及井地联合地震波CT成像方法，包括以下步骤：步骤1：建立井间及井地联合观测系统，所述井间及井地联合观测系统包括设置于井间的井间检波器和炮点，以及设置于地面的地面检波器；步骤2：通过所述井间检波器和地面检波器拾取实际初至波旅行时；步骤3：根据测井资料建立初始速度模型，根据所述初始速度模型计算理论旅行时和射线路径；步骤4：以所述实际初至波旅行时与所述理论旅行时的残差作为目标函数，基于SIRT算法迭代反演速度修正值，对所述初始速度模型进行修正。2.根据权利要求1所述的井间及井地联合地震波CT成像方法，其中，通过滑动时窗能量比法拾取所述实际初至波旅行时。3.根据权利要求1所述的井间及井地联合地震波CT成像方法，其中，基于SIRT算法迭代反演速度修正值，对所述初始速度模型进行修正包括：步骤401：以Sj表示初始模型中的第j个网格的慢度，第k个炮点对应Nk个接收点，其中，j＝1,2,…,M，M为总的网格数量，k＝1,2,…,NUM，NUM为总的炮点数量；步骤402：根据射线追踪方法计算第k个炮点对应的接收点的理论旅行时：其中，表示在第i次迭代中第k个炮点对应的接收点的理论旅行时，i表示迭代次数，i的初始取值为1，anj表示第n条射线在第j个网格内的路径长度；步骤403：计算第k个炮点对应的接收点的实际初至波旅行时Tn和理论旅行时之间的残差：其中，表示在第i次迭代中第k个炮点对应的接收点的实际初至波旅行时Tn和理论旅行时之间的残差；步骤404：计算每个网格的慢度修正值：其中，表示针对第k个炮点在第i次迭代中第j个网格的慢度修正值，Mj表示射线经过第j个网格的次数，即射线密度；步骤405：循环执行步骤402至404，针对所有炮点求出第j个网格的慢度修正值，并计算所有炮点的第j个网格的慢度修正值的和步骤406：对初始速度模型的每个网格的慢度值进行修正：其中，表示在第i次迭代中第j个网格的慢度值；步骤407：重复步骤402至406，对i进行迭代，直到残差小于预定值。4.一种井间及井地...

【专利技术属性】
技术研发人员：须振华，周单，林正良，滕龙，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：北京,11