【答案】应助回帖

» 本帖@通知

@在

» 猜你喜欢

中国矿业报：中山大学曹建劲团队发明纳米微粒找矿技术已经有16人回复
新突破！纳米微粒找矿技术已经有11人回复
地球物理学和空间物理学论文润色/翻译怎么收费? 已经有66人回复
中山大学：非晶微粒携带金属元素迁移和找矿应用已经有10人回复
中山大学曹建劲教授团队：用动植物体内中含金属微粒探测隐伏矿体技术已经有11人回复
地质勘查类的高级工程师评审，需要发表论文已经有1人回复

1楼2019-04-17 10:58:38

diwushuxue

铁虫 (初入文坛)

应助: 2 (幼儿园)
金币: 343.8
红花: 2
帖子: 39
在线: 32.8小时
虫号: 1243949
注册: 2011-03-24
专业: 应用地球物理学

【答案】应助回帖

感谢参与，应助指数 +1

网格多大，你这跟波场分离没关系，是你正演没做好。有可能是频散过于严重

2楼2019-04-17 11:39:28

Ego悄吟

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 114.7
帖子: 31
在线: 4.2小时
虫号: 14884382
注册: 2019-04-15
专业: 结构工程

引用回帖:

2楼: Originally posted by diwushuxue at 2019-04-17 11:39:28
网格多大，你这跟波场分离没关系，是你正演没做好。有可能是频散过于严重

您好，看到您的回帖之后，我回去划细网格了，之前网格大小为10m，计算步长为1.60e-3s,之后网格大小为2.5m，计算步长为0.8e-3s，可是依然出现这个情况。

3D.jpg

水平.jpg

3楼2019-04-18 09:21:33

diwushuxue

铁虫 (初入文坛)

应助: 2 (幼儿园)
金币: 343.8
红花: 2
帖子: 39
在线: 32.8小时
虫号: 1243949
注册: 2011-03-24
专业: 应用地球物理学

引用回帖:

3楼: Originally posted by Ego悄吟 at 2019-04-18 09:21:33
您好，看到您的回帖之后，我回去划细网格了，之前网格大小为10m，计算步长为1.60e-3s,之后网格大小为2.5m，计算步长为0.8e-3s，可是依然出现这个情况。

3D.jpg

水平.jpg
...

还是正演问题，你空间网格10m，时间采样1ms，主频30Hz。再试一下！

4楼2019-04-18 17:34:15

Ego悄吟

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 114.7
帖子: 31
在线: 4.2小时
虫号: 14884382
注册: 2019-04-15
专业: 结构工程

引用回帖:

4楼: Originally posted by diwushuxue at 2019-04-18 17:34:15
还是正演问题，你空间网格10m，时间采样1ms，主频30Hz。再试一下！...

您好，根据您给的建议，我回去试了试，结果还是出现之前的问题。

3D.jpg

水平.jpg

5楼2019-04-19 09:57:57

baobiao007

木虫 (职业作家)

中国特色

应助: 201 (大学生)
金币: 6482.7
散金: 557
红花: 40
帖子: 3050
在线: 1009.9小时
虫号: 505962
注册: 2008-02-18
专业: 应用地球物理学

明显是正演的递推过程没写对。。

我同意叔本华的观点，人们投身艺术和科学领域的强烈愿望之一就是逃离痛苦、残酷和枯燥无味的现实生活，逃离自己飘忽不定的七情六欲的桎梏。--爱因斯坦

6楼2019-04-19 19:21:44

Ego悄吟

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 114.7
帖子: 31
在线: 4.2小时
虫号: 14884382
注册: 2019-04-15
专业: 结构工程

引用回帖:

6楼: Originally posted by baobiao007 at 2019-04-19 19:21:44
明显是正演的递推过程没写对。。

您好，您的意思是差分出错了么？我检查了我的程序，貌似并没有问题。

发自小木虫IOS客户端

7楼2019-04-20 08:39:25

baobiao007

木虫 (职业作家)

中国特色

应助: 201 (大学生)
金币: 6482.7
散金: 557
红花: 40
帖子: 3050
在线: 1009.9小时
虫号: 505962
注册: 2008-02-18
专业: 应用地球物理学

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
Ego悄吟: 金币+5, ★★★★★最佳答案 2019-04-21 09:30:08

送你段例子代码吧，直接运行即可

% Finite-differences in time domain(FDTD) acoustic wave propagation in 2D
% medium

close all;
% Output every ... time steps
IT_DISPLAY = 10;

%% MODEL
% Model dimensions
nx = 401;
nz = 401;
dx = 10; % [m]
dz = 10; % [m]

% Elastic parameters
vp = 3300.0 * ones(nz, nx); % velocity of compressional waves, [m/s]

%% TIME STEPPING
t_total = .55; % [sec] recording duration
dt = 0.7 * min(dx,dz)/max(vp(

); % min grid space / max velocity
nt = round(t_total/dt); % number of time steps
t = [0:nt]*dt;

CFL = max(vp(

) * dt / min(dx,dz);
%% SOURCE
f0 = 10.0;                         % dominant frequency of the wavelet
t0 = 1.20 / f0;                   % half Ricker wavelet excitation time
factor = 1e10;                   % amplitude coefficient
angle_force = 90.0;                % spatial orientation of source
                                 % (not relevant for acoustic case)
jsrc = round(nz/2);                % source location along OZ
isrc = round(nx/2);                % source location along OX

a = pi*pi*f0*f0;
dt2 = dt^2;
source_term = factor * exp(-a*(t-t0).^2);                         % Gaussian

force_x = sin(angle_force * pi / 180) * source_term * dt2 / (dx * dz);

min_wavelengh = 0.5*min(vp(vp>330))/f0;    % shortest wavelength bounded by velocity in the air

%% SUMMARY
fprintf('#################################################\n');
fprintf('2D acoustic FDTD wave propagation in isotripic \nmedium in displacement formulation\n');
fprintf('#################################################\n');
fprintf('Model:\n\t%d x %d\tgrid nz x nx\n\t%.1e x %.1e\t[m] dz x dx\n',nz, nx, dz,dx);
fprintf('\t%.1e x %.1e\t[m] model size\n',nx*dx, nz*dz);
fprintf('\t%.1e...%.1e\t[m/s] vp\n', min(vp(

), max(vp(

));
fprintf('Time:\n\t%.1e\t[sec] total\n\t%.1e\tdt\n\t%d\ttime steps\n',t_total,dt,nt);
fprintf('Source:\n\t%.1e\t[Hz] dominant frequency\n\t%.1f\t[sec] index time\n',f0,t0);
fprintf('Other:\n\t%.1f\tCFL number\n', CFL);
fprintf('\t%.2f\t[m] shortest wavelength\n\t%d, %d\t points-per-wavelength OX, OZ\n', min_wavelengh, floor(min_wavelengh/dx), floor(min_wavelengh/dz));
fprintf('#################################################\n');

%% ALLOCATE MEMORY FOR WAVEFIELD
p3 = zeros(nz+2,nx+2);          % Wavefields at t
p2 = zeros(nz+2,nx+2);          % Wavefields at t-1
p1 = zeros(nz+2,nx+2);          % Wavefields at t-2
% Coefficients for derivatives
co_dxx = 1/dx^2;
co_dzz = 1/dz^2;

%% Loop over TIME
tic;
for it = 1:nt
p3 = zeros(size(p2));
% Second-order derivatives
dp_dxx = co_dxx * (p2(2:end-1,1:end-2) - 2*p2(2:end-1,2:end-1) + p2(2:end-1,3:end));
dp_dzz = co_dzz * (p2(1:end-2,2:end-1) - 2*p2(2:end-1,2:end-1) + p2(3:end,2:end-1));
% U(t) = 2*U(t-1) - U(t-2) + G dt2/rho;
p3(2:end-1,2:end-1) = 2.0*p2(2:end-1,2:end-1) - p1(2:end-1,2:end-1) + (vp.^2).*(dp_dxx + dp_dzz).*dt2;
% Add source term
p3(jsrc, isrc) = p3(jsrc, isrc) + force_x(it);
% Exchange data between t-2 (1), t-1 (2) and t (3) and apply ABS
p1 = p2 .* weights;
p2 = p3 .* weights;
% Output
if mod(it,IT_DISPLAY) == 0
      fprintf('Time step: %d \t %.4f s\n',it, single(t(it)));
      imagesc(p3); colorbar;
      axis equal tight; colormap jet;
      title(['Step = ',num2str(it),'/',num2str(nt),', Time: ',sprintf('%.4f',t(it)),' sec']);
      drawnow;
end
end
toc; disp('End');

赞一下(1人)

8楼2019-04-20 10:36:59

baobiao007

木虫 (职业作家)

中国特色

应助: 201 (大学生)
金币: 6482.7
散金: 557
红花: 40
帖子: 3050
在线: 1009.9小时
虫号: 505962
注册: 2008-02-18
专业: 应用地球物理学

【答案】应助回帖

表情的地方是冒号 :

9楼2019-04-20 10:38:46

Ego悄吟

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 114.7
帖子: 31
在线: 4.2小时
虫号: 14884382
注册: 2019-04-15
专业: 结构工程

引用回帖:

8楼: Originally posted by baobiao007 at 2019-04-20 10:36:59
送你段例子代码吧，直接运行即可

% Finite-differences in time domain(FDTD) acoustic wave propagation in 2D
% medium

close all;
% Output every ... time steps
IT_DISPLAY = 10;

%% MODEL
% M ...

您好，谢谢您的代码，我的模型的问题确实出现在差分上。

发自小木虫IOS客户端

10楼2019-04-21 09:28:27