CN106990438B - 复杂地表区域微测井点位及井深确定方法及装置 - Google Patents
复杂地表区域微测井点位及井深确定方法及装置 Download PDFInfo
- Publication number
- CN106990438B CN106990438B CN201710191079.4A CN201710191079A CN106990438B CN 106990438 B CN106990438 B CN 106990438B CN 201710191079 A CN201710191079 A CN 201710191079A CN 106990438 B CN106990438 B CN 106990438B
- Authority
- CN
- China
- Prior art keywords
- work area
- micro logging
- depth
- well
- surface velocity
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6169—Data from specific type of measurement using well-logging
Abstract
本申请实施例提供了一种复杂地表区域微测井点位及井深确定方法及装置,该方法包括:对工区内指定地震勘探采集数据中的每条测线进行初至走时层析反演,获取每条测线的近地表速度场;基于所述工区内已有的微测井解释成果,标定每个近地表速度场的高速顶界面的速度值;根据每个近地表速度场的高速顶界面的速度值,确定每个近地表速度场所对应测线上各个炮点和检波点的风化层厚度;根据每个近地表速度场所对应测线上各个炮点和检波点的风化层厚度,预测所述工区的风化层厚度;根据所述工区的风化层厚度以及预设的微测井密度,确定所述工区的微测井参数。本申请实施例可以提高表层调查效率,降低表层调查成本。
Description
技术领域
本申请涉及微测井技术,尤其是涉及一种复杂地表区域微测井点位及井深确定方法及装置。
背景技术
在地球物理勘探中,表层调查是野外地震采集过程中的首道工序,也是必不可少的一项基础工作。通过表层调查的解释结果内插后得到的近地表模型,能够为野外静校量的计算、激发井深模型设计以及检波器组合高差设计提供基础数据,同时还可以为叠前深度偏移处理提供准确的近地表速度场。由此可见,表层调查工序直接关系到后续地震采集过程中的激发、接收、现场处理及叠前深度偏移处理等地震勘探的各个方面,是确保野外地震采集效果不可或缺的一个关键环节。
目前表层调查基本以小折射调查为主,但是由于其应用原理的限制,小折射调查只适用于地表平坦、折射界面倾角较小、速度为正序分布的层状介质情况。对于复杂地表区,因该方法的解释精度很低,一般无法满足勘探需求。随着勘探精度要求的不断提高,为了保证近地表模型的精度,在复杂地表区进行微测井已成为目前主要的表层调查方法。
在微测井技术中,微测井参数(例如井深、井点位置)选取十分关键,参数的选取是否科学合理,将直接影响微测井的效果。
目前微测井井深的设计往往是通过对工区低降速带变化规律的了解,根据以往小折射资料的解释成果进行设计。由于在复杂地表区小折射的调查精度很低,致使微测井井深设计出现了较大的偏差:(1)如果在钻机钻井能力极限范围内无法钻至高速层的地方开展微测井调查,则不能获得近表层结构,因而失去了表层调查的意义;(2)如果工区内的风化层厚度很小(比如只有20m),而微测井钻井深度设计很大(比如设计120m),即钻井钻至了不必要的深度,从而增大了表层调查成本;(3)如果工区内的风化层厚度较大(比如120m),微测井钻井深度设计较小(比如微测井井深设计为90m,则不能获得近表层结构,需重新进行微测井施工,从而影响了表层调查的效率。
目前微测井的井点位置布设方法主要是在规则密度的基础上,对地表岩性、速度变化处针对性地加密。但是这种布设方法有一定的弊端,若是在同一地表岩性的较大区域范围内,其风化层厚度超出了钻机钻井能力的极限,如果仍然按照规则密度实施微测井调查,那么这个区域间的每口微测井调查都不能获得所需的高速层速度及其顶界面至地表之间的风化层厚度和速度,从而不能实现表层调查的目的。
由于表层调查是地震采集的首道工序,如果该工序未能完成,则后续工序无法开展。因此,如何合理确定复杂地表区域微测井参数是目前亟待解决的技术问题。
发明内容
本申请实施例的目的在于提供一种复杂地表区域微测井点位及井深确定方法及装置,以提高表层调查效率,降低表层调查成本。
为达到上述目的,一方面,本申请实施例提供了一种复杂地表区域微测井点位及井深确定方法,包括以下步骤:
对工区内指定地震勘探采集数据中的每条测线进行初至走时层析反演,获取每条测线的近地表速度场;
基于所述工区内已有的微测井解释成果,标定每个近地表速度场的高速顶界面的速度值;
根据每个近地表速度场的高速顶界面的速度值,确定每个近地表速度场所对应测线上各个炮点和检波点的风化层厚度;
根据每个近地表速度场所对应测线上各个炮点和检波点的风化层厚度,预测所述工区的风化层厚度;
根据所述工区的风化层厚度以及预设的微测井密度,确定所述工区的微测井参数。
本申请实施例的复杂地表区域微测井点位及井深确定方法,所述基于所述工区内已有的微测井解释成果,标定每个近地表速度场的高速顶界面的速度值,包括:
a、根据所述工区内的已有微测井解释成果,确定该微测井位置处的风化层厚度;
b、选取一个近地表速度场;
c、将选取的近地表速度场对应于所述微测井井点位置处的地表高程值,减去所述微测井井点位置处的风化层厚度,得到所述选取的近地表速度场对应于所述微测井井点位置处的新的高程值;
d、在所述选取的近地表速度场中,确定所述新的高程值处的速度值,并将其作为所述选取的近地表速度场的高速顶界面的速度值;
e、重复上述步骤b~d,直至获取每个近地表速度场的高速顶界面的速度值。
本申请实施例的复杂地表区域微测井点位及井深确定方法,所述根据每个近地表速度场的高速顶界面的速度值,确定每个近地表速度场所对应测线上各个炮点和检波点的风化层厚度,包括:
确定每个近地表速度场所对应测线上各个炮点和检波点的地表高程值;
根据每个近地表速度场的高速顶界面的速度值,确定每个近地表速度场所对应测线上各个炮点和检波点的高速顶界面高程值;
将每个近地表速度场所对应测线上各个炮点和检波点的地表高程值,减去对应位置处的高速顶界面高程值,得到每个近地表速度场所对应测线上各个炮点和检波点的风化层厚度。
本申请实施例的复杂地表区域微测井点位及井深确定方法,所述根据每个近地表速度场所对应测线上各个炮点和检波点的风化层厚度,预测所述工区的风化层厚度,包括:
按照预设的网格参数将所述工区的覆盖范围进行网格化,得到网格化工区;
以每个近地表速度场所对应测线上各个炮点和检波点的风化层厚度作为离散数据;通过预设的内插方法对所述网格化工区进行数据平面内插,预测出所述工区的风化层厚度。
本申请实施例的复杂地表区域微测井点位及井深确定方法,所述预设的内插方法包括克里金内插方法。
本申请实施例的复杂地表区域微测井点位及井深确定方法,所述根据所述工区的风化层厚度以及预设的微测井密度,确定所述工区的微测井参数,包括:
根据预设的微测井密度选定所述工区内各个微测井的井点位置;
根据所述工区的风化层厚度确定所述工区内各个微测井的井点位置处的风化层厚度;
将所述工区内各个微测井的井点位置处的风化层厚度分别加上预设补偿深度值,得到所述工区内每个微测井的井点位置处的井深值。
本申请实施例的复杂地表区域微测井点位及井深确定方法,在所述得到所述工区内每个微测井的井点位置处的井深值之后,还包括:
判断所述工区内每个微测井的井深值是否大于预设的钻井深度极限值;
如果大于,则剔除对应微测井的井点位置。
本申请实施例的复杂地表区域微测井点位及井深确定方法,还包括:
如果所述工区内某个微测井的井深值不大于所述钻井深度极限值,则保留该微测井的井点位置。
本申请实施例的复杂地表区域微测井点位及井深确定方法,所述预设补偿深度值的取值范围包括10米~20米。
另一方面,本申请实施例还提供了一种复杂地表区域微测井点位及井深确定装置,包括:
速度场获取模块,用于对工区内指定地震勘探采集数据中的每条测线进行初至走时层析反演,获取每条测线的近地表速度场;
速度值标定模块,用于基于所述工区内已有的微测井解释成果,标定每个近地表速度场的高速顶界面的速度值;
第一厚度确定模块,用于根据每个近地表速度场的高速顶界面的速度值,确定每个近地表速度场所对应测线上各个炮点和检波点的风化层厚度;
第二厚度确定模块,用于根据每个近地表速度场所对应测线上各个炮点和检波点的风化层厚度,预测所述工区的风化层厚度;
微测井参数确定模块,用于根据所述工区的风化层厚度以及预设的微测井密度,确定所述工区的微测井参数。
由以上本申请实施例提供的技术方案可见,本申请实施例首先对工区内指定地震勘探采集数据中的每条测线进行初至走时层析反演,获取每条测线的近地表速度场;然后基于工区内已有的微测井解释成果,标定每个近地表速度场的高速顶界面的速度值;并根据每个近地表速度场的高速顶界面的速度值预测工区的风化层厚度,最后根据工区的风化层厚度以及预设的微测井密度,确定工区的微测井参数。从而利用工区的测线的近地表速度场来指导工区内微测井的井点位置布设及井深的设计,避免了现有技术在微测井点布设的盲目性,提高了微测井井深设计的准确性,从而达到提高表层调查效率以及降本增效的目的,本申请实施例可适用于任何复杂地表区的表层调查。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。在附图中:
图1为本申请一实施例的复杂地表区域微测井点位及井深确定方法的流程图;
图2为本申请一实施例的三维工区及二维测线位置分布示意图;
图3a为对图2所示主测线04进行初至走时层析反演得到的近地表速度场;
图3b为对图2所示联络线03进行初至走时层析反演得到的近地表速度场;
图4为本申请一实施例中标定微测井时距曲线示意图;
图5为本申请一实施例中,通过已有的微测井解释成果标定图2所示联络测线14的近地表速度场的高速顶界面的地表高程示意图;
图6为本申请一实施例预测的三维工区的风化层厚度平面示意图;
图7为本申请一实施例的三维工区内微测井的井点位置分布示意图;
图8为图7所示三维工区内各个微测井井点位置的井深分布示意图;
图9为本申请一实施例中复杂地表区域微测井点位及井深确定装置的结构框图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
参考图1所示,本申请实施例的复杂地表区域微测井点位及井深确定方法可以包括以下步骤:
步骤S101、对工区内指定地震勘探采集数据中的每条测线进行初至走时层析反演,获取每条测线的近地表速度场。
本申请一实施例中,所述的工区例如可以为三维工区。所述的测线例如可以为以往某次地震勘探采集数据中的二维测线;每条测线上可以布设有若干个炮点和/或检波点。所述对工区内指定地震勘探采集数据中的每条测线进行初至走时层析反演是指,在对工区内每条二维测线进行初至拾取的基础上,对其中每条测线进行初至波走时层析速度反演,从而得到每条测线的近地表速度场。
本申请一些实施例中,初至走时层析反演的网格参数在横向上可以为0.5~4倍道距,纵向上可以为2~20米(m),初至的偏移距范围可以采用1000~6000米(m),一般的通过8~10次迭代反演后,即获得每条测线的(达到一定精度要求)近地表速度场。
在本申请一个示例性实施例中,一个三维工区及二维测线位置分布可以如图2所示;图2中,相对较细的线为三维工区的边界线,这些线围成了一个大致呈长方形的三维工区;相对较粗的线为测线,纵向的测线为主测线(例如主测线04等),其他方向的测线为联络线(例如联络线03、联络线14等)。在对图2所示三维工区内的主测线04进行初至走时层析反演,可获取对应的近地表速度场如图3a所示;在对图2所示三维工区内的联络线03进行初至走时层析反演,可获取对应的近地表速度场如图3b所示。在图3a和图3b中,纵坐标为深度,横坐标为桩号。
步骤S102、基于所述工区内已有的微测井解释成果,标定每个近地表速度场的高速顶界面的速度值。
本申请一实施例中,所述基于所述工区内已有的微测井解释成果,标定每个近地表速度场的高速顶界面的速度值,可以包括如下步骤:
a、根据所述工区内的已有微测井解释成果,确定该微测井位置处的风化层厚度。当然,如果工区内没有微测井可以先实施一口微测井,从而可根据工区内的已有微测井解释成果得到该微测井的井点位置处的风化层厚度。
b、选取一个近地表速度场;所述的选取例如可以为顺序选取或随机选取。
c、将选取的近地表速度场对应于所述微测井井点位置处的地表高程值,减去所述微测井井点位置处的风化层厚度,得到所述选取的近地表速度场对应于所述微测井井点位置处的新的高程值。
在本申请一个示例性实施例中,例如图4所示,微测井时距曲线将地下划分为L1~L5五层状结构,横坐标为垂直时,单位为ms;纵坐标为深度,单位为m;图中的小黑点表示在地下某个深度激发,其产生的垂直时间;小黑点的连线是按照线性拟合获得的层位划分,V表示该层的层速度,H表示该层的厚度。根据微测井的解释成果,所述微测井井点位置处的风化层厚度为10.3+61.6+60.4+21.8=154.1(m),而该微测井的地表高程为1569.8m,在所述微测井井点位置处:地表高程减去风化层厚度值后得到一个新的高程值1569.8-154.1=1415.7(m)。
d、在所述选取的近地表速度场中,确定所述新的高程值处的速度值,并将其作为所述选取的近地表速度场的高速顶界面的速度值。
本申请一实施例中,一般的,近地表速度场的数学表达方式可用(x,y,z,v)来表示,意思为二维坐标位置的某一深度处所对应的速度值,也就是三维位置与速度的对应关系。近地表速度场能够体现不同位置随着深度变化时速度发生的变化,反映了近地表和深层地表条件的变化规律,此时需要一个高速顶界面(也是风化层底界面)来区分近地表和深层地表,而这个界面需要给定具体的速度值才能提取出来。通过微测井的时距曲线可知风化层厚度值,该位置的地表高程与风化层厚度差就是该位置的高速顶界面的高程值,在速度场中会具体对应一个速度值。在确定了近地表速度场的条件下,那么只要给定了具体的位置(x,y,z),则必然能够找到对应的速度值。
在本申请一个示例性实施例中,例如图5所示,在第14条联络测线的近地表速度场中,微测井井点位置处此高程值对应的速度值为1939m/s。
e、重复上述步骤b~d,直至获取每个近地表速度场的高速顶界面的速度值。
步骤S103、根据每个近地表速度场的高速顶界面的速度值,确定每个近地表速度场所对应测线上各个炮点和检波点的风化层厚度。
本申请一实施例中,所述根据每个近地表速度场的高速顶界面的速度值,确定每个近地表速度场所对应测线上各个炮点和检波点的风化层厚度,可以包括如下步骤:
首先,确定每个近地表速度场所对应测线上各个炮点和检波点的地表高程值。
其次,根据每个近地表速度场的高速顶界面的速度值,确定每个近地表速度场所对应测线上各个炮点和检波点的高速顶界面高程值。
然后,将每个近地表速度场所对应测线上各个炮点和检波点的地表高程值,减去对应位置处的高速顶界面高程值,得到每个近地表速度场所对应测线上各个炮点和检波点的风化层厚度。
在本申请一个示例性实施例中,请继续参考图5所示,在第14条联络测线的近地表速度场中,微测井井点位置处此高程值对应的速度值为1939m/s,可以取该速度值的整数1900m/s(严谨来讲应该采用1939m/s而不取整,但是在实际提取高速顶界面的操作过程中,除了给定具体的速度值外,还要给一个速度值的变化范围,而这个变化范围一般采用100m/s,所以才将1939m/s取整为1900m/s)的顶界面做为对应近地表速度场的高速顶界面,则该顶界面在微测井井点位置处的高程值为1418.8m。与微测井的解释成果误差为154.1-151.0=3.1m(2.01%)。可见,其精度是能够满足勘探需求的。
步骤S104、根据每个近地表速度场所对应测线上各个炮点和检波点的风化层厚度,预测所述工区的风化层厚度。
本申请一实施例中,所述根据每个近地表速度场所对应测线上各个炮点和检波点的风化层厚度,预测所述工区的风化层厚度,可以包括如下步骤:
首先,按照预设的网格参数将所述工区的覆盖范围进行网格化,得到网格化工区。其中,网格化可以为均匀网络化,即X、Y方向上均采用250~2000m的网格参数。
其次,以每个近地表速度场所对应测线上各个炮点和检波点的风化层厚度作为离散数据;通过预设的内插方法对所述网格化工区进行数据平面内插,预测出所述工区的风化层厚度。在本申请一些示例性实施例中,所述内插方法例如可以为克里金内插方法、线性差值三角法、反距离加权插值法等。
在本申请一个示例性实施例中,基于上述方法可预测得到一个三维工区的风化层厚度(如图6所示)。
步骤S105、根据所述工区的风化层厚度以及预设的微测井密度,确定所述工区的微测井参数。
本申请一实施例中,所述根据所述工区的风化层厚度以及预设的微测井密度,确定所述工区的微测井参数,可以包括如下步骤:
首先,根据预设的微测井密度选定所述工区内各个微测井的井点位置。
其次,根据所述工区的风化层厚度确定所述工区内各个微测井的井点位置处的风化层厚度。
然后,将所述工区内各个微测井的井点位置处的风化层厚度分别加上预设补偿深度值,得到所述工区内每个微测井的井点位置处的井深值。
需要说明的是,基于微测井表层调查需要穿过风化层进入高速层才能获得高速层的速度。而由于预测的风化层厚度与实际风化层厚度之间可能会存在一定误差,为了能保证拟合高速层的控制点可以在3个以上,需要在预测的风化层厚度的基础上加上预设补偿深度,以保证微测井的解释精度。
在本申请另一个实施例中,在所述得到所述工区内每个微测井的井点位置处的井深值之后,还可以包括:
判断所述工区内每个微测井的井深值是否大于预设的钻井深度极限值;
如果所述工区内某个微测井的井深值大于所述钻井深度极限值,则剔除对应微测井的井点位置;否则,如果所述工区内某个微测井的井深值不大于所述钻井深度极限值,则保留该微测井的井点位置。在本申请一个示例性实施例中,按照上述方法可得到一个三维工区内的每个微测井的井点位置(见图7中各个黑色点)。其中,该三维工区内两个较大的区块内没有微测井分布,是因为这两个区块超出了钻井深度极限值,无法获钻至高速层,因而没有必要在此布设微测井。
在本申请一个示例性实施例中,基于上述方法最终可预测得到一个三维工区内各个微测井井点位置的井深分布(如图8所示)。
在本申请一个实施例中,所述预设补偿深度值的取值范围例如可以为10米~20米。
在本申请一个实施例中,一般的,对于平地,所述钻井深度极限值可以为170m;对于山体等复杂地表,所述钻井深度极限值极限值可以为50m。
虽然上文描述的过程流程包括以特定顺序出现的多个操作,但是,应当清楚了解,这些过程可以包括更多或更少的操作,这些操作可以顺序执行或并行执行(例如使用并行处理器或多线程环境)。
参见图9所示,本申请实施例的复杂地表区域微测井点位及井深确定装置,可以包括:
速度场获取模块91,可以用于对工区内指定地震勘探采集数据中的每条测线进行初至走时层析反演,获取每条测线的近地表速度场;
速度值标定模块92,可以用于基于所述工区内已有的微测井解释成果,标定每个近地表速度场的高速顶界面的速度值;
第一厚度确定模块93,可以用于根据每个近地表速度场的高速顶界面的速度值,确定每个近地表速度场所对应测线上各个炮点和检波点的风化层厚度;
第二厚度确定模块94,可以用于根据每个近地表速度场所对应测线上各个炮点和检波点的风化层厚度,预测所述工区的风化层厚度;
微测井参数确定模块95,可以用于根据所述工区的风化层厚度以及预设的微测井密度,确定所述工区的微测井参数。
本申请实施例的复杂地表区域微测井点位及井深确定装置与上述图1所示的方法实施例对应,因此,有关于本申请实施例的复杂地表区域微测井点位及井深确定装置的细节内容,请参见上述图1所示的方法实施例,在此不再赘述。
为了描述的方便,描述以上装置时以功能分为各种单元分别描述。当然,在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中,计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器,随机存取存储器(RAM)和/或非易失性内存等形式,如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括,但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带,磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。按照本文中的界定,计算机可读介质不包括暂存电脑可读媒体(transitory media),如调制的数据信号和载波。
还需要说明的是,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
本领域技术人员应明白,本申请的实施例可提供为方法、系统或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
以上所述仅为本申请的实施例而已,并不用于限制本申请。对于本领域技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等,均应包含在本申请的权利要求范围之内。
Claims (9)
1.一种复杂地表区域微测井点位及井深确定方法,其特征在于,包括以下步骤:
对工区内指定地震勘探采集数据中的每条测线进行初至走时层析反演,获取每条测线的近地表速度场;
基于所述工区内已有的微测井解释成果,标定每个近地表速度场的高速顶界面的速度值;
根据每个近地表速度场的高速顶界面的速度值,确定每个近地表速度场所对应测线上各个炮点和检波点的风化层厚度;
根据每个近地表速度场所对应测线上各个炮点和检波点的风化层厚度,预测所述工区的风化层厚度;
根据所述工区的风化层厚度以及预设的微测井密度,确定所述工区的微测井参数;其中,
所述基于所述工区内已有的微测井解释成果,标定每个近地表速度场的高速顶界面的速度值,包括:
a、根据所述工区内已有的微测井解释成果,确定该微测井井点位置处的风化层厚度;
b、选取一个近地表速度场;
c、将选取的近地表速度场对应于所述微测井井点位置处的地表高程值,同时减去所述微测井井点位置处的风化层厚度,得到所选取的近地表速度场对应所述微测井井点位置处的新的高程值;
d、在所选取的近地表速度场中,确定所述新的高程值处的速度值,并将其作为所选取的近地表速度场的高速顶界面的速度值;
e、重复上述步骤b~d,直至获取每个近地表速度场的高速顶界面的速度值。
2.根据权利要求1所述的复杂地表区域微测井点位及井深确定方法,其特征在于,所述根据每个近地表速度场的高速顶界面的速度值,确定每个近地表速度场所对应测线上各个炮点和检波点的风化层厚度,包括:
确定每个近地表速度场所对应测线上各个炮点和检波点的地表高程值;
根据每个近地表速度场的高速顶界面的速度值,确定每个近地表速度场所对应测线上各个炮点和检波点的高速顶界面高程值;
将每个近地表速度场所对应测线上各个炮点和检波点的地表高程值,减去对应位置处的高速顶界面高程值,得到每个近地表速度场所对应测线上各个炮点和检波点的风化层厚度。
3.根据权利要求1所述的复杂地表区域微测井点位及井深确定方法,其特征在于,所述根据每个近地表速度场所对应测线上各个炮点和检波点的风化层厚度,预测所述工区的风化层厚度,包括:
按照预设的网格参数将所述工区的覆盖范围进行网格化,得到网格化工区;
以每个近地表速度场所对应测线上各个炮点和检波点的风化层厚度作为离散数据;通过预设的内插方法对所述网格化工区进行数据平面内插,预测出所述工区的风化层厚度。
4.根据权利要求3所述的复杂地表区域微测井点位及井深确定方法,其特征在于,所述预设的内插方法包括克里金内插方法。
5.根据权利要求1所述的复杂地表区域微测井点位及井深确定方法,其特征在于,所述根据所述工区的风化层厚度以及预设的微测井密度,确定所述工区的微测井参数,包括:
根据预设的微测井密度选定所述工区内各个微测井的井点位置;
根据所述工区的风化层厚度确定所述工区内各个微测井的井点位置处的风化层厚度;
将所述工区内各个微测井的井点位置处的风化层厚度分别加上预设补偿深度值,得到所述工区内每个微测井的井点位置处的井深值。
6.根据权利要求5所述的复杂地表区域微测井点位及井深确定方法,其特征在于,在所述得到所述工区内每个微测井的井点位置处的井深值之后,还包括:
判断所述工区内每个微测井的井深值是否大于预设的钻井深度极限值;
如果大于,则剔除对应微测井的井点位置。
7.根据权利要求6所述的复杂地表区域微测井点位及井深确定方法,其特征在于,还包括:
如果所述工区内某个微测井的井深值不大于所述钻井深度极限值,则保留该微测井的井点位置。
8.根据权利要求5所述的复杂地表区域微测井点位及井深确定方法,其特征在于,所述预设补偿深度值的取值范围为10米~20米。
9.一种复杂地表区域微测井点位及井深确定装置,其特征在于,包括:
速度场获取模块,用于对工区内指定地震勘探采集数据中的每条测线进行初至走时层析反演,获取每条测线的近地表速度场;
速度值标定模块,用于基于所述工区内已有的微测井解释成果,标定每个近地表速度场的高速顶界面的速度值;
第一厚度确定模块,用于根据每个近地表速度场的高速顶界面的速度值,确定每个近地表速度场所对应测线上各个炮点和检波点的风化层厚度;
第二厚度确定模块,用于根据每个近地表速度场所对应测线上各个炮点和检波点的风化层厚度,预测所述工区的风化层厚度;
微测井参数确定模块,用于根据所述工区的风化层厚度以及预设的微测井密度,确定所述工区的微测井参数;其中,
所述速度值标定模块基于所述工区内已有的微测井解释成果,标定每个近地表速度场的高速顶界面的速度值,包括:
a、根据所述工区内已有的微测井解释成果,确定该微测井井点位置处的风化层厚度;
b、选取一个近地表速度场;
c、将选取的近地表速度场对应于所述微测井井点位置处的地表高程值,同时减去所述微测井井点位置处的风化层厚度,得到所选取的近地表速度场对应所述微测井井点位置处的新的高程值;
d、在所选取的近地表速度场中,确定所述新的高程值处的速度值,并将其作为所选取的近地表速度场的高速顶界面的速度值;
e、重复上述步骤b~d,直至获取每个近地表速度场的高速顶界面的速度值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710191079.4A CN106990438B (zh) | 2017-03-28 | 2017-03-28 | 复杂地表区域微测井点位及井深确定方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710191079.4A CN106990438B (zh) | 2017-03-28 | 2017-03-28 | 复杂地表区域微测井点位及井深确定方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106990438A CN106990438A (zh) | 2017-07-28 |
CN106990438B true CN106990438B (zh) | 2018-12-28 |
Family
ID=59413430
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710191079.4A Active CN106990438B (zh) | 2017-03-28 | 2017-03-28 | 复杂地表区域微测井点位及井深确定方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106990438B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109581493B (zh) * | 2018-10-22 | 2020-11-06 | 中国石油天然气集团有限公司 | 一种近地表层析模型的构建方法及系统 |
CN111399037B (zh) * | 2019-01-02 | 2022-12-02 | 中国石油天然气集团有限公司 | 高速顶界面提取的方法和装置 |
CN111427084A (zh) * | 2019-01-10 | 2020-07-17 | 中国石油天然气集团有限公司 | 一种基于微测井的近地表速度建模方法及系统 |
CN112052488A (zh) * | 2019-05-20 | 2020-12-08 | 中国石油天然气集团有限公司 | 陆地平台井位位置确定方法及装置 |
CN110501757B (zh) * | 2019-06-27 | 2021-02-26 | 中国石油化工股份有限公司 | 点线面三维逐点井深设计方法 |
CN112764098A (zh) * | 2019-11-04 | 2021-05-07 | 中国石油天然气集团有限公司 | 波场分析系统及方法 |
CN112946755A (zh) * | 2019-12-11 | 2021-06-11 | 中国石油天然气集团有限公司 | 近地表表层的调查方法和装置 |
CN111638555B (zh) * | 2020-05-29 | 2023-02-10 | 中国石油天然气集团有限公司 | 解释微测井的方法及装置 |
CN112346119B (zh) * | 2020-09-16 | 2022-07-05 | 杭州国家水电站大坝安全和应急工程技术中心有限公司 | 一种跨孔地震波层析成像地质解释方法 |
CN112379434A (zh) * | 2020-10-30 | 2021-02-19 | 中国石油天然气集团有限公司 | 适用于沙漠区的层析反演静校正的方法及装置 |
CN112285782B (zh) * | 2020-11-05 | 2023-04-25 | 中国石油天然气集团有限公司 | 近地表地震波吸收衰减调查方法及装置 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102749469A (zh) * | 2012-07-03 | 2012-10-24 | 西安交通大学 | 一种智能手机测量转速的方法 |
US8576663B2 (en) * | 2010-04-30 | 2013-11-05 | Schlumberger Technology Corporation | Multicomponent seismic inversion of VSP data |
CN104133249A (zh) * | 2014-07-31 | 2014-11-05 | 中国石油天然气集团公司 | 一种微测井资料与声波测井资料联合解释的方法及装置 |
CN104316964A (zh) * | 2014-10-29 | 2015-01-28 | 中国石油天然气集团公司 | 一种近地表表层的建模方法及系统 |
CN105277978A (zh) * | 2015-10-15 | 2016-01-27 | 中国石油天然气集团公司 | 一种确定近地表速度模型的方法及装置 |
CN105301638A (zh) * | 2015-10-12 | 2016-02-03 | 中国石油天然气集团公司 | 一种提取风化层底界面的方法和装置 |
WO2016123046A1 (en) * | 2015-01-26 | 2016-08-04 | Schlumberger Technology Corporation | Method for determining petrophysical properties from logging measurements |
-
2017
- 2017-03-28 CN CN201710191079.4A patent/CN106990438B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8576663B2 (en) * | 2010-04-30 | 2013-11-05 | Schlumberger Technology Corporation | Multicomponent seismic inversion of VSP data |
CN102749469A (zh) * | 2012-07-03 | 2012-10-24 | 西安交通大学 | 一种智能手机测量转速的方法 |
CN104133249A (zh) * | 2014-07-31 | 2014-11-05 | 中国石油天然气集团公司 | 一种微测井资料与声波测井资料联合解释的方法及装置 |
CN104316964A (zh) * | 2014-10-29 | 2015-01-28 | 中国石油天然气集团公司 | 一种近地表表层的建模方法及系统 |
WO2016123046A1 (en) * | 2015-01-26 | 2016-08-04 | Schlumberger Technology Corporation | Method for determining petrophysical properties from logging measurements |
CN105301638A (zh) * | 2015-10-12 | 2016-02-03 | 中国石油天然气集团公司 | 一种提取风化层底界面的方法和装置 |
CN105277978A (zh) * | 2015-10-15 | 2016-01-27 | 中国石油天然气集团公司 | 一种确定近地表速度模型的方法及装置 |
Non-Patent Citations (1)
Title |
---|
吕景峰等.复杂地表浅层层析表层调查方法及应用.《石油地球物理勘探》.2012,第335-338页. * |
Also Published As
Publication number | Publication date |
---|---|
CN106990438A (zh) | 2017-07-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106990438B (zh) | 复杂地表区域微测井点位及井深确定方法及装置 | |
Lelievre et al. | A comprehensive study of including structural orientation information in geophysical inversions | |
CN107765298B (zh) | 一种确定近地表速度模型的方法及装置 | |
CN106896406B (zh) | 基于阻抗域属性预测优质储层的方法和装置 | |
CN109870719B (zh) | 一种碳酸盐岩致密薄储层的井位布设方法、装置及系统 | |
CN107817525B (zh) | 一种近地表常速填充方法及装置 | |
CA3023655A1 (en) | Three-dimensional fracture abundance evaluation of subsurface formation based on geomechanical simulation of mechanical properties thereof | |
US10387583B2 (en) | Rotations from gradient directions | |
CN107783186A (zh) | 一种近地表速度建模方法及装置 | |
CN108845350A (zh) | 反演二维速度模型的方法及装置 | |
CN108614293A (zh) | 砂体预测方法及装置 | |
US10113421B2 (en) | Three-dimensional fracture abundance evaluation of subsurface formations | |
US20160377752A1 (en) | Method of Digitally Identifying Structural Traps | |
CN108445532B (zh) | 一种深度域反偏移方法及装置 | |
US20200049844A1 (en) | Computer implemented method for improving a velocity model for seismic imaging | |
CN109597964A (zh) | 理论油藏渗流参数值确定方法、系统、设备、可读介质 | |
US10650107B2 (en) | Three-dimensional subsurface formation evaluation using projection-based area operations | |
US10254441B2 (en) | Method of modelling a subsurface volume | |
CN112684498A (zh) | 一种基于宽方位地震数据的储层裂缝预测方法及系统 | |
US20090043549A1 (en) | Methods, apparatus, and products for seismic ray tracing | |
US20150378039A1 (en) | System and method for seismic imaging of a complex subsurface | |
CN112292714B (zh) | 基于断层辐射的网格隔分 | |
CN107229074A (zh) | 低勘探程度海相探区的变速成图方法和装置 | |
CN112817059A (zh) | 峰峰组内幕储层确定方法及装置 | |
CN113267810B (zh) | 地震勘探全深度速度建模方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |