1.地质碳封存区域大气CO2时空分异特征及影响因子分析方法，其特征在于：将存在废弃或衰竭的油气藏、不可开采煤层并作为CO2的地质碳封存区域，通过考虑常规的地表因素和地质碳封存区域独有的碳泄露因素，分析地质碳封存区域近地面大气CO2时空分异特征和影响因素，其步骤如下：步骤1：采集地质碳封存区域的大气CO2浓度数据、地质构造数据、遥感数据产品、气象观测资料，其中地质构造数据包括地质勘查数据、矿井分布数据，遥感数据产品包括植被覆盖度数据、数字高程模型数据，气象观测资料为风速风向数据；

获取OCO‑2卫星数据作为近地面大气CO2浓度数据，利用ArcGIS软件将OCO‑2卫星数据的原始点带数据重建为0.05°×0.05°空间分辨率的栅格数据作为地质碳封存区域的大气CO2浓度数据，地质构造数据包括地质勘查数据和矿井分布数据，并对其进行数字化和地理配准，其他驱动因子数据通过遥感数据产品和气象观测资料获取；

步骤2：基于数据定性分析地质碳封存区域大气CO2浓度数据的影响因子，影响因子包括地质碳封存区域的CO2泄露影响因子和其他影响因子；其中泄露影响因子包括地质碳封存区域的断层和裂缝分布、井筒分布数据，其他影响因子包括地表地形起伏度、风速风向、植被覆盖度数据；

步骤3：将泄露影响因子和其他影响因子作为候选解释变量，利用空间自相关Global Moran's I工具探析各候选解释变量之间的空间相关性，用Moran's I指数来表示变量之间的空间相关关系，并利用探索性回归工具对候选解释变量的所有可能的变量组合进行测试评估，筛选获得最佳的解释变量组合；

使用ArcGIS软件中的空间自相关工具和探索性回归工具对候选解释变量进行分析处理，以筛选获得最佳的解释变量组合，具体处理步骤如下：首先指定每个普通最小二乘法OLS模型包含的解释变量的最小数量和最大数量，以及

2

Adjusted R、系数p值、方差膨胀因子VIF值、Jarque‑Berap值和空间自相关q值的阈值条2

件，当且仅当找到一个满足可接受的最小Adjusted R 、最大系数p值边界、最大VIF值边界和可接受的最小Jarque Berap值阈值条件的OLS模型时，探索性回归工具会对该模型的残差运行空间自相关Global Moran's I工具，以判断偏低/偏高预计值是否会产生聚集；如果空间自相关q值大于在探索性回归工具的搜索条件中指定的值，即可接受的最小空间自相关q值时，则该模型将被视为合格模型列出于结果输出表；最后选取结果输出表中具有最大2

Adjusted R的合格OLS模型，将其所含解释变量组合作为最佳解释变量组合；

步骤4：将筛选出的最佳解释变量组合作为自变量，将地质碳封存区域的大气CO2浓度作为因变量构建地理加权回归模型，定量分析地质碳封存区域大气CO2的时空分异特征以及最佳解释变量对该区域内大气CO2时空分异的不同影响程度；

将最佳解释变量组合作为自变量，大气CO2浓度作为因变量构建地理加权回归模型，对气象站观测点的观测值进行计算，定量分析地质碳封存区域大气CO2的时空分异特征以及最佳解释变量对该区域内大气CO2时空分异的不同影响程度，具体步骤为：a.确定带宽：采用交叉验证法CV来确定带宽，选取合适的带宽h，使CV(h)达到最小，如下构造CV：其中， 是取带宽h时，去掉观测点(yi；xi1,xi2,…,xin)算得的大气CO2观测值的拟合值，yi表示观测值，xi1,xi2,…,xin表示最佳解释变量，h表示所取带宽，m表示交叉验证数据的个数；

b.确定核函数：根据a所求得带宽确定核函数，采用高斯函数作为核函数，函数表示如下：

其中，h表示带宽，di0是(ui,vi)与(u0,v0)间的距离，(ui,vi)表示样本点坐标，(u0,v0)表示回归点坐标；

c.计算权重：采用高斯函数时，样本点(ui,vi)对回归点(u0,v0)的权重为：

其中，wi(u0,v0)是空间权重，(ui,vi)表示样本点坐标，(u0,v0)表示回归点坐标，h表示带宽；

d.构建回归方程：以其中一个观测点(yi；xi1,xi2,…,xin)为中心，以带宽h为半径确定计算范围，将计算范围里面的所有观测点纳入计算，同时根据核函数确定该观测点的权重；重复以上步骤直至完成所有观测点的权重计算，得到各观测点(yi；xi1,xi2,…,xin)在空间位置为(ui,vi)的大气CO2观测值yi与解释变量xi1,xi2,…,xin的回归方程为：yi＝βi(ui,vi)+β1(ui,vi)xi1+β2(ui,vi)xi2+…+βn(ui,vi)xin+εi2

其中，i＝1,2,…,m，εi独立同分布，εi～N(0,σ)，每个样本点(ui,vi)的回归系数为βk(ui,vi)(k＝1,2,…,n)。

2.根据权利要求1所述的地质碳封存区域大气CO2时空分异特征及影响因子分析方法，其特征在于步骤2中地质碳封存区域泄露影响因子分析，基于OCO‑2卫星数据进行分析，将地质碳封存区域划分成0.05°×0.05°的平面格网，根据断裂带和裂缝的数字化结果，将断裂带和裂缝对应的矢量线数据转为矢量点数据，对照格网填补断裂带和裂缝的点数据，保证断裂带和裂缝经过的格网都有点数据分布，再以OCO‑2卫星数据重建的大气CO2浓度数据为基础数据，将对应的CO2浓度值赋给断裂带生成的点数据，得到断裂带和裂缝处的CO2浓度数据；统计每个与CO2封存有关井筒的地理坐标，包括所有的注入井、封存井，提取对应格网的CO2浓度数据，得到与CO2封存有关井筒处的CO2浓度数据，作为候选解释变量中的泄露影响因子。

3.根据权利要求1所述的地质碳封存区域大气CO2时空分异特征及影响因子分析方法，其特征在于，采用0.05°×0.05°分辨率的数字高程模型DEM分析地形起伏度因子，计算地表地形起伏度的具体公式如下：RFi＝Hmax‑Hmin

其中，RFi为地表地形起伏度，Hmax为0.05°×0.05°格网分析窗口内的最大高程值，Hmin为上述0.05°×0.05°格网分析窗口内的最小高程值；

植被覆盖度因子采用归一化植被指数NDVI反映区域植被覆盖情况，采用0.05°×0.05°分辨率的NDVI数据分析植被覆盖度因子，最终生成目标分辨率栅格数据，每个格网对应一个NDVI值；

风速风向因子分析采用的是各省市的气象站观测资料数据，风向划分为12方向：顺时针将0°和360°位置划分为北、0‑45°为北东、45‑90°为东北、90°为东、90‑135°为东南、135‑

180°为南东、180°为南、180‑225°为南西、225‑270°西南、270°为西、270‑315°为西北、315‑

360°为北西；计算每个气象站的月均风速和主要风向，并根据站点的地理坐标将值对应到每个格网。

4.根据权利要求1所述的地质碳封存区域大气CO2时空分异特征及影响因子分析方法，其特征在于步骤3中地形起伏度因子、风速风向因子、植被覆盖度因子与地质碳封存区域大气CO2的浓度为非线性关系，对地形起伏度因子、风速风向因子、植被覆盖度因子这三个候选解释变量进行对数变换，为避免数量级差异带来的回归参数较小的问题，对所有候选解释变量包括泄露影响因子和其他影响因子进行归一化的处理，具体公式为：其中， 为归一化处理后的值，χi为归一化处理前的值，χmax为待归一化处理数据中的最大值，χmin为待归一化处理数据中的最小值，n为待处理变量数据的个数。