一个计算水文质量变化产生的负荷变形的Matlab代码

地球是一个由大气圈、水圈、岩石圈组成的复杂系统，其内部发生众多的地球动力学过程，比如地震活动、火山喷发、水文负荷等，均涉及到地球在内外力作用下的变形。地球的负荷变形有着相当长的研究历史，在1911年，Love定义了两个无量纲的数h和k来表述地球的垂直位移和引力位变化。而在1912年志田补充了第三个无量纲的数l来描述水平位移。FARRELL在1972年发表的文章 Deformation of the Earth by Surface Loads系统总结了负荷计算的发展历史，通过数值积分计算了更高阶的Love数和Boussinesq近似道出了高阶渐近解，并基于渐近解、圆盘因子、Kummer变换等方法来加速格林函数尤其在近场中计算的收敛速度，首次得到完整的格林函数表，用其来表述位移、重力、倾斜、应变等物理量。

我们知道，地表的质量变化会导致重力变化，这可以由GRACE重力卫星观测得到。而质量的重新分布，则会对地球产生负荷变形。有研究指出，流域水量变化（三峡水库蓄水）、冰川消融、地下水亏损等均会产生符合变形。本小结利用Pan et al.(2000, JGR)的文章Intradecadal Fluctuations and Three-Dimensional Crustal Kinematic Deformation of the Tianshan and Pamir Derived From Multi-Geodetic Imaging里面提供的开源Matlab代码计算CPC土壤水变化对地表产生的垂直、北向、东向位移。原始函数下载地址：水文负荷计算

% Chistrong Wen
% UCAS
% 2023-3-19
% GreensFunction = load('farrell_disp_Green.txt');
grid_EWH_m = (CPC_SOIL.soil(:,:,10)-CPC_SOIL.soil(:,:,9))/1000;
%% 70:75;35:40;
lon = 65:1:90;
lat = 35:1:50;
[lon,lat] = meshgrid(lon,lat);
load.lon = lon;
load.lat = lat;
% for i = 1:16
%     if(i == 1)
%         disp(['total iteration:',num2str(70*40)])
%     end
%     for j = 1:26
%         gps_lon = lon(i,j);
%         gps_lat = lat(i,j);
%         dense_factor = 10;
%         window_radius = 50;
%         [loading_meter_u_global(i,j),loading_meter_n_global(i,j),loading_meter_e_global(i,j)] = get_mascon2loading_scatter_Dense_Global(GreensFunction,grid_EWH_m,gps_lon,gps_lat,dense_factor,window_radius);
%         disp(j)
%     end
% end
load.rg = loading_meter_u_global.*1000;
cpc.lon = CPC_SOIL.lon;
cpc.lat = CPC_SOIL.lat;
cpc.rg  = grid_EWH_m*1000; %mm
subplot(2,2,1)
wzq_plot(cpc)
caxis([-100,100]),xlabel('Lon'),ylabel('Lat'),title('Soil moisture change')
colorbar
xlim([65,90]);ylim([35,50])
subplot(2,2,2)
wzq_plot(load),xlabel('Lon'),ylabel('Lat'),title('vertival deformation')
colorbar
subplot(2,2,3)
load.rg = loading_meter_n_global.*1000;
wzq_plot(load),xlabel('Lon'),ylabel('Lat'),title('Northward deformation')
colorbar
subplot(2,2,4)
load.rg = loading_meter_e_global.*1000;
wzq_plot(load),xlabel('Lon'),ylabel('Lat'),title('Eastward deformation')
colorbar

计算的结果：

一个计算水文质量变化产生的负荷变形的Matlab代码相关推荐

基于BP神经网络的电力负荷预测（Matlab代码实现）
目录 1 概述 2 运行结果 3 完整Matlab代码实现 4 数据 1 概述随着电力工业逐步进入市场化,区域配电网短期电力负荷预测在电力行业中地位越来越重要,精准的短期负荷预测方法对电力系统安全稳 ...
用四阶RungeKutta方程解二阶常微分方程，并计算船舶在规则波中的横摇角（附Matlab代码）
前几天接到师姐分派的任务,让我求解一艘船模的横摇角的时间历程曲线,为后期的减摇控制准备. 1 首先冷静分析一下,原方程如下: 我们要求解的就是theta角和时间t之间的关系曲线,这是一道典型的二阶常微 ...
【负荷预测】基于改进灰狼算法（IGWO）优化的LSSVM进行负荷预测（Matlab代码实现）
基于KPCA 和 STFT 非侵入式负荷监控（Matlab代码实现）
基于机器学习之模型树短期负荷预测（Matlab代码实现）
【负荷预测】长短期负荷预测（Matlab代码实现）
【FFNN负荷预测】基于人工神经网络的空压机负荷预测（Matlab代码实现）
利用Html,Css,javaScript简单写一个计算一天秒的倒计时
利用Html,Css,javaScript简单写一个计算一天还剩多少秒的倒计时页面代码我放下面,不多,就是我那个背景图我没有用精灵图,主要拼接太麻烦了,我直接img了,然后display:none; ...
负荷计算的时候assert失败_负荷计算的入门基础知识，小白入门必备！
负荷计算是电气设计人员必须掌握的一必修课,选导体.设备还有保护开关的选择,电网系统分析,都离不开负荷计算的内容,为此小编特意总结9个负荷计算的入门基础知识,分享给大家,希望能在工作和学习中作为一份参考 ...

一个计算水文质量变化产生的负荷变形的Matlab代码

一个计算水文质量变化产生的负荷变形的Matlab代码相关推荐

最新文章

热门文章