联合地震位错模型和InSAR数据构建2017年九寨沟Mw 6.5地震同震三维形变场

彭颖; 许才军; 刘洋

doi:10.13203/j.whugis20200289

联合地震位错模型和InSAR数据构建2017年九寨沟Mw 6.5地震同震三维形变场

彭颖^1,,
许才军^{1, 2, 3, ,},
刘洋^{1, 2, 3}

1.
武汉大学测绘学院, 湖北武汉, 430079
2.
地球空间环境与大地测量教育部重点实验室, 湖北武汉, 430079
3.
自然资源部地球物理大地测量重点实验室, 湖北武汉, 430079

基金项目:

国家自然科学基金 41721003

国家自然科学基金 41874011

国家重点研发计划 2018YFC1503603

详细信息

作者简介:
彭颖，硕士，研究方向为InSAR高精度三维形变序列测量。pngying@163.com

通讯作者:
许才军，博士、教授。cjxu@sgg.whu.edu.cn

中图分类号: P237
计量
- 文章访问数: 843
- HTML全文浏览量: 202
- PDF下载量: 138
出版历程
- 收稿日期: 2020-06-14
- 网络出版日期: 2022-03-11
- 发布日期: 2022-11-04

Deriving 3D Coseismic Deformation Field of 2017 Jiuzhaigou Earthquake with Elastic Dislocation Model and InSAR Data

PENG Ying^1,,
XU Caijun^{1, 2, 3, ,},
LIU Yang^{1, 2, 3}

1.
School of Geodesy and Geomatics, Whuhan University, Wuhan 430079, China
2.
Key Laboratory of Geospace Environment and Geodesy, Ministry of Education, Wuhan University, Wuhan 430079, China
3.
Key Laboratory of Geophysical Geodesy, Ministry of Natural Resources, Wuhan 430079, China

Funds:

The National Natural Science Foundation of China 41721003

The National Natural Science Foundation of China 41874011

the National Key Research and Development Program of China 2018YFC1503603

More Information

Author Bio:
PENG Ying, master, specializes in 3D deformation field construction with InSAR data. E-mail: pngying@163.com

Corresponding author:
XU Caijun, PhD, professor. E-mail: cjxu@sgg.whu.edu.cn

摘要

摘要: 利用哨兵（Sentinel）-1A卫星升、降轨影像，在地震位错模型约束下获取了2017年九寨沟Mw 6.5地震的高质量三维形变场。首先，利用合成孔径雷达干涉测量技术（interferometric synthetic aperture radar，InSAR）提取九寨沟地震升、降轨同震形变场；然后，通过“两步法”反演获取该地震发震断层的几何参数和分布式滑动模型，以此为约束，采用方差分量估计算法联合解算九寨沟地震三维形变场。结果表明，九寨沟地震同震三维形变场以水平位移为主，垂向形变较弱；南北向形变呈拉张趋势，断层上盘向南、下盘向北滑动，最大位移分别为－19.81 cm和14.38 cm；东西向形变不对称性明显，断层上盘西北部向东水平运动，最大位移为18.37 cm，下盘东南部向西运动，最大位移不足8 cm。将南北、东西向形变与6个全球导航卫星系统（global navigation satellite system，GNSS）台站观测数据进行比较，两者一致性较好且均方根误差较小，分别为1.44 cm和1.77 cm，表明联合升、降轨InSAR观测和地震位错模型约束构建同震三维形变场方法具有较高可行性，显著降低了大地测量数据不足、InSAR观测对南北向形变不敏感等问题的影响。
- 九寨沟地震 /
- 差分干涉测量 /
- 位错模型 /
- 分布式滑动 /
- 三维形变场
Abstract:
Objectives The coseismic deformation field of 2017 Jiuzhaigou earthquake has been investigated by previous numerous studies. However, these studies most focus on the one-dimensional(1D) ground motion along the radar line of sight(LOS), thus resulting in anamorphic coseismic displacements.
Methods We construct high-quality three-dimensional(3D) coseismic deformation field using ascending and descending synthetic aperture radar(SAR) data from Sentinel-1A satellite with constraint of the elastic dislocation model to reflect the real movement patterns of 2017 Jiuzhaigou earthquake. Firstly, we use differential interferometry method to generate LOS coseismic deformation field with SAR images. Then, a two-step inversion algorithm is used to subdivide the fault plane, and obtain the geometric parameters and the optimal fault slip distribution. With constraint of this fault model, we take the variance component estimation approach to reconstruct 3D coseismic deformation field based on interferometric synthetic aperture radar(InSAR) measurements.
Results 3D coseismic deformation results show that the coseismic deformation field is dominated by horizontal displacement. In the north-south(N) direction, the maximum displacement of the hanging wall is －19.81 cm and 14.38 cm in the footwall. In the east-west(E) direction, the maximum displacement of the hanging wall is 18.37 cm for the northwestern, while that of the footwall is －7.84 cm for the Southeastern. In the vertical(U) direction, there is a little uplift in the northern of the fault and the maximum is about 3 cm.
Conclusions The derived north-south and east-west displacements are in good consistency of the global navigation satellite system investigations, indicating that combing InSAR measurements and elastic dislocation model to reconstruction 3D seismic deformation field is feasible and effective, which overcomes the problem of insufficient geodetic data and offers extensive future usage for measuring earthquake deformation.
- Jiuzhaigou earthquake /
- InSAR /
- elastic dislocation model /
- distributed slip /
- three-dimensional deformation field

HTML全文

在直线、平面拟合、空间直角坐标变换、自回归模型求解中,系数矩阵和右端观测量由随机和非随机元素组成,且同一随机元素会在不同的位置出现^[1-4]。常用的变量含误差（errors-in-variables,EIV）模型及相应的整体最小二乘（total least squares,TLS）方法假设系数阵中所有元素含误差。因此,顾及增广误差矩阵元素的随机特性及其相互关系的结构EIV（structured EIV,SEIV）模型和结构加权整体最小二乘（structured weighted TLS,SWTLS）方法应运而生^[5-10]。文献[5]采用极大似然估计方法从含有误差的时间序列数据中确定动态线性系统的参数,被认为是SWTLS问题研究的开端。文献[6]定义了结构整体最小二乘（structured TLS,STLS）这一术语,并将其转换为非线性广义奇异值分解问题求解。文献[7]将结构增广误差矩阵用其中的独立随机元素表示,并命名为约束TLS方法。对于某些列为固定元素而另外一些列为随机元素的混合LS-TLS问题,文献[8]引入正交三角分解计算模型参数。对于误差为非等权相关的情形,文献[9]发展了WLS-WTLS的迭代算法。文献[10]将SEIV模型用部分EIV（partial EIV,PEIV）模型表达,在加权最小二乘准则下推导了参数估计的Gauss-Newton型迭代算法,并进一步给出了参数的一阶近似方差、非线性信赖域区间及偏差。文献[11]给出了PEIV模型的两种迭代算法,其计算效率与系数矩阵中随机元素的数量有关。文献[12]提出了顾及系数阵元素和右端向量相关的广义PEIV模型。文献[13]将系数误差矩阵表示成已知矩阵和独立误差向量的乘积,推导了附有线性和二次约束的STLS问题的迭代算法。文献[14]和文献[15]采用变量投影法将增广系数矩阵表示成仿射结构矩阵与独立随机变量的乘积,然后将STLS问题转换为非线性等式约束优化问题求解。文献[16]研究了含多个右端观测向量的STLS问题。文献[17]提出了结构整体最小范数（structured total least norm,STLN）方法,这一方法可以最小化误差向量的范数,是STLS问题的重要拓展^[17-18]。文献[19-20]研究了STLN问题的快速算法。文献[21]将STLN拓展到系数阵和观测向量具有共同元素的情形。针对自回归（auto-regression, AR）模型这一特定STLS问题,文献[3-4]提出了两种新解法。

在STLS算法推导中,不同文献使用了不同的平差准则,其差异在于是否考虑独立误差的重复次数。如二维坐标转换中,源系统的纵横坐标值在系数矩阵中出现了2次,而目标系统的纵横坐标值在右端向量中只出现了1次。文献[1,7,10,13-15]给出的平差准则未考虑重复次数,而文献[3-4,17,21]均考虑了重复次数。不同的平差准则必定会得到不同的平差结果,其参数估值在统计意义上孰优孰劣,目前尚未给出明确结果。本文从函数模型和数值模拟两方面入手,证明了不考虑独立误差元素重复次数的平差准则能够得到统计意义下更优的解。

1 STLS平差模型及平差准则

EIV模型的函数表达式为^[1]：

y + e_{y} = (A + E_{A}) x

(1a)

式中, $y$ 和 $e_{y}$ 分别表示 $n$ 维观测向量及其误差； $A$ 和 $E_{A}$ 分别表示 $n \times m$ 维系数矩阵及其误差矩阵； $x$ 为 $m$ 维参数向量。误差向量的随机模型为：

e = [\begin{matrix} e_{A} \\ e_{y} \end{matrix}] ~ ([\begin{matrix} 0 \\ 0 \end{matrix}], σ_{0}^{2} Q_{e}) = ([\begin{matrix} 0 \\ 0 \end{matrix}], σ_{0}^{2} [\begin{matrix} Q_{A A} & Q_{A y} \\ Q_{y A} & Q_{y y} \end{matrix}])

(1b)

式中, $e_{A} = v e c (E_{A})$ , $v e c (∙)$ 表示矩阵向量化算子,即将 $n \times m$ 维矩阵的每一列从左至右叠加成一个 $n m$ 维列向量； $σ_{0}^{2}$ 是单位权方差； $Q_{A A}$ 和 $Q_{y y}$ 分别是 $e_{A}$ 和 $e_{y}$ 的对称正定协因数矩阵； $Q_{y A} = Q_{A y}^{T}$ 表示两者的协因数矩阵。当向量e的协因数矩阵 $Q_{e}$ 可逆时,采用WTLS准则求参数的最优估值：

Φ (e) = e^{T} P_{e} e

(2)

式中, $P_{e} = Q_{e}^{- 1}$ 为 $e$ 的权矩阵。当系数矩阵具有某种结构时, $Q_{e}$ 为秩亏矩阵无凯利逆。许多学者将EIV函数模型进行改化,一般是从结构误差矩阵或增广误差矩阵中提取独立误差向量,进一步构造目标函数求解。常用的几种EIV模型修正方法如下：

1） PEIV模型及其平差准则。PEIV模型选取系数阵中独立随机量的真值 $\bar{a}$ 作为待求量,其函数模型为^[10]：

y = (x^{T} \otimes I_{n}) (h + M \bar{a}) + e_{y}

(3a)

a = \bar{a} + e_{a}

(3b)

式中, $a$ 是系数矩阵中随机元素所构成的 $t$ 维列向量； $\bar{a}$ 和 $e_{a}$ 分别是相应的真值和误差向量； $I_{n}$ 表示 $n$ 维单位矩阵； $h$ 是已知的 $n m$ 维常数向量,其元素包含系数阵中的非随机元素； $M$ 是 $n m \times t$ 维常数矩阵,其形式由系数矩阵中非随机元素的个数及元素间的相关性确定。若 $a$ 和 $y$ 相互独立^[10],即 $c o v (e_{a}, e_{y}) = 0$ ,且两者的方差矩阵为 $D (e_{a}) = σ_{0}^{2} ω^{- 1}$ , $D (e_{y}) = σ_{0}^{2} W^{- 1}$ ,其中 $ω$ 和 $W$ 分别为 $e_{a}$ 和 $e_{y}$ 的权矩阵。相应的平差准则为：

Φ (\bar{a}, x) = e_{a}^{T} ω e_{a} + e_{y}^{T} W e_{y}

(4)

可见PEIV的平差准则中没有考虑独立随机误差的重复次数。如果将 $e_{a}$ 和 $e_{y}$ 组合成独立随机误差向量 $γ = {[\begin{matrix} e_{a}^{T} & e_{y}^{T} \end{matrix}]}^{T}$ ,权矩阵 $P_{γ} = d i a g (ω, W) = (d i a g (ω^{- 1}, W^{- 1} {))}^{- 1}$ ,则平差准则式（4）可以写为：

Φ (\bar{a}, x) = γ^{T} P_{γ} γ

(5)

因此PEIV模型的SWTLS解是令独立随机误差的加权平方和最小的参数估值。文献[12]中的GPEIV模型、文献[13]中约束结构SWTLS方法和文献[14-15]中的变量投影法,其本质上都是采用式（5）所示的平差准则。

2） STLN模型及其平差准则。STLN方法是定义一个与参数 $x$ 相关的矩阵 $X$ ,提取系数矩阵中的独立随机误差 $e_{a}$ ,使下式成立^[17]：

X e_{a} = E_{A} x

(6)

联立式（1a）和式（6）,将观测值残差 $r$ 表达成系数阵独立误差 $e_{a}$ 和参数 $x$ 的函数：

r (e_{a}, x) = - e_{y} = y - A x - X e_{a}

(7)

将式（7）在 $e_{a}$ 和 $x$ 的近似值处线性化,舍去二次项后得到线性模型,然后采用下列平差准则^[17]：

Φ (e_{a}, x) = 0.5 e_{y}^{T} e_{y} + 0.5 e_{a}^{T} D_{a}^{2}

(8)

式中, $D_{a} = d i a g (d_{1}, d_{2}, \dots, d_{t})$ , $d_{i}$ （ $i = 1,2, \dots, t$ ）表示 $e_{a}$ 中的第 $i$ 个元素在误差矩阵 $E_{A}$ 中重复出现的次数。尽管STLN 方法没有考虑误差的权值,但可以很方便地将观测误差的权纳入平差准则。仍令 $γ = {[\begin{matrix} e_{a}^{T} & e_{y}^{T} \end{matrix}]}^{T}$ ,定义 $D_{s} = d i a g (D_{a}, I_{n})$ ,由于式（8）中常数项对求极值无影响,则平差准则式（8）等价于：

Φ (x) = γ^{T} D_{s}^{2} γ

(9)

可见,STLN方法考虑了独立误差 $e_{a}$ 的重复次数,且以重复次数的平方将独立误差纳入平差模型。

3）虚拟误差模型及其平差准则。文献[3]以AR模型参数估计为背景,将式（1a）所示的EIV模型在观测值的近似值处线性化。设系数阵真值 $\tilde{A}$ 的近似值为 $A_{0}$ ,改正数为 $∆ A$ ,参数 $x$ 的近似值为 $x_{0}$ ,改正数为 $∆ x$ ,则线性化方程为：

v = A_{0} x_{0} + A_{0} ∆ x + ∆ A x_{0} - y

(10)

式中, $v$ 表示 $y$ 对应的改正数。然后通过矩阵等价变换得到 $∆ A x_{0} = [\begin{matrix} A_{10} & A_{20} \end{matrix}] [\begin{matrix} x_{B} \\ v \end{matrix}]$ , $A_{10}$ 和 $A_{20}$ 表示由 $x_{0}$ 按照一定规则构造的近似矩阵； $x_{B}$ 表示设计矩阵中独立观测值对应的改正数向量。将式（10）化为：

v = {(I_{n} - A_{20})}^{- 1} A_{0} ∆ x + {(I_{n} - A_{20})}^{- 1} A_{10} x_{B} - {(I_{n} - A_{20})}^{- 1} (y - A_{0} x_{0})

(11)

组成虚拟观测值误差方程为：

v_{B} = x_{B}

(12)

联立式（11）和式（12）得到：

v_{g} = B_{g} x_{g} - l_{g}

(13)

式中, $v_{g}$ 为所有随机观测值改正数向量； $B_{g}$ 为等效设计矩阵； $x_{g}$ 为所有待估参数和虚拟参数的改正数； $l_{g}$ 表示 $y - A_{0} x_{0}$ 的线性变换^[3]。采用如下平差准则：

Φ (x) = v_{g}^{T} P_{g} v_{g}

(14)

式中, $v_{g}$ 实际上是 $(E_{A} | e_{y})$ 中的独立误差向量 $γ$ ；P_g=D_gP由两部分组成,其中 $D_{g} = d i a g (d_{1}, d_{2}, \dots, d_{t + n})$ , $d_{i}$ （ $i = 1,2, \dots, t + n$ ,且t+n为独立误差的个数）表示 $v_{g}$ 中的第 $i$ 个元素在 $(E_{A} | e_{y})$ 中的重复次数, $P$ 为观测值向量的权矩阵。可见,虚拟误差模型解AR模型这一典型SWTLS问题考虑了重复次数。与STLN方法的差异在于,它采用的是重复次数本身,而不是STLN中重复次数的平方。

上述3种典型的STLS方法采用了不同的平差准则,不同的平差准则一定会得到不同的平差结果。为了比较不同准则下的平差结果,下面将上述3种典型的平差准则下的解纳入到一个统一的平差模型中。

2 STLS模型的通用表达式及算法推导

EIV模型（式（1a））又可以表示为^[1]：

y - A x - B e = 0

(15)

式中, $B = [\begin{matrix} x^{T} \otimes I_{n} & - I_{n} \end{matrix}]$ 为 $n \times (n m + n)$ 矩阵,其中 $\otimes$ 为Kronecker积符号,定义为 $G \otimes H = [g_{i j} ∙ H]$ ,其中 $G = [g_{i j}]$ ,且 $H$ 为任意矩阵。设 $γ$ 为 $(E_{A} | e_{y})$ 中的 $t$ 个独立随机误差向量,则 $e_{A}$ 和 $e_{y}$ 可以分别表示为：

e_{A} = H_{1} γ

(16a)

e_{y} = H_{2} γ

(16b)

式中, $H_{1}$ 和 $H_{2}$ 分别为 $n m \times t$ 和 $n \times t$ 矩阵,其元素均为常数。联立式（16a）和式（16b）：

e = [\begin{matrix} e_{A} \\ e_{y} \end{matrix}] = [\begin{matrix} H_{1} \\ H_{2} \end{matrix}] γ = H γ

(17)

将式（17）代入式（15）,并令 $G = B H$ ,则EIV模型可以写为：

y - A x - G γ = 0

(18)

独立随机误差的随机模型为：

E (γ) = 0

D (γ) = σ_{0}^{2} Q_{γ} = σ_{0}^{2} P_{γ}^{- 1}

(19)

式中, $Q_{γ}$ 表示 $γ$ 的协因数矩阵； $P_{γ}$ 表示其权矩阵。令 $γ = {[γ_{1} γ_{2} \dots γ_{t}]}^{T}$ ,且 $γ_{i}$ 出现的次数为 $d_{i}$ （ $i = 1,2, \dots, t$ ）。定义 $D = d i a g (k_{1}, k_{2}, \dots, k_{t})$ ,组成综合权矩阵：

P_{s γ} = D P_{γ}

(20)

相应的平差准则为：

ψ (γ) = γ^{T} P_{s γ} γ = m i n

(21)

易知,当 $k_{i} = 1$ 时,式（21）等价于以PEIV模型为代表的准则式（5）；当 $k_{i} = d_{i}^{2}$ 时,式（21）等价于以STLN模型为代表的准则式（9）；当 $k_{i} = d_{i}$ 时,式（21）等价于以虚拟误差模型为代表的准则式（14）。以结构化函数模型（18）为基础,基于平差准则式（21）,构造如下Lagrange目标函数：

φ (x, γ, λ) = γ^{T} P_{s γ} γ - 2 λ^{T} (y - A x - G γ)

(22)

式中, $λ$ 为 $n$ 维Lagrange乘子向量。令目标函数式（22）对各待定量 $x 、 γ 、 λ$ 的偏导数为0,可得：

{\frac{1}{2} \frac{\partial φ (x, γ, λ)}{\partial x}|}_{\hat{x}, \hat{γ}, \hat{λ}} = {(A + {\hat{E}}_{A})}^{T} \hat{λ} = 0

(23a)

{\frac{1}{2} \frac{\partial φ (x, γ, λ)}{\partial γ}|}_{\hat{x}, \hat{γ}, \hat{λ}} = P_{s γ} \hat{γ} + G^{T} \hat{λ} = 0

(23b)

{\frac{1}{2} \frac{\partial φ (x, γ, λ)}{\partial λ}|}_{\hat{x}, \hat{γ}, \hat{λ}} = y - A \hat{x} - G \hat{γ} = 0

(23c)

由式（23b）得到独立误差向量 $γ$ 的估值为：

\hat{γ} = - P_{s γ}^{- 1} G^{T} \hat{λ} = - Q_{s γ} G^{T} \hat{λ}

(24)

式中, $Q_{s γ} = P_{s γ}^{- 1}$ ,表示综合权矩阵的逆； $\hat{x} 、 \hat{γ} 、 \hat{λ}$ 分别为 $x 、 γ 、 λ$ 的估值。将式（24）代入式（23c）并移项可得：

\hat{λ} = - {(G Q_{s γ} G^{T})}^{- 1} (y - A \hat{x})

(25)

将式（25）代入式（23a）可以得到法方程：

\begin{matrix} (A + {\hat{E}}_{A})^{T} (G Q_{s γ} G^{T})^{- 1} A \hat{x} = \\ (A + {\hat{E}}_{A})^{T} (G Q_{s γ} G^{T})^{- 1} y \end{matrix}

(26)

则参数的估值为：

\hat{x} = [(A + {\hat{E}}_{A})^{T} (G Q_{s γ} G^{T})^{- 1} {A]}^{- 1} (A + {\hat{E}}_{A})^{T} (G Q_{s γ} G^{T})^{- 1} y

(27)

由于式（26）法方程矩阵不对称,若在式（26）两端加上 ${(A + {\hat{E}}_{A})}^{T} {(G Q_{s γ} G^{T})}^{- 1} {\hat{E}}_{A} \hat{x}$ ,那么 $\hat{x}$ 可由下式估计：

\hat{x} = [(A + {\hat{E}}_{A})^{T} (G Q_{s γ} G^{T})^{- 1} (A + {\hat{E}}_{A} {)]}^{- 1} (A + {\hat{E}}_{A})^{T} (G Q_{s γ} G^{T})^{- 1} (y + {\hat{E}}_{A} \hat{x})

(28)

根据上述推导过程,可以得出通用模型计算SWTLS问题的迭代步骤为：

1）给定 $A$ 、 $y$ 、 $Q_{γ}$ ,根据 $γ$ 的重复次数和不同平差准则的计入方式确定对角阵 $D$ ,根据 $(E_{A} | e_{y})$ 的结构确定 $H_{1}$ 和 $H_{2}$ 。采用式（20）计算综合权阵 $P_{s γ}$ 和对应的协因数矩阵 $Q_{s γ} = P_{s γ}^{- 1}$ 。计算 $Q_{y y} = H_{2} Q_{γ} H_{2}^{T}$ 和参数的初值 ${\hat{x}}^{(0)} = {\hat{x}}_{L S} = {(A^{T} Q_{y y}^{- 1} A)}^{- 1} A^{T} Q_{y y}^{- 1} y$ ；

2）根据初值 ${\hat{x}}^{(0)}$ 计算矩阵 $B$ 和 $G$ 。采用式（25）和式（24）分别计算 $\hat{λ}$ 和 $\hat{γ}$ ,然后由式（16a）计算 ${\hat{e}}_{A}$ 的值,从而有 ${\hat{E}}_{A} = i v e c ({\hat{e}}_{A})$ , $i v e c (∙)$ 是 $v e c (∙)$ 的逆运算,表示将 $n m 维$ 向量恢复成 $n \times m$ 维矩阵；

3）将 $G$ 、 $Q_{s γ}$ 、 ${\hat{E}}_{A}$ 、 ${\hat{x}}^{(0)}$ 分别代入式（27）或式（28）,计算参数估值 $\hat{x}$ ；

4）如果最后两次的估值 $\hat{x}$ 足够接近于给定的阈值,终止迭代。否则,转向步骤2）。

3 STLS模型平差准则的优化选取

尽管由式（18）和式（19）组成的SWTLS一般模型能通过选取不同的 $D$ 矩阵获得不同准则下的迭代最优解,但无法获得参数的方差或均方误差（mean square error,MSE）等精度评定指标,只能在模拟实验中检验参数精度。由文献[22-23]可知,WTLS并不是一种新的平差方法,仅仅是LS框架下的另外一种平差模型,且经典平差理论中的Gauss-Helmert（GH）模型（附有参数的条件平差模型）是EIV模型的一个特例。因此WTLS问题的解可由非线性GH模型导出。若给定 $x$ 的近似值 $x_{0}$ 和 $γ$ 的近似值 $γ_{0}$ ,式（18）可线性化为：

f (x, γ) = y - A x - G γ = y - A x_{0} - G_{0} γ_{0} + {\frac{\partial f (x, γ)}{\partial x}|}_{x_{0}, γ_{0}} ∆ x + {\frac{\partial f (x, γ)}{\partial γ}|}_{x_{0}, γ_{0}} ∆ γ

(29)

式中, $Δ x = x - x_{0}$ ； $Δ γ = γ - γ_{0}$ ； $G_{0}$ 和 $E_{A}^{0}$ 为近似值 $x_{0}$ 和 $γ_{0}$ 处的 $G$ 和 $E_{A}$ 的值。且有：

{\frac{\partial f (x, γ)}{\partial x}|}_{x_{0}, γ_{0}} = - (A + E_{A}^{0}) = - (A + i v e c (H_{1} γ_{0}))

(30)

{\frac{\partial f (x, γ)}{\partial γ}|}_{x_{0}, γ_{0}} = - G_{0} = - [x_{0}^{T} \otimes I_{n} - I_{n}] H

(31)

将式（30）和式（31）代入式（29）并作适当变换可以得到：

(A + E_{A}^{0}) x + G_{0} γ - (y + E_{A}^{0} x_{0}) = 0

(32)

可见SEIV模型的线性化形式（式（32））是一个标准GH模型,该模型不再体现 $γ$ 的重复次数,实际上重复次数已经由投影矩阵 $H$ 表达了。重复次数体现了EIV模型的结构性,而投影矩阵 $H$ 正是描述上述结构性的量。根据经典LS理论,其平差准则应为：

Φ (γ, x) = γ^{T} P_{γ} γ

(33)

因此,从SWTLS问题的线性化形式证明了其平差准则应满足式（33）。由线性化模型式（32）得到的最小二乘解是在给定 $x_{0}$ 和 $γ_{0}$ 情形下的一个近似解,为了得到更严密的解,将式（32）进一步表达为：

l = (A + E_{A}^{0}) ∆ x + G_{0} ∆ γ + G_{0} γ_{0}

(34)

式中, $l = y - A x_{0}$ 。为了与§2中符号一致,记 $P_{s γ} = I_{t} P_{γ} = P_{γ}$ （ $I_{t}$ 为t维单位矩阵）。根据式（34）及平差准则式（33）,构造如下Lagrange乘子函数：

m i n : φ (∆ x, ∆ γ, λ) = {(γ_{0} + ∆ γ)}^{T} P_{s γ} (γ_{0} + ∆ γ) + 2 λ^{T} (l - (A + E_{A}^{0}) ∆ x - G_{0} ∆ γ - G_{0} γ_{0})

(35)

分别求 $φ (∆ x, ∆ γ)$ 对 $∆ x$ 、 $∆ γ$ 和 $λ$ 的偏导数,并令其值为0,可得：

\frac{1}{2} \frac{\partial φ (∆ x, ∆ γ)}{\partial ∆ x} = - {(A + E_{A}^{0})}^{T} λ = 0

(36)

\frac{1}{2} \frac{\partial φ (∆ x, ∆ γ)}{\partial ∆ γ} = P_{s γ} (γ_{0} + ∆ γ) - G_{0}^{T} λ = 0

(37)

\begin{matrix} \frac{1}{2} \frac{\partial φ (∆ x, ∆ γ)}{\partial λ} = l - (A + E_{A}^{0}) ∆ x - \\ G_{0} ∆ γ - G_{0} γ_{0} = 0 \end{matrix}

(38)

由式（37）可得：

∆ γ = Q_{s γ} G_{0}^{T} λ - γ_{0}

(39)

将式（39）代入式（38）得：

l - (A + E_{A}^{0}) ∆ x - G_{0} Q_{s γ} G_{0}^{T} λ = 0

(40)

由式（40）可得Lagrange乘子向量为：

λ = {(G_{0} Q_{s γ} G_{0}^{T})}^{- 1} (l - (A + E_{A}^{0}) ∆ x)

(41)

式（40）两边同乘以 ${(A + E_{A}^{0})}^{T}$ 并顾及式（36）,则有：

∆ x = [(A + E_{A}^{0})^{T} (G_{0} Q_{s γ} G_{0}^{T})^{- 1} (A + E_{A}^{0} {)]}^{- 1} (A + E_{A}^{0})^{T} (G_{0} Q_{s γ} G_{0}^{T})^{- 1} (y - A x_{0})

(42)

x = x_{0} + Δ x = [(A + E_{A}^{0})^{T} (G_{0} Q_{s γ} G_{0}^{T})^{- 1} (A + E_{A}^{0} {)]}^{- 1} (A + E_{A}^{0})^{T} (G_{0} Q_{s γ} G_{0}^{T})^{- 1} (y + E_{A}^{0} x_{0})

(43)

由GH模型推导得到的参数估计（式（43））与§2中通用模型得到的解（式（28））具有完全相同的形式。GH模型的算法流程如下：

1）给定 $A$ 、 $y$ 、 $Q_{γ}$ ,根据 $(E_{A} | e_{y})$ 的结构确定 $H_{1}$ 、 $H_{2}$ 。计算 $Q_{y y} = H_{2} Q_{γ} H_{2}^{T}$ 和参数的初值 ${\hat{x}}^{(0)} = {\hat{x}}_{L S} = {(A^{T} Q_{y y}^{- 1} A)}^{- 1} A^{T} Q_{y y}^{- 1} y$ ,使用上标（i） $(i = 0,1, 2, \dots, N, N 为迭代次数)$ 对不同迭代次数下的各参数进行区分,则 $γ$ 和 $E_{A}$ 的初值分别为 $γ^{(0)} = 0$ 和 $E_{A}^{(0)} = i v e c (H_{1} γ^{(0)}) = 0$ ；

2）将 ${\hat{x}}^{(i)}$ 代入式（31）计算 $B^{(i)}$ 和 $G^{(i)}$ ,由式（41）计算 $λ^{(i)}$ ,由式（39）计算 $∆ γ^{(i)}$ 并更新 $γ^{(i)} = γ^{(i - 1)} + ∆ γ^{(i)}$ ,由式（30）计算 $E_{A}^{(i)}$ ；

3）由式（42）计算 $Δ x^{(i - 1)}$ ；

4）如果 $Δ x$ 小于给定的阈值,则终止迭代；否则,更新 $x^{(i)} = Δ x^{(i - 1)} + x^{(i - 1)}$ ,转向步骤2）。

注意,若不限定式（35）中 $P_{s γ} = P_{γ}$ ,那么采用线性化的GH模型,在准则式（21）下,即可选取不同类型 $D$ 矩阵的条件下,可以得到和§2中通用模型完全一致的解。采用GH模型进行算法推导有3个方面的意义：（1）证明了SEIV模型的平差准则中不应该再考虑重复矩阵。（2）证明了由通用模型（式（18））和线性化GH模型作为条件,在相同的平差准则下得到的参数估值是一致的。（3）将SEIV模型线性化后,便于采用协方差传播定律计算参数估值的近似精度。

4 STLS模型参数估值的近似方差

根据§2、§3的算法得到SWTLS解以后,可计算独立误差向量的估值 $\hat{γ}$ ,单位权方差计算式为：

{\hat{σ}}_{0}^{2} = \frac{{\hat{γ}}^{T} P_{s γ} \hat{γ}}{n - m}

(44)

SWTLS解的近似协因数矩阵可以根据线性化GH模型的解（式（42））求解。虚拟观测值 $l$ 可以写成：

l = y - A x_{0} = [\begin{matrix} - x_{0}^{T} \otimes I_{n} & I_{n} \end{matrix}] [\begin{matrix} v e c (A) \\ y \end{matrix}] = - B_{0} H η = - G_{0} η

(45)

式中, $η$ 为 $(A | y)$ 中的 $t$ 个独立随机观测向量。根据协因数传播律,向量 $l$ 的协因数矩阵为：

Q_{l l} = G_{0} Q_{s γ} G_{0}^{T}

(46)

实际计算中近似值 $x_{0}$ 可以用SWTLS估值 $\hat{x}$ 代替,相应的残差矩阵近似值 $E_{A}^{0}$ 用对应的估值 ${\tilde{E}}_{A}$ 代替。如果忽略式（42）中 $\tilde{A} = A + E_{A}^{0}$ 这一项中 $A$ 的随机性,那么可由协因数传播律得到SWTLS解的近似协因数矩阵为：

Q_{\hat{x} \hat{x}} = Q_{Δ x Δ x} = {({\tilde{A}}^{T} Q_{l l}^{- 1} \tilde{A})}^{- 1} {\tilde{A}}^{T} Q_{l l}^{- 1} Q_{l l} Q_{l l}^{- 1} \tilde{A} {({\tilde{A}}^{T} Q_{l l}^{- 1} \tilde{A})}^{- 1} = {({\tilde{A}}^{T} Q_{l l}^{- 1} \tilde{A})}^{- 1}

(47)

从而得到SWTLS解的近似方差为：

D_{\hat{x} \hat{x}} = {\hat{σ}}_{0}^{2} {({\tilde{A}}^{T} Q_{l l}^{- 1} \tilde{A})}^{- 1}

(48)

5 数值实验

首先采用二维仿射变换实例验证算法的正确性和可行性。设第 $i$ 个公共点 $P_{i}$ 在源坐标系和目标坐标系中的坐标分别为 $(x_{i}, y_{i})$ 和 $(X_{i}, Y_{i})$ ,仿射变换模型为 $X_{i} = a_{0} + a_{1} x_{i} + a_{2} y_{i}$ , $Y_{i} = b_{0} + b_{1} x_{i} + b_{2} y_{i}$ ,其中 $a_{0}$ 和 $b_{0}$ 为两个坐标系的平移量, $a_{1}$ 、 $a_{2}$ 、 $b_{1}$ 、 $b_{2}$ 分别表示两个坐标系之间的旋转和尺度变换参数。假设有3个或以上的公共点,仿射变换模型可以用EIV模型（式（1））表示,且有：

A = [\begin{matrix} 1 & x_{1} & y_{1} & 0 & 0 & 0 \\ 0 & 0 & 0 & 1 & x_{1} & y_{1} \\ 1 & x_{2} & y_{2} & 0 & 0 & 0 \\ 0 & 0 & 0 & 1 & x_{2} & y_{2} \\ ⋮ & ⋮ & ⋮ & ⋮ & ⋮ & ⋮ \\ 1 & x_{n} & y_{n} & 0 & 0 & 0 \\ 0 & 0 & 0 & 1 & x_{n} & y_{n} \end{matrix}]

x = [\begin{matrix} a_{0} \\ a_{1} \\ a_{2} \\ b_{0} \\ b_{1} \\ b_{2} \end{matrix}]

(49)

给定12个公共点在源坐标系中的坐标分别为（-1,1）、（1.2,-3.0）、（-2.6,3.0）、（3.0,1.5）、（-4.8,-1.0）、（5.2,0.2）、（6.0,5.5）、（-7.2,2.2）、（7.8,-2.0）、（8.5,2.2）、（-9.5,-5.0）、（10.0,-0.8）,参数的真值为 $\tilde{x}$ =[10 4 -2 -10 1 3]^T,根据仿射变换模型计算出12个点在目标系中的真实坐标分别为（4,-8）、（20.8,-17.8）、（-6.4,-3.6）、（19.0,-2.5）、（-7.2,-17.8）、（30.4,-4.2）、（23.0,12.5）、（-23.2,-10.6）、（45.2,-8.2）、（39.6,5.1）、（-18.0,-34.5）、（51.6,-2.4）。对所有坐标值添加方差为 $D_{48} = σ_{0}^{2} I_{48} = I_{48}$ （ $I_{48}$ 为 $48 \times 48$ 的单位矩阵）的随机误差。首先,不考虑系数矩阵中源坐标重复2次而目标坐标重复1次的差异,取重复矩阵 $D = D_{1} = I_{48}$ 时,分别采用本文方法、PEIV模型^[10]和变量投影法^[14]进行求解,得到的参数估计结果如表1所示。

表 1 不同平差方法的SWTLS解（

D = D_{1}

）

Table 1. SWTLS Solutions with Different Adjustment Methods （

D = D_{1}

）

平差方法	${\hat{a}}_{0}$	${\hat{a}}_{1}$	${\hat{a}}_{2}$	${\hat{b}}_{0}$	${\hat{b}}_{1}$	${\hat{b}}_{2}$
PEIV模型^[10]	10.412 055	4.065 602	-2.136 739	-9.641 747	1.116 170	2.912 190
变量投影法^[14]	10.412 055	4.065 602	-2.136 739	-9.641 747	1.116 170	2.912 190
本文方法	10.412 055	4.065 602	-2.136 739	-9.641 747	1.116 170	2.912 190

下载: 导出CSV

| 显示表格

当取 $D = D_{2} = d i a g (4, \dots, 4,1, \dots, 1)$ （4和1的个数都为24）时,采用本文方法和STLN方法^[17]进行求解,得到的参数估计结果如表2所示。

表 2 不同平差方法的SWTLS解（

D = D_{2}

）

Table 2. SWTLS Solutions with Different Adjustment Methods （

D = D_{2}

）

平差方法	${\hat{a}}_{0}$	${\hat{a}}_{1}$	${\hat{a}}_{2}$	${\hat{b}}_{0}$	${\hat{b}}_{1}$	${\hat{b}}_{2}$
STLN方法^[17]	10.418 681	4.056 913	-2.101 082	-9.639 440	1.118 336	2.855 364
本文方法	10.418 681	4.056 913	-2.101 082	-9.639 440	1.118 336	2.855 364

下载: 导出CSV

| 显示表格

从表1可以看出,由于PEIV模型^[10]、变量投影法^[14]均未考虑源系统中坐标的重复次数,两者的平差准则等价于本文提出的通用模型中重复矩阵取单位矩阵时的准则（式（21））,尽管三者对结构误差的处理采用了不同的函数模型,但在相同的准则下都得到了完全一致的平差结果,说明本文提出的通用平差算法是可行有效的。由表2可以看出,STLN方法^[17]顾及了系数阵中随机观测值的重复次数,并且是将观测重复数的平方纳入平差模型,和本文通用模型在同样的平差准则下也得到了一致的平差结果,进一步证明了本文算法的可行性和有效性。

为了验证不同平差准则下SWTLS解的统计性质,模拟系数阵元素和右端项含有公共元素的SEIV模型,系数阵和右端项的真值及误差的结构分别为：

(\tilde{A} | \tilde{y}) = [\begin{matrix} \begin{array}{l} 3.62 \\ - 1.03 \\ 0.75 \\ 4.05 \\ 0.75 \\ - 4.59 \\ - 4.16 \\ 1.25 \\ 2.75 \\ 2.50 \\ - 38.88 \\ 1.55 \\ - 18.89 \\ - 4.50 \\ 5.45 \\ 5.28 \\ 8.25 \\ - 10.85 \\ 4.55 \\ - 16.70 \\ 1.55 \\ 3.20 \\ 3.17 \\ 7.95 \\ 6.00 \end{array} & \begin{array}{l} 1.25 \\ 2.50 \\ 0.20 \\ - 1.72 \\ - 3.21 \\ 0.85 \\ 0.21 \\ 0.55 \\ 0.48 \\ - 3.00 \\ 8.42 \\ 0.85 \\ 17.50 \\ 1.55 \\ - 6.18 \\ 2.19 \\ - 8.01 \\ 4.55 \\ - 2.21 \\ 0.12 \\ 0 \\ 1.32 \\ 0.28 \\ 4.85 \\ - 8.72 \end{array} & \begin{array}{l} 6.75 \\ - 5.36 \\ 2.50 \\ 2.90 \\ 6.75 \\ 1.42 \\ 1.98 \\ 5.50 \\ - 1.80 \\ 15.00 \\ 3.07 \\ 7.20 \\ - 7.93 \\ 2.50 \\ 15.00 \\ 4.06 \\ 17.50 \\ 0.45 \\ 6.25 \\ 24.35 \\ 3.20 \\ 2.49 \\ 0.65 \\ - 10.65 \\ 32.60 \end{array} \end{matrix} |\begin{array}{l} 23.37 \\ 0.75 \\ 6.75 \\ 1.25 \\ - 1.80 \\ 2.50 \\ 0.85 \\ 15.00 \\ 1.55 \\ 17.50 \\ 9.36 \\ 20.20 \\ 52.75 \\ 8.25 \\ 4.55 \\ 24.35 \\ 3.20 \\ 12.80 \\ 6.00 \\ 32.60 \\ 7.95 \\ 14.78 \\ 5.87 \\ 10.90 \\ 27.60 \end{array}]

(50)

(E_{A} | e_{y}) = [\begin{matrix} \begin{matrix} 0 & γ_{4} & γ_{3} \\ 0 & γ_{6} & 0 \\ γ_{2} & 0 & γ_{6} \\ 0 & 0 & 0 \\ γ_{2} & 0 & γ_{3} \\ 0 & γ_{7} & 0 \\ 0 & 0 & 0 \\ γ_{4} & 0 & 0 \\ 0 & 0 & γ_{5} \\ γ_{6} & 0 & γ_{8} \\ 0 & 0 & 0 \\ γ_{9} & γ_{7} & 0 \\ 0 & γ_{10} & 0 \\ 0 & γ_{9} & γ_{6} \\ 0 & 0 & γ_{8} \\ 0 & 0 & 0 \\ γ_{14} & 0 & γ_{10} \\ 0 & γ_{15} & 0 \\ γ_{15} & 0 & 0 \\ 0 & 0 & γ_{16} \\ γ_{9} & 0 & γ_{17} \\ γ_{17} & 0 & 0 \\ 0 & 0 & 0 \\ γ_{21} & 0 & 0 \\ γ_{19} & 0 & γ_{20} \end{matrix} & |\begin{matrix} γ_{1} \\ γ_{2} \\ γ_{3} \\ γ_{4} \\ γ_{5} \\ γ_{6} \\ γ_{7} \\ γ_{8} \\ γ_{9} \\ γ_{10} \\ γ_{11} \\ γ_{12} \\ γ_{13} \\ γ_{14} \\ γ_{15} \\ γ_{16} \\ γ_{17} \\ γ_{18} \\ γ_{19} \\ γ_{20} \\ γ_{21} \\ γ_{22} \\ γ_{23} \\ γ_{24} \\ γ_{25} \end{matrix} \end{matrix}]

(51)

待估参数的真值为 $\tilde{x} = {[\begin{matrix} 1 & 5 & 2 \end{matrix}]}^{T}$ ,独立误差向量 $γ = [γ_{1} γ_{2} \dots γ_{25}]^{T}$ 是零均值独立同分布的随机变量,其方差为 $D_{γ} = σ_{0}^{2} I_{25}$ ,其中 $I_{25}$ 表示25×25阶单位矩阵, $γ_{i}$ （ $i = 1,2, \dots, 25$ ）的重复次数 $d_{i} 可以$ 由 $(E_{A} | e_{y})$ 的结构获得。分别采用3种平差准则计算参数的估值,准则式（21）中对应的 $D$ 矩阵分别取 $D_{1} = I_{25}$ , $D_{2} = d i a g (d_{1}, d_{2}, \dots, d_{25}) = d i a g (1,3, 1,3, \dots, 1,3, 1)$ ,D₃=diag（ $d_{1}^{2}$ , $d_{2}^{2}$ , $\dots$ , $d_{25}^{2}$ ）=diag（1,9,1,9, $\dots$ ,1,9,1）。根据式（44）和式（48）计算单位权方差和参数估值的方差。由于参数真值已知,可以计算均方误差 $M S E (\hat{x}) = {(\hat{x} - \tilde{x})}^{T} (\hat{x} - \tilde{x})$ ,以此来检核3种平差准则下估值的精度,并与近似方差对照。分别取误差水平 $σ_{0}^{2} = 0.25$ 和 $σ_{0}^{2} = 1$ 进行计算。为了体现结果的统计性,将不同误差水平的实验各进行10 000次,将每次运算的结果取平均值,结果分别见表3和表4。

表 3 不同平差准则下的SWTLS结果（

σ_{0}^{2} = 0.25

）

Table 3. SWTLS Results Under Different Adjustment Principles （

σ_{0}^{2} = 0.25

）

统计项	统计量	$D = D_{1}$	$D = D_{2}$	$D = D_{3}$
估计参数值	${\hat{x}}_{1}$	0.999 998	1.000 281	0.999 975
	${\hat{x}}_{2}$	5.002 389	5.007 359	5.009 953
	${\hat{x}}_{3}$	1.999 596	1.999 320	1.998 190
单位权方差	${\hat{σ}}_{0}^{2}$	0.239 737	0.537 349	1.463 591
参数方差估值	${\hat{σ}}_{x_{1}}^{2}$	9.43×10^-6	1.72×10^-5	5.76×10^-5
	${\hat{σ}}_{x_{2}}^{2}$	5.23×10^-4	9.14×10^-4	3.19×10^-3
	${\hat{σ}}_{x_{3}}^{2}$	8.28×10^-6	1.60×10^-5	5.05×10^-5
	${\hat{σ}}_{Σ}^{2}$	5.40×10^-4	9.47×10^-4	3.30×10^-3
MSE	${\hat{m}}_{x_{1}}^{2}$	9.76×10^-6	1.10×10^-5	1.63×10^-5
	${\hat{m}}_{x_{2}}^{2}$	5.44×10^-4	6.21×10^-4	9.40×10^-4
	${\hat{m}}_{x_{3}}^{2}$	8.50×10^-6	9.39×10^-6	1.42×10^-5
	${\hat{m}}_{Σ}^{2}$	5.63×10^-4	6.41×10^-4	9.71×10^-4

下载: 导出CSV

| 显示表格

表 4 不同平差准则下的SWTLS结果（

σ_{0}^{2} = 1

）

Table 4. SWTLS Results Under Different Adjustment Principles （

σ_{0}^{2} = 1

）

统计项	统计量	$D = D_{1}$	$D = D_{2}$	$D = D_{3}$
估计参数值	${\hat{x}}_{1}$	1.000 220	1.003 770	1.007 659
	${\hat{x}}_{2}$	4.997 838	5.022 963	5.049 989
	${\hat{x}}_{3}$	2.000 725	2.004 000	2.007 714
单位权方差	${\hat{σ}}_{0}^{2}$	0.960 341	2.146 473	5.839 121
参数方差估值	${\hat{σ}}_{x_{1}}^{2}$	3.78×10^-5	6.87×10^-5	2.30×10^-4
	${\hat{σ}}_{x_{2}}^{2}$	2.09×10^-3	3.65×10^-3	1.27×10^-2
	${\hat{σ}}_{x_{3}}^{2}$	3.32×10^-5	6.41×10^-5	2.02×10^-4
	${\hat{σ}}_{Σ}^{2}$	2.17×10^-3	3.78×10^-3	1.32×10^-2
MSE	${\hat{m}}_{x_{1}}^{2}$	4.02×10^-5	4.50×10^-5	6.67×10^-5
	${\hat{m}}_{x_{2}}^{2}$	2.21×10^-3	2.53×10^-3	3.85×10^-3
	${\hat{m}}_{x_{3}}^{2}$	3.53×10^-5	3.88×10^-5	5.82×10^-5
	${\hat{m}}_{Σ}^{2}$	2.29×10^-3	2.62×10^-3	3.97×10^-3

下载: 导出CSV

| 显示表格

从表3可以看出,不考虑重复次数的平差准则（采用 $D_{1}$ ）获得的参数估值更接近于真值。参数各分量的MSE最小,从数值上验证了不考虑系数阵误差重复系数的平差准则是最优的。单位权方差的估值略小于真值,是由于式（43）没有考虑非线性模型线性化引入的偏差项。同理,各参数分量的方差均小于对应的MSE,原因是SEIV模型的非线性特性会引入估计偏差,方差加上偏差的平方和等于MSE。但是,对比方差和MSE的差异可知,参数3个分量的方差估值和MSE之间的差异（ $|({\hat{σ}}_{x_{i}}^{2} - {\hat{m}}_{x_{i}}^{2}) / {\hat{m}}_{x_{i}}^{2}|$ ）分别为3.38%、3.86%和2.59%,说明当误差较小时,参数方差估值和MSE的差异很小,可以作为精度评定的指标。因此,本文不再采用非线性最小二乘偏差修正或Monte Carlo模拟的思想求解偏差值^[10,24]。

从上述理论分析可知,将SEIV模型线性化后,其平差准则不应顾及误差重复次数。式（44）和式（48）只适用于 $D$ 为单位阵的情形。从后面两列可以看出,当平差准则中计入误差重复的次数时,单位权方差估值与真值以及参数分量的方差估值与MSE均有较大的差异。当 $D = D_{3}$ 时,单位权方差估值约为真值的6倍,方差估值约为对应MSE的3.5倍,这是由于单位权方差公式中重复计算残差平方和所致,由此进一步证明不能采用这两种平差准则用线性近似方法求单位权中误差及方差。

由表4可知,增大观测误差的方差,能够得到与表3一致的结论,进一步验证了最优平差准则应该选择重复矩阵 $D$ 为单位阵。由式（48）给出的方差估值与MSE在3个分量上的偏差分别为5.97%、5.43%、5.95%,进一步证明了本文的近似精度评定方法在误差较小的情况下是可行有效的。

6 结语

结构EIV模型系数矩阵中的随机元素重复出现的次数是否应计入平差准则以及如何计入平差准则,目前尚未形成定论。本文从模型分析和数值验证两方面入手,证明重复次数不应计入平差准则。主要贡献如下：

1）总结了已有的3种处理SEIV模型的平差准则,指出不同的平差准则会得到不同的平差结果。提出了一种通用的SWTLS平差模型,通过选取不同的综合权矩阵 $P_{s γ}$ 得到的目标函数等价于上述3种不同准则下的目标函数。采用Lagrange乘子法推导了通用模型的解并给出了计算步骤。

2）迭代方法无法给出参数的统计性质,本文将通用模型线性化得到GH模型,从理论上分析了误差重复次数不应计入平差准则的原因。推导了GH模型的算法,证明其与通用模型得到的结果是一致的。根据误差传播律得到了参数的近似方差估值。

3）通过实例验证了本文提出的算法与已有方法结果一致,证明了本文方法可行有效。通过模拟计算证明了不考虑误差重复次数的平差准则得到的解在MSE意义下最优,且参数的近似方差是MSE的良好近似,可以作为精度评定的指标。

致谢: 感谢ESA提供Sentinel-1A卫星数据。

图 1 2017年九寨沟地震构造背景图

Figure 1. Tectonic Setting of the 2017 Jiuzhaigou Earthquake

下载: 全尺寸图片幻灯片

图 2 九寨沟地震InSAR同震形变场

Figure 2. The InSAR Coseismic Interferograms of Jiuzhaigou Earthquake

下载: 全尺寸图片幻灯片

图 3 基于InSAR数据反演的同震滑动分布模型

Figure 3. Fault Slip Distribution from Inversions Constrained by InSAR

下载: 全尺寸图片幻灯片

图 4 2017年九寨沟地震InSAR观测与三维形变解算结果

Figure 4. InSAR Observations and 3D Displacement Result of 2017 Jiuzhaigou Earthquake

下载: 全尺寸图片幻灯片

图 5 九寨沟地震三维形变剖面图

Figure 5. 3D Displacement Profiles of 2017 Jiuzhaigou Earthquake

下载: 全尺寸图片幻灯片

图 6 GPS观测与解算结果一致性比较

Figure 6. Comparisons Between GNSS Observation and Displacement Results Based on InSAR and Fault Model

下载: 全尺寸图片幻灯片

表 1 Sentinel-1A卫星数据基本信息

Table 1 Details of Sentinel-1A Satellite Images

模式	获取时间		时间基线/d	垂直基线/m
模式	主影像	辅影像	时间基线/d	垂直基线/m
升轨	2017-07-30	2017-08-11	12	38.36
降轨	2017-08-06	2017-08-18	12	－70.64

下载: 导出CSV

表 2 不同机构或研究团队发布的九寨沟地震

Table 2 Focal Mechanism Solutions of 2017 Jiuzhaigou

来源	矩震级	地震矩/Nm	走向/(°)	倾角/(°)	滑动角/(°)
USGS	6.51	7.23×10¹⁸	153	84	－33
GCMT	6.5	6.98×10¹⁸	151	79	－8
文献[36]	6.5	6.61×10¹⁸	153	84	－19.5
文献[37]	6.5		153	50	－9
文献[38]	6.5		152	74	8
文献[39]	6.49	6.635×10¹⁸	155	81	－11
本文	6.5	7.025×10¹⁸	153	75	－8.3
注：USGS为美国地质勘探局（United States Geological Survey）；GCMT为哈佛大学全球矩张量解（global centroid moment tensor）

下载: 导出CSV

表 3 GNSS观测到的同震位移^[51]

Table 3 Coseismic Displacements Observed from GNSS^[51]

台站名	经度	纬度	东西向形变/cm	南北向形变/cm
JB33	103.89°E	33.28°N	3.0	9.0
J416	103.80°E	33.32°N	9.0	22.0
H025	103.43°E	32.93°N	0.5	1.0
W014	103.67°E	33.56°N	－0.9	－0.7
W015	103.97°E	33.53°N	2.2	0.4
W037	103.71°E	33.03°N	－4.1	－6.6
SB17	103.65°E	32.81°N	－2.5	－3.0

下载: 导出CSV

参考文献(51)

[1]	Hu J, Li Z W, Ding X L, et al. Resolving Three-Dimensional Surface Displacements from InSAR Measurements: A Review[J]. Earth-Science Reviews, 2014, 133: 1-17 doi: 10.1016/j.earscirev.2014.02.005
[2]	Wright T J, Parsons B E, Lu Z. Toward Mapping Surface Deformation in Three Dimensions Using InSAR[J]. Geophysical Research Letters, 2004, 31(1): L01607
[3]	王永哲, 李志伟, 朱建军, 等. 融合多平台DInSAR数据解算拉奎拉地震三维同震形变场[J]. 武汉大学学报·信息科学版, 2012, 37(7): 859-863 doi: 10.13203/j.whugis2012.07.027 Wang Yongzhe, Li Zhiwei, Zhu Jianjun, et al. Coseismic Three-Dimensional Deformation of L'Aquila Earthquake Derived from Multi-platform DInSAR Data[J]. Geomatics and Information Science of Wuhan University, 2012, 37(7): 859-863 doi: 10.13203/j.whugis2012.07.027
[4]	Fialko Y, Simons M, Agnew D. The Complete (3D) Surface Displacement Field in the Epicentral Area of the 1999 Mw 7.1 Hector Mine Earthquake, California, from Space Geodetic Observations[J]. Geophysical Research Letters, 2001, 28(16): 3063-3066 doi: 10.1029/2001GL013174
[5]	Fialko Y, Sandwell D, Simons M, et al. Three-Dimensional Deformation Caused by the Bam, Iran, Earthquake and the Origin of Shallow Slip Deficit[J]. Nature, 2005, 435(7040): 295-299 doi: 10.1038/nature03425
[6]	Hu J, Li Z W, Ding X L, et al. Two-Dimensional Co-seismic Surface Displacements Field of the Chi-Chi Earthquake Inferred from SAR Image Matching[J]. Sensors, 2008, 8(10): 6484-6495 doi: 10.3390/s8106484
[7]	Werner C, Wegmuller U, Strozzi T, et al. Precision Estimation of Local Offsets Between Pairs of SAR SLCs and detected SAR Images[C]//International Geoscience and Remote Sensing Symposium, Seoul, Korea, 2005
[8]	Bechor N B D, Zebker H A. Measuring Two-Dimensional Movements Using a Single InSAR Pair[J]. Geophysical Research Letters, 2006, 33(16): L16311 doi: 10.1029/2006GL026883
[9]	杨红磊, 彭军还. 基于DInSAR和MAI技术揭示地震三维形变场[J]. 地球物理学进展, 2014, 29(6): 2580-2586 https://www.cnki.com.cn/Article/CJFDTOTAL-DQWJ201406016.htm Yang Honglei, Peng Junhuan. Mapping Three-Dimensional Co-seismic Deformation by Combining Multiple-Aperture Interferometry and Differential Interferometry Synthetic Aperture Radar[J]. Progress in Geophysics, 2014, 29(6): 2580-2586 https://www.cnki.com.cn/Article/CJFDTOTAL-DQWJ201406016.htm
[10]	Jung H S, Yun S H, Jo M J. An Improvement of Multiple-Aperture SAR Interferometry Performance in the Presence of Complex and Large Line-of-Sight Deformation[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2015, 8(4): 1743-1752 doi: 10.1109/JSTARS.2015.2399249
[11]	Gudmundsson S, Sigmundsson F, Carstensen J M. Three-Dimensional Surface Motion Maps Estimated from Combined Interferometric Synthetic Aperture Radar and GPS Data[J]. Journal of Geophysical Research: Solid Earth, 2002, 107(B10): 1-14
[12]	刘斌. InSAR高精度观测地震形变场及其三维重建技术研究[D]. 哈尔滨: 中国地震局工程力学研究所, 2013 Liu Bin. Method for Monitoring Displacement Field of Earthquake with High-Precision and Constructing 3D Coseismic Maps Using InSAR[D]. Harbin: Institute of Engineering Mechanics, China Earthquake Administrtion, 2013
[13]	曹海坤, 赵丽华, 张勤, 等. 利用附加系统误差参数的升降轨InSAR-GPS数据融合方法建立三维形变场[J]. 武汉大学学报·信息科学版, 2018, 43(9): 1362-1368 doi: 10.13203/j.whugis20160461 Cao Haikun, Zhao Lihua, Zhang Qin, et al. Ascending and Descending Orbits InSAR-GPS Data Fusion Method with Additional Systematic Parameters for Three-Dimensional Deformaion Field[J]. Geomatics and Information Science of Wuhan University, 2018, 43(9): 1362-1368 doi: 10.13203/j.whugis20160461
[14]	孙建宝, 梁芳, 徐锡伟, 等. 升降轨道ASAR雷达干涉揭示的巴姆地震(MW6.5)3D同震形变场[J]. 遥感学报, 2006, 10(4): 59-66 https://www.cnki.com.cn/Article/CJFDTOTAL-YGXB200604008.htm Sun Jianbao, Liang Fang, Xu Xiwei, et al. 3D Co-seismic Deformation Field of the Bam Earthquake(Mw 6.5) from Ascending and Descending Pass ASAR Radar Interferometry[J]. Journal of Remote Sensing, 2006, 10(4): 59-66 https://www.cnki.com.cn/Article/CJFDTOTAL-YGXB200604008.htm
[15]	温扬茂, 冯怡婷. 地震破裂模型约束的中国阿里地震三维形变场[J]. 武汉大学学报·信息科学版, 2018, 43(9): 1369-1375 doi: 10.13203/j.whugis20160450 Wen Yangmao, Feng Yiting. Three-Dimensional Deformation Field of Ali Earthquake from InSAR Observations and Earthquake Rupture Model[J]. Geomatics and Information Science of Wuhan University, 2018, 43(9): 1369-1375 doi: 10.13203/j.whugis20160450
[16]	Song X, Jiang Y, Shan X, et al. Deriving 3D Coseismic Deformation Field by Combining GPS and InSAR Data Based on the Elastic Dislocation Model[J]. International Journal of Applied Earth Observation and Geoinformation, 2017, 57: 104-112 doi: 10.1016/j.jag.2016.12.019
[17]	Gan W, Zhang P, Shen Z K, et al. Present-day Crustal Motion Within the Tibetan Plateau Inferred from GPS Measurements[J]. Journal of Geophysical Research Solid Earth, 2007, 112(B8): B08416
[18]	陈长云, 任金卫, 孟国杰, 等. 巴颜喀拉块体北东缘主要断裂现今活动性分析[J]. 大地测量与地球动力学, 2012, 32(3): 27-30 doi: 10.3969/j.issn.1671-5942.2012.03.005 Chen Changyun, Ren Jinwei, Meng Guojie, et al. Analysis of Modern Activity of Major Faults in Northeast Margin of Baryan-Har Block[J]. Journal of Geodesy and Geodynamics, 2012, 32(3): 27-30 doi: 10.3969/j.issn.1671-5942.2012.03.005
[19]	董秀军, 许强, 佘金星, 等. 九寨沟核心景区多源遥感数据地质灾害解译初探[J]. 武汉大学学报·信息科学版, 2020, 45(3): 119-128 doi: 10.13203/j.whugis20190076 Dong Xiujun, Xu Qiang, She Jinxing, et al. Preliminary Study on Interpretation of Geological Hazards in Jiuzhaigou Based on Multi-source Remote Sensing Data[J]. Geomatics and Information Science of Wuhan University, 2020, 45(3): 119-128 doi: 10.13203/j.whugis20190076
[20]	冯万鹏, 李振洪. InSAR资料约束下震源参数的PSO混合算法反演策略[J]. 地球物理学进展, 2010, 25(4): 1189-1196 doi: 10.3969/j.issn.1004-2903.2010.04.007 Feng Wanpeng, Li Zhenhong. A Nover Hybrid PSO/Simplex Algorithm for Determining Earthquake Source Parameters Using InSAR Data[J]. Progress in Geophys, 2010, 25(4): 1189-1196 doi: 10.3969/j.issn.1004-2903.2010.04.007
[21]	温扬茂, 何平, 许才军, 等. 联合Envisat和ALOS卫星影像确定L'Aquila地震震源机制[J]. 地球物理学报, 2012, 55(1): 53-65 doi: 10.6038/j.issn.0001-5733.2012.01.006 Wen Yangmao, He Ping, Xu Caijun, et al. Source Parameters of the 2009 L'Aquila Earthquake, Italy from Envisat and ALOS Satellite SAR Images[J]. Chinese Journal of Geophysics, 2012, 55(1): 53-65 doi: 10.6038/j.issn.0001-5733.2012.01.006
[22]	王振杰, 欧吉坤, 柳林涛. 一种解算病态问题的方法: 两步解法[J]. 武汉大学学报·信息科学版, 2005, 30(9): 821-824 http://ch.whu.edu.cn/article/id/2283 Wang Zhenjie, Ou Jikun, Liu Lintao. A Method for Solving Ill-Posed Problems: Two Step Method[J]. Geomatics and Information Science of Wuhan University, 2005, 30(9): 821-824 http://ch.whu.edu.cn/article/id/2283
[23]	王乐洋, 赵雄, 高华. 大地测量地震断层同震滑动分布反演的两步解法[J]. 武汉大学学报·信息科学版, 2019, 44(9): 1265-1311 doi: 10.13203/j.whugis20170382 Wang Leyang, Zhao Xiong, Gao Hua. A Two-Step Solution Method for the Co-seismic Slip Distribution Inversion of Earthquake Fault in Geodesy[J]. Geomatics and Information Science of Wuhan University, 2019, 44(9): 1265-1311 doi: 10.13203/j.whugis20170382
[24]	Werner C, Wegmuller U, Strozzi T, et al. Gamma SAR Processor and Interferometry Software[C]//The ERS-ENVISAT Symposium, Gothenburg, Sweden, 2000
[25]	廖明生, 林珲. 雷达干涉测量: 原理与信号处理基础[M]. 北京: 测绘出版社, 2003 Liao Mingsheng, Lin Hui. Synthetic Aperture Radar Interferometry: Principle and Signal Processing[M]. Beijing: Sino Maps Press, 2003
[26]	Farr M, Rosen A, Caro E, et al. The Shuttle Radar Topography Mission[J]. Reviews of Geophysics, 2007, 1(2): 361-369
[27]	Goldstein R M, Werner C L. Radar Interferogram Filtering for Geophysical Applications[J]. Geophysical Research Letters, 1998, 25(21): 4035-4038 doi: 10.1029/1998GL900033
[28]	Carballo G F, Fieguth P W. Probabilistic Cost Functions for Network Flow Phase Unwrapping[J]. IEEE Transactions on Geoscience and Remote Sensing, 2000, 38(5): 2192-2201 doi: 10.1109/36.868877
[29]	于勇. 基于网络规划的干涉雷达相位解缠算法研究[D]. 北京: 中国科学院, 2002 Yu Yong. Study on Phase Unwrapping Algorithm based on Network Plan[D]. Beijing: Chinese Academy of Sciences, 2002
[30]	Yu C, Li Z H, Penna N T. Interferometric Synthetic Aperture Radar Atmospheric Correction Using a GPS-Based Iterative Tropospheric Decomposition Model[J]. Remote Sensing of Environment, 2018, 204(1): 109-121
[31]	李鹏, 高梦瑶, 李振洪, 等. 阿尔金断裂带宽幅InSAR对流层延迟估计方法评估[J]. 武汉大学学报·信息科学版, 2020, 45(6): 879-887 doi: 10.13203/j.whugis20190236 Li Peng, Gao Mengyao, Li Zhenhong, et al. Evaluation of Wide-Swath InSAR Tropospheric Delay Estimation Methods over the Altyn Tagth Fault[J]. Geomatics and Information Science of Wuhan University, 2020, 45(6): 879-887 doi: 10.13203/j.whugis20190236
[32]	Okada Y. Surface Deformation to Shear and Tensile Faults in a Half-Space[J]. Bulletin of the Seismological Society of America, 1985, 75(4): 1135-1154 doi: 10.1785/BSSA0750041135
[33]	Okada Y. Internal Deformation due to Shear and Tensile Faults in a Half-Space[J]. Bulletin of the Seismological Society of America, 1992, 82(2): 1018-1040 doi: 10.1785/BSSA0820021018
[34]	Lohman R B, Simons M. Some Thoughts on the Use of InSAR Data to Constrain Models of Surface Deformation: Noise Structure and Data Down Sampling[J]. Geochemistry, Geophysics, Geosystems, 2005, 6(1): 359-361
[35]	Feng W P, Li Z H, Elliott J R, et al. The 2011 Mw 6.8 Burma Earthquake: Fault Constraints Provided by Multiple SAR Techniques[J]. Geophysical Journal International, 2013, 195(1): 650-660 doi: 10.1093/gji/ggt254
[36]	张旭, 冯万鹏, 许力生, 等. 2017年九寨沟Ms 7.0级地震震源反演过程与烈度估计[J]. 地球物理学报, 2017, 60(10): 4105-4116 doi: 10.6038/cjg20171035 Zhang Xu, Feng Wanpeng, Xu Lisheng, et al. The Source-Process Inversion and the Intensity Estimation of the 2017 Ms 7.0 Jiuzhaigou Earthquake[J]. Chinese Journal of Geophysics, 2017, 60(10): 4105-4116 doi: 10.6038/cjg20171035
[37]	单新建, 屈春燕, 龚文瑜, 等. 2017年8月8日四川九寨沟7.0级地震InSAR同震形变场及断层滑动分布反演[J]. 地球物理学报, 2017, 60(12): 4527-4536 doi: 10.6038/cjg20171201 Shan Xinjian, Qu Chunyan, Gong Wenyu, et al. Coseismic Deformation Field of the Jiuzhaigou Ms 7.0 Earthquake from Sentinel-1A InSAR Data and Fault Slip Inversion[J]. Chinese Journal of Geophysics, 2017, 60(12): 4527-4536 doi: 10.6038/cjg20171201
[38]	谢祖军, 郑勇, 姚华建, 等. 2017年九寨沟Ms7.0地震震源性质及发震构造初步分析[J]. 中国科学: 地球科学, 2018, 48(1): 79-92 https://www.cnki.com.cn/Article/CJFDTOTAL-JDXK201801006.htm Xie Zujun, Zheng Yong, Yao Huajian, et al. Preliminary Analysis on the Source Properties and Seismogenic of the 2017 Ms 7.0 Jiuzhaigou Earthquake[J]. Science China: Earth Sciences, 2018, 48(1): 79-92 https://www.cnki.com.cn/Article/CJFDTOTAL-JDXK201801006.htm
[39]	Nie Z, Wang D J, Jia Z, et al. Fault Model of the 2017 Jiuzhaigou Mw 6.5 Earthquake Estimated from Coseismic Deformation Observed Using Global Positioning System and Interferometric Synthetic Aperture Radar Data[J]. Earth, Planets and Space, 2018, 70(1): 55 doi: 10.1186/s40623-018-0826-4
[40]	Burgmann R. Deformation During the 12 November 1999 Duzce, Turkey, Earthquake, from GPS and InSAR Data[J]. Bulletin of the Seismological Society of America, 2002, 92(1): 161-171 doi: 10.1785/0120000834
[41]	温扬茂, 许才军, 刘洋, 等. 升降轨InSAR数据约束下的2007年阿里地震反演分析[J]. 测绘学报, 2015, 44(6): 649-654 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201506009.htm Wen Yangmao, Xu Caijun, Liu Yang, et al. The 2007 Ali Earthquake Inversion from Ascending and Descending InSAR Observations[J]. Acta Geodaetica et Cartographica Sinica, 2015, 44(6): 649-654 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201506009.htm
[42]	崔希璋. 广义测量平差[M]. 武汉: 武汉大学出版社, 2006 Cui Xizhang. Generalized Surveying Adjustment[M]. Wuhan: Wuhan University Press, 2006
[43]	Hu J, Li Z W, Sun Q, et al. Three-Dimensional Surface Displacements from InSAR and GPS Measurements with Variance Component Estimation[J]. IEEE Geoscience and Remote Sensing Letters, 2012, 9(4): 754-758 doi: 10.1109/LGRS.2011.2181154
[44]	房立华, 吴建平, 苏金蓉, 等. 四川九寨沟Ms 7.0地震主震及其余震序列精定位[J]. 科学通报, 2018, 63(7): 649-662 https://www.cnki.com.cn/Article/CJFDTOTAL-KXTB201807007.htm Fang Lihua, Wu Jianping, Su Jinrong, et al. Relocation of Mainshock and Aftershock Sequence of the Ms 7.0 Sichuan Jiuzhaigou Earthquake[J]. Chinese Science Bulletin, 2018, 63(7): 649-662 https://www.cnki.com.cn/Article/CJFDTOTAL-KXTB201807007.htm
[45]	梁建宏, 孙丽, 刘杰. 2017年四川九寨沟Ms 7.0地震及余震精定位研究[J]. 地球物理学报, 2018, 61(5): 2152-2162 https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201805038.htm Liang Jianhong, Sun Li, Liu Jie. A High Precision Relocation Study of the Ms 7.0 Jiuzhaigou Earthquake and the Aftershocks Occurred in 2017[J]. Chinese Journal of Geophysics, 2018, 61(5): 2152-2162 https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201805038.htm
[46]	梁姗姗, 雷建设, 徐志国, 等. 2017年四川九寨沟MS7.0强震的余震重定位及主震震源机制反演[J]. 地球物理学报, 2018, 61(5): 2163-2175 https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201805039.htm Liang Shanshan, Lei Jianshe, Xu Zhiguo, et al. Relocation of Aftershocks of the 2017 Jiuzhaigou, Sichuan, Ms 7.0 Earthquake and Inversion for Focal Mechanism of the Mainshock[J]. Chinese Journal of Geophysics, 2018, 61(5): 2163-2175 https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201805039.htm
[47]	胡晓辉, 盛书中, 万永革, 等. 由余震分布确定2017年九寨沟Ms 7.0地震发震断层面参数[J]. 地震研究, 2019, 42(3): 366-371 doi: 10.3969/j.issn.1000-0666.2019.03.008 Hu Xiaohui, Sheng Shuzhong, Wan Yongge, et al. Fault Plane Parameters of 2017 Jiuzhaigou Ms 7.0 Earthquake Determined by Aftershock Distribution[J]. Journal of Seismological Research, 2019, 42(3): 366-371 doi: 10.3969/j.issn.1000-0666.2019.03.008
[48]	Xie Z J, Zheng Y, Liu C L, et al. An Integrated Analysis on Source Parameters, Seismogenic Structure and Seismic Hazard of the 2014 Ms 6.3 Kangding Earthquake, China[J]. Tectonophysics, 2017, (712/713): 1-9
[49]	Wells D L, Coppersmith J J. New Empirical Relationships Among Magnitude, Rupture Length, Rupture Width, Rupture Area, and Surface Displacement[J]. Bulletin of the Seismological Society of America, 1994, 84(4): 974-1002
[50]	Li Y X, Bürgmann R, Zhao B. Evidence of Fault Immaturity from Shallow Slip Deficit and Lack of Postseismic Deformation of the 2017 Mw 6.5 Jiuzhaigou Earthquake[J]. Bulletin of the Seismological Society of America, 2020, 110(1): 154-165 doi: 10.1785/0120190162
[51]	Liu G, Xiong W, Wang Q, et al. Source Characteristics of the 2017 Ms 7.0 Jiuzhaigou, China, Earthquake and Implications for Recent Seismicity in Eastern Tibet[J]. Journal of Geophysical Research: Solid Earth, 2019, 124: 4895-4915 doi: 10.1029/2018JB016340

施引文献

资源附件(0)

图(6) / 表(3)

计量

文章访问数: 843
HTML全文浏览量: 202
PDF下载量: 138
被引次数: 0

1 STLS平差模型及平差准则
2 STLS模型的通用表达式及算法推导
3 STLS模型平差准则的优化选取
4 STLS模型参数估值的近似方差
5 数值实验
6 结语

联合地震位错模型和InSAR数据构建2017年九寨沟Mw 6.5地震同震三维形变场

作者简介: 彭颖，硕士，研究方向为InSAR高精度三维形变序列测量。pngying@163.com

通讯作者: 许才军，博士、教授。cjxu@sgg.whu.edu.cn

计量

出版历程

Deriving 3D Coseismic Deformation Field of 2017 Jiuzhaigou Earthquake with Elastic Dislocation Model and InSAR Data

Author Bio: PENG Ying, master, specializes in 3D deformation field construction with InSAR data. E-mail: pngying@163.com

Corresponding author: XU Caijun, PhD, professor. E-mail: cjxu@sgg.whu.edu.cn

1 STLS平差模型及平差准则

2 STLS模型的通用表达式及算法推导

3 STLS模型平差准则的优化选取

4 STLS模型参数估值的近似方差

5 数值实验

6 结语

计量

出版历程

目录

1 STLS平差模型及平差准则

2 STLS模型的通用表达式及算法推导

3 STLS模型平差准则的优化选取

4 STLS模型参数估值的近似方差

5 数值实验

6 结语

作者简介:
彭颖，硕士，研究方向为InSAR高精度三维形变序列测量。pngying@163.com

通讯作者:
许才军，博士、教授。cjxu@sgg.whu.edu.cn

Author Bio:
PENG Ying, master, specializes in 3D deformation field construction with InSAR data. E-mail: pngying@163.com

Corresponding author:
XU Caijun, PhD, professor. E-mail: cjxu@sgg.whu.edu.cn