一种油藏开发引起的地面沉降量的确定方法及相关组件

1.本发明涉及油藏开发领域，特别是涉及一种油藏开发引起的地面沉降量的确定方法及相关组件。


背景技术：

2.油藏区域是一个处在地下的地质实体，其中包含有石油、天然气、水以及岩石等各种物质，在对油藏区域进行采油时，油藏区域中的孔隙压力会发生变化，而孔隙压力变化会使得油藏区域的应力同步发生变化，应力发生改变又会导致油藏区域所在的地层产生形变，最终可能会引起油藏区域对应的地面区域的沉降。由于地面沉降可能会对地面上的管道线路以及深井等设备造成严重影响，所以在对油藏区域进行油藏开发时，需要准确地预测本次油藏开发所带来的地面沉降量，以避免对设备造成影响。现有技术在确定地面沉降量时，会通过薄圆形油藏地层压力降低时地面垂直位移的计算表达式来计算油藏区域对应的地面沉降量，或者基于格林函数确定油藏区域的应力发生改变时产生的地面沉降量，但是这些表达式或函数均仅能够在油藏区域为某种特定形状时才能预测出准确结果，例如格林函数仅能够在油藏区域的形状为中心轴对称结构的形状时才能准确预测结果，其不能对任一种形状的油藏区域进行预测，实用性较低且应用范围小。


技术实现要素：

3.本发明的目的是提供一种油藏开发引起的地面沉降量的确定方法及相关组件，能够确定任一种形状的油藏区域的体积模量，进而确定出任一种形状的油藏区域对应的地面区域的沉降量，实用性较好，而且应用范围广。
4.为解决上述技术问题，本发明提供了一种油藏开发引起的地面沉降量的确定方法，包括：
5.在确定进行油藏开发计划的油藏区域后，预测对所述油藏区域进行所述油藏开发计划后的内部孔隙压力；
6.确定所述油藏区域的外部孔隙压力；
7.确定所述油藏区域的表观体积的第一体积模量；
8.确定所述油藏区域内的岩石实体体积的第二体积模量；
9.根据所述第一体积模量、所述第二体积模量以及所述内部孔隙压力分别确定所述油藏区域在预设坐标系中的三个坐标轴方向上的体积力，所述预设坐标系的横轴为x轴，纵轴为y轴，竖轴为z轴；
10.根据所述第一体积模量、所述第二体积模量、所述内部孔隙压力以及所述外部孔隙压力分别确定所述油藏区域的油藏边界在所述预设坐标系中的三个坐标轴方向上的表面力；
11.基于各个所述体积力、各个所述表面力、预设泊松比以及预设剪切模量确定所述油藏区域对应的地面区域的沉降量。
12.优选的，根据所述第一体积模量、所述第二体积模量以及所述内部孔隙压力分别确定所述油藏区域在预设坐标系中的三个坐标轴方向上的体积力，包括：
13.确定所述第一体积模量与所述第二体积模量之间的第一比值；
14.分别确定所述内部孔隙压力在所述油藏区域内的各个坐标点的三个所述坐标轴方向上的偏导数；
15.根据所述第一比值和各个所述偏导数确定所述油藏区域在预设坐标系中的三个坐标轴方向上的体积力。
16.优选的，根据所述第一体积模量、所述第二体积模量、所述内部孔隙压力以及所述外部孔隙压力分别确定所述油藏区域的油藏边界在所述预设坐标系中的三个坐标轴方向上的表面力，包括：
17.确定所述第一体积模量与所述第二体积模量之间的第二比值；
18.确定所述外部孔隙压力与所述内部孔隙压力之间的压力差值；
19.分别确定所述油藏区域的油藏边界在所述预设坐标系的三个坐标轴方向上的单位法向量；
20.根据所述第二比值、所述压力差值和各个所述单位法向量确定所述油藏区域的油藏边界在所述预设坐标系中的三个坐标轴方向上的表面力。
21.优选的，预测对所述油藏区域进行所述油藏开发计划后的内部孔隙压力，包括：
22.确定所述油藏区域中的当前油藏数值；
23.预测对所述油藏区域进行所述油藏开发计划后的剩余油藏数值；
24.根据所述当前油藏数值和所述剩余油藏数值确定所述内部孔隙压力。
25.优选的，基于各个所述体积力、各个所述表面力、预设泊松比以及预设剪切模量确定所述油藏区域对应的地面区域的沉降量，包括：
26.确定所述油藏区域对应的半无限大空间；
27.确定所述油藏区域中的各个坐标点在三个所述坐标轴方向上受到的所述半无限大空间的第一作用力；
28.基于各个所述体积力、各个所述表面力、各个所述作用力、所述预设泊松比以及所述预设剪切模量确定所述油藏区域中的各个坐标点在所述预设坐标系上的垂直位移量；
29.根据各个所述垂直位移量确定所述地面区域的沉降量。
30.优选的，确定所述油藏区域中的各个坐标点在三个所述坐标轴方向上受到的所述半无限大空间的第一作用力，包括：
31.对于任一个所述油藏区域中的坐标点，确定所述坐标点在三个所述坐标轴方向上受到的所述半无限大空间中各个位置的第二作用力；
32.将所述坐标点在同一所述坐标轴方向上受到的所有所述第二作用力的总和作为所述坐标点在所述坐标轴方向上受到的所述第一作用力。
33.本技术还提供一种油藏开发引起的地面沉降量的确定系统，包括：
34.第一压力确定单元，用于在确定进行油藏开发计划的油藏区域后，预测对所述油藏区域进行所述油藏开发计划后的内部孔隙压力；
35.第二压力确定单元，用于确定所述油藏区域的外部孔隙压力；
36.第一模量确定单元，用于确定所述油藏区域的表观体积的第一体积模量；
37.第二模量确定单元，用于确定所述油藏区域内的岩石实体体积的第二体积模量；
38.体积力确定单元，用于根据所述第一体积模量、所述第二体积模量以及所述内部孔隙压力分别确定所述油藏区域在预设坐标系中的三个坐标轴方向上的体积力，所述预设坐标系的横轴为x轴，纵轴为y轴，竖轴为z轴；
39.表面力确定单元，用于根据所述第一体积模量、所述第二体积模量、所述内部孔隙压力以及所述外部孔隙压力分别确定所述油藏区域的油藏边界在所述预设坐标系中的三个坐标轴方向上的表面力；
40.沉降量确定单元，用于基于各个所述体积力、各个所述表面力、预设泊松比以及预设剪切模量确定所述油藏区域对应的地面区域的沉降量。
41.本技术还提供一种油藏开发引起的地面沉降量的确定装置，包括：
42.存储器，用于存储计算机程序；
43.处理器，用于执行所述计算机程序时实现如上述的油藏开发引起的地面沉降量的确定方法的步骤。
44.本技术还提供一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现如上述的油藏开发引起的地面沉降量的确定方法的步骤。
45.本发明提供了一种油藏开发引起的地面沉降量的确定方法及相关组件，涉及油藏开发领域，在确定进行油藏开发计划的油藏区域后，首先确定油藏区域的内部孔隙压力和外部孔隙压力，然后确定油藏区域的表观体积的第一体积模量以及第二体积模量，第二体积模量为油藏区域中成分最多的岩石的体积模量，根据第一体积模量、第二体积模量、内部孔隙压力和外部孔隙压力分别确定油藏区域在预设坐标系中的三个坐标轴方向上的体积力和油藏区域的油藏边界的表面力，最后基于各个体积力、各个表面力、预设泊松比以及预设剪切模量确定油藏区域对应的地面区域的沉降量。由于在确定体积模量时不会受到油藏区域本身形状的影响，所以能够确定任一种形状的油藏区域的体积模量，进而确定出任一种形状的油藏区域对应的地面区域的沉降量，实用性较好，而且应用范围广。
附图说明
46.为了更清楚地说明本发明实施例中的技术方案，下面将对现有技术和实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
47.图1为本技术提供的一种油藏开发引起的地面沉降量的确定方法的流程图；
48.图2为压应力为正时油藏区域的微元体应力的示意图；
49.图3为一种圆盘形油藏区域的示意图；
50.图4为本技术提供的一种地面沉降量结果的示意图；
51.图5为本技术提供的一种半无限大空间内坐标点受力作用的示意图；
52.图6为本技术提供的一种油藏开发引起的地面沉降量的确定系统的结构示意图；
53.图7为本技术提供的一种油藏开发引起的地面沉降量的确定装置的结构示意图。
具体实施方式
54.本发明的核心是提供一种油藏开发引起的地面沉降量的确定方法及相关组件，能够确定任一种形状的油藏区域的体积模量，进而确定出任一种形状的油藏区域对应的地面区域的沉降量，实用性较好，而且应用范围广。
55.为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
56.请参照图1，图1为本技术提供的一种油藏开发引起的地面沉降量的确定方法的流程图，包括：
57.s1：在确定进行油藏开发计划的油藏区域后，预测对油藏区域进行油藏开发计划后的内部孔隙压力；
58.s2：确定油藏区域的外部孔隙压力；
59.s3：确定油藏区域的表观体积的第一体积模量；
60.s4：确定油藏区域内的岩石实体体积的第二体积模量；
61.s5：根据第一体积模量、第二体积模量以及内部孔隙压力分别确定油藏区域在预设坐标系中的三个坐标轴方向上的体积力，预设坐标系的横轴为x轴，纵轴为y轴，竖轴为z轴；
62.s6：根据第一体积模量、第二体积模量、内部孔隙压力以及外部孔隙压力分别确定油藏区域的油藏边界在预设坐标系中的三个坐标轴方向上的表面力；
63.s7：基于各个体积力、各个表面力、预设泊松比以及预设剪切模量确定油藏区域对应的地面区域的沉降量。
64.为了确定出任一种形状的油藏区域在进行油藏开发后引起的地面沉降量，本技术中，在确定对某个油藏区域进行油藏开发时，由于油藏开发计划是已知的，也即在采集该油藏区域内的石油或者天然气时，采集得到的总量或总体积是预先知道的，所以可以根据需要采集的总量如油藏数值来模拟此次采集之后该油藏区域内的内部孔隙压力；而对于油藏区域的外部孔隙压力，其指的是油藏区域外部的地层的孔隙压力，由于在对油藏区域进行油藏开发时对外部的地层产生的影响较小，其孔隙压力变化不大，所以可以将外部孔隙压力视为一个常数；油藏区域的表观体积为包含油藏区域内部空间在内的整个油藏区域的体积，相当于将整个油藏区域视为一个实心的整体时的总体积，而固体实体体积指的是在该油藏区域中所有固体岩石的总体积，例如，当该油藏区域中包含10％石油、15％天然气、25％石灰岩以及50％石英岩时，该油藏区域中的岩石实体体积即为石灰岩和石英岩的体积，可以将其称为岩石骨架或骨架体积。
65.在确定了内部和外部的孔隙压力，第一体积模量和第二体积模量后，根据这些参数来确定油藏区域的体积力和表面力，具体的，由于油藏区域处于地下而且地下环境较为复杂，难以根据油藏区域的实际形状确定油藏区域的体积，此时考虑到作用在油藏区域上的力的大小与油藏区域的体积成正比，所以通过体积力的大小可以确定出油藏区域的体积；而表面力指的是作用到油藏区域内的流体(如原油)表面上的力，表面力与流体表面的面积成正比，当油藏开发采集了一部分的流体时，表面力会发生变化，通过表面力的大小可
以确定出油藏区域在进行油藏开发后的应力变化；所以基于体积力和表面力以及预设的剪切模量，能够确定出该油藏区域在进行油藏开发后对该油藏区域对应的地面区域造成的沉降影响。
66.更具体的，由于油藏区域的组成复杂且其会受到来自各个方向的力，可以将油藏区域本身视为一个各向同性弹性介质，在油藏区域中，岩石骨架作用在油藏区域的表观体积上的平均骨架应力σ等于各个在坐标点上的实际骨架应力与1减去骨架岩石的孔隙度的差之间的乘积之间的平均值，请参照图2，图2为压应力为正时油藏区域的微元体应力的示意图，在对油藏区域中某个坐标点(x，y，z)进行受力分析时，由于其会受到来自各个方向的各个点上的力的作用，所以可以将该点视为一个微元体，在压应力为正时该坐标点的应力平衡方程为：
[0067][0068][0069][0070]
其中，fi为该点受到的总外力，σ
ij
为该点在对应方向上的应力，其中i和j均可以是x或y或z。此外，该点上的柯西应变为：
[0071][0072][0073][0074]
其中，ε
ij
为该点在对应方向上的应变量，ui为该点在对应方向上的经过油藏开发后的位移量，i和j均可以是x或y或z。基于该点上的骨架应力以及柯西应变，可以确定出各向同性弹性介质的线性本构方程，为：
[0075][0076]
其中，e
p
油藏区域的表观弹性模量，v
p
为油藏区域的表观泊松比，g为预设的剪切模量。当对油藏区域进行油藏开发时，油藏区域内部的可流动区域的孔隙压力会发生变化，而
应力与孔隙压力为正相关；同时，由于岩石孔喉中的流体渗流速度极慢，可以将这些流体视为静态的流体，也即这些流体的受力平衡，所以岩石骨架给孔喉中流体的作用力与该流体的压差提供的力的大小相等且方向相反，而该作用力的反作用力会作用到岩石骨架上，从而会引起岩石骨架的应力变化，而流体对岩石骨架中作用的体积力为：
[0077][0078]
其中，f为体积力，φ为孔隙度，p为孔隙压力，l为指向的某个方向，l可以是x或y或z，将式4代入式1可以得到油藏区域的骨架应力平衡公式：
[0079][0080][0081][0082]
考虑到岩石骨架中的孔隙或者孔喉被流体充满，可以在岩石骨架外部的不同方向均匀施加一个对岩石骨架的应力t，而t＝p，可以将此时的岩石骨架视为一个无孔的完整固体结构，此时可以确定岩石骨架的应变为：
[0083][0084]
其中，e0为岩石骨架材料本身的弹性模量，v0为岩石骨架材料本身的泊松比，为岩石骨架在对应方向上的第一应变量，i和j均可以是x或y或z，然后在对岩石骨架外部施加总应力，有：
[0085]
t＝σ
t,ij-pδ
ij
＝σ
ij-(1-φ)pδ
ij
ꢀꢀꢀ
(式7)
[0086]
其中，σ
t，ij
为对应方向上的总应力，δ
ij
为克罗内克符号，i和j均可以是x或y或z，这时岩石骨架相当于一个孔隙度为φ但是孔隙中无流体的多孔固体，此时岩石骨架的应变为：
[0087][0088][0089][0090]
[0091]
其中，为岩石骨架在对应方向上的第二应变量，i和j均可以是x或y或z，根据线弹性假设，应用叠加原理可以得到岩石骨架的总应变为：
[0092][0093]
其中，为岩石骨架的总应变，基于此，可得：
[0094][0095][0096][0097][0098]
其中，k
p
为油藏区域的表观体积模量也即第一体积模量，k0为岩石骨架材料本身的体积模量也即第二体积模量，将上述公式表示为应力解的形式有：
[0099][0100][0101][0102][0103]
其中，g
p
为油藏区域的表观剪切模量，然后，将上述式2、式5以及式11相结合，得到由位移表示的变形微分方程组：
[0104][0105]
由于上述公式与弹性力学位移法求解的lame-navier方程形式一致，而lame-navier方程中的体积力项等于上述公式中的孔隙压力梯度项(1—k
p
/k0)。所以油气藏开发过程中的骨架应力变化所遵循的力学规律等价于受一定体积力的连续介质应力变化，其等
效体积力为：
[0106][0107][0108][0109]
其中，f
′i为油藏区域在对应方向上的等效体积力，其是一种虚拟力但其是一种能够等效于实际体积力的参数，i可以是x或y或z，在油藏区域的封闭边界处，也即油藏区域中流体能够进行渗流的边界处，该边界两侧的法向总应力平衡，确定一个该边界的法向单位向量n＝(n
x
，ny，nz)，n
x
表示为该边界往x方向上的单位向量，ny表示为该边界往y方向上的单位向量，nz表示为该边界往z方向上的单位向量，可得：
[0110][0111][0112][0113]
其中，为流体外侧的对应方向上的应力，为流体内侧的对应方向上的应力，p
out
为油藏区域外侧的孔隙压力，p
in
为油藏区域内侧的孔隙压力，结合式11与式14，可得：
[0114][0115][0116][0117]
其中，为流体外侧的对应方向上的应变，为流体内侧的对应方向上的应变。由于上述公式的等号左侧部分与常规的无孔介质的平面静力平衡方程的等号左侧部分一致，而常规的无孔介质的平面静力平衡方程的等号右侧部分为0，可见封闭边界面引起的地层变形相当于在油藏区域内部对应为置施加的表面力，其表面力为：
[0118]
[0119][0120][0121]
其中，为油藏区域在对应方向上的等效表面力，i可以是x或y或z，通过上述公式，能够确定出油藏区域的在各个方向上的体积力和油藏边界的表面力，以便于后续通过等效体积力f
′i、等效表面力和预设的剪切模量g来确定油藏区域对应的地面区域的沉降量。请参照图3和图4，图3为一种圆盘形油藏区域的示意图，图4为本技术提供的一种地面沉降量结果的示意图，该圆盘形油藏区域纵深也即厚度为10m，与地面之间的距离为3000m，半径为1000m，利用上述公式对图3中的圆盘形油藏区域，若预设的剪切模量为30gpa，预设泊松比为0.25，k
p
为40gpa，k0为500gpa，若对该油藏区域进行油藏开发后，油藏区域中各点的压力均下降10mpa，经过上述计算后可得到如图3中的地面沉降量结果，图3中的横轴为地面上的某个地点与油藏区域圆心对应的地点之间的距离，单位为米，纵轴为沉降量，单位为毫米，可见，油藏区域圆心对应的地点的沉降量约为0.08毫米，距离该圆心地点5000米远的地点的沉降量约为0.015毫米。
[0122]
综上，在确定进行油藏开发计划的油藏区域后，首先确定油藏区域的内部孔隙压力和外部孔隙压力，然后确定油藏区域的表观体积的第一体积模量以及第二体积模量，第二体积模量为油藏区域中成分最多的岩石的体积模量，根据第一体积模量、第二体积模量、内部孔隙压力和外部孔隙压力分别确定油藏区域在预设坐标系中的三个坐标轴方向上的体积力和油藏区域的油藏边界的表面力，最后基于各个体积力、各个表面力、预设泊松比以及预设剪切模量确定油藏区域对应的地面区域的沉降量。根据对体积力与表面力的分布式表征能够确定任一种形状的油藏区域的体积模量，进而确定出任一种形状的油藏区域对应的地面区域的沉降量，实用性较好，而且应用范围广。
[0123]
在上述实施例的基础上：
[0124]
作为一种优选的实施例，根据第一体积模量、第二体积模量以及内部孔隙压力分别确定油藏区域在预设坐标系中的三个坐标轴方向上的体积力，包括：
[0125]
确定第一体积模量与第二体积模量之间的第一比值；
[0126]
分别确定内部孔隙压力在油藏区域内的各个坐标点的三个坐标轴方向上的偏导数；
[0127]
根据第一比值和各个偏导数确定油藏区域在预设坐标系中的三个坐标轴方向上的体积力。
[0128]
在确定油藏区域在各个坐标轴方向上的体积力时，有：
[0129][0130]
[0131][0132]
其中，f
′i为i方向上的等效体积力，i可以是x或y或z，k
p
为油藏区域的表观体积模量也即第一体积模量，k0为岩石骨架材料本身的体积模量也即第二体积模量，p为内部孔隙压力。
[0133]
作为一种优选的实施例，根据第一体积模量、第二体积模量、内部孔隙压力以及外部孔隙压力分别确定油藏区域的油藏边界在预设坐标系中的三个坐标轴方向上的表面力，包括：
[0134]
确定第一体积模量与第二体积模量之间的第二比值；
[0135]
确定外部孔隙压力与内部孔隙压力之间的压力差值；
[0136]
分别确定油藏区域的油藏边界在预设坐标系的三个坐标轴方向上的单位法向量；
[0137]
根据第二比值、压力差值和各个单位法向量确定油藏区域的油藏边界在预设坐标系中的三个坐标轴方向上的表面力。
[0138]
在确定表面力时，有：
[0139][0140][0141][0142]
其中，为i方向上的等效表面力，i可以是x或y或z，k
p
为油藏区域的表观体积模量也即第一体积模量，k0为岩石骨架材料本身的体积模量也即第二体积模量，p
out
为外部孔隙压力，p
in
为内部孔隙压力也即p。
[0143]
作为一种优选的实施例，预测对油藏区域进行油藏开发计划后的内部孔隙压力，包括：
[0144]
确定油藏区域中的当前油藏数值；
[0145]
预测对油藏区域进行油藏开发计划后的剩余油藏数值；
[0146]
根据当前油藏数值和剩余油藏数值确定内部孔隙压力。
[0147]
为了确定内部孔隙压力，本技术中，考虑到油藏开发时会对油藏区域中的油藏数值产生影响，而油藏数值发生变化时内部孔隙压力也会随之产生变化，可见油藏数值与内部孔隙压力之间有对应关系，根据其对应关系可以简单的确定出内部空隙压力。
[0148]
作为一种优选的实施例，基于各个体积力、各个表面力、预设泊松比以及预设剪切模量确定油藏区域对应的地面区域的沉降量，包括：
[0149]
确定油藏区域对应的半无限大空间；
[0150]
确定油藏区域中的各个坐标点在三个坐标轴方向上受到的半无限大空间的第一作用力；
[0151]
基于各个体积力、各个表面力、各个作用力、预设泊松比以及预设剪切模量确定油藏区域中的各个坐标点在预设坐标系上的垂直位移量；
[0152]
根据各个垂直位移量确定地面区域的沉降量。
[0153]
为了确定出油藏区域的沉降量，本技术中，由于实际的油藏区域是一片面积很大但是纵深有限的区域，相当于油藏区域是一个垂直于地面的方向上的范围无限而平行与地面的方向上的范围无限，所以可以将油藏区域视为一个半无限大的空间，通过映射确定出油藏区域的半无限大空间，再利用半无限大空间受到作用力时的基本解来确定各个坐标点的沉降量，而对于每个坐标点，均有：
[0154][0155]
其中，uz(x,y,z)为该坐标点的沉降量，ω为油藏区域整体的范围，г为油藏区域中的油藏边界，为坐标点在i方向上受到的第一作用力。具体的，对于半无限大空间内的任意一点的垂向位移为：
[0156][0157]
其中，f为半无限大空间内沿坐标系横轴方向上存在的力也即某一无限远方向上存在的力，π为圆周率，uz为该点的垂向位移，x，y，z分别为该点的在预设坐标系上的坐标值，v为油藏区域的泊松比，g为预设剪切模量，x
′
，y
′
，z
′
分别为发力点，也即给(x，y，z)这个坐标点作力的点，r1和r2分别为：
[0158][0159][0160]
基于此，可以确定不同方向上的位移量，公式为：
[0161][0162][0163][0164]
其中，为平行与预设坐标系竖轴方向上的位移量，为预设坐标系竖轴方向较偏向于横轴方向上的位移量，为预设坐标系竖轴方向较偏向于纵轴方向上的位移量。
[0165]
作为一种优选的实施例，确定油藏区域中的各个坐标点在三个坐标轴方向上受到的半无限大空间的第一作用力，包括：
[0166]
对于任一个油藏区域中的坐标点，确定坐标点在三个坐标轴方向上受到的半无限大空间中各个位置的第二作用力；
[0167]
将坐标点在同一坐标轴方向上受到的所有第二作用力的总和作为坐标点在坐标轴方向上受到的第一作用力。
[0168]
在确定坐标点受到的第一作用力时，考虑到该坐标点会受到各个方向给予的第二作用力，请参照图5，图5为本技术提供的一种半无限大空间内坐标点受力作用的示意图，而不同方向上的不同距离给予该坐标点的第二作用力也不同，相当于在同一方向上有多个发力点向该坐标点发力，所以需要确定出坐标点在各个坐标轴方向上收到的第二作用力的总和，其总和即为坐标点在该坐标轴方向上的全部受力，也即第一作用力。
[0169]
请参照图6，图6为本技术提供的一种油藏开发引起的地面沉降量的确定系统的结构示意图，包括：
[0170]
第一压力确定单元21，用于在确定进行油藏开发计划的油藏区域后，预测对油藏区域进行油藏开发计划后的内部孔隙压力；
[0171]
第二压力确定单元22，用于确定油藏区域的外部孔隙压力；
[0172]
第一模量确定单元23，用于确定油藏区域的表观体积的第一体积模量；
[0173]
第二模量确定单元24，用于确定油藏区域内的岩石实体体积的第二体积模量；
[0174]
体积力确定单元25，用于根据第一体积模量、第二体积模量以及内部孔隙压力分别确定油藏区域在预设坐标系中的三个坐标轴方向上的体积力，预设坐标系的横轴为x轴，纵轴为y轴，竖轴为z轴；
[0175]
表面力确定单元26，用于根据第一体积模量、第二体积模量、内部孔隙压力以及外部孔隙压力分别确定油藏区域的油藏边界在预设坐标系中的三个坐标轴方向上的表面力；
[0176]
沉降量确定单元27，用于基于各个体积力、各个表面力、预设泊松比以及预设剪切模量确定油藏区域对应的地面区域的沉降量。
[0177]
对于本技术提供的一种油藏开发引起的地面沉降量的确定系统的详细介绍，请参照上述油藏开发引起的地面沉降量的确定方法的实施例，本技术在此不再赘述。
[0178]
请参照图7，图7为本技术提供的一种油藏开发引起的地面沉降量的确定装置的结构示意图，包括：
[0179]
存储器31，用于存储计算机程序；
[0180]
处理器32，用于执行计算机程序时实现如上述的地面沉降量的确定方法的步骤。
[0181]
对于本技术提供的一种油藏开发引起的地面沉降量的确定装置的详细介绍，请参照上述油藏开发引起的地面沉降量的确定方法的实施例，本技术在此不再赘述。
[0182]
本技术还提供一种计算机可读存储介质，计算机可读存储介质上存储有计算机程序，计算机程序被处理器执行时实现如上述的地面沉降量的确定方法的步骤。
[0183]
对于本技术提供的一种计算机可读存储介质的详细介绍，请参照上述油藏开发引起的地面沉降量的确定方法的实施例，本技术在此不再赘述。
[0184]
还需要说明的是，在本说明书中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个
……”
限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
[0185]
对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。
对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其他实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。