版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、声 波 加 p ml 边 界 条 件 雷克子波 因计算机资源有限,不便在太大的区域做 PML的计算,一般只在边界上取有限 宽度的区域作为 PML计算区域。根据相关文献的研究, PML边界区域最少长度应为 半个波长 6 。本文综合考虑了效果与开销等因素,选取了边界上50 层作为 PML的 计算区域。常规计算区域与 PML边界区域的如图 2-4 所示。 ax Px t xPx K vx x Pz vz az z zPz K z t z P vx vz ax xP K t x z xz 2 vP lo g 2 L3 衰减系数 式中 L为 PML层的厚度, 为层内的点距 PML与非 PML的边界的距离
2、, vP 为纵 vx1 P 波速度, R 那么在 PML边界区域内,对于式( 2-13) tx xvxx 即为理论反 tx 射系数,一般取 0.001 较为合适, 为 方向的空间步长。 vx t xvx ,可看作为在常规的计算方程基础上,减去一项进行 PML的 阻尼修正项 因本文中只考虑各项同性介质中的地震波传播规律,故可做 假设。在此利用一下三个假设: k 1 vx 1 i,j vx i, j 2 2 k 1 vz 1 vz i 2, j 2 k1 1 k1 P 2 P i , j 2 i,j x z i k1 2 1 i,j 2 k1 2 1 2,j k1 Pik,j1 vx k1 2 1
3、 i,j 2 k1 2 1 2,j 因为以上三个近似精度均为时间方向上的近似,且时间精度均为二阶精度,因 交错网格技术的时间精度为二阶,故以上近似不影响本式的计算精度 故可得: vx k1 2 1 i,j 2 1 1 0.5 0.5t vx k1 2 1 i,j 2 k i,j xm am Pi,j m 1 Pi ,kj m 1 i ,j m 1 同理: vz 1 k 2 1 2,j 1 1 0.5 2-18a ) 0.5 k t vz 1 2 1 2,j 1 k 1 i,j 2 zm k am Pi m,j 1 Pi km 1,j i m 1,j 2-18b ) Pik,j i,j 1 1
4、0.5 t 0.5 t P ik,j i,j k i,j N am z m 1 vz tN xam xm1 k1 2 1 m ,j 2 k1 2 x1 i,j m vz k1 2 1 i m , j 2 k1 vz i ,j2m 1 z i ,j m 12 2-18c ) 此为在 PML边界区域内的弹性声波应力 - 速度方程组。 % %交错网格- 非均匀介质二维声波方程 (一阶压力- 速度) 、2阶时间差分、 2阶 空间差分精度 %加上边界 % close all ;clear,clc tic tt=-0.06:dtt:0.06; %* dtt=0.0001; 震源为 Ricker 子波 fm
5、=30; A=0.01; wave=A*(1-2*(pi*fm*tt).2).*exp(-(pi*fm*tt).2); plot(wave),title( 震源子波 -Ricker 子波 ); %* % 模型参数设置 dz=5;% 纵向网格大小,单位m dx=5;% 横向网格大小,单位m dt=0.0001;% 时间步长,单位 s T=0.5;% 波动传播时间,单位s wave(round(T/dt)=0; % 将子波后面部分补零 % % 研究区域 % z=-750:dz:750; x=-1000:dz:1000; pml=50;% 吸收层的网格数 plx=pml*dx;% 上下吸收层的厚度
6、plz=pml*dz;% 左右吸收层的厚度 z=-750-plz:dz:750+plz; x=-1000-plx:dx:1000+plx; % 采样区间 n=length(z); m=length(x);% 采样点数 z0=round(n/2); x0=round(m/2);% 震源位置 Vmax=0; % 纵波最大速度 %Setting Velocity xt=-1000-plx:dx/2:1000+plx; % 速度与密度采样区间 nt=length(zt); mt=length(xt); % 速度与密度采样点数目 V=zeros(n,m);% 介质速度 ,m/s d=zeros(nt,m
7、t);% 介质密度 ,kg/m3 %均匀介质模型 for i=1:n for k=1:m V(i,k)=2.0e3; end end for i=1:n for k=1:m d(2*i,2*k)=2.3e3; end end % % %层状介质模型 % % for i=1:n % % for k=1:m % % if i Vmax Vmax=V(i,k); end end end %* 衰减系数 % ddx 、ddz 即,x方向和 z方向的衰减系数 R=1e-6; % 理论反射系数 ddx=zeros(n,m); ddz=zeros(n,m); for i=1:n for k=1:m % 区域
8、 1 if i=1 z=i-(n-pml); ddx(i,k)=-log(R)*3*Vmax*x2/(2*plx2); ddz(i,k)=-log(R)*3*Vmax*z2/(2*plz2); % 区域 2 elseif ipml z=0; ddx(i,k)=- log(R)*3*Vmax*x2/(2*plx2);ddz(i,k)=0; end end end % figure(1),imagesc(ddz),title(z方向衰减系数 ); % figure(2),imagesc(ddx),title(x方向衰减系数 ); %* p0=zeros(n,m); p1=zeros(n,m); p
9、x0=zeros(n,m); px1=zeros(n,m); pz0=zeros(n,m); pz1=zeros(n,m); K=zeros(n,m); %* 波场模拟 Vx1=zeros(nt,mt); Vx0=zeros(nt,mt); Vz1=zeros(nt,mt); Vz0=zeros(nt,mt); for t=dt:dt:T p0(z0,x0)=dt*V(z0,x0)2*wave(round(t/dt); for i=2:n-1 for k=2:m-1 K(i,k)=d(2*i,2*k)*V(i,k)2; Vz1(2*i+1,2*k)=(1- 0.5*dt*ddz(i,k)*Vz
10、0(2*i+1,2*k)-dt*(p0(i+1,k)- p0(i,k)/(d(2*i+1,2*k)*dz)/(1+0.5*dt*ddz(i,k); Vx1(2*i,2*k+1)=(1- 0.5*dt*ddx(i,k)*Vx0(2*i,2*k+1)-dt*(p0(i,k+1)- p0(i,k)/(d(2*i,2*k+1)*dx)/(1+0.5*dt*ddx(i,k); pz1(i,k)=(1-0.5*dt*ddz(i,k)*pz0(i,k)- K(i,k)*dt*(Vz1(2*i+1,2*k)-Vz1(2*i- 1,2*k)/dz)/(1+0.5*dt*ddz(i,k); px1(i,k)=(1-0.5*dt*ddx(i,k)*px0(i,k)- K(i,k)*dt*(Vx1(2*i,2*k+1)-Vx1(2*i,2*k- 1)/dx)/(1+0.5*dt*ddx(i,k); p1(i,k)=px1(i,k)+pz1(i,k); end % figure(2),imagesc(pz1);title(z % figure(3),imagesc(px1);title(x 分量 ); 分量 ); figure(1),imagesc(p1);title( figure(2),imagesc(p);title( %* Toc end p0=p1; pz0=pz1;px0=px1; V
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 中学手机使用规范管理实施方案
- XX学校智慧校园应用方案
- 汽车配件采购与维修方案
- 采矿工程师职业技能提升总结
- 适合小学生的预防艾滋病教育
- 外卖骑手安全意识的培训
- 文化艺术团体团建活动合同
- 2024-2025学年高中地理选择性必修1中图中华地图版(单元测试)试题合集
- 10KV铁塔施工应急预案方案
- 幼儿园新入园免疫接种制度解析
- 大学生青年志愿者协会退部制度
- 口腔诊所物品管理与设备维修制度
- 2022-203学年(中职)《餐饮服务与管理》试题3试卷带答案
- 清华大学工程热力学课件 第1章01
- 2022-2023学年山东省东营市实验中学物理高一第一学期期中达标测试试题(含解析)
- 中式面点师初级(五级)教学计划、大纲
- 小学综合实践二年级上册第三单元《神奇的影子》教材分析及全部教案
- 小学语文人教六年级上册(统编2023年更新)第六单元-《浪淘沙》导学单
- 测试计划-10篇模板
- 交联电缆运行时的热特性参数设计计算
- TSG-R0005-2022《移动式压力容器安全技术监察规程》(2022版)
评论
0/150
提交评论