利用SARscape对日本填海造陆和天然气开采进行地表形变监测

        日本千叶市,是日本南部重要的工业港市。位于西部的浦安市是一个典型的"填海造田"城市,东南部的东金区有一片天然气开采区域,本文利用SARscape,用干涉叠加的方法,即PS和SBAS,对这两个区域进行地表沉降现象的监测。

1、监测区域

日本千叶市:

图1 日本千叶市监测区范围

2、数据情况:

2006-2010年,26期ALOS PALSAR-1数据。

2006-2010年,32期ENVISAT ASAR数据。

两组数据,分别是L波段的PALSAR数据和C波段的ASAR数据,覆盖相同的区域。数据量都大于20景,满足PS和SBAS的要求。

3、监测结果与分析

对两组数据,整个区域分别进行SBAS处理,形变速率图上可看出,两个结果是一致的,然而,PALSAR在空间的连续性上更好,ASAR则有较多的无效值,这也说明L波段的数据相干性要比C波段的高。不过,对已知的形变区域,两个结果表现出很高的吻合度。

图2 SBAS处理得到的平均形变速率(蓝色:-15mm/年,红色5mm/年),左图为PALSAR的结果,右图为ASAR的结果

在整体监测的结果图上,监测到了两个形变区域:一个是位于西部的一个垃圾填埋场,另一个是位于东南部的天然气开采区域,这两个区域也是我们已知的形变区域。下面分别对这两个已知的形变区域再做局部的分析。

1、垃圾填埋(填海造陆)区域:下图为千叶县浦安市在1950年、1975年、1980年的填海面积:

图3 千叶浦安市的填海面积

对这个区域,分别用PS和SBAS的方法进行SAR时序分析,得到的结果如下图所示:

图4 浦安地区2006-2010年的平均位移速速率图,26期PLASAR的结果(上),32景ASAR数据的结果(下),左边是SBAS结果,右边是PS结果

分析PALSAR和ASAR数据的PS和SBAS结果,均表现为:在城市区域,PS的点密度很高,和SBAS得到的结果一样,而PS结果的空间分辨率要比SBAS的高。而分析同一个点在监测时间段内的形变过程,如下图,绿色为SBAS的结果,红色为PS的结果,无论是哪一种数据,SBAS得到的结果更为平滑。

图5 同一个点在时间序列上的形变过程(绿色为SBAS结果,红色为PS结果,上图为PALSAR结果,下图为ASAR结果)

  1. 南部东金区域的天然气开采区

下图为天然气开采公司(KNG)提供的天然气资源储存和开采图以及相应的地质调查图。在调查范围5年内,绝大部分区域都有沉降发生,最大达到了20cm。这是由天然气开采再注入水造成的地表形变,

图6 南部东金天然气开采区域

图7 5年的东金区地表调查图

对该局部区域分别用PALSAR数据和ASAR数据做PS和SBAS处理,结果如下:

图8 东南地区2006-2010年的平均位移速速率图mm/年

26期PLASAR的结果(上),32景ASAR数据的结果(下),左边是SBAS结果,右边是PS结果。对比两个数据的PS和SBAS的结果,PALSAR数据在非城市区域得到的PS点密度和SBAS的一致,而ASAR数据,由于在天然气开采之后注入大量的水,PS得到了很少的一些PS点,不过SBAS能反应出整个区域的形变情况。

4、结果验证

用从GSI地理网得到的实测的GPS数据进行结果验证。下图是稳定的GPS点的位置分布:

图9 稳定的GPS点的位置

对GPS点测量到的值绘制在时间-位移坐标系,可以看到GPS的位移随着时间有所波动,整体来看表现出一个线性的变化,对实测值进行线性拟合,得到红色的这条趋势线。如下图所示:

图10 一个稳定的GPS点的位移实测值曲线(黑色线)及拟合结果(红色线)

将拟合的值作为GPS点的实测值,和SBAS得到的结果进行求差,两者之差的绝对值小于等于5mm/年。

图11 PALSAR的SBAS结果和GPS点拟合值的差

图12 ASAR的SBAS结果和GPS点拟合值的差

将GPS实测值的拟合值和SBAS得到的年平均速率之差进行统计分析,得到的结果如下表所示:

表1:GPS实测拟合值和SBAS得到的平均速率之差的统计值

GPS-PALSAR

GPS-ASAR

平均差

-0.28mm/年

-1.1mm/年

标准差

1.9mm/年

1.8mm/年

通过与实测值对比分析,SAR干涉叠加技术可以得到较高的形变测量精度。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:/a/538190.html

如若内容造成侵权/违法违规/事实不符,请联系我们进行投诉反馈qq邮箱809451989@qq.com,一经查实,立即删除!

相关文章

36-代码测试(上):如何编写Go语言单元测试和性能测试用例?

每种语言通常都有自己的测试包/模块,Go语言也不例外。在Go中,我们可以通过testing包对代码进行单元测试和性能测试。 如何测试 Go 代码? Go语言有自带的测试框架testing,可以用来实现单元测试(T类型)和性…

Point-Nerf复现及解析

Point-Nerf复现及解析 鸣谢:同组的李xx师兄博士(交流思路)、辰昶仪器的狗哥等人(帮忙down资源) 0.0我自己的复现工程0.1相关库介绍0.1.1 pytorch0.1.2 h5py0.1.3 Scikit-Image0.1.4 imageio0.1.5 scipy0.1.6 Matplotlib0.1.7 fonttools 0.2…

JAVA的学习日记DAY6

文章目录 数组例子数组的使用数组的注意事项和细节练习数组赋值机制数组拷贝数组反转数组添加 排序冒泡排序 查找多维数组 - 二维数组二维数组的使用二维数组的遍历杨辉三角二维数组的使用细节和注意事项练习 开始每日一更!得加快速度了! 数组 数组可以…

16. 网络编程(1)

Hi,大家好!从本节开始我们学习网络编程相关的知识。基于TCP服务器和客户端实现流程框架。 本节目录: 网络编程在软件开发中具有相当重要的作用,涉及到各方各面: 网络通信: Linux系统作为一个多用户、多任务的操作系统,网络通信是其重要的功能之一。通过网络编程,可以实现…

稀碎从零算法笔记Day46-LeetCode:互质树

这几天有点懈怠了 题型:树、DFS、BSF、数学 链接:1766. 互质树 - 力扣(LeetCode) 来源:LeetCode 题目描述 给你一个 n 个节点的树(也就是一个无环连通无向图),节点编号从 0 到 …

5款好用又免费的UI设计软件

之前我们分享了五款好用的制作原型的工具,制作完了原型,就要对界面进行优化,这个时候就是 UI 设计师的任务了,UI 设计软件对于设计师们来说是很重要的,UI 设计工具是否好用直接影响到最后结果的好坏,那么就…

[lesson20]初始化列表的使用

初始化列表的使用 类成员的初始化 C中提供了初始化列表对成员变量进行初始化 语法规则 注意事项 成员的初始化顺序与成员的声明顺序相同成员的初始化顺序与初始化列表中的位置无关初始化列表先于构造函数的函数体执行 类中的const成员 类中的const成员会被分配空间的类中…

Baichuan-7B-chat WebDemo 部署调用

Baichuan-7B-chat WebDemo 部署调用 Baichuan2 介绍 Baichuan 2 是百川智能推出的新一代开源大语言模型,采用 2.6 万亿 Tokens 的高质量语料训练。在多个权威的中文、英文和多语言的通用、领域 benchmark 上取得同尺寸最佳的效果。 环境准备 在autodl平台中租一…

MySQL排序原理与优化方法(9/16)

order by排序优化 MySQL排序策略 内存临时表 or 磁盘临时表? **内存临时表排序:**在MySQL中,使用InnoDB引擎执行排序操作时,当处理的数据量较小,可以在内存中完成排序时,MySQL会优先使用内存进行排序操作…

【LeetCode】动态规划类题目详解

所有题目均来自于LeetCode,刷题代码使用的Python3版本 动态规划 问题分类 如果某一个问题有重叠的子问题,则使用动态规划进行求解是最有效的。 动态规划中每一个状态一定是由上一个状态推导出来的,这一点区别于贪心算法 动态规划五部曲 确…

Qt控件---按钮类型

文章目录 QPushButton(普通按钮)QRadioButton(单选按钮)按钮分组 QCheckBox(复选按钮) QPushButton(普通按钮) 属性说明text按钮中的⽂本icon按钮中的图标iconSize按钮中图标的尺寸…

每日Bug汇总--Day05

Bug汇总—Day05 一、项目运行报错 二、项目运行Bug 1、**问题描述:**前端将从后台查询的数据作为参数进行get请求,参数为空 原因分析: 这种写法可能只支全局的参数调用方法的传参响应 代码实现 if (this.jishiName) {this.$http({url…

K8S node节点配置

1.开始操作之前要先关闭防火墙,SELinux,swap分区 关闭防火墙 sudo systemctl stop firewalld关闭SELinux sudo setenforce 0 # 临时关闭 sudo sed -i s/^SELINUXenforcing$/SELINUXper…

数据结构-----Lambda表达式

文章目录 1 背景1.1 Lambda表达式的语法1.2 函数式接口 2 Lambda表达式的基本使用2.1 语法精简 3 变量捕获3.1 匿名内部类3.2 匿名内部类的变量捕获3.3 Lambda的变量捕获 4 Lambda在集合当中的使用4.1 Collection接口4.2 List接口4.3 Map接口 HashMap 的 forEach() 5 总结 1 背…

56、巴利亚多利德大学、马德里卡洛斯三世研究所:EEG-Inception-多时间尺度与空间卷积巧妙交叉堆叠,终达SOTA!

本次讲解一下于2020年发表在IEEE TRANSACTIONS ON NEURAL SYSTEMS AND REHABILITATION ENGINEERING上的专门处理EEG信号的EEG-Inception模型,该模型与EEGNet、EEG-ITNet、EEGNex、EEGFBCNet等模型均是专门处理EEG的SOTA。 我看到有很多同学刚入门,不太会…

C++学习知识

C知识小菜单: 备赛蓝桥杯过程中的一些小知识积累,持续更新中! 文章目录 C知识小菜单:1.小数取整:2.小数点后保留几位:3.数字占几位字符:4. 求x 的 y 次幂(次方)5. 求平方…

Spring Security——13,认证成功失败注销成功处理器

认证成功&&失败&&注销成功处理器 说明:一、认证成功处理器1.1 自定义成功处理器1.2 配置自定义成功处理器 二、认证失败处理器2.1 自定义失败处理器2.2 配置自定义失败处理器 三、登出成功处理器3.1 自定义登出处理器3.2 配置登出处理器 四、完结撒…

NineData创始人CEO叶正盛受邀参加『数据技术嘉年华』的技术大会

4月13日,NineData 创始人&CEO叶正盛受邀参加第13届『数据技术嘉年华』的技术大会。将和数据领域的技术爱好者一起相聚,并分享《NineData在10000公里跨云数据库间实时数据复制技术原理与实践》主题内容。 分享嘉宾 叶正盛,NineData CEO …

node后端上传文件到本地指定文件夹

实现 第一步,引入依赖 const fs require(fs) const multer require(multer) 第二步,先设置一个上传守卫,用于初步拦截异常请求 /*** 上传守卫* param req* param res* param next*/ function uploadFile (req, res, next) {// dest 值…

【刷题】图论——最小生成树:Prim、Kruskal【模板】

假设有n个点m条边。 Prim适用于邻接矩阵存的稠密图,时间复杂度是 O ( n 2 ) O(n^2) O(n2),可用堆优化成 O ( n l o g n ) O(nlogn) O(nlogn)。 Kruskal适用于稀疏图,n个点m条边,时间复杂度是 m l o g ( m ) mlog(m) mlog(m)。 Pr…