双星物语2存档(别拿豆包不当干粮，剑走偏锋的又耐玩的游戏走一波吧~~)
双星物语2存档文章列表:
- 1、别拿豆包不当干粮，剑走偏锋的又耐玩的游戏走一波吧~~
- 2、65周年 | 李志伟：InSAR部分地学参数反演
- 3、助力企业复工，交运董家口公司运送职工2000余名
- 4、16天一闪的它暗藏了多少信息？
- 5、星座周运指南4月2日—4月8日
别拿豆包不当干粮，剑走偏锋的又耐玩的游戏走一波吧~~
俗话说得好：别拿豆包不当干粮。在现在游戏动不动就十几或几十个G的今天，一些游戏则剑走偏锋，容量不大但是绝对可玩性极强。大有一种：别看我小我又技巧的架势。现在就让大家一起看看这些小游戏吧。
《饥荒》 这款由Klei Entertainment开发的一款动作冒险类求生游戏，玩家将以一名科学家的身份出现在神秘的异世界中，玩家则要开动大脑想尽办法在这个世界生存下来，并逃出这个异世界。
点评:上手简单，物品制作系统十分有趣，但单一的任务会让人感觉无聊。
《模拟山羊》
一款由Coffee Stain Studios开发制作的动作冒险类（ACT）游戏，这是一只山羊的故事，纳尼？破坏东西可以兑换积分，从此你就可以向其他同伴吹一辈子了，老子是这个世上最羊逼的山羊！
点评：当我第一次打开这游戏时心里的第一念头竟然是：这到底是什么鬼？这山羊也太厉害了吧。 游戏十分有趣有各种有趣的bug等你去发掘，楼主最喜欢在里面空中飞翔和日母羊了，那真是特酸爽
点评∶ 玩法新颖、自由度高、操作充满魔性，真的是 中二山羊欢乐多。
《以撒的结合：重生》
是一款动作角色扮演游戏。是以撒的的重做增强版，降低了难度，加了许多新的道具，更加灵活的电脑AI，怪物的造型难度也都有些改变。全新的手工绘制的像素风格的游戏画面，高度抛光的视觉效果，再加上全新的配乐和音频设计，这款新作带给玩家的体验将直接提升多个档次。
说实话这个游戏的剧情也太厉害了吧，以撒的妈妈要杀死这个以撒，面对母亲的屠刀，以撒展开一场精彩逃亡与反抗，游戏中以撒都是一副泪眼汪汪的模样，还挺萌的，武器就是自己的眼泪。这个游戏的怪物千奇百怪 ，极尽变态之能事，从大便、蛆虫、骷髅，一直到心脏、肉块甚至睾丸等等，都成为了能走会蹦的怪物，其中最让我恶心的是以撒妈妈的大毛腿，让我感叹这制作组也太没节操了吧。
点评： 整个游戏是哥特式风格，随机的地图，让游戏时时刻刻充满惊喜。道具千变万化，敌人造型也很有创意。游戏也有一定的难度，可算是手残党的克星。
火炬之光2》
是一部ARPG游戏，是《火炬之光》的续作游戏，由位于美国西雅图的Runic Games开发和《火炬之光1》代相比，《火炬之光2》在整体感觉上更像是一次冒险，你将看到广阔的室外世界，从岩石纵横的荒地到黄沙漫漫的沙漠，从天然洞穴到神秘地牢。所有的这些都让新作比原作更为吸引人，也更为丰富，充满惊喜。游戏中一共有4种职业，玩家可更具自己的喜好和战斗风格来选择，在装备方面，玩家可以通过多种方法强化自己的装备比如宝石、附魔等等，另外，寻找罕见宝石、套装和稀有武器也是玩家不断游戏的动力之一。如果你能找到与当前技能选择十分搭配的武器，那么你会发现战斗真的非常吸引人。不然也不会被IGN评为9.1的高分。
点评：硕果仅存的RPG之光，可算得上是浓缩版暗黑破坏神，配置也极其亲民，没玩过的朋友可以去体会下，你会感觉到刷刷刷果真是奇趣无穷
《闪克2》
（Shank2）是2010年横版动作游戏黑马《闪克》的续作，基于一代的创作灵感与基础，二代更进一步的增强了操作与画面，游戏采用以2D为主 部分的3D效果，音乐多以摇滚乐为主，操作更加灵活，在战斗上，玩家可以通过战斗而获取敌人的武器来使用，节奏十分爽快和血腥。武器道具也多种多样，的确是为横版ACT力作。
点评：横版中的鬼泣，其简约的动漫暴力美学，让玩家大呼过瘾。如果你是一个ACT暴力血腥斗士，那么一定要玩玩闪客,相信它一定不会让你失望。
《风卷残云》
是为一款3D动作过关游戏，以硬派武侠风格为基调、以主要水墨风格为主要特效渲染方式，让玩家轻松驾驭流畅多变的功夫技巧，勇闯重重难关。在游戏中玩家扮演一位被流放边疆的王子–罗劲，由于国家的连年战事，罗被封为将军并召回故土帮助平定纷争。期中会遇到各种各样的武林高手，学习到各种招数，提升技能，招募帮手，内容丰富且充满挑战！
点评：作为国产仅有的几个act游戏，真的是挺不错了，游戏上手简单，但是操作起来，你会发现难，真他妈的难。里面的小兵个个都战斗力爆表，对于我这种手残党来说简直是噩梦。游戏虽然有一定的难度，但打击感和连招系统做得不错，大赞。
《瘦长鬼影：降临》
是Blue Isle打造的一款恐怖FPS作品，本作使用Unity引擎制作，游戏将给予玩家身临其境的恐怖体验。或许某一天，在一个漆黑的午夜或是晴朗的白天，一个没有五官，身穿西服扎领带的神秘男子忽然出现在你的视野里。无论何时何地，或是在你的梦里，都会遇到这个人。他会时常出现在你的世界，逐渐摧残与折磨你的生活和神经，让你几近无法喘息。
直到有一天..........
你，消失了.......
听起来是不是很恐怖，玩起来那就更酸爽了。当初玩逃生都没吓到我，居然被这个游戏吓得毛骨悚然，游戏玩了第一关，心脏实在是受不了，就删了。这个游戏只有游侠网有汉化版，喜欢恐怖类游戏的朋友可以去尝试下，玩的时候一定要好好注意身后的鬼影
点评：恐怖爱好者的福音，胆小者慎入！
《骑马与砍杀》
骑马与砍杀系列是由土耳其的taleworlds开发制作，由paradox代理发行，该游戏最初仅由一对夫妇开发制作，虽然早先的版本十分粗糙，但是仍受到颇多好评，游戏的世界为架空的中古世界，玩家在以、游戏中可以扮演骑士、商人、盗匪等多种角色，甚至可以割据一方，直至统一卡拉迪亚大陆。内容丰富，自由度极高，同时众多制作精良的MOD也是该游戏的亮点之一。
这是一个非常真实的游戏，游戏系统包括了动作/战术模式和战略模式，在M&B中，玩家的游戏历程完全由玩家自行决定，你可以做一位除暴安良、为国王英勇作战的骑士领主，也可以做一个云游四方的商人，还可以辅佐某位被驱逐的大贵族帮他夺回他的王国，甚至可以劫掠村庄、攻击商队当一把“十恶不赦”的强盗。听起来是不是非常心动，成为战争领主，在中世纪的世界里面横行霸道，你值得尝试！
评价：代入感强，兵种的熟练度和武器的属性使得战斗非常接近于现实，游戏内容丰富而又充实，每次战斗都让人热血沸腾，电脑硬盘必备游戏
《双星物语2》
是Falcom研发的一款日式经典ARPG（动作角色扮演游戏），2008年9月25日在日本上市，2009年6月28日经北京娱乐通汉化在国内发行。《双星物语2》与前作一样采用了双人合作冒险模式，玩家可随时切换操作两位主角。男主角是“年轻的寻宝猎人”拉格那，擅长武器为链剑。女主角则是“吸血鬼族公主”爱尔温，善用魔法。游戏以天空中的浮游大陆为舞台，前作中吃食物来提升等级和携带宠物要素都被保留。还会有新的宠物登场亮相。
说起双星物语，印象最为深刻的就是清新亮丽的画风、那跟在主角身后的猫猫狗狗以及吃食物升级的特殊设定，刺激的BOSS战、搞怪的竞技场、个性化的计步器和怒击指示器。总而言之，《双星物语2》是一款时下难得的单机精品，不论你是否玩过一代，都会立刻被其水墨版纯净的画面所吸引。在单机游戏如此萧条的时代，暂时放下网游PK的血腥与杀戮，享受一次单纯的游戏乐趣，相信你也会在不经意间喜欢上它。
点评：很不错的小清新游戏，萌萌的画风，有趣的设定，总会让人不经意的爱上这个游戏。
《疯狂骑手》
是一款Techland开发Ubisoft发行的超爽快越野摩托车竞速游戏。玩家可以在全球45个丰富多彩的赛道上进行角逐，还可以做出各种各样精彩的空中特技动作。为了取得比赛的胜利，玩家必须合理地运用加速技能，还要充分利用赛道上隐藏的捷径。
赛车游戏都是以第一个冲过终点线为目的。而本作的亮点在于，融合了极限运动的动作元素，在高速行驶的过程中，会有很多悬崖需要飞跃，在空中就可以做各种动作，也会解锁很多成就。在飞跃的过程中，也会有很多意想不到的奖励，例如飞跃圆圈获取能量、落地进入靶心等等极限运动的小细节。极大的加深了挑战度与刺激感。
点评：作为一个1g多的赛车游戏，画面还是挺不错的，关卡众多，游戏玩法类似于跑跑卡丁车。
《雷曼：起源》
育碧曾经辉煌的《雷曼》系列近几年已经演变成疯兔那样的恶搞作品，忽略了正统的雷曼主角的动作冒险游戏，《雷曼：起源》以雷曼本身为主角的全新2D动作游戏。本作立足于回复系列传统的动作乐趣，该作包含超过60个2D关卡，游戏的故事模式支持最大四人合作通关，在游戏中玩家可以解锁新的技能，并使用新技能重新挑战之前的关卡以重新发现新的道路和隐藏道具等。
这款游戏可以说是当代最具艺术性的游戏之一，有着完美的想象力和永无止境的笑料。它给人的感觉就好像是：艺术总监Michel Ancel 用一张小小的餐巾纸，异想天开的画出了整个世界，并使其不断的扩张。游戏给玩家注入了一种生气，崩塌破碎的场景被茂密葱翠的场景所取代，搞笑的角色们嬉笑怒骂。尽管那幻想般的画面初看给人感觉一片混乱，但是它十分吸引玩家的眼球。
评测：游戏的背景音乐对于场景和气氛都烘托作用非常完美，画面实在太漂亮啦。还有超棒的挑战与成就元素，游戏有着极佳的控制手感。这个游戏我大概玩了好几个月，现在还没通关，楼主真是太失败了
《兽人必须死2》
是由Robot Entertainment开发的一款塔防类(TD)游戏。本作的剧情紧接着1的剧情发生，当主人公完成所有防御的第二天，发现自己突然又加入了与疯狂兽人们的战斗。新作引入了合作模式，玩家可邀请好友来一起猎杀兽人，支持双人合作。游戏也将包含新的战役模式，新的敌人等。
新作包含诸多新元素，例如重新设计的魔法兽人必须死2书，可建造更多的塔防堡垒，新的升级系统能让猎杀兽人更有趣。重新设计的魔法书可以让玩家能够有更多的选择去建造防御设施。升级了骷髅头的获取系统允许玩家使用自己最喜欢的方式去把兽人们烧死、碾死、砸死、肢解等等。
评测：游戏画面比之一代进步不少，不再像一代那样满屏狗牙了。结合了动作类型游戏的优点以及塔防游戏的策略性使得整个游戏的可玩性大大提高 ，陷阱(超过50种)可尽情虐杀兽人，这听起来挺不错，不是吗？
这些你都玩过多少？不玩可是错过很多有趣的玩意哦，经典就称不上了
65周年 | 李志伟：InSAR部分地学参数反演
本文内容来源于《测绘学报》2022年第7期（审图号GS京（2022）0495号）
I
nSAR部分地学参数反演
李志伟1
, 许文斌1,胡俊1, 冯光财1,杨泽发1,李佳1, 张恒2,陈琦2,朱建军1,王琪洁1,赵蓉3,段梦3
1. 中南大学地球科学与信息物理学院, 湖南 长沙 410083;
2. 中国资源卫星应用中心, 北京 100094;
3. 中南林业科技大学土木工程学院, 湖南 长沙 410018
基金项目：国家杰出青年科学基金(41925016);国家自然科学基金(42030112);湖南省教育厅资助科研项目(21A000)
摘要：自然因素或人类活动可使地球表层或内部应力发生变化, 进而产生灾害事件。获取灾害事件与灾害过程的关键地学参数对于准确理解灾变过程、科学解释灾变机制、正确拟定应对策略具有重要意义。InSAR技术已广泛应用于自然因素或人类活动导致的灾害事件与灾变过程的参数反演。本文首先介绍了InSAR卫星的发展和地表形变监测的基本原理; 然后分析了InSAR地学参数反演在地震、火山活动、地下水抽取、矿山开采、冻土冻融、冰川运动、地下流体迁移等各类潜在致灾事件中的研究现状; 最后分析了InSAR地学参数反演存在的主要挑战和问题。
关键词：星载SAR InSAR 地表形变 地学参数反演
引文格式：李志伟, 许文斌, 胡俊, 等. InSAR部分地学参数反演[J]. 测绘学报，2022，51(7)：1458-1475. DOI: 10.11947/j.AGCS.2022.20220156
LI Zhiwei, XU Wenbin, HU Jun, et al. Partial geoscience parameters inversion from InSAR observation[J]. Acta Geodaetica et Cartographica Sinica, 2022, 51(7): 1458-1475. DOI: 10.11947/j.AGCS.2022.20220156
阅读全文：http://xb.chinasmp.com/article/2022/1001-1595/20220728.htm
引 言
自然因素或人类活动可使地球表层或内部应力发生变化，进而导致地震、火山活动、滑坡、泥石流、基础设施损毁等一系列地质灾害。在地质灾害预警与防控体系中，获取灾害事件与灾害过程的关键地学参数是至关重要的。然而，受限于技术手段或观测空间，部分关键地学参数难以直接测量(如地震震源参数)。地质灾害事件的第一表征往往是地表形变，通过卫星遥感获取的地表形变数据间接反演地学参数是目前获取关键地学参数的主要途径之一。通过反演获得的关键地学参数对于准确理解灾变过程、科学解释灾变机制、正确拟定应对策略具有重要意义。受观测成本高、工作范围小、空间分辨率低等因素的限制，传统GNSS、水准等测量手段监测的地表形变通常是局部且离散的，难以宏观、精细地反映地表真实形变场，限制了地学参数反演的精度和可靠性。
星载合成孔径雷达干涉测量(interferometric synthetic aperture radar，InSAR)是一种主动微波遥感技术。自1989年，文献[1]首次提出差分InSAR(differential InSAR，D-InSAR)技术以来，InSAR已成为地表形变监测的重要工具之一[2-6]。相较于传统大地测量手段，InSAR在形变监测上具有全天时、全天候、空间上连续覆盖、分辨率高、数据处理效率高等优势。此外，InSAR可利用存档影像回溯地表历史形变(最早可到1991年)，从而为已发生的形变事件参数反演提供观测资料基础。鉴于此，InSAR已被广泛应用于自然因素或人类活动导致的灾害事件的参数反演。
为了更系统地阐述InSAR地学参数反演领域的研究进展，本文梳理分析了目前该领域的主要研究现状和挑战。
1 InSAR卫星及测量原理
1.1 InSAR卫星发展现状
作为一种新兴的对地观测手段，InSAR技术具有全天时、全天候、大范围、高精度、低成本的优势，在地震、火山、滑坡等地质灾害的快速响应，城市、矿区、人工建构筑物形变的长期监测，以及冰川运动、冻土冻融、地下水开采等活动的动态监测中得到了广泛应用。随着硬件技术的日益发展，几十年间InSAR观测平台实现了从空基向星基的跨越，多元化成像方式、多极化信号收发及高时空分辨率观测已成为SAR卫星发展的新趋势。本节将从在轨SAR卫星和即将发射的SAR卫星两个方面，对星载InSAR平台的发展历程进行介绍。
1.1.1 在轨SAR卫星
21世纪，航天工业的蓬勃发展极大地推动了SAR卫星技术的发展和进步。SAR卫星系统朝着广域覆盖、星座联合、超高分辨率、多模式、全极化的方向发展。由表 1可知，所有在轨SAR卫星均采用了多种成像模式，以此满足高分辨率或广域覆盖的监测需求。特别地，为满足广域覆盖需求，除小型SAR卫星ICEYE-X1/X2外，其他在轨卫星均具备扫描式成像模式，特别是Sentinel-1A/B宽幅成像模式采用的TOPS技术，能更好地满足常态化的广域地表观测需要。此外，为提高SAR卫星的时间分辨率，实现对突发事件的快速响应，很多SAR卫星采用了多星联合的星座运行模式。例如，RadarSat-2在与RCM串联组成星座后，最快4 d即可实现地面重访观测；在欧洲等重点监测区域，COSMO-SkyMed四星星座可实现24 h内的快速重访；TerraSAR-X/TanDEM-X星座采用单发双收模式，时间基线几乎为0，最大限度地降低了时间去相关和大气对干涉相位的影响。
表 1 在轨SAR卫星系统及主要参数
Tab. 1 Main parameters of the orbiting SAR satellite systems
1.1.2 即将发射的SAR卫星
当今时代，SAR卫星观测刈幅越来越宽，重访周期越来越短，全极化、简缩极化已成为多数SAR系统的标准配置，亚米级分辨率已实现常规化。星载SAR平台的巨大潜力和广阔前景引起各国广泛关注，政府、军方甚至商业遥感公司均在计划投资运营新的SAR卫星。除传统X、C、L波段的SAR卫星外，P和S波段的SAR卫星的研发也在持续深入(见表 2)。欧空局预计于2022年发射的BIOMASS卫星有望成为全球第1颗P波段SAR卫星；美国NASA与印度合作研发的NISAR将成为全球第1颗L和S双频段SAR卫星。在NISAR与德国宇航局的下一代TanDEM-L任务中，SAR凝视成像模式有望得到进一步发展。承继ALOS-2卫星、搭载新一代PALSAR-3传感器的ALOS-4，将具备更高的观测性能。
表 2 即将发射的SAR卫星系统
Tab. 2 SAR satellite systems to be launched
1.2 InSAR高程测量原理
DEM获取是InSAR技术的一个重要应用，通过两次在不同位置对同一地面区域进行观测，InSAR可以解算大范围的地表高程[7-10]。图 1为InSAR高程测量的成像几何。其中，A1、A2分别为两次成像的位置，B为空间基线，B‖、B⊥分别为平行基线和垂直基线，ρ1、ρ2为雷达天线到地面点的距离，α为基线倾角，θ为地面点的入射角，θ0为参考椭球面上点的入射角，δθ=θ-θ0，h为该地面点的椭球高，H为卫星高度。
图 1 InSAR高程测量成像几何
Fig. 1 Imaging geometry of InSAR elevation measurements
图选项
设Ψ1为地面点在A1处的观测相位，Ψ2为同一地面点在A2处的观测相位，则干涉相位Φinf由平地相位(Φflat)、地形相位(Φtopo)、形变相位(Φdef)、大气相位(Φatm)和噪声相位(Φnoi组成，可以表示为
(1)
在InSAR高程测量中，需满足两次成像期间没有地表形变发生、大气状态一致(如双天线或双星绕飞观测)，或地表形变、大气相位可以被很好地扣除或抑制的前提，根据成像几何关系，高程h可以近似表示为
(2)
式中，λ为雷达波长。通过式(2)可以计算地表相对高程，其中地形相位可以通过相位去平地处理并解缠得到。
2000年，美国国家航空航天局采用航天飞机携带C波段和X波段雷达干涉仪在11 d内对全球80%的地表进行了至少一次观测，获取了第1个InSAR全球DEM产品SRTM(shuttle radar topography mission)，提供30 m和90 m两种分辨率产品，相对高程精度约为10 m[9]。单轨双天线的设计使得相位测量中的环境误差被极大削弱。其生产的全球DEM数据已有多个版本，部分开放获取。
2007和2010年，德国宇航局与空客防务公司发射了TerraSAR-X和TanDEM-X卫星。两颗卫星参数配置一致，以250~500 m的基线距离进行双星绕飞，可实现交轨或沿轨干涉测量。该任务陆续发布了12、30、90 m分辨率的3种全球DEM产品，绝对高程精度可达10 m，相对高程精度可达2 m，很好地展示了星载InSAR获取高精度DEM的可行性和优越性[10]。
1.3 InSAR形变测量原理
1.3.1 D-InSAR基本原理
在重轨InSAR系统中，假设获取了两景覆盖同一地区的单视复数影像，通过配准和干涉生成一幅干涉图[11]，其干涉相位则同时包含地表高程和形变信息，可以表示为[12]
(3)
式中，dlos代表LOS(line of sight)向的形变量；等号右侧分别表示平地相位分量、地形相位分量、形变相位分量、大气相位分量、噪声分量及整周相位，其中k为整周模糊度。
D-InSAR技术是对形变相位分量除外的其他相位分量进行估算并去除，或通过滤波等抑制至很低水平以获取dlos。地形相位分量和平地相位分量可以利用外部DEM产品提供的较高精度高程h和准确估计的基线B从干涉相位去除[13]；通过多视处理和滤波可以抑制失相干噪声φnoise；通过二维空间解缠方法可以估计整周模糊度k[14]；大气相位分量可以通过外部数据辅助或滤波方法进行改正[15]。
1.3.2 MT-InSAR基本原理
为了克服D-InSAR技术的局限性，MT-InSAR(multi-temporal InSAR)技术应运而生，凭借能够有效抑制失相干噪声、大气延迟、高程残余误差及提供时间序列形变等优点，一系列时序算法相继被提出，如永久散射体干涉法(persistent scatterer，PS)[16-18]、短基线集方法(small baseline subsets，SBAS)[19]、部分时域相干点法(temporarily coherent point, TCP)[20]、干涉点目标分析法(interferometric point target analysis, IPTA)[21]及分布式散射体干涉(SqueeSAR)[22]等。MT-InSAR技术对选取的反射相位稳定点进行建模和时序分析，能够更加精确地获取地表形变速率和时间序列。其中，相干点目标选取、形变参数估计及大气误差改正是MT-InSAR的3大关键技术环节。
相干点目标选取在MT-InSAR技术中占据非常重要的地位[23-24]，其质量和密度直接影响形变参数的提取精度。目前主流的方法主要有相干系数阈值法、振幅离差指数法、相位稳定性分析法、偏移量估计法等。选取相干点目标后，将对相干点的形变参数进行估计，MT-InSAR形变参数解算的关键就是如何解算或避免相位模糊度[25-26]，目前已有的方法主要包括最小二乘法、解空间搜索法、LAMBDA方法、三维解缠法、模糊度探测法等。随着MT-InSAR技术的发展，大气延迟误差改正的精度决定着InSAR形变监测的精度，目前主流方法主要包括采用外部数据辅助法(如GNSS、MODIS、MERIS及数值天气模型等提供的大气数据[15])和时空滤波法[19]进行改正。
2 InSAR参数反演
2.1 地震震源参数反演
地震震源参数是指地震发生时震源处的一些特征量(如空间位置、强度等)或震源物理过程的一些物理量(如发震断层的几何特性、地震位错的空间展布、震源时间函数等)。震源参数是制定地震应急救援策略和开展救援工作的关键震情信息，也是研究发震断层机制及其破裂特征、分析区域应力触发与危险性评估的重要基础。随着卫星遥感技术的日新月异及处理能力的不断提高，震源参数反演由基于单一观测数据的反演逐渐发展为基于多源观测数据的联合反演，大大推动了震源参数反演进程。水准测量、GNSS、光学、InSAR等大地测量形变数据是震源参数反演的重要数据源。利用InSAR数据进行地震震源参数反演源于文献[5]，利用两景ERS-1 SAR影像对1992年美国加州Landers Mw7.3级地震的同震形变进行研究。随着InSAR技术在同震形变监测方面的理论与方法日趋成熟，国内外越来越多的学者开始利用InSAR地表形变反演震源参数。文献[27]首次利用二轨D-InSAR技术提取的视线向同震形变场为数据约束，基于Okada模型[28]反演得到了Landers地震的震源机制解。文献[29]率先在国内利用InSAR数据反演了1998年张北-尚义Mw6.2级地震的震源机制解。迄今，InSAR技术已成功应用于诸多同震、震后、震间地表形变监测与震源参数反演案例中[30-34]。
大地测量震源参数反演的方法主要包括解析法与数值法两大类。其中解析法一般适用于反演线性问题和简单的非线性问题，震源参数反演研究主要集中在数值模拟算法方面。对于由断层错动而引发的构造地震而言，震源参数即为断层参数。通过对大地测量数据、地震波资料等进行地球物理反演，即可获得更具细节的震源描述，包括发震断层的长度、宽度、深度、倾角、走向角、震源位置等几何参数，以及滑动角和滑动量等运动参数。其中，断层几何参数的非线性反演通常采用模拟退火算法、遗传算法、神经网络法、粒子群算法、贝叶斯法等搜索算法[35-36]，断层运动参数的线性反演通常可以采用附加光滑约束的非负约束最小二乘法[37]、总体最小二乘法等进行求解。图 2为2021年青海玛多Mw7.4级地震的震源反演结果[34]。利用Sentinel-1升降轨SAR数据获取了该地震的视线向和距离向同震形变场，采用四叉树降采样算法对形变数据进行了降采样，构建了变走向及倾角的6段断层几何，引入了二阶拉普拉斯平滑算子进行平滑约束，同时采用L曲线法和格网搜索法分别确定了最优平滑因子和6个断层分段的最优倾角，最后采用非负最小二乘方法求解了该地震的最优同震滑动分布。
图 2 青海玛多地震的震源反演结果
Fig. 2 The earthquake source inversion results of Maduo earthquake
图选项
InSAR震源参数反演中，简化的断层几何、断层面的离散化方法与平滑因子的选取、约束数据的多源性、数据定权方法与降采样方法及位错模型的选取(如弹性半空间位错模型与球体分层位错模型)等，是制约地震震源参数反演精度的重要因素。另外，Okada位错模型是基于各向同性、均匀半空间体假设的，是对真实地球空间结构的粗糙近似，并未考虑地壳横向和纵向分层结构，这也会导致断层滑移参数估计不可避免地存在偏差。如何准确确定不同观测数据集之间的相对权比与平滑因子、发展顾及地壳横向分层结构的地震震源参数反演仍然值得进一步研究。
2.2 火山活动参数反演
InSAR可以监测火山浅部岩浆运动、岩浆补给和亏损引起的地表变形。通过地表形变和模型反演得到的岩浆囊参数有助于识别火山喷发的前兆，跟踪火山喷发的演变，能够深化对火山岩浆活动的理解，有利于减轻火山活动及伴随的次生灾害带来的危害[3]。Mogi模型是最早用于解释火山形变的几何模型，由日本火山学家Mogi于1914年提出并用于Sakurajima火山的重复水准观测建模[38]。文献[6]首次利用InSAR技术发现了埃特纳火山的收缩现象，并利用Mogi模型反演了岩浆囊的位置和体积变化。Mogi模型假设在地壳的弹性半空间中嵌入一个球形源(源半径α远小于源深度d)，球形源内静水压力变化产生的地表位移和体积变化为
(4)
式中，u、v、w为地表(x，y，0)处的位移；球形源的腔体中心为(0，0，-d)；
，为腔体中心到自由表面的距离；ν为泊松比，一般设置为0.25；ΔP表示腔体内部的压力变化；G为剪切模量。ΔVcavity实际是压力增量对弹性半空间所做功的度量，只考虑了周围弹性半空间的力学性质，并不能在真正意义上等同于可压缩流体(如岩浆)的注入体积。随后，文献[28]给出了一套完整的有限矩形位错模型(Okada模型)的解析表达式。文献[39]推导了有限长椭球源(Yang模型)的精确解析解，并成功用于基拉韦厄火山的参数反演。基于上述Mogi、Okada及Yang这3种基本模型，又发展出了扁圆盘状、封闭(或开放)的管道，以及具有更高自由度的复合位错模型等模型[30, 40]。这些简单、有效的火山活动解析模型被广泛应用于板块主动分裂边界的火山运动系统的研究，包括岩浆囊参数、岩脉入侵过程、岩浆喷发时空分布、火山喷发与断层运动的耦合关系[41-43]。图 3为2014年冰岛Bardarbunga火山岩墙入侵与断层滑动模型[42]、2017年埃塞俄比亚Erta Ale火山喷发后引起的基台形变模型及该火山浅层内部结果示意图[41]。
图 3 火山活动参数反演[41, 43]
Fig. 3 Volcanic activity inversion results[41, 43]
图选项
然而，解析模型是对研究对象进行高度抽象后的建模，因此揭示的规律也相对模糊。数值模型为更准确地理解火山事件的构造过程提供了可能性。数值模型可以将岩浆房建模为有限大小的任意形状及由不同力学性能的岩层组成的复合火山，还可以引入各种类型的断层，使结果更逼真。在数值模型中，容纳岩浆囊的地壳被分为许多离散单元，采用适当的本构方程(描述岩石的流变性质)，结合质量守恒定律(连续性方程)和动量守恒定律(动量方程)建立微分方程，通过联立代数方程近似求解。目前最常用的数值模型求解方法是有限元法和边界元法。此外，构建多种地球物理手段结合的火山监测系统有助于更清晰地呈现岩浆系统的构造形态和运动状态；综合岩石学、地球化学和地球物理的观点处理有关岩浆管道系统的问题也有利于加深对岩浆系统演化的认识[44]。
2.3 地下水文学参数反演
地下水超采导致的地面沉降、地裂缝、海水侵入等灾害，影响工农业生产和城市可持续发展，造成供水、粮食、生态安全等一系列问题。对地下水资源进行科学管理和开展相关防灾减灾工作，需准确估计含水层的参数。InSAR能够获取高分辨率地表形变场，建立地表形变与地下水变化的函数模型，探测地下水变化的物理特征及反演相应的含水层参数，为地下水文学应用提供了崭新的视角[45-46]。文献[47]首次利用InSAR技术和含水层压缩模型证明了加州羚羊谷地下水开采与地表沉降间的耦合关系。随后，在拉斯维加斯、内华达州等具有良好观测条件的区域相继开展了含水层地表形变及参数反演的InSAR研究[48]。在我国华北平原、汾渭盆地、长江三角洲等重点区域，InSAR技术被用于分析地下水开采引起的地表沉降的时空分布特征和形变规律[49]。利用InSAR观测的季节/年际形变数据和一定数量的地下水测井数据，结合土力学、应力应变等理论模型，可识别和反演一系列地下水文学参数[50-53]：①含水层的边界和内部特征；②含水层储量(含水层厚度、储量)；③含水层贮水参数(贮水度、弹性/非弹性骨架贮水系数)；④含水层响应状态(渗透系数、滞后系数)；⑤地下水载荷应力变化；⑥水头及采水量估计。这些参数可应用于：①提取含水层空间-时间分布的异质性，弥补水位测井网络的分布；②校正已有的地下水流动模型，用于地下水流动模拟及地表沉降预测；③评估地下水储量变化，及其对构造活动和地震风险的影响；④精化区域含水层和岩石圈结构。
如前文所述，利用InSAR开展地下水文学参数的反演研究，可以获得高空间分辨率的地下水文学参数，有效弥补地面观测数据不足难以估计空间上详细的水文参数的问题，对探测含水层物理结构、校正地下水文学模型及开展地面沉降模拟具有良好的应用前景。然而，利用InSAR观测反演地下水参数仍存在以下两大问题。
(1) 承压含水层模型简化对参数反演的影响。目前，以InSAR地表形变作为观测值反演水文学参数主要基于Terzaghi有效应力原理所导出的线弹性模型[41]。在数值求解中一般采用“两步法”，即将三维地下水流动和一维垂向固结过程进行分步计算(如MODFLOW-SUB方法)[40]。文献[54]研究表明，这一类简化土体骨架水平位移模型，在含水层边界或采水井周边等存在水平形变的区域会造成明显误差，并且利用线弹性模型进行参数反演时，往往忽略含水层参数随时间的变化特性(如渗透系数随有效应力增加而减小)，造成反演结果无法准确反映含水层的动态变化。除了线弹性形变，含水层土体也可能存在塑性和黏弹性等非线性形变(如弱透水层的塑形形变)[55]，如何准确结合InSAR观测结果与不同本构关系，是反演地下水参数的关键问题。
(2) 不同地质结构对参数反演的影响。除了地下水储量变化对承压含水层带来的直接影响，断层结构与不透水层/半透水层的夹层结构等低渗透率地质单元的分布会降低含水层系统的地表响应，从而减弱InSAR对含水层形变的捕获能力，增加参数反演的不确定性[40, 44]。文献[44]利用InSAR资料估计了边界断层的平均水平传导率，然而由于含水层与断层结构的非均一性，在断层两侧存在明显的拟合残差。文献[56]利用贝叶斯估计方法和最大后验概率的准则，基于单个测井数据和InSAR观测资料，成功反演了含水层渗透率的侧向空间分布，为识别含水层中不同地质结构提供了新的思路。但采用贝叶斯估计和最大后验概率的方法计算效率不高，且未对多变量进行联合反演，限制了对含水层中不同地质单元影响的评估与模型的准确性。
如何融合InSAR水平形变(如含水层边缘、采水井周边)与不同的含水层变形机制进行地下水文学参数反演，获取含水层的动态响应、应力场各向异性等参数，以及如何分离和重建微小的地表响应信号，联合反演获得不同场景下水文参数的时空分布，仍值得进一步研究。
2.4 地下采空区几何参数反演
地下矿产资源开采后留下的采空区容易引发一系列灾害事故。例如，采空区上覆岩体因失去支撑容易产生移动形变，从而引发地表塌陷、基础设施破坏等地质灾害；上覆岩体破裂后，地表或地下水容易通过裂缝流入废弃采空区并形成“老窑水”，一旦“老窑水”在未知情况下被导通，极易引发矿井突水等重大安全事故。因此，精确探测或反演地下采空区几何参数(如位置、深度、长度、宽度、厚度、倾角、方位角等)，对于减少采空区灾害事故至关重要。传统方法大都以物探等接触式手段为主[57]，虽精度较高，但工作范围小，且费时费力。
矿山地表移动形变与地下采空区几何参数密切相关[58]。因此，若能精确监测地表移动形变，理论上便可反演地下采空区几何参数。如前文所述，InSAR可低成本、大范围、高空间分辨率地获取地表形变值，其在地下采空区几何参数反演中具有较大的应用潜力。近年来，国内外学者在该领域展开了探索，并发展了多种算法。根据反演依据不同，算法大致可分为两大类。
(1) 联合InSAR观测与形变几何特征反演法。2013年，有学者首次以InSAR视线向形变观测值为基础，根据煤矿区地表沉降盆地中心位置、沉降梯度变化等形变几何特征反演了大范围地下采空区的中心位置，但该研究未能反演深度、长度、宽度等其他关键几何参数。为克服该局限，文献[59—62]提出引入开采沉陷盆地半径、拐点偏移距、开采影响传播角等变形几何特征参数反演地下采空区深度、长度、宽度等其他几何参数。该方法计算工作量小，但因其仅利用了部分(如中心点、走向或倾向剖面等)InSAR观测值，导致采空区几何参数反演精度对InSAR观测误差较为敏感。
(2) 联合InSAR观测与理论模型反演法。2018年，文献[63]引入概率积分法(我国最常用的开采沉陷理论预计模型之一)，建立了采空区几何参数与InSAR视线向形变的函数模型，并通过非线性搜索率先实现了基于InSAR观测值的地下采空区几何参数定量反演。图 4为安徽某矿区实测采空区几何参数(蓝色矩形)与InSAR反演的几何参数(红色矩形)。由图 4可看出InSAR实测形变与基于反演参数和概率积分法模型预测的形变的差异。2021年，文献[64]引入反向学习法提高了采空区参数反演的效率和精度。相较于联合形变几何特征方法，联合理论模型反演法能充分利用InSAR整体观测(对误差敏感性较前者低)，且理论基础更坚实；但因该方法涉及大量非线性搜索，反演效率不佳。
图 4 安徽某矿区实测采空区几何参数与InSAR反演的几何参数[63]
Fig. 4 The in-site measured goaf parameters and the InSAR inverted parameters in a mining area in Anhui province[63]
图选项
总体而言，虽然两类算法在一定程度推动了采空区几何参数InSAR反演研究，但目前的反演精度都不高。此外，现有研究大都集中于煤矿规则采空区(如长壁采煤)，而甚少关注非规则工作面参数反演。因此，面向更复杂的开采场景，并提出全新算法以提升采空区参数反演精度，是未来的重要研究方向。
2.5 冻土活动层厚度反演
随着全球气候变暖，多年冻土已经发生明显且持续的退化，对碳循环、水循环、生物多样性、工程安全等产生了重大影响。多年冻土区活动层是多年冻土与大气圈层进行物质交换的纽带，它随着气温的变化冬季发生冻结、夏季发生融化，引发地表产生周期性形变[65-70]。活动层厚度是冻土的重要物理参数，其变化直接反映多年冻土的存继及其退化或变化情况[65, 70]。大面积获取和反演活动层厚度对冻土科学研究、冻土区工程建设和防灾具有非常重要的意义。传统获取活动层厚度的方法主要有实地测量方法、解析模型反演法[70-73]。随着遥感技术的发展，基于遥感数据的深度学习方法[73-74]也可被用于估计冻土区活动层厚度，但受不同条件的制约。近年来，具有高分辨率、高精度、大范围、长时段优势的InSAR技术可以准确获取冻土区地表形变，为大范围、高精度反演冻土物理参数(如融化深度、活动层厚度)提供了观测基础。国内外学者通过建立冻土区InSAR地表形变与活动层厚度之间的函数关系，估计多年冻土活动层厚度，相关研究大致可以分为以下两类。
(1) 基于活动层变化与InSAR形变观测量之间的函数关系。文献[66]采用时序InSAR技术探究了冻土区地表形变与活动层厚度之间的关系，其假设冻土区地表沉降形变是由活动层融化引起的，并基于水质量守恒定律，建立了形变和活动层厚度之间的积分函数模型，通过逆运算可以求得研究区域的活动层厚度。但该模型水质量守恒的假设在实际中并不能满足。文献[68, 75]简化文献[66]模型中的土壤孔隙、土壤湿度饱和度，用于反演青藏高原多年冻土区的活动层厚度。文献[69]建立了活动层土壤水分含量模型，得到不同深度土壤水分含量的分布，将其代入文献[66]反演模型中估算活动层厚度；但该方法在建立冻土土壤含水量沿土壤深度方向的分布模型时，依赖实测的土壤含水量数据，影响了活动层厚度反演结果的精度。除此之外，文献[76]在文献[66]模型基础之上，采用L波段InSAR结果和P波段极化SAR数据联合估计土壤湿度和活动层厚度，解算过程实际为求解联合目标函数最小值问题。该方法减少了土壤湿度信息不准确对活动层厚度估计的影响，但目前P波段数据获取困难。文献[67]将工程应用中采用的土壤分层模型引入活动层厚度反演研究中，提出了一种基于InSAR形变与土壤分层模型反演活动层厚度的方法，该模型假设冻土区地表融化形变由两部分组成，即冻土融化产生的下沉量和压缩变形，在反演过程中土壤的热参数采用经验值，并且假设活动层下多年冻土刚性不可压缩，因而导致反演精度对这些假设较为敏感。
(2) 基于土壤一维热传导方程。文献[68]提出了基于InSAR形变及土壤一维热传导方程反演活动层厚度的方法。该方法首先利用时序InSAR获取多年冻土区的地表时间序列形变数据，分析冻土时序形变与温度之间的时间延迟关系；然后根据土壤一维热传导方程推导并建立活动层厚度与延迟时间之间的函数关系；最后利用冻土区融化沉降形变曲线计算地表温度和融化深度达到最大时的延迟时间，进而可以估算多年冻土区土壤的最大融化深度，即活动层厚度。该方法间接建立了InSAR观测值和冻土活动层厚度之间的函数关系，反演精度对InSAR形变结果的变化趋势依赖性较高。
基于InSAR技术反演活动层厚度具有全天时、全天候、大范围、高空间分辨率等优势，在冻土物理参数反演研究中具有很好的应用前景。但通过建立活动层厚度与冻土区InSAR形变量之间的函数模型反演活动层厚度，受InSAR形变结果为相对形变量的影响，反演结果多为相对活动层厚度，其结果的精度受形变测量精度的影响。建立活动层厚度与冻土形变时间延迟(相对温度)之间的函数关系反演活动层厚度的方法，可以得到绝对的活动层厚度，但其精度受模型采用的土壤热参数经验值精度及只考虑土壤热变化的影响。实际中，冻土变化机理与过程复杂，如何准确地建立InSAR观测量和活动层厚度的函数关系，是限制InSAR反演冻土相关物理参数的关键问题。在未来研究中，利用InSAR技术监测获取准确的活动层厚度需要全面考虑实际冻土变化的机理，实现高精度InSAR冻土物理参数可靠反演方法仍需深入研究。
2.6 冰川跃动参数反演
近年来，高亚洲发生的一系列与冰川跃动相关的灾害引起了全世界对气候变暖背景下冰川动态的关注。冰川跃动是由冰川内部压力状态发生突破性变化而触发的[77]。在不同的内外因素作用下，跃动会使冰川发生多种几何形态变化，包括冰面高程变化、流速变化、边界变化、表面形态特征变化[78]。观测这些几何参数的时空演化是反演冰川跃动持续时间、转移物质当量、触发机理等信息的关键所在，而这些跃动信息又是灾害防范和应急救援的科学基础。冰川边界、表面形态特征变化可以通过光学遥感分类解译的方法获取，方法相对成熟。观测的难点在于冰川的表面高程变化和流速变化。利用InSAR技术估计冰川表面高程变化主要有两种方法：一是利用多轨SAR数据估计冰川三维形变，然后结合冰面局部坡度估计冰川表面高程变化[77-78]；二是利用InSAR技术生成多期冰面DEM，对多期DEM进行配准后再进行差分[79-83]。DEM差分法依赖于InSAR DEM的生产能力，但应用更广。由于冰面散射特征变化较快，SAR冰川流速变化估计一般采用偏移量跟踪技术(offest tracking)，即影像匹配技术，包括斑点匹配(speckle tracking)[80]和强度匹配(intensity tracking)[81]。以SAR影像匹配技术获取的多轨多时相像素偏移量为基础观测值，通过时序形变处理和三维形变模型可以反演冰川三维时序流速[82-83]。
为提高冰川跃动参数反演的可信度，学者们一般综合利用InSAR和其他手段(光学遥感和激光测高等)的多项观测结果。文献[84]率先利用SAR影像匹配估计的冰川流速变化和光学影像解译的冰川形态变化，反演喀喇昆仑地区冰川跃动机理，认为该地区冰川跃动是由内部温度环境的改变引起的。此后，文献[85]利用同样手段获取的冰川流速变化和形态变化，反演了西昆仑山地区冰川跃动机理，认为该地区冰川跃动行为受到了热力学条件改变和内部水文系统季节性变化的共同作用。文献[86]结合SAR影像匹配和光学影像匹配获取的冰川二维流速序列和InSAR DEM差分技术估计的冰川多时段表面高程变化，反演了帕米尔高原Bivachny冰川近期跃动物质转移量、跃动体移动速度及跃动触发机制，认为该冰川跃动由内部水文条件变化引起。文献[83]通过基于匹配偏移量的SAR多维小基线集技术估计喀喇昆仑山Hispar冰川的三维时序流速，同时利用InSAR DEM差分技术估计该冰川的多时段表面高程变化，再结合光学影像匹配技术估计的冰川二维时序流速、光学遥感解译获取的冰川表面形态变化和冰川边界变化，反演了该冰川近期跃动的起止时间、跃动体分布范围、跃动体移动速度、跃动中物质转移量及跃动触发机制，结果表明该冰川跃动主要由冰川底部水压突破极限而触发，北部支流先于主干开始跃动。需要注意的是，光学遥感的星下点成像方式相比SAR的斜距成像方式受地形影响更小，在地形特别复杂的地区应该将两种手段的优势互补，以实现长时序高精度冰川流速和高程变化观测，更好地服务于冰川跃动参数反演和跃动触发机制解译。
2.7 大气水汽参数反演
水汽是地球大气中变化最快的组分，在大气辐射平衡、能量转换、云的形成，以及降雨、气候研究等过程中发挥着重要的作用[87]。高分辨率精确获取大气水汽的空间分布对大气动力学的研究也十分重要。InSAR技术作为一种新的大气水汽估计手段，可以弥补传统气象学水汽监测方法成本高、空间分辨率低、受天气影响较大等不足。文献[88]率先发现InSAR干涉图可以提供空间分辨率高至米级的高精度差分大气水汽，从而开启了InSAR气象学研究。更进一步，文献[89]推导了常用SAR传感器反演差分水汽的理论精度，计算结果表明，该理论精度可达1 mm。然而，InSAR获取的是主、辅SAR影像成像时刻对应的水汽差分值，难以直接应用于中尺度气象学和数值天气预报等。
文献[17]通过时空域滤波的方法首次将大气水汽作为副产品从InSAR干涉图中估计出来。虽然时空域滤波的方法被广泛应用，然而目前还没有研究可以提供统一的标准以实现不同试验条件下最优滤波器和滤波窗口的自适应选择。随后，文献[90]创新性地利用最小二乘配置方法将大气水汽作为参数，与地表形变同时进行估计，然而该方法仅适用于稳定且渐进的地表形变。在忽略地表形变影响的前提下，文献[91]通过假设湍流大气在时序上的统计均值为零，获取了莱茵河上游地区17景ASAR影像成像时刻的非差分水汽分布结果。虽然这种附加假设条件的方法便于计算，但其需要大量SAR数据确保统计意义，且反演的水汽结果精度受极端天气的影响较大。为了克服这些限制，文献[92]以时序水汽均值不变作为约束条件，获取了沙特阿拉伯地区的5景SAR影像成像时刻对应的非差分水汽，试验结果如图 5所示；同时证明了利用时序水汽均值为零的约束条件反演的水汽结果在时序上整体存在一个固定的偏差，且该偏差等于时序上的水汽均值。虽然时序水汽均值不变的约束条件可以确保反演的时序水汽的整体精度，然而由于大部分地区难以在反演之前获取研究区域时序上的水汽均值，因此一定程度上限制了该方法的应用。除此之外，融合各种外部气象数据(如GPS、MERIS)和气象模型(如WRF模型)的InSAR非差分水汽反演方法被不断发展和优化[93-95]。然而，目前外部气象数据的空间分辨率普遍较低，且受云、雨等条件的影响较大，该方法反演的非差分水汽精度直接取决于输入数据的质量。
图 5 基于时序水汽均值不变约束条件下的吉达地区的非差分水汽反演结果[92]
Fig. 5 Undifferenced water vapor inversion results in Jeddah region with the constraint of time series water vapor mean being constant[92]
图选项
尽管InSAR技术已经成功反演了特定试验区的大气水汽，但是如何发展一种通用的、高精度的水汽反演方法仍然是InSAR在气象学研究中亟须解决的瓶颈问题之一。随着SAR卫星的不断发展，可用于水汽反演的InSAR干涉图越来越多，如何在充分利用海量InSAR数据的基础上，结合数理统计的知识较好地解决时序InSAR非差分水汽反演中的秩亏问题，是未来重要的研究方向之一。同时，伴随着获取的非差分水汽数据的不断丰富，如何将这些高精度、高分辨率的水汽数据与气象资料同化，提高气象预报特别是降水预报的准确性，是未来InSAR气象学研究的重要方向。
2.8 地下流体体积变化参数反演
地下流体(石油、天然气、地热等)的承载、卸载、运移及开采等会引发显著的地表形变，从而导致地面沉降、坍陷等地质灾害，已经成为影响区域经济和社会可持续发展的重要因素。国内外学者为了准确理解地下流体灾变过程并科学解释地下流体灾变机制，将地下流体与地表作用机理进行耦合研究，构建出一系列物理解析模型，将地表形变与模型相结合，以反演地下流体几何参数、体积变化参数及静力压强参数，已取得较好的应用效果。其中，地下流体的体积变化可以解译地下流体活动、评估地下流体活动造成的灾害严重程度等，是地下流体反演中最重要的参数之一。
反演地下流体参数需要精确获取地下流体活动导致的地表形变，而InSAR技术能低成本、大范围、高空间分辨率地获取地表面域形变值，因此在地下流体参数反演中具有巨大应用潜力。国内外学者在地下流体体积变化参数反演领域构建了多种模型，大致可分为如下两大类。
(1) 基于确定几何形状的地下流体体积反演。基于地下流体的先验信息(勘探资料、地质构造等)，将地下流体真实形状近似为一个确定的几何形状，从而进行单一源的体积反演，获得地下流体体积变化。这种反演方法称为基于确定几何形状的地下流体体积反演。简化为简单的几何形状，降低了参数反演的复杂程度。此类反演模型包括Mogi模型、Okada模型、CDM(compound dislocation model)模型、Yang's模型和扁椭球模型等。1995年文献[6]首次利用InSAR技术和Mogi模型发现了埃特纳火山的收缩现象，并获得了埃特纳火山岩浆囊位置及体积变化。2017年，文献[40]发展了一种自由度更高的CDM，并应用于卡尔布科火山，基于InSAR观测实现了地表三维形变及岩浆囊体积变化反演。2021年，文献[96]利用InSAR技术监测呼图壁储气库注/采气过程，并引入CDM模型实现了地下储气库的体积变化反演和地表三维形变估计，为我国战略储气库的安全监测和体积反演提供了科学参考。这些模型在进行地下流体体积反演时依赖于先验的源形状，因此，假定的与真实的地下流体几何形状是否一致是地下流体反演的关键。然而，实际中地下流体几何形状一般是不规则的，难以获得准确的反演结果。
(2) 基于非确定几何形状的地下流体体积反演。在无须先验假定地下流体几何形状的前提下，基于弹性半空间理论将地下流体储层分为若干个均质单元(如平行六面体、柱体)，先对每个单元分别进行地下流体体积反演，最后通过体积可叠加原理获得整个地下流体体积变化。这种反演方法被称为基于非确定几何形状的地下流体体积反演，可以减少由于先验信息不足导致的反演精度不够的问题。此类反演模型包括柱体模型(Geertsma's模型)、分布式点源模型等。2002年，文献[97]首次将InSAR数据用于分布式点源模型的地下流体体积参数反演中，获得了Coso地热田的地下流体体积变化，InSAR数据的使用极大地改善了地下流体参数反演的准确性。随后，文献[98—99]利用该模型与InSAR数据分别反演了加利福尼亚州威尔明顿油田的体积变化和与黄石火山温泉的体积变化。2011年，文献[100]提出了一种将分布式点源模型与椭球体模型相结合的自由几何体建模方案，该方案联合InSAR数据和重力数据，实现了同时反演火山口下方岩浆囊的体积变化分布与重力压力变化。2016年，文献[101]利用InSAR数据、分布式点源模型及PCA技术，实现了基拉韦厄火山口长时序体积变化，以及火山喷发期和休眠期不同深度的体积源分布。2017年，文献[102—103]在分布式点源模型的基础上，提出了基于地下流体模型约束的三维形变模型，将体积变化参数与地表三维形变联合求解，并成功反演了涩北气田的动态体积变化。
以上反演模型皆为高度非线性模型，在实际应用中需假定地下流体几何参数(如地下流体深度、厚度等)，使非线性问题转化为线性问题，这些简化在地下流体反演参数中会导致一定的不确定性。因此，通过利用非线性反演算法(模拟退火法、粒子群算法、遗传算法等)搜索得到模型的最优反演参数将是未来研究的重点。其次，地下流体活动与断层活动导致的地表形变常常出现混叠现象，如何从形变结果中区分断层活动地表形变与地下流体活动引起的地表形变，从而获得更准确的地下流体体积变化，未来也需要进一步研究。以上基于弹性半空间理论的解析模型是高度抽象化后的模型，因此揭示出的地下流体体积变化规律相对不够全面，而数值模型(有限元法和边界元法)可将地下流体形状建模为有限大小的任意形状并结合不同的力学性质，采用适当的本构方程(描述地下流体的流变性质)联立求解地下流体参数。因此，如何高效准确地建立地下流体反演的数值模型，也是未来进一步研究的方向。
3 InSAR参数反演面临的问题和挑战
3.1 InSAR观测误差与单一观测几何制约参数反演精度
3.1.1 低相干区InSAR观测降低参数反演精度
利用InSAR技术可以获得高空间分辨率的地表形变信息，可为与地表形变相关的潜在致灾过程参数反演提供可靠的数据源。然而受时空去相干、观测误差等因素的限制，一些InSAR误差源很难被准确地抑制或去除，从而制约了参数反演的精度。断层区域的植被、水域覆盖往往会导致InSAR相干性降低，严重的失相干会导致InSAR技术无法获取有效的地表形变监测结果。受两次SAR成像时刻电离层和对流层的温度、相对湿度、气压等大气物理参数发生时空变化的影响，InSAR差分干涉相位信息中残留了大气相位，很难对其进行消除，严重制约了InSAR技术监测地表形变的精度和后续参数反演精度。其次，卫星轨道误差、DEM误差、数据解缠误差及数据处理过程中引入的误差等也会影响地表形变获取的精度和参数反演精度。
3.1.2 InSAR一维形变导致部分反演参数不敏感
除地学模型的准确和可靠外，观测值能否全面描述地表形变的关键特征也是地学参数精确反演的关键。受SAR斜视成像的影响，单一平台或轨道InSAR仅能监测地表真实三维形变在雷达视线方向的一维投影，难以反映地表真实形变特征。因此，仅利用InSAR一维形变观测值反演地学参数可能出现部分参数对观测值的不敏感现象(即观测值难以有效约束相关参数)，降低参数反演的精度和可靠性。例如，在目前SAR卫星近极地轨道场景下，地表南北方向形变分量对InSAR一维视线向形变贡献较小。这就意味着，与南北向形变密切相关的模型参数对InSAR观测不敏感，该模型参数反演的误差会较大，可靠性不高。鉴于此，在InSAR数据允许的前提下，充分利用多视角、多平台及多轨道InSAR形变观测值，有助于提高地学参数反演精度和可靠性。
3.1.3 InSAR几何条件相近引起的病态问题
不同卫星成像几何条件相近，会大幅降低观测的图形强度。此外，同一卫星获取的InSAR形变观测数据，也不是完全独立的。这些因素会造成地学参数反演中法方程存在一定的病态问题，导致反演的地学参数不稳定、不精确。目前，主要通过附加约束(如正则化等)稳定待求参数克服病态问题。不同地学反演模型，参数数量及参数间可能存在的相关性也会造成模型求解方程病态，需要通过线性化或有效的非线性估计方法减少模型待估参数，从而减弱病态问题对参数反演精度的影响。针对大地测量地球物理反演中的病态问题，文献[104—105]提出了很多解决方法(如谱修正迭代法、正则化法、岭估计、主成分分析、贝叶斯方法等)，在满足结果精度的要求下，评估病态问题对特定地学模型的影响程度将深化对特定病态问题的认识。
3.2 InSAR观测与地学参数之间模型关系不准制约参数反演
InSAR观测得到的是地表三维形变沿LOS向的一维投影。地面位移通常发生在三维立体空间内，在某些情况下，如沿南北向的走滑、与LOS向近乎垂直的滑坡等，InSAR观测值可能并不能很好地捕捉这些地表形变，也不能很好地作为形变观测值约束地学模型的参数反演，从而使得既有InSAR观测与模型地学参数之间的关系不准确或约束性弱。另外，由于目前SAR卫星的重返周期仍然较大，部分潜在的致灾事件或地球物理现象，如地震、滑坡、冰崩等，是在极短的时间内发生的，且之后还会出现一些残余的运动或环境因素变化。如果上述灾害事件或地球物理现象发生在InSAR两次观测之间，目前的InSAR形变结果很难区分出主要变形与残余运动和环境因素变化的影响。直接将InSAR观测用于震源参数等模型参数反演过程中，将会导致输入数据中存在非主要形变事件和环境因素变化的形变污染，影响模型参数反演的准确性。此外，对于不同类型的灾害事件或地球物理现象，如果对其发生形变的物理机制认知不准确甚至错误，而运用不合适的物理机制或假设进行建模，将严重影响地学参数与InSAR形变观测值之间的关系。这时，很有可能反演出一组伪地学参数，得到不准确甚至错误的物理解释和结论。
3.3 海量观测数导致的反演效率低
InSAR定期重返和高空间分辨率形变监测优势能为地学参数反演提供海量观测数据。例如，面对10 km2的研究区域，3 m空间分辨率的单个InSAR干涉对理想条件下可获得约11万个地表形变观测值。随着时间序列和多平台/多轨道InSAR数据集的加入，观测值数量将呈指数级增长。然而，目前大多地学模型为非线性，且现有方法基本通过非线性搜索反演地学参数。由于非线性反演算法的核心思想是通过解空间搜索寻找最优化地学参数，因此，其面对海量InSAR观测值时反演效率较低。目前，通过四叉树降采样等手段减少观测值数量[106]或使用更高效的反演算法[36]，是克服该问题的途径之一，但是如何降采样获得有代表性的观测样本，如何平衡反演精度和效率，仍是需要研究的问题。
3.4 多次地球活动反演的参数可区分性
受限于当前SAR卫星较长的重访周期(即低时间分辨率)。两景SAR影像差分得到的InSAR形变场往往包含了影像获取时间内所有地球活动所造成的形变[21]。例如，InSAR同震形变中往往叠加了较大余震和震后早期形变(如余滑、孔隙弹性回弹、非弹性形变等)信号，造成基于InSAR反演的同震震源参数具有一定的不确定性。因此，如何有效区分多次地球活动反演的参数，是限制InSAR技术在地学领域应用和推广的瓶颈之一。为了克服上述难题，除了寄希望于未来高时间分辨率的SAR卫星数据之外，部分学者开始引入特定地球活动形变模型(如将余滑模型引入同震滑动反演模型中)，进而构建更为复杂的参数模型以求逼近真实地球物理过程[107]，为后续深入研究指明了方向。然而，特定模型的引入并不能解释所有相关地球活动事件(如引入余滑模型并不能解释较大余震、孔隙弹性回弹等活动造成的形变)，而且联合模型的复杂性(如强非线性等)增强，造成了模型求解稳定性差、耗时量大等。此外，发展基于特定模型有效融合高频GNSS和InSAR联合模型，将为区分多次地球活动反演参数提供有效解决途径。
3.5 InSAR与其他大地测量观测联合反演问题
随着卫星遥感技术的日新月异、计算机数据传输能力与处理能力的不断提高，震源参数反演由基于单一观测数据的反演逐渐发展为基于多源观测数据(水准、GNSS、光学、InSAR等数据)的联合反演，大大推动和加速了震源参数反演进程。联合多源大地测量数据不仅可以提高震源参数解的稳健性，而且具有更好的数据分辨率，进而提高对反演参数的约束能力。InSAR数据具有较高的空间分辨率，可以为浅源地震的反演提供紧约束。GNSS数据具有较高的时间分辨率，可以弥补InSAR技术低时间分辨率的不足。光学影像匹配技术可以监测得到大范围的东西向和南北向水平地表形变场，可以较好地弥补InSAR视线向形变对南北向不敏感和GNSS数据空间上点位稀疏的不足。针对地震产生的大范围、大梯度复杂形变场，综合利用光学影像匹配、InSAR、GNSS等空间对地观测技术各自的优势，对快速恢复地震同震形变场及准确获取震源参数至关重要。利用多源大地测量观测数据联合反演震源参数时，不同观测数据在目标函数中的相对权重决定了各类数据对反演结果的贡献，从而对参数估值的最优性质产生影响。如何准确确定不同观测数据集之间的相对权比是提高震源参数反演精度的关键因素之一。
3.6 InSAR地学参数反演精度评定问题
使用有限的地表观测数据约束复杂的地学参数是非常难的反演问题。常规的InSAR地学参数反演主要包含非线性全局最优估计方法和线性最小二乘方法，这样求解得到的地学参数往往是一组全局最优解。任何反演方法的差异都可能导致截然不同的模型参数解。只有充分考虑反演过程中的不确定性，才能更加客观、准确地反演得到震源参数，更好地解译其机制。随着InSAR数据的类别、模式及获取平台的增多，如何有效考虑InSAR大气误差、卫星轨道误差、DEM误差、数据解缠误差，以及数据处理过程中引入的误差、参数先验信息、系数矩阵误差等重要信息，从而定性、定量地分析和解释这些误差对反演参数不确定性的影响，是InSAR地学参数反演需要进一步解决的难题之一。为了克服上述难题，可以考虑将数据误差引入随机模型，联合贝叶斯反演方法和蒙特卡罗马尔科夫链中的抽样方法评估高维参数后验概率分布，解决传统反演方法难以定性、定量评估参数误差分布的问题，并提供参数误差分布的置信区间。
4 结语
目前，InSAR已被广泛应用于自然因素或人类活动导致的灾害事件和灾变过程的关键地学参数反演。本文系统介绍了国内外学者在地震震源参数、火山活动参数、地下水文学参数、地下采空区几何参数、冻土活动层厚度、冰川物质平衡变化、冰川跃动参数、大气水汽参数及地下流体体积变化参数等地学参数反演领域的成功应用和经验方法。但是，顾及星载InSAR固有的成像模式、有限的观测条件、亟待发展的反演算法等因素，InSAR地学参数反演仍面临着观测误差复杂与单一观测几何制约反演精度的问题、观测量与地学参数之间模型关系不准确的问题，以及海量观测数据导致的反演效率低、多次地球活动引发的反演参数的不确定性、InSAR与其他大地测量观测联合反演、InSAR地学参数反演精度评定的问题。
作者简介
第一作者简介：李志伟(1974-), 男, 教授, 研究方向为InSAR理论与应用。
E-mail: zwli@csu.edu.cn
初审：张艳玲
复审：宋启凡
终审：金 君
资讯
助力企业复工，交运董家口公司运送职工2000余名
大众网·海报新闻2月11日讯(记者 蒋甜 通讯员 杨辉丽)2月10日，交运集团董家口公司服务单位双星集团全面复工，截止到昨晚12点，发出车辆39个车次，运送双星集团职工2000余名。
据悉，为妥善做好复工期的疫情防控工作，董家口公司召开工作会议。交运集团董家口公司为每辆车配备了消毒液、电子体温计、一次性口罩。在发车前，公司机务处对车辆安全技术状况、车容车貌等进行检查，确保车况良好;驾驶人员对参运车辆的方向盘、地板、门窗、扶手等关键位置进行重点消毒;对车辆进行通风换气，时间不低于10分钟，确保车箱内空气流通。
每名驾驶员在发车前进行体温检测，体温正常者方能上岗。除此之外，公司对每名职工的同住人员每天进行体温检测管控，对职工同住人员信息进行登记存档，织密防控网络。同时，加强乘客上车前的体温检测，如遇发热症状的乘客第一时间移交卫生健康部门。
16天一闪的它暗藏了多少信息？
快速射电暴是宇宙中一种非常剧烈的射电闪光现象。随着探测手段的提升以及大型射电装置的投入使用，人们不断地探索着快速射电暴蕴含的丰富的天体物理信息。近日，《自然》期刊发表了一篇关于周期性重复快速射电暴的报道，其约为16天的周期特征将为我们揭示快速射电暴的起源迈进巨大的一步。本文基于张冰教授在《自然》发表的评论文章整理而成。
图片来源：pixabay
快速射电暴（Fast Radio Bursts， FRB）是一种神秘的来自银河系外的射电天文现象。爆发的持续时间仅为几个毫秒，却可在这极短的时间内显示出极高的亮度。根据估计，FRB在几个毫秒释放的能量，相当于太阳燃烧一天的能量总和。
2007年，Lorimer团队从澳大利亚Parkes望远镜历史的存档数据中找到了人类历史上第一个被观测到的FRB事件：FRB 010724。随后，Keane、Thornton等团队陆续发现了其他类似的事件，这标志着FRB的研究正式成为时域天文学的一个重要新方向。
每个FRB发射出的电磁波在宇宙中传播时，通常都会受到星际介质的影响，高频电磁信号总是比低频信号更早到达接收机，这一现象称为电磁波的色散。色散量的大小由介质的等离子体密度和传播的路径决定。自第一例FRB发现以来，不论是观测的色散量超出银河系预期的色散量贡献，还是精确到宿主星系的定位，都表明FRB是一些来自银河系外的辐射过程，其起源一直是天文界的热点话题。
起先，人们认为FRB起源于一些极端的灾难性事件，如双致密天体并合、超新星爆发等。2016年，第一例重复暴事件FRB 121102的发现引起了天文学界的热议。重复暴是指能够在同一个位置产生两次以上FRB事件的源。重复暴的发现使人们意识到FRB绝对不是简单的灾难性事件所能解释的。重复暴的研究尤为重要，它们可以在天空的同一个地方重复地产生。由于FRB持续的时间极短，科学家很难对其进行精确定的空间定位。但重复暴的发生在不断地为我们提供机会，一次不行两次，两次不行多试几次，总是能确定暴发源的位置。
2017年，科学家们实现了对重复暴FRB 121102的精确定位，并找到了其距离我们30亿光年左右的宿主星系。距离这个FRB源120个光年以内的距离，存在一个未知的射电持续源。随着技术手段的不断提升，目前已有9个FRB定位到其宿主星系，暗示着FRB宇宙学的时代或将到来。2018年，在FRB 121102的偏振测量中，科学家发现这个FRB源可能处于一个极高的磁场环境中，这是目前其他FRB源所不具有的特征。尽管人们在重复暴宿主星系以及磁场环境的研究上做了许多努力，FRB的起源依然成谜。
2019年，多个新的重复暴源被发现，天文学者们自此确信重复暴不再是某种个例。迄今为止，已经有超过100个FRB源被发现，而其中有20多个源可以重复地产生爆发。绝大多数重复暴中可以发现一些奇怪的子脉冲结构，这些子脉冲具有不同的中心频率，并在不同时刻抵达探测器，这种现象称为时间-频率漂移。2020年4月份，这一现象在一个银河系内的磁星射电暴发中被发现，暗示着磁星或许与FRB之间有着千丝万缕的联系。
不过。不论磁星也好、中子星也罢，它们都是宇宙中高速旋转的强磁场致密天体。理论上这类天体是产生FRB的绝佳场所，特别是脉冲星这种高速旋转的中子星，其射电辐射束总是周期性地扫过地球视线方向。人们自然联想到，旋转的中子星或者磁星作为FRB的起源，极有可能存在某种周期性特征。
近日，加拿大氢强度测绘实验（CHIME，图1）发布了FRB 180916.J0158 65存在大约16天为重复周期的观测结果，这为人类揭示FRB的起源迈出了巨大的一步！
FRB 180916.J0158 65是CHIME发现的最早的重复暴之一。它也是迄今为止发现的距离我们最近的FRB（红移为0.034，距离我们4.8亿光年），其宿主星系是一个旋涡星系。
该星系拥有较高的恒星形成率，这与另一个重复暴FRB 121102的宿主星系类似，却不同于其他一些非重复的低恒星形成率宿主星系，这或许说明在FRB王国中的确存在着重复暴和非重复暴两类成员。
图1：CHIME望远镜台址。这些与传统意义上的望远镜截然不同的“铁架子”其实是CHIME的天线，图片来源：http://chime.phas.ubc.ca
CHIME望远镜位于加拿大西南部，由4个100米长的柱形天线组成，观测频段为400-800 MHz。CHIME望远镜拥有极大的视场，平均每天可探测到1-10个FRB事件，是捕捉FRB的超级能手。由于FRB 180916.J0158 65经常落入CHIME望远镜的视野，因此每天都会对其进行长时间的自动监控。自2018年9月16日到2020年2月4日以来，共观测到它的38次爆发。此外欧洲超长基线干涉测量网络（EVN）以及德国的Effelsberg望远镜也参与其中，EVN此次也探测到若干爆发，而Effelsberg给出探测的流量上限约为太阳辐射的千万分之一。令人惊讶的是，这些都爆发呈现了16.35天的周期特征。每个时段的活动窗口约为5天，该窗口中的大多数暴发事件都集中在大约0.6天的时间中（图2）。
图2：FRB 180916.J0158 65的周期特征。灰色为5天的射电活动窗口，三角形表示不同望远镜观测到的FRB信号，图片来源：引自The CHIME/FRB Collaboration et al。 2020
在如此长时间尺度的天体物理现象中寻找周期并不是件容易的事，尤其是只有几十个观测样本。这需要仔细分析观测数据来寻找活跃的时间窗口，而这一任务通常会因假定的规则爆发周期未知而变得非常复杂。另外，具有宽功率谱的随机过程，即红噪声，其功率密度平滑地增加到低频。红噪声对周期的寻找也是一个棘手的难题。如果忽略红噪声的影响，随机的变化也可能产生看似周期间隔的行为，如先前对某些类星体的观测（图3）。
图3：类星体PG 1302−102的光变曲线，图片来源：引自Vaughan et al。 2016
对于周期特征的解释，一种说法认为FRB起源于轨道周期为16天的中子星双星系统[13，14]。当中子星与其大质量伴星靠的比较近时，双星系统进入活跃状态。根据轨道周期计算，这需要双星之间的距离大约为四分之一个天文单位（日地距离，约1.5亿公里），然而在这个距离下中子星与伴星的相互作用很难达到产生FRB的要求。我们考虑FRB辐射来自于中子星自身的活动，伴星强大的星风会遮挡来自中子星的射电辐射，如果中子星也存在强力的星风，那么它将扫开周围的恒星星风物质，打开一个射电窗口，使FRB有机会扫过地球。这个窗口以16天为间隔，周期性地扫过我们的视线（图4a）。
图4：周期性FRB的两种图像，图片来源：引自Zhang 2020。（a）中子星双星系统；中子星（Neutron star）的星风将扫开大质量伴星星风的抛射物质，形成一个射电透明的观测窗口，使FRB得以逃脱。随着双星系统的运动，这个射电窗口扫过地球视线的周期即为双星系统的轨道周期。（b）中子星的进动；中子星的射电辐射束非常窄。当具有一定椭率的中子星自由进动或者受外界影响受迫进动时，其辐射束以进动周期等间隔性地扫过地球视线方向。
中子星的进动周期是另一种可能的解释。“周期性”一词对于中子星而言，往往自然让人联想到自转周期，然而实际上中子星具有极高的自转速度（周期大小为1毫秒到1秒），16天的长度似乎不太能解释为中子星的自转周期。由于中子星自身存在一定的椭率，高速旋转的中子星可产生周期为几十天的自由进动。又或者在伴星的引力作用下，中子星同样可产生类似周期大小的受迫进动。这里在几何上要求FRB产生于中子星的磁层，极端相对论电子使大部分辐射流量被集束在一个很窄的光锥内。这些射电束就像陀螺仪一样，随着中子星的进动，周期性地进入地球的视野（图4b）。
总之，FRB周期的探测，可谓是一项艰巨而精细的工程，它对于揭示FRB的起源有着重大意义。同时，FRB的周期活跃性也为后续跟踪观测提供指导。当然，是否存在比16天更短的公约数周期问题依然值得探索。我们不禁思考，周期性特征究竟是不是所有重复暴共同具有的属性。非重复暴本身是一些重复率极低的事件，还是本征的一次性事件。这需要我们在未来使用大型观测设备投入更多的时间跟精力进行探索。相信揭秘FRB真相的这一天已经不远了。
星座周运指南4月2日—4月8日
今天，金星换座进入金牛座。而持续的，火土双星于本周合相于摩羯座并刑克逆行水星。这让这一周既充满困难，又不乏收获。
勤劳的人总是会有收成，畏惧困难的人总会离成功很远。
请参考太阳和上升星座
白羊座
拼命三郎是不是说的就是现在在为你的事业而努力拼搏的白羊座呢？这可能是你事业上的瓶颈期，也可能是一个大的挑战期，这不会吓到你。相反，这反倒激起了白羊座的求胜心，开始变得比往常更加有恒常心，来对待你所拼搏的这份事业。财富有了起色会让你感到富足，满意。而最重要的是，你要重新整理你的理财保险或资源，可能会有新的发现。
金牛座
很有魅力的时刻来了，不像狮子座，这段时间沉默低调的牛儿不需卖力就自然的风情万种，魅力无限。这是金牛座相亲约会的好时节。有什么比和自己的心上人到花树下野餐更加浪漫呢？高等教育考试、评级可能会给你很多压力，让你投入很多精力，但仍然感到不是很自信，但是相信这是属于你的时刻，如果是面试，你高升的魅力值会帮到你哦。
双子座
工作上得到的一些利益和嘉奖在这段时间容易出现被掠夺、失去价值、被攻击的可能。尤其一些让你莫名其妙的事情都可能发生。双子们还是平日里说话做事要再谨慎一些，不然，不知不觉间自己已经得罪了人，埋下了日后遭到报复的种子，自己却还不自知。你或许真的该好好思考一下，谁才是你真正的友人了。另外，双子也有出现意外的梦中情人的状态。
巨蟹座
婚姻的压力可能对巨蟹来说是分外痛苦的，在你做了一些展现自己才能，在事业上大干一场的计划之后，现在或许更应该你来思考，到底自己可以为这个社会贡献出什么？这份贡献的价值和你对家人的重要性来说，到底哪一个对你来说更重要，更具备价值。单身巨蟹有在朋友之中或者通过朋友介绍遇见自己心上人的可能。也许就是这一次聚餐和郊游？！把握好机会哦。
狮子座
狮子座需要在这段时间注意自己的身体状况，这一点除了普通的理疗和检查，保养之外，很重要的一点是去查除自己潜意识里的罪恶。这对狮子来说有点困难，因为让狮子座说自己是错的，是不好的真的还不如让他们去死。但是很多时候，有一些情况就会这样发生，当你拒绝去发现和改进的时候，他们已经在你的身躯里病变了。
处女座
来自远方（也可能是法律亲属，比如继父母等）的问候会让你感觉幸福，甚至有时候，他会以一种馈赠的形式出现。当然这也是你度假的好时机。不过，说到放松，或许心理咨询是更加适合你此时此刻状态的。再或者是和自己的敌人做一些推心置腹的交谈，化解你们之间的隔阂，你要记住的是，处女座最忌讳的是把所有的标准都添加在了别人和那个梦想出来的自己身上。而对自己的真实视而不见。损人不利己。
天秤座
因为你善巧、温和又高情商的社交，你会因为你的谦逊和善于妥协在这段时间收获到一些回报。另外，最近你的家庭情况可能会给你造成一些压力。父母的身体还好吗？要记得在繁忙生活中去探望，去观察，不要等到一些情况发生了再去补救。已婚天秤可能因为一些私情的爆发而引发婚姻的不稳定和伤害。你们的沟通需要更加深入，而不是把旧账都扯出来不断互相攻击。
天蝎座
你的亲密关系令你感到一些困扰和艰涩。甚至有可能出现分手的苗头。但是如果你是单身天蝎，那么则可能突然间出现一些对你而言非常有魅力的人，你满载活力却也谨慎周密地开始靠近。这种感觉让你有一点激动，但也并不让你感到畏惧。你会倾向于隐藏自己，而更多时候，你会在情感面前发现自己令自己惊诧的一面。
射手座
这段时间持续的艰辛总算是换来了一些收益。这令你感到满足。你会继续这样勤奋地工作，赚钱，养家，实现自我价值。不过适当的也可以放松一些，在你的生活中你可以适当地添置一些柔软亲和的内容成分，比如定期到茶楼约三两好友高谈阔论。这并不一定是在浪费时间，相反，他们也许会帮助你发现自己其他的潜能。你不仅需要死守一些东西，你还需要一些灵感。
摩羯座
成为健身达人并不是坏事，毕竟在春暖花开的时候如同村上春树一般在花树间慢跑，对于摩羯来说是非常享受的事情。多和家人在一起，思维一下你的家庭结构和分工，有时候间或性地调整一下未尝不是对家庭的好事。周末了，安排一次和家人的郊游，到青山绿水间好好放松。这是你需要的。这种行为有利于打开你的财脉。
水瓶座
一些奇怪的事情可能发生在你的生活中，交流上要注意谨慎和严密，如果感觉最近火气有点旺，那么要多做呼吸静心。交流上的严谨是为了减少爆发一些意外攻击以及暗算的发生。而思考上的严谨则是需要你把你手上的资源重新梳理一下，看看是否还有什么漏掉的价值存档没有开发。做了很不好的噩梦，记得要忏悔。
双鱼座
在社交层面，一向非常如鱼得水的双鱼座这段时间可能会感到束手束脚、不大顺利。有时候，冲突也是在所难免。这个时候你最重要的就是要反思一下，自己的价值观是否出了什么问题。出差、游历都是一种可能，而发生法律事务则成为另一种极大的可能性。你需要让自己在一种田园诗般的环境中放松一下。
作者：新月文化