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

本文主要是介绍利用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干涉叠加技术可以得到较高的形变测量精度。

这篇关于利用SARscape对日本填海造陆和天然气开采进行地表形变监测的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/898974

相关文章

Pandas使用AdaBoost进行分类的实现

《Pandas使用AdaBoost进行分类的实现》Pandas和AdaBoost分类算法,可以高效地进行数据预处理和分类任务,本文主要介绍了Pandas使用AdaBoost进行分类的实现,具有一定的参... 目录什么是 AdaBoost?使用 AdaBoost 的步骤安装必要的库步骤一:数据准备步骤二:模型

使用Pandas进行均值填充的实现

《使用Pandas进行均值填充的实现》缺失数据(NaN值)是一个常见的问题,我们可以通过多种方法来处理缺失数据,其中一种常用的方法是均值填充,本文主要介绍了使用Pandas进行均值填充的实现,感兴趣的... 目录什么是均值填充?为什么选择均值填充?均值填充的步骤实际代码示例总结在数据分析和处理过程中,缺失数

QT进行CSV文件初始化与读写操作

《QT进行CSV文件初始化与读写操作》这篇文章主要为大家详细介绍了在QT环境中如何进行CSV文件的初始化、写入和读取操作,本文为大家整理了相关的操作的多种方法,希望对大家有所帮助... 目录前言一、CSV文件初始化二、CSV写入三、CSV读取四、QT 逐行读取csv文件五、Qt如何将数据保存成CSV文件前言

通过Spring层面进行事务回滚的实现

《通过Spring层面进行事务回滚的实现》本文主要介绍了通过Spring层面进行事务回滚的实现,包括声明式事务和编程式事务,具有一定的参考价值,感兴趣的可以了解一下... 目录声明式事务回滚:1. 基础注解配置2. 指定回滚异常类型3. ​不回滚特殊场景编程式事务回滚:1. ​使用 TransactionT

Java中使用Hutool进行AES加密解密的方法举例

《Java中使用Hutool进行AES加密解密的方法举例》AES是一种对称加密,所谓对称加密就是加密与解密使用的秘钥是一个,下面:本文主要介绍Java中使用Hutool进行AES加密解密的相关资料... 目录前言一、Hutool简介与引入1.1 Hutool简介1.2 引入Hutool二、AES加密解密基础

SpringSecurity6.0 如何通过JWTtoken进行认证授权

《SpringSecurity6.0如何通过JWTtoken进行认证授权》:本文主要介绍SpringSecurity6.0通过JWTtoken进行认证授权的过程,本文给大家介绍的非常详细,感兴趣... 目录项目依赖认证UserDetailService生成JWT token权限控制小结之前写过一个文章,从S

使用Jackson进行JSON生成与解析的新手指南

《使用Jackson进行JSON生成与解析的新手指南》这篇文章主要为大家详细介绍了如何使用Jackson进行JSON生成与解析处理,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起学习一下... 目录1. 核心依赖2. 基础用法2.1 对象转 jsON(序列化)2.2 JSON 转对象(反序列化)3.

C#使用SQLite进行大数据量高效处理的代码示例

《C#使用SQLite进行大数据量高效处理的代码示例》在软件开发中,高效处理大数据量是一个常见且具有挑战性的任务,SQLite因其零配置、嵌入式、跨平台的特性,成为许多开发者的首选数据库,本文将深入探... 目录前言准备工作数据实体核心技术批量插入:从乌龟到猎豹的蜕变分页查询:加载百万数据异步处理:拒绝界面

Python使用自带的base64库进行base64编码和解码

《Python使用自带的base64库进行base64编码和解码》在Python中,处理数据的编码和解码是数据传输和存储中非常普遍的需求,其中,Base64是一种常用的编码方案,本文我将详细介绍如何使... 目录引言使用python的base64库进行编码和解码编码函数解码函数Base64编码的应用场景注意

Java进行文件格式校验的方案详解

《Java进行文件格式校验的方案详解》这篇文章主要为大家详细介绍了Java中进行文件格式校验的相关方案,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起学习一下... 目录一、背景异常现象原因排查用户的无心之过二、解决方案Magandroidic Number判断主流检测库对比Tika的使用区分zip