反演二维速度模型的方法及装置.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN108845350A(43)申请公布日2018.11.20(21)申请号201810214398.7(22)申请日2018.03.15(71)申请人中国石油天然气集团有限公司地址100007北京市东城区东直门北大街9号申请人中国石油集团东方地球物理勘探有限责任公司(72)发明人张文波詹毅李合群陈宇峰朱旭光于永来(74)专利代理机构北京三友知识产权代理有限公司11127代理人李辉刘飞(51)Int.Cl.G01V1/28(2006.01)G01V1/30(2006.01)权利要求书4页说明书13页附图5页(54)发明名称反演二维速度模型的方法及装置(57)摘要本申请提供了一种反演二维速度模型的方法及装置，该方法包括：确定炮点到各个接收点的观测直达波纵波旅行时、观测反射波纵波旅行时及理论直达波纵波旅行时；构造观测直达波纵波旅行时与理论直达波纵波旅行时的残差目标函数并确定其极小值，以获取井旁初始速度模型中各个反射层的最优层慢度；根据最优层慢度确定每个反射界面的纵波速度模型；网格化井旁二维地质模型；确定网格化后的井旁二维地质模型中，每个反射界面对应的接收点上的理论反射波纵波旅行时；根据观测直达波纵波旅行时、理论直达波纵波旅行时、观测反射波纵波旅行时和理论反射波纵波旅行时，获取以多项式描述的反射界面和纵波速度。本申请可获得更加精确的VSP地震成像数据。CN108845350ACN108845350A权利要求书1/4页1.一种反演二维速度模型的方法，其特征在于，包括：确定炮点到各个接收点的观测直达波纵波旅行时、观测反射波纵波旅行时及理论直达波纵波旅行时；构造所述观测直达波纵波旅行时与所述理论直达波纵波旅行时的残差目标函数；确定所述残差目标函数的极小值，以获取预设的井旁初始速度模型中各个反射层的最优层慢度；根据所述最优层慢度确定每个反射界面的纵波速度模型；网格化预设的井旁二维地质模型，将所述纵波速度模型填充到相应的网格位置上，并根据炮点和接收点的空间坐标，将炮点和接收点映射到相应的网格位置上；确定网格化后的井旁二维地质模型中，每个反射界面对应的接收点上的理论反射波纵波旅行时；根据所述观测直达波纵波旅行时、所述理论直达波纵波旅行时、所述观测反射波纵波旅行时和所述理论反射波纵波旅行时，获取以多项式描述的反射界面和纵波速度。2.如权利要求1所述的反演二维速度模型的方法，其特征在于，在所述获取以多项式描述的反射界面和纵波速度之后，还包括：对所述以多项式描述的反射界面和纵波速度进行优化。3.如权利要求1所述的反演二维速度模型的方法，其特征在于，所述确定炮点到各个接收点的理论直达波纵波旅行时，包括：确定井旁初始速度模型的反射界面深度；根据公式确定所述井旁初始速度模型的慢度；根据公式确定炮点到各个接收点的理论直达波纵波旅行时；其中，T(z)为深度z处的理论直达波纵波旅行时，αi为第i个反射界面的慢度，Hi和Hi-1分别为第i个和第i-1个反射界面的深度，和分别为第i个和第i-1个反射界面对应深度的接收点的直达波纵波旅行时间，T0为炮点到第1个接收点的直达波纵波旅行时，M为井旁初始速度模型中反射界面的数量。4.如权利要求1所述的反演二维速度模型的方法，其特征在于，所述残差目标函数包括：其中，Φ为残差目标函数，T(hj)为深度hj处的理论直达波纵波旅行时，为到达第j个接收点的观测直达波纵波旅行时，hj为第j个接收点的深度，αi为第i个反射界面的慢度，Hi为第i个反射界面的深度，T0为炮点到第1个接收点的直达波纵波旅行时，N为接收点数量，M为井旁初始速度模型中反射界面的数量。2CN108845350A权利要求书2/4页5.如权利要求1所述的反演二维速度模型的方法，其特征在于，所述根据所述最优层慢度确定每个反射界面的纵波速度模型，包括：根据公式确定直达波到达所述井旁初始速度模型中每个反射界面上的直达波纵波旅行时；根据公式确定每个反射界面的纵波速度模型；其中，Tk(Hk)为深度Hk处的理论直达波纵波旅行时，T0为炮点到第1个接收点的直达波纵波旅行时，M为井旁初始速度模型中反射界面的数量，αi为第i个反射界面的慢度，Hi和Hi-1分别为第i个和第i-1个反射界面的深度，Hk为第k个反射界面的深度，Vi为第i个反射界面的纵波速度模型，Ti和Ti-1分别为炮点到第i个和第i-1个反射界面对应深度的接收点的直达波纵波旅行时。6.如权利要求1所述的反演二维速度模型的方法，其特征在于，所述确定网格化后的井旁二维地质模型中，每个反射界面对应的接收点上的理论反射波纵波旅行时，包括：根据公式确定网格化后的井旁二维地质模型中射线与反射界面的交点；以所述交点作为接收点，确定网格化后的井旁二维地质模型中炮点到反