1.一种利用相邻地层互相关约束的地层面模型修正方法，其特征在于，包括：获取观测样本数据和待修正地层界面的已有模型；

根据所述观测样本分别统计所述待修正地层界面和辅助地层界面的高程边际特性函数，其中，所述辅助地层界面为所述待修正地层界面的相邻地层界面；

估计所述待修正地层界面的估计概率密度函数，作为对应的高程先验概率密度；

根据所述辅助地层界面在待修正位置上的数据精度确定辅助地层数，并根据所述辅助地层数、所述待修正地层界面的高程边际特性函数、所述辅助地层界面的高程边际特性函数以及相邻地层界面的高程变量的Copula函数，确定联合分布函数；

根据所述待修正地层界面的高程边际特性函数和所述联合分布函数，确定所述待修正地层界面至少一个坐标点对应的似然函数；

根据所述至少一个坐标点对应的似然函数，对所述待修正地层界面的高程先验概率密度进行贝叶斯更新并进行归一化处理，确定所述至少一个坐标点的高程后验概率密度；

根据所述至少一个坐标点的高程后验概率密度，确定对应的模型修正值，对所述待修正地层界面的已有模型中的高程值进行更新。

2.根据权利要求1所述的利用相邻地层互相关约束的地层面模型修正方法，其特征在于，所述根据所述观测样本分别统计待修正地层界面和辅助地层界面的高程边际特性函数包括：

分别将所述待修正地层界面和所述辅助地层界面的高程值范围划分为多个等间距区间；

分别统计所述待修正地层界面和所述辅助地层界面的观测样本中的高程值落入不同的等间距区间的频率数据；

根据所述频率数据，分别确定所述辅助地层界面关于高程值的第一边际概率密度函数和所述待修正地层界面关于高程值的第二边际概率密度函数；

根据所述频率数据，分别确定所述辅助地层界面的第一累积概率密度函数和所述待修正地层界面的第二累积概率密度函数；

将所述第一累积概率密度函数变换为[0，1]区间内均匀分布下的分位数，确定所述辅助地层界面对应的第一分位数，将所述第二累积概率密度函数变换为[0，1]区间内均匀分布下的分位数，确定所述待修正地层界面对应的第二分位数。

3.根据权利要求2所述的利用相邻地层互相关约束的地层面模型修正方法，其特征在于，所述估计所述待修正地层界面的估计概率密度函数，作为对应的高程先验概率密度包括：

获取所述待修正地层界面的地层接触样本数据；

通过克里金方法，确定所述至少一个坐标点对应的克里金方差；

根据所述待修正地层界面的已有模型，确定所述至少一个坐标点对应的期望高程；

将所述克里金方差作为正态分布的方差、所述期望高程作为正态分布的期望，建立正态分布函数，其中，所述正态分布函数作为所述待修正地层界面的高程值在所述至少一个坐标点下的估计概率密度函数；

将所述待修正地层界面的估计概率密度函数作为对应的所述高程先验概率密度。

4.根据权利要求3所述的利用相邻地层互相关约束的地层面模型修正方法，其特征在于，所述辅助地层界面包括位于所述待修正地层界面上侧的第一辅助地层界面、位于所述待修正地层界面下侧的第二辅助地层界面，所述根据所述辅助地层界面在待修正位置上的数据精度确定辅助地层数包括：

分别获取所述第一辅助地层界面的已有模型和所述第二辅助地层界面的已有模型；

通过模型不确定性评价方法，分别确定所述第一辅助地层界面的已有模型中的第一高程误差、所述第二辅助地层界面的已有模型中的第二高程误差；

根据所述第一高程误差和所述第二高程误差，分别判断所述第一辅助地层界面的已有模型和所述第二辅助地层界面的已有模型是否满足数据精度条件；

若其中之一满足，则将满足所述数据精度条件的所述第一辅助地层界面或所述第二辅助地层界面作为所述辅助地层界面，所述辅助地层数为一层；

若都满足，则所述第一辅助地层界面和所述第二辅助地层界面都作为所述辅助地层界面，所述辅助地层数为两层。

5.根据权利要求4所述的利用相邻地层互相关约束的地层面模型修正方法，其特征在于，所述根据所述辅助地层数、所述待修正地层界面的高程边际特性函数、所述辅助地层界面的高程边际特性函数以及相邻地层界面的高程变量的Copula函数，确定联合分布函数包括：

若所述辅助地层数为一层，则根据所述第一边际概率密度函数、所述第二边际概率密度函数、以及所述相邻地层界面的高程变量的Copula函数，确定所述联合分布函数，通过以下公式得到所述联合分布函数：

fA,B(ZA,ZB)＝fA(ZA)·fB(ZB)·cA,B(uA,uB)其中，fA,B(ZA,ZB)表示所述辅助地层数为一层时的所述联合分布函数，fA(ZA)表示所述第一边际概率密度函数，fB(ZB)表示所述第二边际概率密度函数，cA,B(uA,uB)表示所述相邻地层界面的高程变量的Copula函数，uA为所述第一分位数，uB为所述第二分位数。

6.根据权利要求5所述的利用相邻地层互相关约束的地层面模型修正方法，其特征在于，所述根据所述待修正地层界面的高程边际特性函数和所述联合分布函数，确定所述待修正地层界面至少一个坐标点对应的似然函数包括：若所述辅助地层数为一层，则在所述第一辅助地层界面或所述第二辅助地层界面的已有模型中，确定所述至少一个坐标点对应的模型高程值；

根据所述模型高程值、所述第二边际概率密度函数和所述联合分布函数，通过以下公式确定所述似然函数：

其中，L(ZB|ZA＝ZA)、fA|B(ZA＝ZA|ZB)表示所述待修正地层界面在所述至少一个坐标点处对应的所述似然函数，ZA表示所述辅助地层界面对应的高程变量，ZB表示所述待修正地层界面对应的高程变量，VA表示所述模型高程值，fA，B(ZA＝VA，ZB)表示所述辅助地层数为一层时的所述联合分布函数，fB(ZB)表示所述第二边际概率密度函数。

7.根据权利要求4所述的利用相邻地层互相关约束的地层面模型修正方法，其特征在于，所述第一边际概率密度函数包括所述第一辅助地层界面的边际概率密度函数、所述第二辅助地层界面的边际概率密度函数，所述根据所述辅助地层数、所述待修正地层界面的高程边际特性函数、所述辅助地层界面的高程边际特性函数以及相邻地层界面的高程变量的Copula函数，确定联合分布函数包括：若所述辅助地层数为两层，则根据第二边际概率密度函数、所述第一辅助地层界面的边际概率密度函数、所述第二辅助地层界面的边际概率密度函数、以及所述相邻地层界面的高程变量的Copula函数，确定所述联合分布函数，通过以下公式得到所述联合分布函数：fA,B,C(ZA,ZB,ZC)＝fA(ZA)·fB(ZB)·fC(ZC)·cA,B(uA,uB)·cB,C(uB,uC)·cA,C|B(uA|B,uC|B)其中，fA，B，C(ZA，ZB，ZC)表示所述辅助地层数为两层时的所述联合分布函数，fB(ZB)表示所述第二边际概率密度函数，fA(ZA)表示所述第一辅助地层界面的边际概率密度函数，fC(ZC)表示所述第二辅助地层界面的边际概率密度函数，cA，B(uA，uB)、cB，C(uB，uC)以及cA，CB(uAB，uC|B)表示所述相邻地层界面的高程变量的Copula函数，uA为所述第一辅助地层界面对应的所述第一分位数，uB为所述第二辅助地层界面对应的所述第一分位数，uC为所述第二分位数。

8.根据权利要求7所述的利用相邻地层互相关约束的地层面模型修正方法，其特征在于，所述根据所述待修正地层界面的高程边际特性函数和所述联合分布函数，确定所述待修正地层界面至少一个坐标点对应的似然函数包括：若所述辅助地层数为两层，则在所述第一辅助地层界面和所述第二辅助地层界面的已有模型中，确定所述至少一个坐标点对应的模型高程值；

根据所述模型高程值、所述第二边际概率密度函数和所述联合分布函数，通过以下公式确定所述似然函数：

其中，L(ZB|ZA＝VA，ZC＝VC)、fA，C|B(ZA＝VA，ZC＝VC|ZB)表示所述待修正地层界面在所述至少一个坐标点处对应的所述似然函数，ZA表示所述第一辅助地层界面对应的高程变量，ZB表示所述待修正地层界面对应的高程变量，ZC表示所述第二辅助地层界面对应的高程变量，VA表示所述第一辅助地层界面对应的所述模型高程值，VC表示所述第二辅助地层界面对应的所述模型高程值，fA，B，C(ZA，ZB，ZC)表示所述辅助地层数为两层时的所述联合分布函数，fB(ZB)表示所述第二边际概率密度函数。

9.根据权利要求4所述的利用相邻地层互相关约束的地层面模型修正方法，其特征在于，所述根据所述至少一个坐标点对应的似然函数，对所述待修正地层界面的高程先验概率密度进行贝叶斯更新并进行归一化处理，确定所述至少一个坐标点的高程后验概率密度包括：

在所述至少一个坐标点处，根据所述似然函数和所述高程先验概率密度的乘积，对所述高程先验概率密度进行贝叶斯更新：对更新后的高程先验概率密度进行归一化处理。

10.根据权利要求1所述的利用相邻地层互相关约束的地层面模型修正方法，所述根据所述至少一个坐标点的高程后验概率密度，确定对应的模型修正值，对所述待修正地层界面的已有模型中的高程值进行更新包括：根据所述至少一个坐标点的高程后验概率密度，确定所述高程后验概率密度对应的条件期望值；

将所述条件期望值作为所述至少一个坐标点在所述待修正地层界面的已有模型中对应的模型修正值；

根据所述模型修正值，对所述至少一个坐标点在所述待修正地层界面的已有模型中的高程值进行更新。