基于相位差谱合成时频非平稳地震波

基于相位差谱合成时频非平稳地震波 第 卷 第 期地 理 与 地 理 信 息 科 学 ,:,,,, ，:,, ,,, 年 月 ,:,,,,,:，,,,,,,,:,,,,:，,,,,:,,,,,,,, ,,,,,,,:,, 基于相位空间相关性分析的 技 术 ,,,,,,, 在地面沉降监测中的研究与实践 ,,,,，，，刘 欢 欢范 景 辉陈 建 平郭 小 方 ,，,,，中国地质博物馆 北 京 中国国土资源航空物探遥感中心 北 京 ,,,:::,,,,,:::,,，,中国地质大学北京市国土资源信息开发研究实验室 北 京 ,,,:::,, 摘 要 ,技术克服了图像失相干和大气影响等限制因素 ，扩 展 了 在提取大面积微小形变 上 的 应 用 ， ,,,,,,, ,,,,,, 但 基 于亮度离散度和先验形变模型的 有一定的局限性 。 基于相位空间相关性分析的 技 术 通 过,,,,,,, ,,,,,,, 分析相位各个组分的空间相关性 ，完 成 的 选 择 和 形 变 信 息 提 取 ，可以在没有亮散射体的情况下选取 点 和 ,,,,提 测 ，，。 取各种地 表 关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf 形变信息 且不受先验的形变时间 模 型 限 制 可应用于各种形变研究 该文以天津地区地面沉降监 实 践 为 例 ，描述了该方法的关键技术和工作流程 ，利 用 幅 数 据 ，得到了较好的监测结果 。,,，，,,,,,,,, ,,, 关 键 词 相位空间相关性 地 面 沉 降,,,,,,,,,,,,中 图 分 类 号 文 献 标 识 码 文 章 编 号 ,,,,,,,,,,,,,:,:,,:,,:,,::,,,:, 。地面沉降 在我国许多地区均有出现 ,,,,,,数 据 预 处 理 , ，、、对 地 表 微 小 形 变 很 敏 感 且 具 有 全 天 候 全 天 时 广 、，覆 盖 高 分 辨 率 等 优 点 但 由 于 两 次 观 测 条 件 变 化 造 实 验 数 据 是 年 月 至 年 月 的 ,::,,:,::,,,,, ,，景 天 津 地 区 数 据 极 化 成 像 成 的 干 涉 图 像 时 空 失 相 干 以 及 大 气 不 均 匀 产 生 的 大 ，，,,,,, ,,,, ,, ,。，模 式 为根 据 最 小 化 干 涉 图 的 非 相 干 之 和 以 气 相 位 延 迟 等 严 重 限 制 了 其 应 用 范 围 和 监 测 精 ,,, 。及度 ， 技 术 通 过 选 取 点 并 分 析 这 些 呈 类 点 散 最 大化干涉图的相干性之和的 主 影像选择 原 则 以,,,,,,, ,, ，年 月 日 的 为 主 影 像 生 成 幅 差 分 ,::,,,,,，: ,,，射 体 的 像 素 的 相 位 来 反 演 形 变 信 息 保 证 了 相 位 ，，干 涉 图 研 究 范 围 为 行 列 幅 辅 影 质 ,,:::,,::,,,, 。量 通 过建立形变模 型和对这些像素的时序 滤 波 可像与主影像的时间和空 间 垂直 基线 如 表 所 示 , ，去 除 大 气 相 位 等 误 差 克 服 了 图 像 失 相 干 ,,,,,, 表 幅辅影像与主影像的时间和空间垂直基线 ,,,,, ,, ,,,, ,,,:,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, ,,,,,,,,,,, ,,,,时 间 基 线 空 间 垂 直 基 线 编 号 成 像 时 间 。 ,,和 大 气 影 响 等 限 制 因 素但上 述方法大多通过 ,,,,,,,,,,,,,,,,,,,,,,,,, ,::,,,:,,, ,,,, ,,, ，分 析 像 素 的 亮 度 值 及 其 变 化 来 选 取 点 在 城 市 地 ,, ,,::,,,,,,,,,,,,,, ,,::,,:,,,:,,,:,,,，，区 比 较 适 用 应 用 于 自 然 地 貌 时 会 因 为 亮 度 值 或 亮 ,,::,,:,,:,,,,,,,,,。 ， 度 变 化 的 阈值太低 而不能得到可 靠 的 结 果 再 者 ,,::,,:,,:,,,,:,,: ,,::,,:,,,,,,:,,,,,通 常 需 要 用 形 变 的 时 间 模 型 识 别 点 并 去 除 其 干 ,, ,,::,,:,,,,,,:,,,，,扰 信 息 会 产 生 如 下 问 题 一 是 模 型 通 常 被 假 设 为 匀 ,,::,,:,,,,,,,,,, ,,::,,,,,,:,:,,,,,，,速 或 周 期 性 等 不 适 用 于 不 规 律 形 变 二 是 通 常 地 表 ,:,::,,:,,,,,,:,,,,，随 时 间 的 形 变 规 律 本 来 就 是 研 究 目 的 之 一 先 验 的 ,,,::,,:,,:,,,:,,,, ,,,::,,:,,,,,,,,,:。时 间 模 型 显 然 会 有 限 制 ,,,::,,,,,,:,,:,, ，天 津 在 地 质 构 造 上 易 发 生 地 面 沉 降 虽 然 自,前 期 的 主 要 工 作 流 程 是 读 取 主 辅 影 像 的 ,， ，年 加 强 了 沉 降 控制和监控体 系 建设 取 得 了 一 ,,,, , ，精 轨 数 据 完 成 辅 影 像 到 主 影 像 的 粗 配 准 和 精 ，。 :,定 成 果 但 沉 降 形 势 仍 比 较 严 峻 本 文 以 天 津 地 ，，配 准 读 取 高 程 数 据 重 采 样 幅 辅 ,,,, ,，, ,, ，区 地 面 沉 降 为 例 通过 分析差分干涉图上各像 素 的 影 。道 计 算 出 的平地相位 及 地 形 相位分量和地理 编 码 位 组 分 及各组 分 的 空 间 相 关 性 完 成 点 的 选 相 ,,，， 像 对主辅影像共轭计算 生 成 干 涉 图 去 除 用 精 密 轨，在 对 辅 影 像 到 主 影 像 进 行 子 像 素 级 精 配 准 时 考 ，取几 ，乎适用于所 有 地 表 情 况且 不 受 形 变 模 型 限 虑 ，，制能收 获 稿 日得 期 各 , 个 成 像 时, 期 修 的 订 形 日 期变 , 信 息类 似 的 方 法 已 ,:,,,,,,,,,:,,,:,,:,,,“”基 金 项 目 国家公益性基金项目 地质灾害雷达干涉测量预警系统的研究与 设计 领导形象设计圆作业设计ao工艺污水处理厂设计附属工程施工组织设计清扫机器人结构设计 ,::,,,:,,，,,:,,被 应 用 于 监测火山及青藏高原北麓地区的冻土,，。,,,,，，，，作 者 简 介 刘 欢 欢 女 回 族 硕 士 研 究 生 助 理 工 程 师 研究方向为遥感与地理信息系统应用 ,,,,,,，,,,,,,,,,,,,,,,,,:,,,,,,。形变 地 理 与 地 理 信 息 科 学 第 页第 卷 ,,,, :, 到 有 些 辅 影 像 与 主 影 像 的 相 关 性 太 低 会 降 低 配 准 精 :,,,, ,,，,,,，，，，，，γ ,?,,,,,,,,,,,,Δ,:,,,, , : ,,,，度 依 据辅影像与主 影像空间垂直基线值的 正 负 将 :,，式 中 是 可 用 干 涉 图 的 个 数 是 高 程 误 差 : ，，,:,, Δ,，其 分 为 两 组 分别对 两组辅影像按空间垂直 基 线 的 的 估 计 值 。 ，，Δ,:,,,，绝 对 值 从 低 到 高 排 序 各组中先进行与主影 像 空 间 ，计 算 点一定邻域内的相位均值 需 要 其 他 ,,， 基 线 最 近 的 三 幅 辅 影 像 与 主 影 像 的 配 准 然 后 依 次,点 的 分 布 情 况 在 没 有 任 何 点分布信息的条 ,,,, ，进 行 某 一 幅 辅 影 像 与 其 后 的 幅 辅 影 像 的 配 准 不 , ，，件 下 可以对所有分布点 进 行 迭 代 计 算 最 终 识 别 出 ，足 幅 时 递 减 直至 对与主影像空间垂直基 线 最 远 , 。，点 虽然在计算平 均 相 位 及 估算 和时 ,,,，，,:,,γ, , ，的 两 幅 辅 影 像 配 准 结 束 最终用加权最小二 乘 法 反 ，候选点的相位 可 能 受 较 强 的 噪 声 干 扰 使 得,,′。 演 得 到 辅 影 像 到 主 影 像 的 配 准 系 数 因 原 影 像 为 降 ，， ,,,，，轨 数 据 对 差 分 干 涉 图 做 距 离 向 的 镜 像 处 理 基 线 较 ，，不 可 忽 略 但 从统计角度分析 较 大 的 像 素 更 有 可 γ, 、短 相 干性较好的差 分干涉相位图上有较清 晰 的 相 。，能 是 点 每 次 迭 代 计 算 中 只 选 取 较 高 值 ,,γ, ，位 条 纹 可 从 中 大 概 推 测 出 形 变 的 程 度 和 趋 势 分 的 ，的 像 素 被 去 掉 剩 余 候 选 点 的值 也 会 逐 ′，，, ,,,γ, 布 而 相 干 性 较差的差分 干涉相位则受干 扰 比较 ，，候 选 点 来 估 算 和随 着 迭 代 次 数 的 增 加 ，，渐 降,,:,,γ, ,， ，。严 重 条 纹 不 清 晰 不 能 直 观 反 映 地 表 形 变 信 息 低 ，低 受的 影 响 逐 渐 变 小 而 主 要 取 决 于，，，，γ′, ,,,,,点 的 选 择 , ,,。 , ,,当的 均 方 根 小 于 某 一 给 定 阈 值 经 验 值 时 :,:,γ, 通过分析差 分 干 涉 相 位 的 组 成 及 各 组 分 的 空 间 ，。 达 到 收 敛 要 求 迭 代 终 止 此 次 试 验 选 取 亮 度 离 散，，相关性得到一个相位的稳定性 判 定 标 准 量通 过 计 γ ，度 高 于 的 点 作 为 候 选 点 进 行 迭 代 计 算 迭 :,,,,, 算其概率密 度 函 数 及 亮 度 离 散 度 和 值 的 相 关 性 约 γ 。代 终 止 至的 均 方 根 小 于 :,::,γ, ，。。束确定的阈值最后对选取出的 点进行精简 ,,γ 点 的 选 择,,, ,,点 的 判 定 标 准 量 ,,, ,,γ由于随机相位地面 目 标 点 也 有 一 定 的 概 率 具 有 对 于 差 分 干 涉 图 序 列 中 任 一 幅上 的 第个 像 , , ,,,,,,:,, ，较 高 的确 定的 一 个 原 则 就 是 使 真 正 点 ,,γγ，,素 其 相 位 值 可 以 表 示 为 ，最 大 化 地 保 留 下 来 同 时 非 点 被 误 选 的 概 率 也 要 ,,,, ，，，，，，，，，，，，,,，，，， ,,,,,,,,,,,:,,,,,, ,,:,,,,,,, 。,,小于某一合理 的 阈 值 的 概 率 密 度 函 数 ,,式等 号 右 边 依 次 是 像 素 视 线 向 的 形 变 ,γ, ,γ, ,, 相 ,,为对 应 的 随 机 相位像素的概 率 密度函数 ,,,,,,,,,, γ, ,, γ, ,,,,， ,γ, α,, γ, α,,, γ, 、、、位 高 程 误 差 相 位 大 气 相 位 延 迟 轨 道 误 差 相 位 ，,,，,,,,假 设 当则 有 时 和 点 的 概 率 密 度 函 数 的 加 权 和 ::,,,,γγ, ,,,, γ, ?,,, , :,, :,, ,、、,。 和 噪 声 相 位 散 射 体 噪 声 热 噪 声 配 准 误 差 等 ,,,,,,, ,,,,,γ, α ,, γ??: : ,笔 者 ,, γ,,, 分 别 在 式 等号左侧采用 试 验 数 据 的 候 选 ,,, ，,认 为越 小 像 素 越 有 可 能 为 点 认 为 形 变 ，，，，。 ,,,, 点 右 侧采用大量模拟 随 机 点 计 算 出 非 ,,,,, α、 相 位 大气相位延迟 和轨道误差的相位在特 定 范 围,点 被误选的概率为 , ，内 是 空 间 相 关 的 高 程误差和噪声误差在此 范 围 内 ,,,,,,,α,, γ, ?,,,,,,γ,, ,, ,, ,, ,γ,,γ, ,,,,,,?γ。，不 相 干 且 均 值 为 假 设 其 他 点 位 置 已 确 定 :,,,, ,，，，，，，，，，，，，,，，，， ,,,,,,,,,,,:,,,,,, ,,:,,,,,,,,,,,,,:,, γ,。，，出这 样 将和带 入 上 式 就 可 求 γ, α ，对 像 素 周 围 一 定范围 内像素的相位求均值 则 , ,,,,,将 式和 式相 减 可 得 ,,由 于 亮 度 离 散 度 和 相 位 稳 定 性 在 判 断 点 的 ,, ,, ,,,，,′,，，，，，，，，有 ,,,,,,,,,,,:,,,,,,Δ :，， ,,,， 过程中有一定的相关性 根 据 亮 度 离 散 度 对 ,, ,，,,,、式 中 是 形 变 相 位 大 气 延 迟 相 位 和 轨 道 误 ′，，,,,,,,,,,:,, :差 相，,,，候 选 点 分 组 分别计算每组的和 选 ， ,α,, γ 。 取 各 组 中高 于 阈 值 的 点 考 虑 到 像 素 的 值 一 。 位 与 其 均 值 分 别 求 差 值 的 总 和 γ, γ, :，般 会 随 着的 增 加 而 降 低 用最 小二乘法拟合出 , 由 于 高 程 误 差 与 垂 直 基 线 呈 正 比 ,，,, ,, ,,,,，，，，，， Δ,:,,,? ,,,:,,, ,二者之间的线性关系 ,,, 将 其 代 入 式 有 ,,,,,,,:,, :,, ,:,,,，,, γ ,,,,,,,，,′，，，，，，，，，，，,,,,,,,,,,? ,,,:,,,,,,,， ,, ,,,,,,:,,:。的 像 素 为 点 最 后 选 择 ,，式 中 比 例 常 数 是 唯 一 与 基 线 相 关 的 量 可 利，,,,, ,γ,，，,, ,:,,, ，。 用 所 有 干 涉 图 用 最 小 二 乘 法 估 计 得 到 点 的 精 简 ,,, ,, 根据上述相位稳定 性 测 度 和 亮 度 离 散 度 选 择 出，,,，假 设很 小 可 以 基 于 式 得 到 判 定 ′，，,,,, 像 素 ,，。， 相 位 稳 定 性 的 标 准 的 点不全是真正的 点 还 需 精 简 首 先 由 于,,,, ,刘 欢 欢 等 基于相位空间相关性分析的 技术在地面沉降监测中的研究与实践 ,,,,,,, 第 页第 期 ,,, ，,,强 散 射 体 经 常 会 影 响 到 周 围 像 素 的 散 射 特 性 并 差主要是空间非相干的 视 线 角 误 差 高 程 误 差 和 与 使 ，主 影 像 有 关 的 空 间 非 相 干 误 差 会 使 得 相 邻 点 ,, ，，邻 近 像 素 也 表 现 为 特 性 对 于 这 样 的 点 群 ,,,,间 ，，的 相 位 差 大 于大 大 降 低 了 相 位 解 缠 的 精 度 需 π只 ,，需 保 留 高值 的 点 再 者 对 于 一 些 稳 定 性 随 γ, 。要 在相位解缠前去除 其 中 空 间 非 相 干 的 视 线 角 误 ，改 变 的 像 素 利 用 方 法 迭 ,::,,,，,,，,,,,,,,时 间 ,,的 估 计 值 进 行 高 程 误 差 校 利 用 ，，差 ′Δ,:,,, 代 计 ，算 其的 变 化 剔 除 标 准 方 差 超 过 某 一 阈 值 γ, ，正 由 于 的 估 算 不 完 全 可 靠 带 来 的 空 间 相 干 ,，，,:,,, ,,为 的 点 最 后 去 除 被 误 判 为 的 噪 声 ,一 般 :,,,,误 差 和 其 ,、他 空 间 相 干 误 差 项 如 大 气 轨 道 误 差 ，。点 共 选 出 个 点 ,,,,,,,, 点间引入的空间 相 干 差 分 误 差 很 ,，等 在 邻 近 ,,形变信息的提取 , ,,，。小 在 解 缠 中 可 暂 时 忽 略 经 过 空间非相干误差 校 空 间 非 相 干 的 误 差 校 正 ,,, 正 的 点 的 差 ，,, 分 干 涉 图 序 列 各 个 时 期 的 干 涉 条 点 的 差 分 干 涉 相 位 中 较 大 的 空 间 非 相 干 误,, ,,,。性 较 差 的 像 对 更 清 晰 图 ,纹 包 括 原 本 相 干 , 图 经过空间非相干误差校正的 点的差分干涉图序列 , ,, ,,,, ,,,,,,,,,,,,,,,,,,,,,,:,,,,,,,,,,,,,,,,,,,:,,,,,,:,:,,:,,,,,,,,,:,,,,,,,,,,,:,, ,,,,, ，相 位 解 缠 的 实 现 滤 波进行估算和去除 这些误差相位一般分为 与 主 ,,, ,,。 影 像 有 关 的 误 差 和 与 辅 影 像 有 关 的 误 差 空 ,:，，该 解 缠 工 作 属 于 三 维 相 位 解 缠 问 题 被 分 解 ,,间 相 干 的 视 线 角 误 差 多由 于高程数据精度 为 时 间 维 的 相 位 解 缠 和 二 维 空 间 上 的 离 散 不 规 则 ,:，, ，分 布 不 够和高程数据向雷达坐标系 空 间转换时 造 成 需 。、。点 相 位 解 缠 由 于 点 在 二 维 空 间 中 呈 离 散 的 单 独 估 算 其 他 空 间 上 的 非 相 干 误 差 项 用 噪 声 模 型 ,, 不 。，，估 算 规 则 格 网 分 布 在 进 行 二 维 空 间 相 位 解 缠 时 首 先 对 因 为 会 随空间垂直基 线 增 大 而 呈 比 例 增 ,:，, ，点 进 行 格 网 化 和 格 网 化 滤 波 用 最 邻 近 距 离 法 ,,，，大 而垂直基线越大干涉 相 位 中 噪 声 成 分 越 多 另 外 进 ，行 格 网 插 值 然后利 用统计网络费用流算法 对 格 。 对 相 位 质 量 也 有 较 大 影 响 因 此 在 完 成 第 一 ,:， 网 迭 代 校 正 空 间 相 干 的 误 差 项 ,,, ，次相位解缠和误差 相 位 估 算 之 后 将 和 ,:，, ,:， ，化 数 据 进 行 相 位 解 缠 最 后 进 行 格 网 内 的 相 位 解 ，相 位 解 缠 之 后 空 间 相 干 的 误 差 项 仍 会 对 形 ，，。 去 除 再进行相位解缠 会 提 高 相 位 解 缠 的 精 度 重 ， 缠 ，、、变 相 位 产 生 干 扰 其 中 对 于 大 气 相 位 延 迟 轨 道 误 复 进 行 估 算 去 除 误 差 相 位 和 相 位 解 缠 的 迭 代 过 ，程 。 得 到 每 个 点 的 真 实 相 位 ,,，，差 相 位 等 空 间 相 干 的 误 差 组 分 可利用其时间不 直 至相邻两次迭代计算的结果 非 常接近 得 到 更 稳 。 定和精确的解缠结果 相 干 性 ，与 形 变 相 位 不 同 的 特 征 综合利用时间和空 间 上 的 地 理 与 地 理 信 息 科 学 第 页第 卷 ,,,, ，分 别 计 算 年 和 年 点 的 沉 降 速 率 ,::,,::,,,天津地面沉降监测结果与验证 , ， 并 转绘到地理坐标系下 可以看到天津地区的 沉 降 雷 达 视 线 向 的 累 积 沉 降 量 ,,, ,，,。 速 率 和 趋 势 分 布 变 化 较 大 图 见 封 在 ,,,::, , 雷达视线向的累积形变量与其解缠相位呈正比， 年 的 形 变 速 率 图 上 呈小面积沉 降中心离散分布 ,::: ,,,λ, ,，，到 年 市 区 以 西 的 沉 降 得 到 控 制 一 ,, ,::, ,,:,, ,,,,,,,,, ,π 些 ，沉 降 中 心 消 失 而市区南部的沉 降 中心 和 沉 降 , , 在 得 到 幅 差 分 干 涉 相 位 图 的 相 位 解 缠 结 ,,,,，。 更 加 严 重 沉 降 中 心 扩 大 成 沉 降 区 由 于 点 ,,,，，果 后 即 可 算 得 雷 达 视 线 向 的 形 变 量 此 次 研 究 范 ，在 空 间 上 呈 离 散 分 布 分 别 对 其 年 和 年 ,::, ,::, 围，平 均 沉 降 速 率 进 行 线 性 差 值 得 到 该 区 的 ，,,，, ，的 视 线 角 均 在 左 右 用 雷 达 视 线 向 的 形 变 量 做 沉 ,:? 年 。 降 趋 势 的 比 较 分 析 影 响 不 大 以 天 津 西 北 角 某 相 对 ，，均 沉 降 速 率 场 又 以 各 自 的 沉 降 最 低 值 为 起 点 生 ，沉 降 量 最 小 的 个 点 群 为 稳 定 参 考 点 以 ,,,,,:: 成沉降中心的形变分析 ,,, , ,，,。 了 间 隔 的 形 变 等 高 线图 见 封 ， ,:,, ,,以 主 影 像 成 像 时 期 为 基 准 对 各 沉 降 中 心 周 围的辅 影像与主 影像的差 分 干涉图的解缠结 ,,,,:,, 范 围 内 约 个 点 进 行 形 变 统 :,:,:,,,，,,, ,, :，果 为 基 准 计 算 各 个 点 在 之 后 各 成 像 时 期 的 相 ,,,,。 计图 各 组 点 各 个时期的形变 量 分布较集 ,,, ，位 解 缠 结 果 与 基 准 值 的 差 值 可 得 到 ,::,,,,,:，，中 每 组 点 的形变规律相似 由 于形变量与形变 ,,到, ，相 位 线 性 相 关 验 证 了 这 些 点 形 变 相 位 的 空 间 相 ,,。 对 应 辅 影 像成像时期 的雷达视线向的 累 积形 变 量 ,,，，，从 图 见 封 上 可 见 天津主城区的沉降很小 市 关 性 较 好 故 可 利 用 形变相位 的空间相关性 识 别 ,,,, 、。，，区 南 部 在 之 后 发 生 了 大 面 积 大 幅 度 点 距主影像较近时期 的 形 变 量 绝 对 值 较 小 空 间 ,::,,:,,:, ，，，。 相 关 性 好 数 值 分 布 集 中 而 距 主 影 像 成 像 时 间 较 的 沉 降 其 他 较 小 区 域 沉 降 中 心 离 散 分 布 在 最 、，，长 形 变 值 较 大 的 尽 管 空 间 相 关 性 较 好 但 由 于 形 一 期 图 上 标 注 了 形 变 参 考 点 和 个 形 变 严 重 后 , ，。变 绝 对 值 较 大 数 值 分 布 相 对 较 分 散 由 于 该 方 法 的 区 ,、、，域 天 津 市 区 以 东 的 胜 芳 王 庆 坨 西 青 区 市 ，不 受形变模型的限制 监测的形变规律更符合 实 际 ，区 以 南 的 静 海 县 杨 成 庄 乡 靠近市区的北辰区与东 ，。情 况 更利于探讨沉降成因 并 因 地 制 宜地 治理沉降 、。 丽 区 交 界 附 近 西 青 区 大 寺 镇 及 津 南 区 由 于 本 ，假 设 的 形 变 模 型 的 限 制 图 理论上能方 法 不 受 , 反映 天 津 各 个 时 期 真 实 的 雷 达 视 线 向 的 地 表 形 。变 年 和 年 雷 达 视 线 向 的 形 变 速 率,,, ,::,,::, 图 沉降中心及参考点的雷达视线向形变时序统计 , ,,,, ,,,,:,,,,:,,,,,:,,,,,,,,:,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,:,,,, ,,, ,， 结 果 验 证 的 杨 柳 青 的 沉 降 量 为 发 展 速 度 比 当,,, ,:,:, ,,: ,，《》年老沉降中心还要快 而 本 次 试 验 中 主 影 像 的 雷 根 据 天津 市 地面沉降调查 与 监 测 所 绘 制 的 达 —，年 地 面 沉 降 速 率 等 值 线 图 统 计 市 区 地 ,::,,::,，、视 线 角 约 年 年 这 些 沉 降 中 心 雷 达 视 ,:?,::,,::,,，面 沉 降 已 控 制 在塘沽区沉降量由,:,,,,,,:: :，线向的沉降速 率 约 为 多 可 以推测其垂直方 ,:,,,,，、降 至 大 港 区 汉 沽 区 地 面 沉 降 ,,,,,,:,,,:，向上的速率约为 多 考 虑 到 监测时段不完全 ,:,,,也 不 同 程 度 得 到 控 制 天津市区西南部东丽 区 军 粮 ，。一 样 可以认为两个结果基 本 吻 合 、、、、城 津 南 区 咸 水 沽 西 青 区 杨 柳 青 武 清 区 杨 村 静 海 “本文沉降中心的分 布 和 天 津 市 区 沉 降 很 。 县 城 关 至 大 邱 庄 一 带 形 成 了 新 的 沉 降 中 心 本 研 ”，小 的 结 论与已有研究结果基本一致 沉 降 的 趋 。 ，体 现 了 上 述 结 果 其 中 津 南 区 的 咸 水 沽 和 西 青 究,，，,,,,,区 势 分 布 和 沉 降 量 也 与 范 景 辉 等的 监 测 结 果 接 ,刘 欢 欢 等 基于相位空间相关性分析的 技术在地面沉降监测中的研究与实践 ,,,,,,, 第 页第 期 ,,, ,,,,,参 考 文 献 ,点 与 董 克 刚 等的研究结果较一致 宝 坻 城 关 以 南 ,, ，，,,,,，,,，,,,,,,,,, :,:::, ,,,,，,,,,,,,,,,,,，,,,,,, ， 广 大 平 原 均 有 不 同 程 度 的 地 面 沉 降 南 部 和 滨 海 地,,,,，,,，:,,,，,,,，，， ,，,,,,,,,:,,:, ,,:,,,,,,,,,, ,,,:,,, ，。区 尤 为 严 重 且 与 河 北 省 沉 降 区 连 成 一 片 ,::,，，,,,,,,,,,,,::,,,,,,,:,, ,,,,年 ,, 津 南 区 沉 降 与 任 立 生 等测 得 的 年 沉 降 量 范 景 辉 基于相干目标的 技 术地表形变监测研究与应 ,,,,,,,, ,,, ,,，用 中国科学院遥感应用研究所 ,,,::,, ,:,, ，，，范 景 辉 郭 华 东 郭 小 方 等 基 于相干目标的干涉图叠加方法 ,,，、、接 近 而 沉 降 中 心 的监 测结果与母景琴 ,,, ,,,,,，，,,,监测天津地区地面沉降 遥 感 学 报 ,,,::,,,,,,,,,,,, ,,,,。 的 监 测 结 果 接 近据 天 津 地 区 年 垂 ,, ，，，范景辉李梅郭 小 方等基 于 方 法 和 数 据 监 测 ,,,,,,::, ,,,,,,,,, ,,,, ,,，，,,,天津地面沉降的试验研究国土资源遥感,,,::,,,,,,,,, ，直 形 变 速 率 图 天津中部地区的杨柳青 咸 水 沽 一, ,,，，葛 大 庆 王 艳 范 景 辉 地 表 形 变 监测方法及关键问题 ,,,,,,,, ,， 带 都 是 沉 降 速 率 较 大 的 中 心 地 带 其 中 沉 降 中 心 南分 析 ,,国 土 资 源 遥 感 ，,,, ,,,::,,,,,,,, ， ,, ，，，侧 南 移 至 陈 官 庄 以 南 天津地区南部的沉降 速 率 增葛 大 庆 王 艳 郭 小 方 等 基于相干点的多基线 技 术 ,,,,,,,, , 与地表形变监测 ,,遥 感 学 报 ，，,,, ,,,::,,,,,,:,,,,,,，大 咸 水 沽 大 苏 庄 的 显 著 沉 降 带 向 北 萎 缩 退 出 大 ,,, ，傅 文 学 田 庆 久 技 术 及 其 在 地 表 形 变 监 测中的应用现状与 ,,,,,苏 庄 地 带 天津西部 的杨柳青后期比前期的 沉 降 速 ,,，，,,,发 展 地 球 科 学 进 展 ,,,::,,,,,,,,,,,,,, , ,,，，，率 有 所 减 少 南部的 陈官屯地段和东北部的 宁 河 地 薛 跃 明 郭 华 东 王 长 林 等 基 于 技术的矿区地表形 ,,,,,,,, , ,变 检 测 研 究 ,,遥 感 应 用 ，,,, ,,,,::,,,,,,,,，段 沉 降 速 率 增 大 在 宁 河 还 形 成 了 局 部 沉 降 中 心,, ，，汤 益 先 张 红 王 超 基 于 永 久 散射 体 雷 达干涉测量的苏州地区 ,,,,。 。上 述 结 果 与 本 文 的 监 测 结 果 基 本 一 致 ,,，，,,,沉 降 研 究 自 然 科 学 进 展 ,,,::,,,,,:,,,,:,:, 结 语 , ,, ,:,，,,，,,，,,,，,,,,,,,,,,,,,，,,,,，,,,,，,,，:,,,，,:，:，,,,, ,,,,,,:，,,,,:,,,,,,,,,,,,:,,,,,:,,:,,,,,,,,,:，,,,,:,, 技 术 方 法 通 过 对 相 位 组 分 的 空 间 相 ,,,,,,, ,,， ,,,,,,,:，,,,,,,，,,,,::,,,，，关 性 分 析 几 乎 适 用 于 所 有 地 表 点 的 选 取 且 不 ,, ,, ，，，:,,,，,,，, ，,，，，,,，,,,，,,,，,，，,,,,,,::,,,，,,,,, ó ,,,,,,,,，,,,:,,,，,,,,,,，,,,,,,，,,,,,:,,,,:,,,,,,,,, ，受 形 变 模 型 限 制 可 以获得各个成像时期的 形 变 信 ，,, ,,,,,,,,，:,,,,,,,，,,,,,,,,,,，,,，:,,,，,,,，,,,，,,,:,，，息 提 高 了 时 间 分 辨 率 对 形 变 规 律 及 成 因 的 分 析 ，,,, ,,,,,,,,:,,,,,,,,,,, 有,, ,,谢 酬 基于永久散射 体 的青藏高原冻土形变研究 中 国 科 ,,,,, ,，学 院 研 究 生 院中国科学院遥感应用研究所 。 ,::,, 重 要 的 指 导 意 义 地 面 沉 降 监 测 结 果 的 精 度 受 ,,, ，,,姚 国 清 母 景 琴 技 术 在 地 面沉降监测中的应用 ,,,,,,,,, ,,,、, 点 的 质 量 误差相位 的估计和相位解缠算法 的 稳 地 学 前 缘 ，，,,, ,::,,,,,,,,,,,,，健 性 等 影 响 如何优化 参数设置以得到密度和 质 量 ,, ，，，,,董 克 刚 周 俊 于 强 等 天 津市地面沉 降的特征及危害 地 ,,,,, ，，,,,质灾害与环境保护 ,::,,,,,,,,:, ，较 高 的 点 集 更 符 合 实 际 地 进 行 误 差 估 计 和 都 ,, ,, ，，，任 立 生 黄 立 人 周 俊 等 天 津 沉降监测的十年试验结果 ,,,,,,去 ，,除 还 需 进 一 步 研 究 本 文 相 位 解 缠 部 分 仍 未 做 到 ,,测 绘 科 学 ，，,,, ,,,::,,,,,,,,:, ,, ，母 景 琴 姚 国 清 干 涉 雷 达 时 间序 列 分 析方法在地面沉降监测 ,,,，真 正 的 三 维 空 间 相 位 解 缠 所得地表形变量是 雷 达 ,,，,,,中 的 应 用 国 土 资 源 遥 感 ,,,::,,,,,,:, ，，视 线 向 的 要 获 得 垂 直 高 度 上 的 形 变 量 除 了 要 考 虑 ,, ，，郭 良 迁 薄 万 举 陈 宇 坤 天 津 地区的垂直形变与构造活动研 ,,,，、视 线 角 的 变 化 还 要 考 虑 各 种 地 形 造 成 的 叠 掩 倒 置 ,,，，,,,究 大地测量与地球动力学 ,,,::,,,,,,,, ，。和 盲 区 等 问 题 这 有 待 今 后 研 究 解 决 ,,,,,,,,,,,,,,,,:,,,,,,,,,,,,,,,,:,,,:,,,,,,,,,,,,,,,,,,,:,,,,,,,,,,,,,:,,,,,,,:,,,,,,,, ,,,,,,,,，，，，,, ，,,,,,,,,,,，,,,,,,,:，，，,,,,,,,,,: ,,,:,,,, ,,,,,，, ,,,,,,,:,:,,,,,,,,,,:,,,,,,,,,,,:::,,,,,,，, ,,,,,,,,,,:,,:,,,,,,,,,,,, ,,,:,,,,,,,,,,,,,,:,,,,, , ,,,:,,,,,,,,,,,:::,,,,, ,,,,，，, ,,,,,,,,,,:,,,,,,,, ,,,,,,,,,,,,,,,,,,,,,,,,,:,,,,,,,,,,,,,,:,,:,,,,,,,,,,,,,,:::,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, ,,，,,,,,,,,,,,,，,,,,,,,:,,,,,,,:,,:,,,,,:,,,,,,,,,,,,:,,:,,:,,,,,,:,,,,,,,,,,:,,,,,,,,,,,，,,,,,,, ,,,,,，,,,,,,,,,,,,,,,,，,:,,,,,,,,,，,,,,，,,,,,,:,,,,，,,:,,:，，,,,,,:,,,,,,,,,,,:,,,,,:,,,,,,,,:,,,,, ，，，,,,,,,,,,,,,,,,,,:,,,,，,,,,,,,,,,,,,,,,,,,,,,,,，,,,,:,,,,,,,,，,,,:,,,,,:,,,,,,,,,,:,,，,,,,，,,,,,,,,,, ,,,,,,,，,，,,,,,:,,,,,,,,,,,,,,,,,,,,,,,，,,:,,,,,,,,,,,,,,,,,,,,,，,,，,:,,,:，,:,,，,,,,,,,，,,,,,,,,,,，:，,,,,,,,,,,,, ，，，,,,:，,,,:,,,,,，,，,,,，,,,:，,,,,,,,,,,,,,,,,,,，,,:,,,,,,,,,,,,，,,,,,,,,,,,，,,,,:，,,,,,,,,,，:,,,,,,,,, ,,,,,,,,,,,,,,,,,,,,,,:，，,,,,,:,:,,,,,，,,，:，,,,,,,,,,,,:，，:,,,,,,,,，,,,,,,,,:,,,,,,,,,:,,,,,,,，，,,,,,,,:,,，,:， ，,,,,,,,,,，:,,,,,:,,,,:，,,,,，,,,,:,,,,,,,,,,:，,,,,:,,,,,,,,,, ,,,,:,,,,:,，,,,，,,,,,,,,,:，,,,,:,,:,,,,,,,,,,,,,,,,, ，,,,:，,,,,:,,,,:，,,,,:,,:，,,,,,,,,,,,,,,,,:,,,,,,,,,,,,，:,,,,,,,,,,，,,:,,,,:,,,,,,,,,,,,,,,,,,,,,,,:,,,, ,,,,,,，，，,,,,,,:,,,,,,,:，,,,,:,,,,,,,,,,,,,,,,,,,,,,:，,,,,,,,,:,,,,,,,,,:，，,,:,,,,,，,,，,,,,,,,,:,,,,,,,,,,,:,,,,, ,,,,,，,，,,,,,,,,,:，,,,,:,,,,:，,,,,:,,，,，,,,,:,,,,,,,,,,，,,,,,,,，，,,,,, ,,,,,，:,,:,,，,,,，:, :,,,::,,,,,,,,, ，，,:,,,,::,,:，,,,,,,,，,,,,,, ,，,,,,,，:,，,,,,,:，,,,,,,,,,,,,:,,,,,,,,,:,:,,,:，,,,,,,,,,,,,,,,,,,,, ,,,,,,,,，，，,,，,,,,,,,，,,,,,,,，,,，:，,,,,，,,,,:，,,,,:,,,,,,,,,,,,,,,,,,,,:,,,,,,,:，,,,,:,,,,,,,,,:,,,,,,,,,:，,,,,,,,,,, ，,:，,,,,:,,,,,,,,,,,,,，,,,,,,,,,,,,，,,,,,,:，,,,,:,,,,:,,,,,，:,,，,:,,，，,,,,,,,，:,,,,:,,,,,,,,,,,,,,,: ,,,,,,,,,,,, ,,,,:,,,,,:,,,,,:，,,,,:,,,,,,,,,,,,,，,,,,,，:,,,,,,,,,,,,,,:,,,,,,,,,,,,, ,,,,,,,,:,,,,,,,,,,,,,,,,,,,,,,,, ,,,,,,,,,,,:，，,,,,,:,,,,,,,,,,,,