最新要闻
- 每日焦点!生态河湖 诗画江城
- 2.4万元!苹果首款MR头显Vision Pro发布:单眼像素超越4K电视
- 真颠覆XR!苹果Vision Pro技术、体验几乎完美 但是一根线毁所有|世界快报
- 苹果Mac三大新品国行价格汇总:巅峰顶配9.7万元-环球速看
- 8450米的珠峰上舍命救人 女子承诺1万美元只肯给4000美元:贵吗?
- 1100万年后 决定性别的Y染色体将消失:人类会怎么样?
- 旅游景区清洁工人工作情况影像记录带定位时间时间水印相机
- 全球动态:父母去世了,为子女买的保险能用遗产续费吗?
- 天天讯息:高中教材选修和必修有什么区别_高中必修是什么意思 选修呢
- 油毡纸多少钱一卷(油毡纸) 环球播资讯
- 当前短讯!不负如来不负卿txt下载完整版_不负如来不负卿txt下载
- UFOU首获投资,升降桌售价上万元有人望而却步,孟振亮曾连续创业
- 每日短讯:管理能力有哪些方面_能力有哪些方面
- 阿里云宣布举办首届AI黑客马拉松挑战赛 最终获胜团队将获得3万奖金
- 现实题材电影《八角笼中》发布一“骗”真心版预告 由王宝强与肖央主演
- 人造电子皮肤首次亮相 每层主干均为长分子链组成
广告
手机
iphone11大小尺寸是多少?苹果iPhone11和iPhone13的区别是什么?
警方通报辅警执法直播中被撞飞:犯罪嫌疑人已投案
- iphone11大小尺寸是多少?苹果iPhone11和iPhone13的区别是什么?
- 警方通报辅警执法直播中被撞飞:犯罪嫌疑人已投案
- 男子被关545天申国赔:获赔18万多 驳回精神抚慰金
- 3天内26名本土感染者,辽宁确诊人数已超安徽
- 广西柳州一男子因纠纷杀害三人后自首
- 洱海坠机4名机组人员被批准为烈士 数千干部群众悼念
家电
焦点观察:MRST绘制三维网格图
MRST绘制三维网格图
Matlab储层模拟工具箱(The MATLAB Reservoir Simulation Toolbox, MRST)是一款用于储层建模的免费、开源的软件,主要由 SINTEF Digital 数学与控制论系的计算地球科学小组开发。
更多介绍MRST (sintef.no)
下载地址Download (sintef.no)
(资料图片仅供参考)
在下载完工具箱之后,将其加入库路径,随后开始绘图。
首先是针对规则网格的绘制与展示:
1 规则网格
Step1用cartGrid可以创建一个空的笛卡尔网格,如下图(a)所示;
G = cartGrid([3, 3, 3], [1, 1, 1].*meter);%创建一个3*3*3的网格plotGrid(G), view([-40, 32]), xlabel("x"), ylabel("y"), zlabel("Depth")
Step2对空网格进行赋值;
在这一步中我们对这3×3×3=27个网格进行赋值,然后成图,如图中(b)
G = cartGrid([3, 3, 3], [1, 1, 1].*meter);Values = linspace(1,27,27).";%将数据转成与对应网格数相同的n×1的序列即可plotCellData(G, Values)%通过plotCellData显示view([-40, 32]), xlabel("x"), ylabel("y"), zlabel("Depth"), colorbar
如果我们需要的网格就是这类规则的方格,那成图过程简单概括也就是“建网格,贴数值”。
2 不规则网格
Step1用cartGrid可以创建一个空的笛卡尔网格,抠除掉一个无效的网格,使之成为透明的,如下图(a)所示;
%同样还是先创建一个3*3*3的网格% 创建 3x3x3 的笛卡尔网格G = cartGrid([3, 3, 3]);actnum = ones(G.cells.num, 1);actnum(21) = 0;% 移除不可见的网格G = removeCells(G, actnum == 0);plotGrid(G), view([-40, 32]), xlabel("x"), ylabel("y"), zlabel("Depth")
Step2对空网格进行赋值;
%% 抠网格贴颜色%同样还是先创建一个3*3*3的网格% 创建 3x3x3 的笛卡尔网格G = cartGrid([3, 3, 3]);actnum = ones(G.cells.num, 1);actnum(21) = 0;% 移除不可见的网格G = removeCells(G, actnum == 0);% 为剩余网格生成随机数值values = rand(G.cells.num, 1);plotCellData(G, values);%通过plotCellData显示view([-40, 32]), xlabel("x"), ylabel("y"), zlabel("Depth"), colorbar
将得到这样一个结果:
3 任意网格
当我们需要对一个任意的三维数据,进行3D网格展示时,执行的操作可概括为,创建一个范围足够的规则网格,随后从规则的网格里抠除无值部分。具体到如切片数据时,仅需依据Line、CDP、Time范围对三维网格赋值、重排,坐标排序格式如下,包含X,Y,Z和指定位置处的Value,有值时Actnum为1,否则为0:
完整代码如下:
filename = "Demo.slice";%切片文件[line,cdp,X,Y,Time,value]=readslicedata(filename);%读取切片函数showslice3D(line,cdp,Time,value);%切片三维可视化
其中包含两个子函数
(1)读取切片函数readslicedata
function [line,cdp,X,Y,Time,value]=readslicedata(filename) fid=fopen(filename,"rt"); if fid<0 warndlg("打开文件失败!"); return; end FormatString=repmat("%f",1,7); dataoutput=textscan(fid,FormatString,"HeaderLines",1); line=dataoutput{2}; cdp=dataoutput{3}; X=dataoutput{4}; Y=dataoutput{5}; Time=dataoutput{6}; value=dataoutput{7}; fclose(fid);end
(2)切片三维可视化函数showslice3D
function showslice3D(line,cdp,Time,value)%汇总切片信息num = size(line,1);sliceInfo = zeros(num,5);sliceInfo(:,1)=line;sliceInfo(:,2)=cdp;sliceInfo(:,3)=int64(Time);sliceInfo(:,4)=value;sliceInfo(:,5)=1;%获取线道号、时间最大、最小范围minLine=min(line);maxLine=max(line);minCDP=min(cdp);maxCDP=max(cdp);minTime=min(Time);maxTime=max(Time);%重组切片坐标[X, Y, Z] = meshgrid(minLine:maxLine, minCDP:maxCDP, minTime:maxTime);COO = [X(:), Y(:), Z(:)];% 数组转表格COOtable = array2table(COO, "VariableNames", {"X", "Y", "Z"});Slicetable = array2table(sliceInfo, "VariableNames", {"X", "Y", "Z", "Value","Valid"});% 连接两表joinedTable = outerjoin(COOtable, Slicetable,"Keys", {"X", "Y", "Z"});selectedTable = joinedTable(:, {"X_COOtable", "Y_COOtable", "Z_COOtable", "Value","Valid"});selectedTable = sortrows(selectedTable, [3, 2, 1]); %按z->y->x的顺序重新排序data=table2array(selectedTable);values = data(~isnan(data(:,4)),4); %只取有值的部分data(isnan(data)) = 0; %对空值填0actnums = data(:,5);Xnum = maxLine-minLine+1;Ynum = maxCDP-minCDP+1;Znum = maxTime-minTime+1;G = cartGrid([Xnum, Ynum, Znum]);% 移除不可见的网格G = removeCells(G, actnums == 0);plotCellData(G, values,"EdgeColor", "none");%通过plotCellData显示view([-40, 32]), xlabel("Line"), ylabel("CDP"), zlabel("Time"), colorbarend
最后结果如下:
关键词:
焦点观察:MRST绘制三维网格图
每日焦点!生态河湖 诗画江城
2.4万元!苹果首款MR头显Vision Pro发布:单眼像素超越4K电视
真颠覆XR!苹果Vision Pro技术、体验几乎完美 但是一根线毁所有|世界快报
苹果Mac三大新品国行价格汇总:巅峰顶配9.7万元-环球速看
8450米的珠峰上舍命救人 女子承诺1万美元只肯给4000美元:贵吗?
1100万年后 决定性别的Y染色体将消失:人类会怎么样?
旅游景区清洁工人工作情况影像记录带定位时间时间水印相机
全球动态:父母去世了,为子女买的保险能用遗产续费吗?
天天讯息:高中教材选修和必修有什么区别_高中必修是什么意思 选修呢
油毡纸多少钱一卷(油毡纸) 环球播资讯
当前短讯!不负如来不负卿txt下载完整版_不负如来不负卿txt下载
UFOU首获投资,升降桌售价上万元有人望而却步,孟振亮曾连续创业
每日短讯:管理能力有哪些方面_能力有哪些方面
阿里云宣布举办首届AI黑客马拉松挑战赛 最终获胜团队将获得3万奖金
现实题材电影《八角笼中》发布一“骗”真心版预告 由王宝强与肖央主演
人造电子皮肤首次亮相 每层主干均为长分子链组成
美国一男子连吃100天麦当劳后减重53斤 网友困惑垃圾食品是否需重新定义
巴奴火锅反向推出“天价土豆”套餐 套餐几乎囊括巴奴所有热门菜品
日本地下偶像女子团体成员灌男粉丝喝洗洁精 看其吐泡泡后捧腹大笑
在线访问量直线飙升 OpenAI已跻身全球20家流量最高网站榜单
2023年夏日游戏节发布预热宣传片 将于6月9日进行直播
比亚迪上调部分车型的保养费 幅度达到50%左右
清华大学成立秀钟书院:今年开始招生 目标宏大_当前热文
语音系统用上全自研算法!李想:95%以上控制张嘴就能搞定-天天亮点
浙江大学学生被录用为“场务及驱鸟员” 网友惊呼太卷了
创美药业(02289)将于7月31日派发末期股息每股0.45元
hackthebox sniper medium
致乘风破浪的青春|胡千行:求是创新,多维发展-消息
网上出现假冒北京市人社局的网站!官方发布风险提示_头条焦点
警惕!身上有这7个特征的人 容易得肺癌
男子嘴贴保鲜膜忍不住亲亲宝宝引热议 医生科普:被老人、家人亲吻会遭殃
太卷了!浙大学霸应聘机场驱鸟员:要求英语四级425分 环球微头条
文心一言 VS 讯飞星火 VS chatgpt (32)-- 算法导论5.2 4题_全球看点
环球微动态丨【解决方法】网络设备使用CLI命令行模式进入Telnet登录,如交换机,路由器
厦门大学有哪些专业 厦门大学专业排名
天天热头条丨信泰航空拟对全资子公司信泰克增资1000万
5月国内汽车召回数量排名:美系特斯拉遥遥领先 全球滚动
珠峰为救遇险女子放弃登顶当事人发声:停止网暴_环球视讯
OpenAI跻身全球前20大网站:月活用户近9亿
购物追剧两不误!腾讯视频会员年卡+京东PLUS年卡降价:到手138元
排海计划坚持不变!福岛核电站鱼体放射性物质超标180倍 环球头条
当前动态:2023年南疆冬小麦成熟期及复播玉米播种期预报金十期货6月5日讯,根据目前南疆冬小麦发育进程和6月天气气候趋势预测,预计2023年南疆冬小麦成熟期大部
全国一体化算力算网调度平台正式发布
比iPhone 14更薄的机身塞进4600mAh!vivo S17本周开卖:2499起
AMD闪电发起价格战!RTX 4060 Ti这下太难过了 环球报道
每日精选:NVIDIA黄仁勋已返回美国 此前被传本周前往上海
Midjourney最强平替来了!广告大片分分钟直出:网友直呼太逼真
上千块的吹风机到底是不是智商税?看完秒懂
消息!2023年618淘宝好价节官方立减规则是什么?是否计入最低价?
React学习时,自己拟定的一则小案例(table表格组件,含编辑)
投资者权衡美联储利率政策 美债收益率周一全线走高
Intel Arc独立显卡第一次上水冷!这是何必呢?_环球速看料
上千师生操场烧烤缓解考前压力 老师:学校自费 图个好寓意 环球时讯
损失过亿!唯品会崩溃遭P0级故障:负责人被免职
国产显卡第一家!摩尔线程支持DX11:这下稳了
Edge独占结束!微软确认BingChat将向第三方浏览器推出
知乎,在短故事赛道跑出一个爆款
今日快看!Python第三方库批量下载到本地,并离线批量安装第三方库
AI版女网红“半藏森林”上线,服务项目让人意想不到 世界热门
全球观天下!在win2016服务器上安装Loadrunner 11版本
商品日报(6月5日):商品市场走势分化 原油系领涨金属全线回落
iPhone 15系列将于本月在郑州富士康量产:备货量巨大
特斯拉车主高速违停飞无人机:风景太美忍不住|全球观速讯
腾讯天美联合开发:育碧中国背景《刺客信条》新图公布
环球头条:iPhone 15将全系采用三星M12屏幕:Pro版依然独享ProMotion动态高刷
高人一等!男子开路虎插队洗车 我就猖狂了:嚣张到网友狂吐槽
最美贵州在等你_对于最美贵州在等你简单介绍|世界动态
东铁线是火车还是地铁(东铁线)
API NEWS | Money Lover爆出潜在API漏洞_世界速看
探究GIS地图在城市规划、环境管理和农业领域的应用
【全球速看料】druid活跃线程数量持续增长问题
世界快报:数据的完整性是企业数字化转型的基础
近万条儿童科普知识大全ACCESS\EXCEL数据库
三缸豪车 雷克萨斯最入门SUV车型LBX官图发布:用上1.5L混动-世界微速讯
腾讯FPS王炸大作!《无畏契约》国服客户端预载:6月8日终测 焦点
夏季风扇和空调到底谁更健康?看完明白了
2023年5月新能源车批发销量出炉:比亚迪能打10个长城
半藏森林克隆人暂下线:还在内测 名额已满_焦点简讯
Java虚拟线程
9千多少儿百科全书百科知识ACCESS\EXCEL数据库|天天看热讯
环球速讯:【解决方法】loopback口添加同目的网段的地址后,无法用默认路由ping通该网段地址
ppt应该保存什么格式
【环球财经】日经股指涨破32000点 为1990年7月以来首次
每日播报!高中最特别的时刻!高考前最后一个晚自习还记得多少?
哪吒CEO张勇:哪吒S是最漂亮的B级车 还是100万内最好的轿跑 热门
迪士尼一边涨价一边裁员:门票卖到799元
女子酒驾被查 喝了几瓶菠萝啤 以为是饮料 天天新资讯
天天热头条丨《小美人鱼》反派人类公主凡妮莎剧照出炉:比美人鱼颜值高、超甜
ST华铁(000976)6月5日主力资金净买入651.43万元 全球热闻
ChatGLM-6B int4的本地部署与初步测试
Ambient Mesh:Istio 数据面新模式
某OA 11.10 未授权任意文件上传
判断非String对象是否为null,小伙竟然用StringUtils.isEmpty(obj+"")
全球负收益率债券规模反弹至近2万亿美元
西部矿业: 目前上半年尚未结束,公司暂时没有收到锂资源公司分红的消息_每日快讯
今日黄金td行情分析(2023年6月5日)
前沿资讯!净化器原厂/第三方滤芯对比:差距高下立判
网传北京将出台摩托车新规 京A禁止过户 多方回应
《暗黑4》玩家希望新角色继承地图数据:重新跑图太无聊 实时焦点