Test

这篇文章创建所有需要的Categories

Some contents in the article.

When I was a boy, my grandfather told me how a German taxi driver, Franz Bussman, found a brother who
was thought to have been killed twenty years before.

$ \hat{\mathbf{Z}}^{\dagger} \mathbf{J}=\mathbf{J}\left(\hat{\mathbf{Z}}^{\prime}\right)^{-1} $

If the math formula can be shown correctly.

二级题目测试

是否在右边栏显示2级标题？—— 会的！Categories

Code

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
%% Information
% Title: AESM1590-EX4
% Time: 30-11-2022
% Name: Boyu Gao + Miaoyang Yuan
% Group: 7

%% Clear
clear all;
clear variables;
clc;


%% Read data
% Note: geo(3) = 60: number of shots ns. Finally, we will plot all 60 shots. 
filename_data = 'slopesurf.segy';
[Data_ex4, H_CSP, geo] = segyread(filename_data);
[Data_CMP, H_CSP_CMP] = sortdata(Data_ex4, H_CSP, 5); 


%% Velocity parameters
% midpoint
midpnt = [1200, 2000, 2600];
% when cmp = 1200
t_1 = [74, 227, 296].*(geo(1)*10^(-3));
v_1 = [2010, 2010, 1800];
% when cmp = 2000
t_2 = [215,262,305,384].*(geo(1)*10^(-3));
v_2 = [3500,2000,2500,3900];
%when cmp = 2600
t_3 = [262,296,348,384].*(geo(1)*10^(-3));
v_3 = [2000,2200,3900,3900];
t = {t_1,t_2,t_3};
v = {v_1,v_2,v_3};
tvpicks(:,:,1) = [74, 227,296,296;
                  215,262,305,384;
                  262,296,348,384].*(geo(1)*10^(-3));  % add the same number of last one in t_1 to make it fit

tvpicks(:,:,2) = [2010,2010,1800,1800;
                  3500,2000,2500,3900;
                  2000,2200,3900,3900];                % add the same number of last one in v_1 to make it fit
smute = 0.05;  % 0.5 / 0.05


%% NMO
for mid_tmp = 1:3
    [Data_CMP_out,H_CMP_out] = selectcmp(Data_CMP, H_CSP_CMP, midpnt(mid_tmp));
    t_matrix = cell2mat( t(mid_tmp) );
    v_matrix = cell2mat( v(mid_tmp) );
    nmo_vlog(Data_CMP_out, H_CMP_out, geo, t_matrix, v_matrix);
end


%% Find CMP locations + Generate velocity model  
[gather_positions, gather_folds] = analysefold(H_CSP, 5);
vmodel = generatevmodel2(midpnt, tvpicks, gather_positions, geo);


%% Stacking the NMO-ed data
[Data_CMP_out,H_CMP_out] = selectcmp(Data_CMP, H_CSP_CMP, midpnt(1));
gather = nmo_vlog(Data_CMP_out,H_CMP_out,geo,t_matrix,v_matrix,smute);
gather_stack = stackplot(gather,geo);
title('without smute');


%% NMO to all the data and stack them without smute
figure();
[cmp_positions, cmp_folds] = analysefold(H_CSP_CMP, 5);
zosection = nmo_stack(Data_CMP, H_CSP_CMP, cmp_positions, cmp_folds, geo, vmodel, 0);
imageseis(zosection);
xlabel(' the spacing between CMPs ');
title('without smute');


%% NMO to all the data and stack them with smute
figure();
zosection = nmo_stack(Data_CMP,H_CSP_CMP, cmp_positions, cmp_folds, geo, vmodel, smute);
imageseis(zosection);
xlabel(' the spacing between CMPs ');
title('with smute');