生成储层的地下结构的速度模型和密度模型的制作方法

文档序号:28315330发布日期:2022-01-01 01:27阅读:205来源:国知局
生成储层的地下结构的速度模型和密度模型的制作方法
生成储层的地下结构的速度模型和密度模型
1.要求优先权
2.本技术要求于2019年3月21日递交的美国专利申请no.16/360,174的优先权,其全部内容通过引用并入本文。
技术领域
3.本公开涉及生成储层的地下结构的速度模型和密度模型。


背景技术:

4.开发用于储层的地下结构的准确速度模型和密度模型是储层地质分析的重要方面。沉积岩随着时间通过沉积形成。岩石层由不同的材料或在不同的条件下形成。因此,每一层岩石可能具有不同的属性。速度模型绘制出地理区域中的岩石层和地震波(例如,p波或s波)将在每一层中传播的预期速度。密度模型绘制出地理区域中岩层的密度分布。


技术实现要素:

5.本公开描述了包括可以用于生成速度模型和密度模型的计算机实施方法、计算机程序产品和计算机系统的方法和系统。在一些实施方式中,用于生成速度模型和密度模型的计算机实现的方法包括由数据处理装置接收油气储层的一组地震数据;由数据处理装置设置初始速度模型和初始密度模型;由所述数据处理装置基于所述一组地震数据生成所述油气储层的波场;由数据处理装置选择空间方向;由数据处理装置基于波场产生所选择的空间方向的速度梯度和反射率梯度;以及由数据处理装置使用选择的空间方向的速度梯度和反射率梯度更新速度模型和密度模型。
6.该方面的其他实施方式包括相应的计算机系统、装置和记录在一个或多个计算机存储设备上的计算机程序,它们均被配置为执行方法的动作。一个或多个计算机的系统可以被配置为通过在系统上安装的、在操作中导致系统执行动作的软件、固件、硬件或软件、固件或硬件的组合来执行特定的操作或动作。一个或多个计算机程序可以被配置为通过包括指令来执行特定操作或动作,所述指令在被数据处理装置执行时使得该装置执行动作。
7.在具体实施方式、附图和权利要求书中阐述本说明书的主题的一个或多个实施方式的细节。通过具体实施方式、权利要求书和附图,所述主题的其他特征、方面和优点将变得明显。
附图说明
8.图1示出了根据实施方式的示例波场传播模型。
9.图2示出了根据实施方式的基于在一个方向上的速度梯度和反射率梯度来更新速度模型和密度模型的示例过程。
10.图3示出了根据实施方式的基于迭代频带生成速度模型和密度模型的示例过程。
11.图4示出了根据实施方式的示例速度模型和密度模型生成过程。
12.图5a和图5b示出了根据实施方式的合成速度模型和合成密度模型。
13.图6a和图6b示出了根据实施方式的初始速度模型和初始密度模型。
14.图7a和图7b示出了根据实施方式的基于本公开中讨论的过程获得的速度模型和密度模型。
15.图8示出了根据实施方式的炮点道集的示例图表。
16.图9是示出了根据一些实施方式的用于提供与如本公开中描述的算法、方法、功能、处理、流程和过程相关联的计算功能的示例计算机系统的框图。
17.各附图中相似的附图标记和标志指示相似的要素。
具体实施方式
18.以下具体实施方式描述了用于生成速度模型和密度模型的技术。在不背离本公开的范围的情况下,对所公开的实施方式的各种修改、改变和置换可以被实现并且对本领域普通技术人员而言将显而易见,并且本公开中所定义的一般原理可适用于其他实施方式和应用。在一些情况下,可以省略对于获得所描述的主题的理解不必要的细节,以便不会用不必要的细节混淆一个或多个所描述的实施方式,并且因为这样的细节在本领域普通技术人员的技能范围内。本公开并非意在限于所描述的或示出的实施方式,而应赋予与本公开中描述的原理和特征一致的最宽范围。
19.本公开总体上描述了用于生成速度模型和密度模型的方法和系统,包括计算机实现的方法、计算机程序产品和计算机系统。速度模型和密度模型基于地震数据生成。在一些情况下,地震信号可以通过振源设备发送到源位置处的土地的地下。地震信号的示例包括声学信号。地震信号传播通过地下,并且可以被置于接收器位置的接收器设备接收。在一些情况下,源设备、接收器设备或其组合可以置于地表上。信号可以向下传播,直到到达反射结构并向上朝向地表反射。因为信号通过地下被折射、反射,因此所接收信号的特性包含地下结构的信息。
20.在海上勘测中,气枪和水听器可以分别用作源设备和接收器设备。在采集期间,地震源从气枪阵列中爆炸。被反射和被折射的信号由水听器的拖缆获取。在陆地采集中,可以使用炸药作为爆炸源,并且采用地震检波器作为接收器设备。在另一示例中,振动卡车可以用作振源设备。也可以使用产生和接收地震信号的其他设备。
21.在一些情况下,接收的信号数据可以包括在接收器位置处由接收器设备接收的接收器信号数据的时域样本。例如,在地震采集期间,源(炸药、振动卡车、气枪阵列等)被激活,来自地下地质边界的反射/折射/透射由位于地表上的接收器设备记录。顺序或同时针对每个炮点道重复进行这种类型的采集,直到针对该勘测区域采集了所有地震数据为止。这些采集的地震数据被称为接收器信号数据。在一些情况下,采集的地震数据在现场被收集、被传输到办公室(经由计算机网络、物理网络或其组合存储和传输)、以及被用作计算设备的输入以生成速度模型和密度模型并生成地下结构的图像。接收的信号的数据可以被称为地震数据或测量数据。
22.在一些情况下,可以收集包括采样的接收信号的地震数据以构建公共图像集(cig)或公共反射点(crp)。在地球物理学分析中,cig是指具有固定表面位置的整个图像的子集,该固定表面位置用于测量在固定图像点的局部图像之间的变化。如果部分图像是反
射角的函数,则相应的cig被称为公共角集(cag)或角域公共图像集(adcig)。
23.cag可以用于通过角度分析、迁移速度质量控制和迁移速度分析来确定幅度变化。将地震数据映射到包含地下中的反射边界的图像的过程被称为迁移。基于地下区域的速度模型和密度模型开发迁移过程。当迁移速度正确时,来自不同入射角的图像聚焦在相同的深度上,从而在cag上产生平坦的事件。相反,当速度有误差时,cag中的事件变得非平坦。除了速度模型,密度模型也可以用于迁移过程。使用正确的速度模型和密度模型,可以基于收集的地震数据生成子结构的图像。
24.地下区域的解空间可以分为多层。层之间的界面可以通过与传播角无关的标量反射率来表征。可以在层之间定义速度。图1示出了根据实施方式的示例波场传播模型100。尽管模型100被示为2维(2d),但是模型100可以通过在3d向量中对声音信号传播方向进行建模来表示3维(3d)地下结构。在图1中,ξ表示层之间的界面。然而,图1只示出了三个子结构界面:对于接口110(ξ=ξ
m

i
)、接口120(ξ=ξ
m
)和接口130(ξ=ξ
m+i
),可以以类似的方式对附加的子结构界面进行建模。
25.对于界面120(ξ=ξ
m
),有六个不同的波场需要考虑:前向入射波场p
+

m
)、后向入射波场p


m
)、前向出射波场q
+

m
)、后向出射波场q


m
)、前向源波场s
+

m
)和后向源波场s


m
)。前向出射波场q
+

m
)和后向出射波场q


m
)可以通过使用双向波动方程和波场分离技术使用波场传播模拟来产生。波场分离技术的示例包括双希尔伯特变换和f

k滤波。这些波场可以被表示为具有相同大小的向量。在傅立叶域中,对于给定频率,这六个波场通过波场传播、传输和反射效应连接。以下是双向波动方程:
[0026][0027][0028]
p
+

m
)==w(ξ
m
,ξ
m
‑1)*q
+

m
‑1)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(3)
[0029]
p


m
)=w(ξ
m
,ξ
m+1
)*q


m+1
)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(4)
[0030]
其中,表示前向方向上界面120处的角度相关透射矩阵、表示前向方向上界面120处的角度相关反射矩阵、表示后向方向上界面120处的角度相关透射矩阵、表示后向方向上界面120处的角度相关反射矩阵、w(ξ
m
,ξ
m
‑1)表示由界面120和110之间的单向波动方程构造的前向传播矩阵、以及(ξ
m
,ξ
m+1
)表示由界面120和130之间的单向波动方程构造的前向传播矩阵。注意,矩阵t
±
、r
±
和w是方阵,其列大小与方程(1)至(4)中其他向量的大小相同。为了减少由角度相关的透射矩阵t
±
和反射矩阵r
±
引入的解空间,可以使用角度无关的系数将方程(1)和(2)改写为以下方程:
[0031]
q
+

m
)=s
+

m
)+t
+

m
).*p
+

m
)+r
+

m
).*p


m
)
ꢀꢀ
(5)
[0032]
q


m
)=s


m
)+t


m
).*p


m
)+r


m
).*p
+

m
)
ꢀꢀ
(6)
[0033]
其中,t
+

m
)和r
+

m
)分别是前向方向上界面120处的角度无关的透射和反射系数的向量,t


m
)和r


m
)分别是后向方向上界面120处的角度无关的透射和反射系数的向量,以及符号“.*”表示两个向量的逐元素乘法。在方程(5)和(6)中,向量t
±
和r
±
的大小与方程(1)至方程(4)中的其他向量相同。
[0034]
t
+

m
)、t


m
)和r
+

m
)与r


m
)相关,如以下方程所示:
[0035]
r
+

m
)=

r


m
)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(7)
[0036]
t
+

m
)=i+r


m
)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(8)
[0037]
t


m
)=i

r


m
)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(9)
[0038]
基于方程(5)至方程(9),的梯度(表示界面120(ξ=ξ
m
)上第i个位置处的前向反射率)可以用以下方程表示:
[0039][0040]
其中,表示的梯度,re{a}表示取变量a的实部,(a)
h
表示取变量a的厄米共轭,和分别表示向量p


m
)和p
+

m
)中的第i个元素。
[0041]
v
m,i
的梯度,表示层m中第i个位置的速度,可以用以下方程表示:
[0042][0043]
其中,g(v
m,i
)表示v
m,i
的梯度。
[0044]
使用局部一维波场传播假设,方程(11)可以被你简化为:
[0045][0046]
其中,和分别表示向量q


m
)和q
+

m
)中的第i个元素,和分别表示向量λ
+

m
)和λ


m
)中的第i个元素。λ
+

m
)和λ


m
)分别表示界面120(ξ=ξ
m
)处的前向出射波场和后向出射波场。a
+

m
)和a


m
)可以通过时间反向传播和地表数据残差的前向波场和后向波场分离来产生。
[0047]
如前所述,速度和反射率的梯度可以基于方程(10)和(12)计算。此外,可以基于速度模型和密度模型对解空间进行参数化。入射反射率r可以用以下方程表示:
[0048][0049]
其中ρ
m
和v
m
分别表示界面120(ξ=ξ
m
)处的密度和速度,而ρ
m+1
和v
m+1
分别表示界面130(ξ=ξ
m+1
)处的密度和速度。因此,如果ρ
m
和v
m
是固定的,反射率的梯度与ρ
m+1
和v
m+1
变化的趋势相匹配。
[0050]
前面讨论的速度梯度和密度梯度是从沿多个方向传播的波场获得的。例如,在具有x、y和z空间方向的3d空间中,可以分别在x、y和z方向执行方程(5)至(9)中的波场分离。因此,可以获得每个方向的速度梯度和反射率梯度。可以分别使用每个方向的速度梯度和反射率梯度来更新速度模型和密度模型。图2示出了根据实施方式的基于在一个方向上的
速度梯度和反射率梯度来更新速度模型和密度模型的示例过程200。
[0051]
在202,使用方程(12)计算速度梯度。然后速度梯度在范围[

1,1]内被归一化。归一化的速度梯度被表示为ngv。
[0052]
在204,使用方程(10)计算反射率梯度。反射率梯度也在范围[

1,1]内被归一化。归一化的反射率梯度被表示为ngr。
[0053]
在210,归一化的速度梯度(被表示为ngv)用于更新界面120(ξ=ξ
m
)上方的地下层的速度模型,例如界面110(ξ=ξ
m
‑1)和界面120(ξ=ξ
m
)之间的地下层。在所示示例中,使用了获得步长搜索过程。在212,初始速度更新步长dv被设置为dv0。在214,确定速度更新步长是否可以减少表面数据残差。基于当前速度模型v
current
生成新的速度模型v
updated
,通过使用以下方程:
[0054]
v
updated
=v
current
+dv*ngv
ꢀꢀꢀꢀꢀꢀ
(14)
[0055]
表面数据残差(被表示为ds)表示测量数据与模拟数据之间的差异。测量数据(被表示为d)表示在接收器设备处接收的信号。基于一组给定的模型获得模拟数据(被表示为dq),该模型包括速度模型和密度模型。因此,ds=d

dq。将使用新速度模型v
updated
计算的表面数据残差与使用当前速度模型v
current
计算的表面数据残差进行比较。如果基于新速度模型获得的表面数据残差等于或大于基于先前速度模型获得的表面数据残差,则搜索过程进行到216,其中使用线搜索方案调整dv。搜索过程继续进行到214,其中调整的dv用于确定基于调整的dv获得的新速度模型是否会减少表面数据残差。搜索过程一直持续,直到找到一个dv值以使新的速度模型减少表面数据残差。在218,使用找到的dv值获得更新的速度模型。
[0056]
在220,归一化的放射率梯度(被表示为ngr)用于更新界面120(ξ=ξ
m
)上方的地下层的密度模型,例如界面110(ξ=ξ
m

i
)和界面120(ξ=ξ
m
)之间的地下层。在所示示例中,使用了自适应步长搜索过程。在222,初始反射率更新步长dρ被设置为dρ0。在224,确定反射率更新步长是否可以减少表面数据残差。基于当前密度模型ρ
current
生成新的密度模型ρ
updated
,通过使用以下方程:
[0057]
ρ
updated
=ρ
current
+dρ*ngr
ꢀꢀꢀꢀꢀꢀ
(15)
[0058]
将使用新反射率模型p
updated
计算的表面数据残差与使用当前反射率模型p
current
计算的表面数据残差进行比较。如果基于新反射率模型获得的表面数据残差等于或大于基于先前反射率模型获得的表面数据残差,则搜索过程进行到226,其中使用线搜索方案调整dρ。搜索过程继续进行到224,其中调整的dρ用于确定基于调整的dρ获得的新反射率模型是否会减少表面数据残差。搜索过程一直持续,直到找到一个dρ值以使新的反射率模型减少表面数据残差。在228,使用找到的dρ值获得更新的反射率模型。
[0059]
在230,归一化的速度梯度用于更新速度模型以用于反演的下一次迭代。可以对界面120(ξ=ξ
m
)下方的地下层(例如界面120(ξ=ξ
m
)和界面130(ξ=ξ
m+1
)之间的地下层)执行反演的下一次迭代。这里可以使用与步骤210中使用的相同的自适应步长进搜索算法。
[0060]
在240,归一化的反射率梯度用于更新密度模型以用于反演的下一次迭代。可以对界面120(ξ=ξ
m
)下方的地下层(例如界面120(ξ=ξ
m
)和界面130(ξ=ξ
m+i
)之间的地下层)执行反演的下一次迭代。这里可以使用与步骤220中使用的相同的自适应步长进搜索算法。
[0061]
可以基于其他方向的速度梯度和反射率梯度重复过程200中的步骤。过程200中的
步骤可以重复多次迭代以更新解空间。
[0062]
速度模型和密度模型可以使用可调频带的迭代过程进行更新。测量数据的低频部分最初用于更新速度模型和密度模型。使用迭代过程以包括测量数据的附加频率分量,直到使用了测量数据的全部带宽。图3示出了根据实施方式的基于迭代频带生成速度模型和密度模型的示例过程300。
[0063]
在302,接收地震信号的测量数据。测量数据也可以被称为地震数据。在304,设置初始频率f。在一些实施方式中,可以选择小于5赫兹(hz)的频率作为初始频率f。例如,初始频率f可以被设置为2hz。在306,确定测量数据中的最大频率。最大频率可以通过执行傅立叶变换将时域数据变换到频域数据来确定。通过基于频域数据选择能够满足配置的最小信噪比标准的最大频率来确定最大频率。将测量数据中的最大频率与初始频率f进行比较。如果最大频率大于或等于f,则过程300继续进行到310,其中使用低通滤波器以f作为上限频率对测量数据进行滤波。因此,滤波后的测量数据中的频率低于或等于f。
[0064]
在312,当前速度模型和密度模型被设置为初始速度模型和密度模型。在一些情况下,可以基于配置的默认模型选择初始速度模型和密度模型,例如线性梯度模型或均匀背景模型。备选地或附加地,可以基于附加的先验信息来选择或修改初始速度模型和密度模型。在314,双向波动方程与当前速度模型和密度模型一起用于模拟地下区域中的波场。在316,选择一个方向用于波场分离。在具有x、y和z空间方向的3d空间中,所选方向可以是x、y或z方向中的一个。在318,所选方向应用于先前讨论的波分离技术以获得所选方向的速度和反射率梯度。
[0065]
在320,基于先前讨论的过程200使用所选方向的速度梯度和反射率梯度来更新速度模型和密度模型。在322,当前速度模型和密度模型被设置为更新的模型。在330,确定是否存在其他未选择的空间方向。如果还有其他未选择的方向,则过程300返回到步骤312。重复步骤312至330,直到处理完所有三个空间方向。在处理完所有三个空间方向之后,过程300进行到332,其中通过将初始频率值增加df来更新频率f。在一些情况下,df值可以由用户或系统管理员配置。过程300返回到306,其中是否比较更新的频率f和测量数据中的最大频率。过程300继续进行直到更新的频率f大于测量数据的最大频率。此时,过程300进行到308,其中速度模型和密度模型被输出为解模型。
[0066]
根据本公开中讨论的过程获得的速度模型和密度模型提供了地下结构的准确特征。因此,它们可以用于提高地下结构的成像分辨率。在地震勘探中,这些质量提高的地下图像可以提供对地下结构更好的了解。因此,根据本公开中讨论的过程获得的速度模型和密度模型可以用于指导油气勘探或生产操作,并在生产过程中提高钻井的成功率和减少地质灾害。
[0067]
图4示出了根据实施方式的示例速度模型和密度模型生成过程400。为了清楚地呈现,下文的描述在图1至图3和图5a至图9的上下文下总体地描述处理400。然而,应该理解的是,处理400可以例如通过任何其他合适的系统、环境、软件和硬件、或者系统、环境、软件和硬件的适当组合来执行。在一些情况下,处理400可以在计算机集群、超级计算机或任何其他计算设备或计算设备的集合上执行。在一些实现中,处理400的各个步骤可以并行、组合、循环或以任意顺序运行。
[0068]
在402,数据处理装置接收所述油气储层的一组地震数据。在404,数据处理装置设
置初始速度模型和初始密度模型。在406,数据处理装置基于所述一组地震数据产生油气储层的波场。在408,数据处理装置选择空间方向。在410,数据处理装置基于波场产生选择的空间方向的速度梯度和反射率梯度。在412,数据处理装置使用选择的空间方向的速度梯度和反射率梯度更新速度模型和密度模型。
[0069]
使用基于二维(2d)合成示例的模拟来说明本公开中讨论的建模过程的效果。在模拟中,震源深度为5m,并且主频为10赫兹(hz)的ricker小波用作震源小波。在横向上,震源位置从100m开始到2300m结束,间距为50m,这导致45个炮点道集。所有炮点道集的接收器几何形状相同:接收器深度为10m,接收器横向覆盖范围[0,2400m],间距为25m。反演计划中的最大反演频率被设置为30hz。在过程300中,初始频率f被设置为f=2hz,并且最大频率是30hz。因为这是一个二维样本,所以在反演过程中选择两个方向(z方向和x方向)来更新模型。图5a和图5b示出了根据实施方式的合成示例的速度模型510和密度模型520。图6a和图6b示出了根据实施方式的初始速度模型610和初始密度模型620。图7a和图7b示出了根据实施方式的基于本公开中讨论的过程获得的反演速度模型710和反演密度模型720。
[0070]
图8示出了根据实施方式的炮点道集的示例图表800。该图表包括部分810中的炮点道集(表示测量数据)和部分820中的炮点道集(表示基于获得的速度模型和密度模型的模拟数据)。获得的速度模型和密度模型提供了准确的结果。因此,部分830中的炮点道集(表示数据残差)是弱的。
[0071]
图9是示出了根据一些实施方式的用于提供与如本公开中描述的算法、方法、功能、处理、流程和过程相关联的计算功能的示例计算机系统900的框图。所示系统900包括与网络930耦接的计算机902。
[0072]
所描述的图示仅是所描述的主题的一种可能的实现,并且不旨在将本公开限制于该单个描述的实现。本领域的普通技术人员将认识到,所描述的部件可以以与本公开一致的替代方式连接、组合或使用。
[0073]
网络930促进计算机902与其他组件(例如,获得针对位置的观测数据并将观测数据发送到计算机902的组件)之间的通信。网络930可以是无线或有线网络。网络930还可以是存储器管道、硬件连接或组件之间的任何内部或外部通信路径。
[0074]
计算机902包括被配置为执行本公开中所描述的方法的计算系统。例如,计算机902可以用于实现图2

4所示的过程。在一些情况下,所述方法的算法可以用可执行计算代码(例如,c/c++可执行代码)来实现。在一些情况下,计算机902可以包括运行批量应用的独立的linux系统。在一些情况下,计算机902可以包括移动或个人计算机。
[0075]
计算机902包含任何计算设备,例如服务器、台式计算机、膝上型/笔记本计算机、个人数据助理(pda)、平板计算设备或这些设备内的一个或多个处理器,该计算设备包括物理实例、虚拟实例或两者。计算机902可以包括输入设备,例如小键盘、键盘、触摸屏、麦克风、语音识别设备和可以接受用户信息的其他设备。计算机902可以包括可传达与计算机902的操作相关联的信息的输出设备。该信息可以包括数字数据、视觉数据或音频信息,或信息的组合。该信息可以在图形用户界面(ui)(或gui)中呈现。
[0076]
计算机902可以用作用于执行本公开中描述的主题的计算机系统的客户端、网络组件、服务器、数据库、持久性或组件。所示出的计算机902可通信地与网络930耦接。在一些实现中,计算机902的一个或多个组件可以被配置为在包括基于云计算的环境、本地环境、
全局环境、或者环境的组合在内的不同环境中操作。
[0077]
从高层面来看,计算机902是可操作用于接收、发送、处理、存储和管理与所描述的主题相关联的数据和信息的电子计算设备。根据一些实施方式,计算机902还可以包括或可通信地耦接到应用服务器、电子邮件服务器、web服务器、缓存服务器、流数据服务器、或服务器的组合。
[0078]
计算机902可以通过网络930从客户端应用程序(例如,在另一台计算机902上执行)接收请求。计算机902可以通过使用软件应用程序处理接收的请求来响应接收的请求。还可以从内部用户(例如,从命令控制台)、外部(或第三方)、自动化应用程序、实体、个人、系统和计算机向计算机902发送请求。
[0079]
计算机902的组件中的每个组件可以使用系统总线903通信。在一些实施方式中,计算机902的任何或所有组件(包括硬件或软件组件),可以通过系统总线903彼此连接或与接口904(或两者的组合)连接。接口可以使用应用程序编程接口(api)912、服务层913或api912和服务层913的组合。api912可以包括针对例程、数据结构和对象类的规范。api912可以独立于或依赖于计算机语言。api912可以指完整的接口、单个功能或一组api。
[0080]
服务层913可以向计算机902和可通信地耦接到计算机902的其他组件(无论是否被示出)提供软件服务。计算机902的功能对于使用该服务层的所有服务消费者可以是可访问的。软件服务(例如由服务层913提供的软件服务)可以通过定义的接口提供可重用的、定义的业务功能。例如,接口可以是以java、c++或以可扩展标记语言(xml)格式提供数据的语言所编写的软件。尽管被示出为计算机902的集成组件,但是在备选实施方式中,可以将api912或服务层913示出为相对于计算机902的其他组件或可通信地耦接到计算机902的其他组件的独立组件。此外,在不脱离本公开的范围的情况下,api912和/或服务层913的任意或所有部分可以被实现为另一软件模块、企业应用或硬件模块的子模块或副模块。
[0081]
计算机902包括接口904。尽管在图9中被示为单个接口904,但是可以根据计算机902的特定需要、期望或特定实现和所描述的功能而使用两个或更多个接口904。接口904由计算机902可以用于与连接到网络930的分布式环境中的其它系统(无论是否示出)通信。通常,接口904可以包括或被实现为使用以软件或硬件(或软件和硬件的组合)编码的逻辑,该编码可操作用于与网络930通信。更具体地,接口904可以包括支持与通信相关联的一种或多种通信协议的软件。因此,网络930或接口的硬件可以用于在所示计算机902之内和之外传送物理信号。
[0082]
计算机902包括处理器905。尽管在图9中被示为单个处理器905,但是可以根据计算机902的特定需要、期望或特定实现和所描述的功能而使用两个或更多个处理器905。一般地,处理器905可以执行指令并操纵数据以执行计算机902的操作,该操作包括使用如本公开中所描述的任何算法、方法、功能、处理、流程和过程。
[0083]
计算机902进一步包括存储器906,其可以保存用于计算机902或可以连接到网络930的组件的组合(无论是否被示出)的数据。存储器906可以存储与本公开一致的任意数据。在一些实现中,根据计算机902的特定需要、期望或特定实现和所描述的功能,存储器906可以是两个或更多个不同类型的存储器的组合(例如,半导体和磁存储器的组合)。尽管在图9中被示出为单个存储器906,根据计算机902的特定需要、期望或特定实现和所描述的功能,可以使用相似或不同类型的两个或更多个存储器906(相同类型、不同类型或类型组
合)。尽管存储器906被示出为计算机902的集成组件,但是在备选实施方式中,存储器906可以在计算机902的外部。
[0084]
应用程序908可以是根据计算机902和所描述的功能性的特定需要、期望或特定实施方式提供功能性的算法软件引擎。例如,应用908可以用作一个或多个组件、模块、应用等。此外,尽管被示为单个应用908,但是应用908可以被实现为计算机902上的多个应用908。另外,尽管被示出为与计算机902集成在一起,但是在备选实施方式中,应用908可以在计算机902的外部。
[0085]
计算机902还可以包括电源914。电源914可以包括可以被配置为用户或非用户可更换的可再充电或不可再充电电池。在一些实现中,电源914可以包括电力转换和管理电路(包括再充电、备用和电力管理功能)。在一些实现中,电源914可以包括电源插头,以允许计算机902插入墙上插座或电源以例如为计算机902供电或为可再充电电池充电。
[0086]
可以存在与包含计算机902的计算机系统相关联或在其外部的任意数量的计算机902,每个计算机902通过网络930进行通信。此外,在不脱离本公开的范围的情况下,术语“客户端”、“用户”和其他适当的术语可以适当地互换使用。此外,本公开包含许多用户可以使用一个计算机902和一个用户可以使用多个计算机902。
[0087]
所描述的主题的实施方式可以单独或组合地包括一个或多个特征。
[0088]
例如,在第一实施方式中,用于生成速度模型和密度模型以获得油气储层中地下结构的图像的计算机实现方法包括:由数据处理装置接收油气储层的一组地震数据;由数据处理装置设置初始速度模型和初始密度模型;由所述数据处理装置基于所述一组地震数据生成所述油气储层的波场;由数据处理装置选择空间方向;由数据处理装置基于波场产生所选择的空间方向的速度梯度和反射率梯度;以及由数据处理装置使用选择的空间方向的速度梯度和反射率梯度更新速度模型和密度模型。
[0089]
前述和其他所述实施方式可以各自可选地以单独或组合的方式包括以下特征中的一个或多个:
[0090]
第一特征,可与一般实施方式和以下任何特征结合,其中操作进一步包括:选择第二空间方向;由所述数据处理装置产生所述第二空间方向的第二速度梯度和第二反射率梯度;以及使用所述第二速度梯度和所述第二反射率梯度来更新所述速度模型和所述密度模型。
[0091]
第二特征,可与前面或以下的任何特征结合,其中更新速度模型和密度模型包括:使用第一迭代过程来更新所述速度模型,其中,所述第一迭代过程包括第一多次迭代,所述第一多次迭代中的每一次包括生成速度更新步长值、基于所述速度更新步长值和所述速度梯度生成更新的速度模型,并且在所述第一多次迭代中的每一次中,所述初始速度模型被在先前迭代中生成的所述更新的速度模型替换;以及使用第二迭代过程来更新所述密度模型,其中,所述第二迭代过程包括第二多次迭代,所述第二多次迭代中的每一次包括生成密度更新步长值、基于所述密度更新步长值和所述反射率梯度生成更新的密度模型,并且在所述第二多次迭代中的每一次中,所述初始密度模型被在先前迭代中生成的所述更新的密度模型替换;以及
[0092]
第三特征,可与前面或以下特征中的任一个结合,其中,当所述更新速度模型没有减少第一表面数据残差时所述第一迭代过程终止,并且当所述更新密度模型没有减少第二
表面数据残差时所述第二迭代过程终止。
[0093]
第四特征,可与前面或以下特征中的任一个结合,其中该方法进一步包括:使用基于初始频率的低通滤波器对所述一组地震数据进行滤波;以及其中,所述油气储层的所述波场通过使用所述滤波的一组地震数据生成。
[0094]
第五特征,可与前述或以下特征中的任一个结合,其中该方法进一步包括:使用迭代过程更新所述速度模型和所述密度模型,其中,所述迭代过程包括多次迭代,所述多次迭代中的每一次包括生成频率步长值、基于所述频率步长值生成更新的频率、以及基于所述更新的频率对所述一组地震数据进行滤波以生成更新的速度梯度和更新的反射率梯度,以及使用所述更新的速度梯度和所述更新的反射率梯度更新所述速度模型和所述密度模型。
[0095]
第六特征,可与前面或以下的特征中的任一个组合,其中,当更新频率大于一组地震数据中的最大频率时,迭代过程终止。
[0096]
在第二实施方式中,一种存储指令的非暂时性计算机可读介质,所述指令在被执行时使计算机系统执行操作,所述操作包括:由数据处理装置接收油气储层的一组地震数据;由数据处理装置设置初始速度模型和初始密度模型;由所述数据处理装置基于所述一组地震数据生成所述油气储层的波场;由数据处理装置选择空间方向;由数据处理装置基于波场产生所选择的空间方向的速度梯度和反射率梯度;以及由数据处理装置使用选择的空间方向的速度梯度和反射率梯度更新速度模型和密度模型。
[0097]
前述和其他所述实施方式可以各自可选地以单独或组合的方式包括以下特征中的一个或多个:
[0098]
第一特征,可与一般实施方式和以下任何特征结合,其中操作进一步包括:选择第二空间方向;由所述数据处理装置产生所述第二空间方向的第二速度梯度和第二反射率梯度;以及使用所述第二速度梯度和所述第二反射率梯度来更新所述速度模型和所述密度模型。
[0099]
第二特征,可与前面或以下的任何特征结合,其中更新速度模型和密度模型包括:使用第一迭代过程来更新所述速度模型,其中,所述第一迭代过程包括第一多次迭代,所述第一多次迭代中的每一次包括生成速度更新步长值、基于所述速度更新步长值和所述速度梯度生成更新的速度模型,并且在所述第一多次迭代中的每一次中,所述初始速度模型被在先前迭代中生成的所述更新的速度模型替换;以及使用第二迭代过程来更新所述密度模型,其中,所述第二迭代过程包括第二多次迭代,所述第二多次迭代中的每一次包括生成密度更新步长值、基于所述密度更新步长值和所述反射率梯度生成更新的密度模型,并且在所述第二多次迭代中的每一次中,所述初始密度模型被在先前迭代中生成的所述更新的密度模型替换;以及
[0100]
第三特征,可与前面或以下特征中的任一个结合,其中,当所述更新速度模型没有减少第一表面数据残差时所述第一迭代过程终止,并且当所述更新密度模型没有减少第二表面数据残差时所述第二迭代过程终止。
[0101]
第四特征,可与前面或以下特征中的任一个结合,其中所述操作进一步包括:使用基于初始频率的低通滤波器对所述一组地震数据进行滤波;以及其中,所述油气储层的所述波场通过使用所述滤波的一组地震数据生成。
[0102]
第五特征,可与前面或以下特征中的任一个结合,其中所述操作进一步包括:使用
迭代过程更新所述速度模型和所述密度模型,其中,所述迭代过程包括多次迭代,所述多次迭代中的每一次包括生成频率步长值、基于所述频率步长值生成更新的频率、以及基于所述更新的频率对所述一组地震数据进行滤波以生成更新的速度梯度和更新的反射率梯度,以及使用所述更新的速度梯度和所述更新的反射率梯度更新所述速度模型和所述密度模型。
[0103]
第六特征,可与前面或以下的特征中的任一个组合,其中,当更新频率大于一组地震数据中的最大频率时,迭代过程终止。
[0104]
在第三实施方式中,一种设备包括:至少一个硬件处理器;以及非暂时性计算机可读存储介质,耦接到至少一个硬件处理器并存储用于由至少一个硬件处理器执行的编程指令,其中编程指令在被执行时,使至少一个硬件处理器执行包括以下各项的操作:接收油气储层的一组地震数据;设置初始速度模型和初始密度模型;基于所述一组地震数据生成所述油气储层的波场;选择空间方向;基于波场产生所选择的空间方向的速度梯度和反射率梯度;以及使用所选择的空间方向的所述速度梯度和所述反射率梯度来更新速度模型和密度模型。
[0105]
前述和其他所述实施方式可以各自可选地以单独或组合的方式包括以下特征中的一个或多个:
[0106]
第一特征,可与一般实施方式和以下任何特征结合,其中操作进一步包括:选择第二空间方向;产生所述第二空间方向的第二速度梯度和第二反射率梯度;以及使用所述第二速度梯度和所述第二反射率梯度来更新所述速度模型和所述密度模型。
[0107]
第二特征,可与前面或以下的任何特征结合,其中更新速度模型和密度模型包括:使用第一迭代过程来更新所述速度模型,其中,所述第一迭代过程包括第一多次迭代,所述第一多次迭代中的每一次包括生成速度更新步长值、基于所述速度更新步长值和所述速度梯度生成更新的速度模型,并且在所述第一多次迭代中的每一次中,所述初始速度模型被在先前迭代中生成的所述更新的速度模型替换;以及使用第二迭代过程来更新所述密度模型,其中,所述第二迭代过程包括第二多次迭代,所述第二多次迭代中的每一次包括生成密度更新步长值、基于所述密度更新步长值和所述反射率梯度生成更新的密度模型,并且在所述第二多次迭代中的每一次中,所述初始密度模型被在先前迭代中生成的所述更新的密度模型替换;以及
[0108]
第三特征,可与前面或以下特征中的任一个结合,其中,当所述更新速度模型没有减少第一表面数据残差时所述第一迭代过程终止,并且当所述更新密度模型没有减少第二表面数据残差时所述第二迭代过程终止。
[0109]
第四特征,可与前面或以下特征中的任一个结合,其中所述操作进一步包括:使用基于初始频率的低通滤波器对所述一组地震数据进行滤波;以及其中,所述油气储层的所述波场通过使用所述滤波的一组地震数据生成。
[0110]
第五特征,可与前面或以下特征中的任一个结合,其中所述操作进一步包括:使用迭代过程更新所述速度模型和所述密度模型,其中,所述迭代过程包括多次迭代,所述多次迭代中的每一次包括生成频率步长值、基于所述频率步长值生成更新的频率、以及基于所述更新的频率对所述一组地震数据进行滤波以生成更新的速度梯度和更新的反射率梯度,以及使用所述更新的速度梯度和所述更新的反射率梯度更新所述速度模型和所述密度模
型。
[0111]
第六特征,可与前面或以下的特征中的任一个组合,其中,当更新频率大于一组地震数据中的最大频率时,迭代过程终止。
[0112]
在一些实施方式中,所描述的方法可被配置为向计算机实现的控制器、数据库或其他计算机实现的系统发送消息、指令或其他通信,以动态地启动对另一计算机实现的系统的控制,控制或使得另一计算机实现的系统执行计算机实现的操作。例如,可以传输基于数据的操作、操作、输出或与gui的交互,以使与计算机、数据库、网络或其他基于计算机的系统相关联的操作执行存储效率、数据检索或与本公开一致的其他操作。在另一示例中,与任何所示的gui进行交互可以自动地导致从gui传输的一个或多个指令触发对数据的请求、数据的存储、数据的分析或与本公开一致的其他操作。
[0113]
在一些情况下,传输的指令可以导致关于有形的真实世界的计算设备或其他设备的控制、操作、修改、增强或其他操作。例如,所描述的gui可以发送请求以减慢或加速计算机数据库磁/光盘驱动器,激活/停用计算系统,使网络接口设备禁用、节制或增加通过网络连接允许的数据带宽,或发出可听/可视警报(例如,机械警报/发光设备)作为对于与所描述的方法相关联的计算系统的结果、行为、确定或分析的通知,或者作为与所描述的方法相关联的计算系统的交互。
[0114]
在一些实施方式中,所描述的方法的输出可以用于动态地影响、指导、控制、影响或管理与烃生产、分析和回收有关的或用于与本公开一致的其它目的有形的现实世界的设备。例如,与处理的地震数据有关的数据可以用于提高生成的地震/结构图像的质量或用于其他分析/预测过程。作为另一个示例,与处理的地震数据有关的数据可以用于修改井眼轨迹、增加/降低油气钻探的速度或停止/启动油气钻探;激活/停用警报(例如,视觉、听觉或语音警报),或影响炼油厂或泵送操作(例如,停止、重启、加速或减少)。其他示例可以包括在检测到地下障碍物时警告(例如,用视觉、听觉或语音警报)地质导向和定向钻井人员。在一些实施方式中,所描述的方法可以被集成为动态的计算机实现的控制系统的一部分,从而用于控制、影响与本公开一致的与烃相关的任何设备或其他有形的现实世界的设备,或与与本公开一致的与烃相关的任何设备或其他有形的现实世界的设备一同使用。
[0115]
在本说明书中描述的主题和功能操作的实施方式可以在数字电子电路中、在有形地被实现的计算机软件或固件中、在计算机硬件中实施,包括在本说明书中公开的结构及其结构等同物、或它们中的一个或多个的组合中实施。所描述的主题的软件实现可以被实现为一个或多个计算机程序。每个计算机程序可以包括一个或多个计算机程序指令的模块,该模块被编码在有形的、非暂时的、计算机可读的计算机存储介质上,用于由数据处理装置执行或控制数据处理装置的操作。备选地或附加地,程序指令可以被编码在人工生成的传播信号中/之上。例如,该信号可以是机器生成的电、光或电磁信号,其被生成以对信息进行编码以传输到合适的接收器装备以供数据处理装置执行。计算机存储介质可以是机器可读存储设备、机器可读存储基板、随机或串行存取存储器设备、或计算机存储介质的组合。
[0116]
术语“数据处理装置”、“计算机”和“电子计算机设备”(或本领域普通技术人员所理解的等同物)是指数据处理硬件。例如,数据处理装置可以包括用于处理数据的所有种类的装置、设备和机器,示例性地包括可编程处理器、计算机、或多个处理器或计算机。该装置
还可以包括专用逻辑电路,该专用逻辑电路包括例如中央处理单元(cpu)、现场可编程门阵列(fpga)或专用集成电路(asic)。在一些实施方式中,数据处理装置或专用逻辑电路(或数据处理装置或专用逻辑电路的组合)可以基于硬件或基于软件(或基于硬件和基于软件的组合)。可选地,装置可以包括为计算机程序创建执行环境的代码,例如,构成处理器固件、协议栈、数据库管理系统、操作系统或者执行环境的组合的代码。本公开考虑具有或不具有常规操作系统(例如linux、unix、windows、macos、android、或ios)的数据处理装置的使用。
[0117]
也可以被称为或描述为程序、软件、软件应用程序、模块、软件模块、脚本或代码的计算机程序可以以任何形式的编程语言编写。编程语言可以包括例如编译语言、解释语言、声明性语言或过程语言。程序可以以任何形式部署,包括作为独立程序、模块、组件、子例程或在计算环境中使用的单元。计算机程序可以(但不是必须)与文件系统中的文件相对应。程序可以被存储在保存其他程序或数据(例如,存储在标记语言文档中的一个或多个脚本)的文件的一部分中、被存储在专用于所讨论的程序的单个文件中、或者被存储在多个协同文件中,该协同文件存储一个或多个模块、子程序或代码部分的文件。计算机程序可以被部署为在一个计算机上或者在位于例如一个站点或分布在通过通信网络互联的多个站点的多个计算机上执行。尽管各图中所出示的程序的部分可以被示出为通过各种对象、方法或过程实施各种特征和功能的各个模块,但是在适当时程序可以替代地包括多个子模块、第三方服务、组件、库等。相反,各种组件的特征和功能在适当时可以被组合成单个组件。可以统计地、动态地或者统计地且动态地确定用于进行计算确定的阈值。
[0118]
本说明书中描述的方法、过程或逻辑流可以由一个或多个可编程计算机来执行,所述一个或多个可编程计算机执行一个或多个计算机程序以通过操作输入数据并且生成输出来执行功能。方法、处理或逻辑流也可以由专用逻辑电路(例如cpu、fpga或asic)来执行,并且装置也可以实现为专用逻辑电路(例如cpu、fpga或asic)。
[0119]
适合执行计算机程序的计算机可以基于一个或多个通用和专用微处理器以及其他种类的cpu。计算机的的元件是用于执行指令的cpu和用于存储指令和数据的一个或更多个存储器设备。通常,cpu可以从存储器接收指令和数据并写入存储器。计算机还可以包括用于存储数据的一个或多个大容量存储设备,或者可操作地耦接到用于存储数据一个或多个大容量存储设备;在一些实施方式中,计算机可以从大容量存储设备接收数据和将数据传输到大容量存储设备,大容量存储设备包括例如磁盘、磁光盘或光盘。此外,计算机可以嵌入在另一设备中,例如,移动电话、个人数字助理(pda)、移动音频或视频播放器、游戏机、全球定位系统(gps)接收器或者便携式存储设备(例如,通用串行总线(usb)闪存驱动器)。
[0120]
适合存储计算机程序指令和数据的计算机可读介质(适当的暂时性或非暂时性的)可以包括所有形式的永久/非永久和易失性/非易失性存储器器、介质和存储器设备。计算机可读介质可以包括例如半导体存储器设备,例如随机存取存储器(ram)、只读存储器(rom)、相变存储器(pram)、静态随机存取存储器(sram)、动态随机存取存储器(dram)、可擦可编程只读存储器(eprom)、电可擦可编程只读存储器(eeprom)和闪存器件。计算机可读介质还可以包括例如磁性设备,例如磁带、盒式磁带、盒式磁带和内部/可移动磁盘。计算机可读介质还可以包括磁光盘和光学存储器设备和技术,其包括例如数字视频盘(dvd)、cdrom、dvd+/

r、dvd

ram、dvd

rom、hd

dvd和bluray。存储器可以存储各种对象或数据,其包括缓存、类、框架、应用程序、模块、备份数据、作业、网页、网页模板、数据结构、数据库表、存储库
和动态信息。存储器中存储的对象和数据的类型可以包括参数、变量、算法、指令、规则、约束和引用。此外,存储器可以包括日志、策略、安全或访问数据以及报告文件。处理器和存储器可以由专用逻辑电路来补充或者并入到专用逻辑电路中。
[0121]
本公开中描述的主题的实施方式可以在具有显示设备的计算机上实现,该显示设备用于提供与用户的交互,该交互包括向用户显示信息(和从用户接收输入)。显示设备的类型可以包括例如阴极射线管(crt)、液晶显示器(lcd)、发光二极管(led)和等离子监视器。显示设备可以包括键盘和包括例如鼠标、轨迹球或触控板的定点设备。还可以通过使用触摸屏向计算机提供用户输入,例如具有压力灵敏度的平板计算机表面或使用电容或电感应的多点触摸屏。其他种类的设备可以用于提供与用户的交互,该交互包括接收用户反馈,该反馈包括例如感官反馈,该感官反馈包括视觉反馈、听觉反馈或触觉反馈。可以以声音、语音或触觉输入的形式接收来自用户的输入。此外,计算机可以通过向用户使用的设备发送文档或者从该设备接收文档来与用户交互。例如,计算机可以响应于从用户客户端设备上的web浏览器接收到的请求而向所述web浏览器发送网页。
[0122]
术语“图形用户界面”或“gui”可以以单数或复数形式使用,以描述一个或更多个图形用户界面以及特定图形用户界面的每一次显示。因此,gui可以表示任意图形用户界面,包括但不限于web浏览器、触摸屏或处理信息并且有效地向用户呈现信息结果的命令行界面(cli)。通常,gui可以包括多个ui元素,其中一些或全部与web浏览器相关联,诸如交互式字段、下拉列表和按钮。这些和其他ui元素可以与web浏览器的功能相关或表示web浏览器的功能。
[0123]
本说明书中描述的主题的实施方式可以在包括后端组件(例如,作为数据服务器)或包括中间件组件(例如,应用服务器)的计算系统中实现。此外,计算系统可以包括前端组件,例如,具有图形用户界面或web浏览器中的一者或两者的客户端计算机,通过该客户端计算机用户可以与计算机进行交互。系统的组件可以通过通信网络中有线或无线数字数据通信(或数据通信的组合)的介质或任意形式(例如通信网络)互相连接。通信网络的示例包括局域网(lan)、无线电接入网络(ran)、城域网(man)、广域网(wan)、全球微波接入互操作性(wimax)、无线局域网网络(wlan)(例如,使用802.11a/b/g/n或802.20或协议的组合)、互联网的全部或一部分、或一个或多个位置处的任意其它通信系统(或通信网络的组合)。网络可以在网络地址之间传输例如互联网协议(ip)分组、帧中继帧、异步传输模式(atm)小区、语音、视频、数据或通信类型的组合。
[0124]
计算机系统可以包括客户端和服务器。客户端和服务器一般相互远离并且通常通过通信网络进行交互。客户端和服务器的关系可以通过在相应计算机上运行并且相互具有客户端

服务器关系的计算机程序来产生。
[0125]
集群文件系统可以是可从多个服务器访问以进行读取和更新的任何文件系统类型。由于交换文件系统的锁定可以在应用层完成,锁定或一致性跟踪可能不是必需的。此外,unicode数据文件可以不同于非unicode数据文件。
[0126]
尽管本说明书包含许多特定的实现细节,但是这些细节不应被解释为对可被要求保护的范围的限制,而是对具体示例特有的特征的描述。在单个实施方式中,还可以组合实现本说明书中在独立实施方式的上下文中描述的特定特征。反之,在单个实施方式的上下文中描述的各种特征也可以在多个实施方式中分开地或以任何合适的子组合来实施。此
外,尽管可能将先前描述的特征描述为以某些组合起作用并且甚至最初要求如此保护,但是在一些情况下,来自所要求保护的组合的一个或多个特征可以从组合中删除,并且所要求保护的组合可以指子组合或子组合的变化。
[0127]
已经描述了本主题的特定实施方式。对于本领域技术人员显而易见的是,所描述的实施方式的其他实施方式、改变和置换在下文的权利要求的范围内。尽管在附图或权利要求中以特定顺序描述了操作,但这不应被理解为:为了实现希望的结果,要求按所示出的特定顺序或按相继的顺序来执行这些操作,或者要求执行所有所示出的操作(一些操作可以被认为是可选的)。在一些情况下,多任务或并行处理(或者多任务和并行处理的组合)可以是有利的并且视情况来执行。
[0128]
此外,在前述的实现中的各种系统模块和组件的分离或集成不应被理解为在所有实施方式中要求这样的分离或集成,并且应该理解的是,所描述的程序组件和系统一般可以一起集成在单个软件产品中或封装为多个软件产品。
[0129]
因此,前述示例实施方式不限定或限制本公开。在不脱离本公开的精神和范围的情况下,还可以存在其他改变、替换和变化。
[0130]
此外,任何要求保护的实施方式被认为适用于至少一种计算机实施的方法;存储用于执行计算机实现的方法的计算机可读指令的非暂时性计算机可读介质;以及计算机系统,该系统包括与硬件处理器可互操作地耦接的计算机存储器,所述硬件处理器被配置为执行计算机实现的方法或存储在非暂时性计算机可读介质上的指令。
当前第1页1 2 
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1