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

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

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

倒计时4天!百度Create AI开发者大会“大模型与深度学习技术”论坛亮点抢鲜看!

作为人工智能的核心基础技术,深度学习具有很强的通用性,大模型技术在深度学习的基础上,通过构建更加庞大神经网络模型和应用transformer等更加领先的算法,使模型的处理能力产生质的飞跃。飞桨(PaddlePaddle&#xff09…

MySQL分区表(14/16)

分区表 基本概述 分区表是数据库中一种用于优化大型表数据管理和查询性能的技术。它将一个表的数据根据特定的规则或条件分割成多个部分,每个部分称为一个分区。每个分区可以独立于其他分区进行存储、管理和查询,这样可以提高数据处理的效率&#xff0…

VS Code中“@“符号如何自动补全导入路径

一、下载 Path Intellisense 插件 二、打开设置,在扩展中选择该插件,点击setting.json 三、添加配置: "":"${workspaceRoot}/src" 如图: 四、在项目src目录中新建jsconfig.json文件 (一定要是src目…

动态规划(背包问题)

一:动态规划概述: 动态规划实际上是一种将原本的 大 方面的问题转化为许许多多的 小方面 的一种应用, 在一定程度上避免数据的重复, 并且能够将数据以自己希望的方式进行存储, 用来解决多阶段的数学问题, 从而提高算法的效率 在算法当中, 动态规划主要包括有: 递推, 线性DP 记忆…

【Java核心技术】第3章 Java的基本程序设计结构

1 数据类型 Java一共有8种数据类型: 4种整型 类型存储需求int4字节short2字节long8字节byte1字节 2种浮点型 类型存储需求float4字节double8字节 1种字符型 1种布尔型 2 变量声明 2.1 局部类型推断 如果可以从变量的初始值推断变量类型,只需要使用…

【数组】5螺旋矩阵

这里写自定义目录标题 一、题目二、解题精髓-循环不变量三、代码 一、题目 给定⼀个正整数 n,⽣成⼀个包含 1 到 n^2 所有元素,且元素按顺时针顺序螺旋排列的正⽅形矩阵。 示例: 输⼊: 3 输出: [ [ 1, 2, 3 ], [ 8, 9, 4 ], [ 7, 6, 5 ] ] 二、解题精髓…

JVM参数列表

-client :设置JVM使用client模式,特点启动较快(神机不明显(I5/8G/SSD)) -server :设置JVM使用server模式。64位JDK默认启动该模式 -agentlib:libname[options] :用于加载本地的lib -agentlib:hprof :用于获取JVM的运行情况 -agentpath:pathnamep[options] :加载制定路径的本…

Day:007(1) | Python爬虫:高效数据抓取的编程技术(scrapy框架使用)

Scrapy的介绍 Scrapy 是一个用于抓取网站和提取结构化数据的应用程序框架,可用于各种有用的应用程序,如数据挖掘、信息处理或历史存档。 尽管 Scrapy 最初是为网络抓取而设计的,但它也可用于使用API提取数据或用作通用网络爬虫。 Scrapy的优势…

【Nacos】Nacos最新版的安装、配置过程记录和踩坑分享

Nacos是什么?有什么功能?大家可以自行联网(推荐 https://cn.bing.com/)搜索,这里就不做介绍了。 简单的看了下官网,安装最新版的Nacos(v2.3.2)需要使用到JDK(1.8.0&…

应急响应-战前反制主机HIDSElkeid蜜罐系统HFish

知识点 战前-反制-平台部署其他更多项目: https://github.com/birdhan/SecurityProduct HIDS:主机入侵检测系统,通常会有一个服务器承担服务端角色,其他主机就是客户端角色,客户端加入到服务端的检测范围里&#xff…

滑动窗口用法

文章目录 1. 长度最小的子数组(模板)2. 无重复字符的最长字串3. 最小覆盖字串4. 加油站5. 替换字串得到平衡字符串 1. 长度最小的子数组(模板) 题目分析 直接用步骤分析示例1,[]表示窗口,min_length表示满…

软件测试级别和对应要求

软件测试级别指的是将软件测试活动按照不同的开发阶段和测试目的进行分类,形成不同层次的测试过程。 分级依据,根据V模型 单元测试:这是软件生命周期中的第一个测试级别,主要针对软件的最小单元模块进行,例如类、函数…

Robotstudio2024中从备份文件恢复和创建工作站的具体方法演示

Robotstudio2024中从备份文件恢复和创建工作站的具体方法演示 如下图所示,打开Robotstudio2024软件,有需要的可以从以下链接获取: ABB机器人编程仿真软件RobotStudio 2024.1-链接baiduyun 点击“新建”—工作站—创建, 如下图所示,点击“ABB模型库”,选择自己使用的机器…

Proxmox VE qm 方式一键创建Linux虚拟机

前言 实现qm 方式一键创建Linux虚拟机,提高效率。 qm 一键创建Linux 虚拟机 以下实现在线下载镜像,创建虚拟机,安装系统需要自己手动安装哦,如果想实现全自动安装系统,建议部署自己的内网pxe server 系统参考各参数…

EditPlus来啦(免费使用!)

hello,我是小索奇 今天推荐一款编辑器,是索奇学习JavaSE时入手滴,非常好用哈,小索奇还是通过老杜-杜老师入手滴,相信很多人也是通过老杜认识嘞,来寻找破解版或者准备入手这个间接使用的编辑器~ EditPlus是…

HarmonyOS开发实例:【菜单app】

简介 分布式菜单demo 模拟的是多人聚餐点菜的场景,不需要扫码关注公众号等一系列操作,通过分布式数据库可以方便每个人可及时查看到订单详情,数量,总额等;效果如下 demo效果 工程目录 完整的项目结构目录如下 ├…

如何编译OpenHarmony自带APP

作者:王石 概述 OpenHarmony 的主干代码是开源社区的重要学习资源,对于想进行应用开发和熟悉 OpenHarmony 能力的同学主干代码是非常重要的资源,在主干代码的 applications 目录里聚集了很多原生的应用实现,那么如何编译这些代码…

广州南沙番禺联想SR530服务器主板传感器故障维修

今日分享一例广州市南沙区联想ThinkSystem SR530服务器sensor sysbrd vol故障问题维修案例; 服务器型号是:Lenovo thinksystem sr530 g6服务器 服务器所在位置:广东省广州市南沙区 服务器故障问题:机房异常停电,来电后…

【Entity Framework】EF连接字符串和模型

【Entity Framework】EF连接字符串和模型 文章目录 【Entity Framework】EF连接字符串和模型一、概述二、使用 Code First 按约定创建连接三、使用 Code First 和指定的数据库名称按约定创建连接四、将 Code First 与 app.config/web.config 文件中的连接字符串结合使用五、将 …

请核对您的姓名、证件号码、有效期和年限是否选择正确,请勿使用挂失过的身份证

问题 请核对您的姓名、证件号码、有效期和年限是否选择正确,请勿使用挂失过的身份证 详细问题 笔者在专利业务办理系统进行新用户注册,注册时间为晚上大概22:00以后。注册时已核对姓名、证件号码、有效期和年限,已确保正确&…