数字信号处理实验报告
实验一、地震子波波形显示及一维地震记录合成
一、实验目的
1、认识地震子波（以雷克子波为例），对子波有直观的认识。
2、利用线性褶积公式合成一维地震记录。
二、实验内容
1、 雷克子波：
()()t f e
t w m t f m πγπ2cos 22/2-=（零相位子波）、 ()()t f e t w m t f m πγ
π2sin 22/2-=（最小相位子波）， 其中m f 代表子波的中心频率，γ代表子波宽度,随着γ的增大，子波能量后移，当γ=7时，最小相位子波可视为混合相位子波，这里取m f = 25 Hz ，γ= 3;
2、 根据公式编程实现零相位子波、最小相位子波的波形显示；
3、 设计反射系数)(n r (n=500)，其中0.1)100(=r ，7.0)200(-=r ，5.0)300(=r ，4.0)400(=r ，6.0)500(=r ，其它为0;
4、 应用褶积公式∑=-=
*=N
m m n w m r n w n r n f 1)()()()()(合成一维地震记录，并图形显示； 5、 根据所学知识对实验结果进行分析。
三、实验结果：
1、零相位子波：
（1）程序源代码：
%编写零相位子波
t=0.002;
r=3;
fm=25;
for n=1:51
w(n)=exp(-(2*pi*fm/r)^2*(t*n)^2)*cos(2*pi*fm*t*n); end
plot(w)
xlabel('n')
ylabel('w')
title('零相位子波')
（2）图像：
2、最小相位子波：
（1）程序源代码：
%最小相位子波
t=0.002;
r=3;
fm=25;
for n=1:51
w(n)=exp(-(2*pi*fm/r)^2*(t*n)^2)*sin(2*pi*fm*t*n);
end
plot(w)
xlabel('n');ylabel('w');
title('最小相位子波')
（2）图像：
3、对比 不同时的波形图
（1）程序：
t=0.002;
r=3;
fm=25;
for n=1:51
w1(n)=exp(-(2*pi*fm/r)^2*(t*n)^2)*cos(2*pi*fm*t*n);
end
r=4;
for n=1:51
w2(n)=exp(-(2*pi*fm/r)^2*(t*n)^2)*cos(2*pi*fm*t*n);
end
r=5;
for n=1:51
w3(n)=exp(-(2*pi*fm/r)^2*(t*n)^2)*cos(2*pi*fm*t*n);
end
subplot(1,3,1),plot(w1);
axis([0,55,-0.7,1]);xlabel('n');title('r=3时');
subplot(1,3,2),plot(w2);
axis([0,55,-0.7,1]);xlabel('n');title('r=4时');
subplot(1,3,3),plot(w3);
axis([0,55,-0.7,1]);xlabel('n');title('r=5时');
（2）图像：
（3）分析：γ代表子波宽度，随着γ的增大，子波能量后移。
4、一维地震记录：
（1）零相位子波程序：
t=0.002;
r=3;
fm=25;
for n=1:51
w(n)=exp(-(2*pi*fm/r)^2*(t*n)^2)*cos(2*pi*fm*t*n);
end
%设置反射系数
r=zeros(500);
r(100)=1.0;
r(200)=-0.7;
r(300)=0.5;
r(400)=0.4;
r(500)=0.6;
%编写褶积公式
f=zeros(1,550);
for n=1:550
for m=1:500
if(1<=(n-m)&&(n-m)<=51)
f(n)=f(n)+r(m)*w(n-m);
end
end
end
plot(f)
（2）零相位子波图像：
（3）最小相位子波程序：
t=0.002;
r=3;
fm=25;
for n=1:51
w(n)=exp(-(2*pi*fm/r)^2*(t*n)^2)*sin(2*pi*fm*t*n);
end
%设置反射系数
r=zeros(1,500);
r(100)=1.0;
r(200)=-0.7;
r(300)=0.5;
r(400)=0.4;
r(500)=0.6;
%编写褶积公式
f=zeros(1,550);
for n=1:550
for m=1:500
if (1<=(n-m)&&(n-m)<=51)
f(n)=f(n)+r(m)*w(n-m);
end
end
end
plot(f)
(4)最小相位子波图像：
（5）对比
录： -0.6-0.4-0.20
0.2
0.40.60.8
最小相位子波的一维地震记录
250300350400450500
零相位子波的地震记录要滞后。最小相位子波的能量要稍小于零相位子波的能量。
程序：
t=0.002;
r=3;
fm=25;
for n=1:51
w1(n)=exp(-(2*pi*fm/r)^2*(t*n)^2)*cos(2*pi*fm*t*n);
w2(n)=exp(-(2*pi*fm/r)^2*(t*n)^2)*sin(2*pi*fm*t*n);
end
r=zeros(1,500);
r(100)=1.0;
r(200)=-0.7;
r(300)=0.5;
r(400)=0.4;
r(500)=0.6;
f1=zeros(1,550);
f2=zeros(1,550);
for n=1:550
for m=1:500
if(0<(n-m)&&(n-m)<=51)
f1(n)=f1(n)+r(m)*w1(n-m);
f2(n)=f2(n)+r(m)*w2(n-m);
end
end
end
plot(f1)
hold
plot(f2,'r')
title('对比最小相位和零相位子波的一维地震记录')
四、结果分析：
1、离散雷克子波，取采样间隔一般是2毫秒或者1毫秒，在这里我取了2毫秒。采样点数取51个。
2、γ代表子波宽度，随着γ的增大，子波能量后移。
3、褶积公式：∑
=-
=
*
=
N
m
m
n
w
m
r
n
w
n
r
n
f
1
)
(
)
(
)
(
)
(
)
(中，r的长度为M=500，w的长度为N=51，则f的长度为N+M-1=550。所以计算过程中应该用二重循环，外层是n从1到550，内层是m从1到500.同时循环过程要满足1<=n-m<=51。
实验二、带通滤波及频谱分析
一、实验内容
对复合频率信号进行频谱分析，并根据其振幅谱设计带陷滤波器，滤掉某些频率成分。
二、实验步骤
1、设计某一信号)(t x ，包含多种频率成分（可用雷克子波()()t f e t w m t f m πγπ2cos 2
2/2-=）； 2、将)(t x 离散，并应用fft 进行频谱分析，绘出振幅谱；
3、分析振幅谱有什么特点，在频率域设计带陷滤波器（可加斜坡），以消除某频段（大于62.5Hz ）的频率成分，并显示滤波后的振幅谱。（要求绘出滤波器图形）
4、将滤波后的信号反变换回时间域，并绘出信号曲线，观察其与原信号的差别。
5、根据所学知识对实验结果进行分析。
三、实验过程
以主频为25Hz 的雷克子波为例。
设置的子波长度为200，滤波器长度为51。
（1） 不加斜坡的滤波器
分析：关于N/2=128对称。但是吉卜斯现象波动明显。
程序如下：
clf
t=0.002;
r=3;
fm=25;
for n=1:200
x(n)=exp(-(2*pi*fm/r)^2*(t*n)^2)*cos(2*pi*fm*t*n);
end
figure(1)
plot(x)
title('主频为25Hz 的雷克子波')
for n=1:256
for n=1:200
x(n)=exp(-(2*pi*fm/r)^2*(t*n)^2)*cos(2*pi*fm*t*n);
end
for n=201:256
x(n)=0;
end
end
s=fft(x);
real_s=real(s);
imag_s=imag(s);
p=sqrt(real_s.^2+imag_s.^2);
figure(2)
plot(p)
title('快速傅立叶变换后的振幅谱') ylabel('振幅')
fx=62.5; %消除大于62.5Hz的成分
N=256;
ff=1/(N*t);
k1=fx/ff;
k2=N-k1;
for n=1:256
for n=1:k1
r(n)=1;
end
for n=k1-1:k2
r(n)=0;
end
for n=k2+1:N
r(n)=1;
end
xx(n)=r(n)*p(n);
end
figure(3)
plot(r);
title('滤波器图形')
xxx=abs(xx);
figure(4)
plot(xxx);
title('滤波后的振幅谱')
ylabel('振幅')
s_s=ifft(xxx);
figure(5)
plot(real(s_s))
title('滤波后的信号曲线')
（2）加斜坡的滤波器
程序：除滤波器处的程序不同外，其余相同，不做赘述。下面附上加斜坡的滤波器的程序如
下：
for n=1:256
for n=1:k1-5
r(n)=1;
end
for n=k1-4:k1-1
r(n)=8-k1/4;
end
for n=k1:k2
r(n)=0;
end
for n=k2:k2+4
r(n)=k2/4-56;
end
for n=k2+5:N
r(n)=1;
end
xx(n)=r(n)*p(n);
end
分析，在0和1过渡时没有直接跳跃，而是设计了一个线性函数，即加了一个斜坡过渡。
吉卜斯效应：当频率特性曲线是不连续函数而对滤波因子去有限项时，导致对应的频率
特征曲线是一波动的曲线，频率域滤波算子反生畸变。
（3）对比两种滤波器
（4）对比两种滤波后的信号曲线
（5）对比不同斜率的斜坡的滤波作用：
clc; %s1=input('请输入文件名: ','s'); fid=fopen('yy-10.txt','r'); c1=fscanf(fid,'%f'); N=length(c1); for i=2:3:N k=(i-2)/3+1; deltt(k)=c1(i); vv(k)=1000000/deltt(k); rr(k)=0.31*vv(k)^(1/4); Z(k)=vv(k)*rr(k); end n1=N/3; dp=360.2:0.2:2303 for k=1:n1-1 R(k)=(Z(k+1)-Z(k))/(Z(k+1)+Z(k)); end figure(9); plot(dp,R); %============================================================= %对反射系数序列进行低通滤波 %============================================================== r1=fft(R); r1(1001:8716)=0.0; figure(10); plot(abs(r1)); r2=ifft(r1); R1=real(r2); figure(11); plot(dp,R1); for i=1:n1-1 if(abs(R1(i))<0.01) R1(i)=0.0; end end figure(12); plot(dp,R1); f=30; wl=50; t=-wl:wl; deltt=0.002; b=(1-2*(pi*f*t*deltt).^2).*exp(-(pi*f*t*deltt).^2); figure(1);
数字信号处理实验报告 实验一、地震子波波形显示及一维地震记录合成、实验目的 1、认识地震子波(以雷克子波为例) ，对子波有直观的认识。 2、利用线性褶积公式合成一维地震记录。 、实验内容 1、雷克子波： w t e 2 fm / t cos 2 f m t (零相位子波)、 w t e 2 fm / t sin 2 f m t (最小相位子波)， 其中f m 代表子波的中心频率，代表子波宽度, 随着的增大，子波能量后移，当=7 时，最小相位子波可视为混合相位子波，这里取f m = 25 Hz ，= 3; 2、根据公式编程实现零相位子波、最小相位子波的波形显示； 3、设计反射系数r (n) (n=500) ，其中r (100) 1.0 ，r(200) 0.7 ，r(300) 0.5 ， r(400) 0.4 ，r (500) 0.6 ，其它为0; N 4、应用褶积公式f (n) r(n) w(n) r(m)w(n m) 合成一维地震记录，并图 m1 形显示； 5、根据所学知识对实验结果进行分析
三、实验结果： 1 、零相位子波： (1)程序源代码： %编写零相位子波 t=0.002; r=3; fm=25; for n=1:51 w(n )=exp(-(2*pi*fm/rF2*(t* nF2)*cos(2*pi*fm*t* n); end plot(w) xlabel( 'n' ) ylabel( 'w' ) title( ' 零相位子波' ) (2)图像： 2、最小相位子波：
1）程序源代码： %最小相位子波 t=0.002; r=3; fm=25; for n=1:51 w( n)=exp(-(2*pi*fm/rF2*(t* nF2)*si n(2*pi*fm*t* n); end plot(w) xlabel( 'n' );ylabel( 'w' ); title( '最小相位子波' ) (2)图像： 3、对比不同时的波形图 (1)程序： r=3;
版主你好，我有个问题想问问你，就是我在往discover中输入井位文件的时候（用的Microsoft Excel格式），点开wellbase模块，一输入数据就会显示打开数据库失败（当然是英文，我用中文表达的，呵呵），特征显示就是Microsoft Excel表格刚一出现黄框框的时候就会弹出个对话框说打开数据库失败！我也请教了许多人，该修改的地方都改好了，例如discover的服务是不是打开了，病毒软件是不是关了，等等，全部不起作用，我也不知道是什么原因，难道是这个软件和我的电脑不兼容？！？那也太郁闷了把~！ 希望版主能尽早给小弟我一个回复！先在这里谢了!~ 当你再出现如此错误的时候，请打开运行--〉odbcad32,第一个选项卡，选中gxdb,配置，connect test,看一下是否连通，如果不连通，那么选择database选项卡，看一下server name 是否为GGX_+n你的计算机名，network 中，tcp/ip中写入“ip=127.0.0.1",重新connect test,如果还是不能连通，确定你在安装discovery的时候，选择了所有用户都能使用，或者只限自己则没有更换登陆账户， 如果还是不行，请打开excel，设置一下安全性，使之能运行gxdb.xla宏 如果还是不行，我挂了 各位大侠，怎么在ＰＲＩＺＭ中怎么把岩心测试数据加入到测井曲线上谢谢了！ curves菜单中--〉core，选择show core curves,然后再某一线性道中添加岩性离散数据点，然后用aera fill 我是一个Dicovery软件初学者，现在遇到了一个比较棘手的问题。我目前研究的工区逆断层发育，有些井是过断层的，出现地层重复，在分层数据表中，同一地层名字下有两个深度数值，请问这种情况下，分层数据如何加载？ 在数据文件中增加一列，例如well2存在底层重复，则 wellname obs formati onname well2 ES1/default/1 ES1/defautl/topMD well2 Es1/default/2 ES1/default/topMD 我合并两个工区时遇到错误，两个有重叠的工区，线道号不连续，怎么合并到一起呢 如果重合的部分都拥有相同的线道号，例如seis1线号1-100，seis2线号50-200，并且道间距、旋转角度、采样间隔等都相同，也就是说，这两块地震体原来为一块统一的大地震体分割而成，那么加载时候如同线道连续的两块地震体 请教一个问题:如何在DISCOVERY中制作油藏剖面图? s首先在atlas中制作各层构造图，然后再xsection中的isomap选项喀中，把各个构造图选中，其余的和底层对比的设置相同‘（除了岩性填充）（、如果构造面上存在断层，其显示效果并不好，而如果手工绘制断层，那只是个示意图罢了） 请问： 我的断层文件由蓝马输出，好多条断层使用相同的断层名字，所以输入discovery中后，断层在剖面上显示比较乱。我想对输入数据进行重新命名，我的问题就是，怎么把一条合并的断层分开，并使之存在各自独立的名字？ 谢谢 fault manager中，interpolate 选择no 请教两个问题： 1、做连井对比时，能否加入岩性柱子？ 2、作连井是不是在选模版的时候不同井测井曲线的名称必须一致？因为电阻率测井有些井是Rt、Rs，有些井是M1R1等 1 可以加入岩性竹子如果从外界倒入的岩性曲线可以直接加入，如果在discovery解释模型中的岩性曲线，在log/wells选项卡的ude中选择你的解释模型 2、不是必须一致 例如你的电阻绿曲线有好多名字，那么就在default settings中设置别名，然后再模板中选择rt 3 跨带工区就是坐标系统选择莫卡特 mercator 投影系统
地震勘探原理实验一 地震子波波形显示及一维地震合成记录 姓名： 学号：专业：地球物理勘察技术与工程 级 一、实验目的 1. 认识子波，对子波的波形有直观的认识。（名词：零相位子波，混合相位 子波，最小相位子波；了解子波的分辨率与频宽的关系；） 2. 利用褶积公式合成一维地震记录。 二、实验步骤 1. 雷克子波 ()() ))(21(22 t f e t r m t f m ππ-=- 零相位子波 ())2sin() ln(222 t f e t w m n t f m π-= （最小相位子波） n= m1/m2为最大波 峰m1和最大波谷m2之比 ()())2cos(log *22xw t f e t w m m t f m +=-π 钟型子波 xw 为初相m 为时间域主波峰与次波峰之比 w(t)=exp(-2*Fm^2*t^2*ln(n))*sin(T-2*pi*Fm*t) n=m1/m2 最大相位子波 （最大相位子波请同学们自己查找相关文献完成，非必须完成）
其中 f代表子波的中心频率， t =i*dt，dt为时间采样间隔，i为时间 m 离散点序号; 这里可以为 f = 10，25，40，100 Hz等，采样间隔dt=0.002 m 秒，i为0～256; 2.根据公式编程实现不同频率的零相位子波的波形显示； 不同中心频率的零相位子波图 f = 25: m f = 100: m 3.其地质模型为：
设计反射系数)(n r (n=512)，n 为地层深度，其中0.1)100(=r ，为第一层介质深度；7.0)200(-=r ，为第二层介质深度；5.0)300(=r ，为第三层介质深度；4.0)400(=r ，为第四层介质深度；6.0)450(=r ，为第五层介质深度；其它为0。 地震波在介质中传播，当到达介质分界面时，发生反射和透射，反射波被检波器接受，生成地震记录。反射系数表示地震波在两层介质分界面的能量重新分配，如r(100)=1.0，表示地震波入射到分界面时，只有一种波，反射纵波（或反射横波）。反射系数不为1.0时，表示当地震波入射到分界面时，产生两种反射波。反射系数为正，表示反射波相位与入射波相位相差2π；反射系数为负，表示反射波相位与入射波相位相差π，存在半波损失。 4. 应用褶积公式∑=-=*=N m m n w m r n w n r n f 1)()()()()(合成一维地震记录，并图 形显示； 应用褶积公式求f （n ）的程序为： #include
#include % 地震合成记录 % 日期：07.07.19 % clc clear reply = input('请输入层数n(Default=5):','s'); %层数为n ifisempty(reply) n = 5; else n = sscanf(reply,'%f',[1 1]); end reply = input... ('请输入各层速度、密度及层厚(Defaul=[600 1000 1500 2000 2500;1500 1800 2000 2500 3000;500 700 400 300]):','s'); ifisempty(reply) V = [600 1000 1500 2000 2500]; dens = [1500 1800 2000 2500 3000]; %速度和密度v和den h = [500 700 400 300]; else clear a; a = sscanf(reply,'%f',[3 n]); V = a(1,:); dens = a(2,:); h = a(3,:); end % % 计算反射系数R % forilayer = 1:n-1 z1(ilayer) = V(ilayer) * dens(ilayer); z2(ilayer) = V(ilayer+1) * dens(ilayer+1); %各层反射系数R R(ilayer) = (z2-z1) / (z2+z1); end % % 计算各反射界面所对应的时间tlength % tlength(1) = 2*h(1)/V(1); forilayer = 2:n-1 tlength(ilayer) = tlength(ilayer-1) + 2*h(ilayer)/V(ilayer); end reply = input('请输入Ricker子波的频率f和采样间隔dt(Defalt=40 0.004):','s'); ifisempty(reply) f = 40; %子波频率f和采样间隔dt dt = 0.004; else clear a; a = sscanf(reply,'%f',[2 1]); f = a(1); dt = a(2); end % % 计算各反射界面所对应的采样点数nR % nsample = floor(tlength(n-1)/dt); forilayer = 1:n-1 nR(ilayer) = floor(tlength(ilayer)/dt); end % % 形成反射系数序列RR % RR(1:2*nsample) = 0;%？这个地方反射系数的长度应该是nsample/2 forilayer = 1:n-1 RR(nR(ilayer)) = R(ilayer); %只有在有界面的地方反射系数才有值end %subplot(2,2,1); stem(RR); title('反射系数序列'); % % 形成一个Ricker子波wavelet % wavelet = ricker(f,dt); fori = 1:length(wavelet); Landmark 快捷键 Q键：“+”指针，主要用于显示位置。若要取消，可连续点击三次Q键即可。点击第一次Q键，指针变为”|”;点击第二次Q键，指针变为”-“;点击第三次指针，即可取消剖面上的”+”指针。 W键：隐藏解释的层位。若要显示，可重新刷新即可。 Y键：沿某一层位拉平地震数据体（注意：只针对当前的剖面）；如果要取消拉平的话，再次点击Y键or Ctrl+Y。 A键:auto-dip,层位解释时，自动倾角追踪方式。 P键: 层位解释时，手动追层位。 F:断层解释快捷键。 H:层位解释快捷键。 K键：层位解释时，各种追踪方式的参数选择快捷键。 L键：激活层位。 Z键：地震剖面显示局部放大（选中Z键后，选”+”剖面局部放大，选”-“剖面局部缩小） C键：层位解释时，追踪模式为相关模式(correlation). V键：剖面显示时，剖面的CDP大小方向反转。 B键：记录在seiswork中的各个操作步骤的标签。 / 键：在解释断层和层位时，计算断层的断距。 Ctrl z键：可快速地选择纵向的时间范围。 Copy键：复制terminal窗口中的内容（命令或者路径等） Paste键：粘贴copy键所拷贝的东西到所需要的地方。 补充：i键：在层位解释状态下可激活层位。 注意：键盘上的快捷键功能可以改变， 在seiswork主菜单下> Defaul t s >expert keys下 设置键盘快捷键，将键盘字母与需要的功能选中，应用一下即可； 如果不需要键盘快捷键功能，只需要选中键盘字母，unset即可取消快捷径功能。 当所有的快捷功能都设置完后，在expert keys窗口下保存一下，存成一个后缀为.xks的文件。下次打开地震工区，进入在seiswork主菜单下> Defaul t s >expert keys下，重新加载.xks文件，即可实现键盘指定的快捷功能。 合成记录“二次标定法则”及其在地震解释中的应用 栗宝鹃①董春梅②宋亚民③张木森④刘斌⑤ （①②③④中国石油大学（华东）地球科学与技术学院，山东青岛，266580；⑤中石化胜利 油田河口采油厂，山东东营，257200） 摘要：人工合成记录是联系测井和地震的桥梁，精确的合成记录标定在精细地震解释中起着至关重要的作用。影响合成记录的因素主要有反射系数和地震子波，反射系数由声波数值计算得到，人为干涉的作用不大，为此，子波是影响合成记录标定的重要因素。由于子波旁瓣较多，导致合成记录波形改变、波组增多，影响标定的准确性。本文提出一种合成记录的“二次标定法则”，该方法先用与地震数据相同频率的标准的雷克子波进行初次标定，标定准确主要波组；然后用井旁道子波进行二次标定，对初次标定的结果做细微调整，由此得到准确的合成记录标定结果。 关键词：合成记录子波地震解释 1 引言 地震解释中，地震反射同相轴的地质层位和岩性意义是通过合成记录标定得到的。合成记录是联系地震和测井的桥梁，在正确进行地震解释中起着至关重要的作用。 合成记录在地震解释中的作用体现在以下方面：（1）精确的合成记录标定可以建立测井和地震之间的准确对应关系[1]，由此可以根据合成记录来确定要追踪的层位和辅助进行层位的闭合解释；（2）地震子波是有极性的，由地震子波的极性很容易定义地震剖面的极性[2]。因此，根据合成记录标定过程中提取的子波，可以进行地震剖面的极性判断。 2 合成记录原理 制作人工合成记录，首先利用声波和密度测井资料求取反射系数序列，然后将反射系数序列与地震子波褶积获得[3]。具体公式如下： (t)(t)*r(t) s w =（1）其中，(t) s为合成地震记录，r(t)为反射系数序列，(t) w为地震子波。 由公式（1）可以看出，合成记录的好坏与反射系数和子波有关。 一个界面的反射系数是由上下两层的波阻抗得到的，其表达式为： 1 1 2 2 1 1 2 2 v v v v R ρ ρ ρ ρ + - =(2) 式中，R为反射系数， 2 1 ρ ρ,为上、下 两层的密度， 2 1 v v,为上、下两层的速度。 由公式（2）可知，反射系数由两层的速度和密度参数共同决定的。通过声波曲线可以得到速度值。 本文部分内容来自网络整理，本司不为其真实性负责，如有异议或侵权请及时联系，本司将立即删除！ == 本文为word格式，下载后可方便编辑和修改！ == 蓝马R5000软件学习个人总结 软件竞技之LandmarkR5000初学之我见 我在这里要分享的是我对LMR5000学习的一点心得和经验，首先声明，R5000 我用的不是很好，常用的一些模块和操作还行，更高级一点就不是那么熟悉了，只能算作是入门吧，因为201X上的有些东西都还没搞明白，然后到了R5000， 一些数据管理和操作的界面都发生了改变，特别是操作界面有了很大的变化。 下面简要谈谈自己的一些心得体会，和大家分享，有什么不妥或者错误的地方，还望大家批评指正，同时也希望能够给初学者一点启示。 1、对Landmark的基本操作要有所了解，比如数据的加载，这里面包括建立投影系统、建 立OpenWorks数据库、建立井的资料等等，这些都是一些最基本的东西，这些 要是有问题，那以后的工作都是在错误的基础来完成了。（曾见过有人竟然把 投影系统就给弄错了的）这样做的工作简直就是稀里糊涂的。 2、基础资料要整理好，特别是R5000对数据管理更加严格了，这里说的基础 资料最主要的 还是井的资料，比如井的坐标、钻井分层、时深表、井轨迹、测井曲线等。这 些都是基础数据，要把这些东西整理好了，以后干活就会减少很多的麻烦了。 3、了解R5000相对201X的一些变化，R5000对OpenWorks/SeisWorks工区数 据环境做了 很多变化，最明显的是：没有了地震工区的概念,取而代之的是解释工区；2D 地震层位存储于Oracle中；所有地震数据和层位在OpenWorkss中管理；解释 数据和地震数据都有了版本号；对数据的大小的限制取消了。 4、简要了解R5000对数据的管理，以前版本的地震数据及解释数据等很多东 西都放在了 project下面，现在有了改变，地震数据和层位数据放在了所建的survey目录下，而其它的信息都放在project下面，如色棒、格式定义等，在SWDATA目录下，而且在project下面每个工区都对应着你所使用过的模块的信息及操作。 人工合成地震记录作业-CAL-FENGHAI.-(YICAI)-Company One1 人工合成地震记录程序设计 （一）、人工合成地震记录原理： 地震记录上看到的反射波波形是地震子波在地下各反射界面上发生反射时形成的。反射波的振幅有大有小（决定于界面反射系数的绝对值）、极性有正有负（取决于反射系数的正负）、到达时间有先有后（取决于反射界面的深度）的地震反射子波叠加的结果。 如果地震子波的波形用S (t )表示，地震剖面的反射系数为双程垂直反射时间t 的函数，用R (t )表示，那么反射波地震记录形成的物理过程在数学上就可以用S (t )的R (t )的褶积表示，即某一时刻的反射波地震记录f (t )是: )()()(t R t S t f *= 其离散形式为： ))(()()(1 t m n R t m S t n f M m ?-??=?∑= 如果大地为多层介质，在地面记录长度内可接收的反射波地震记录为： ))(()()(11t m n R t m S t n f M m N n ?-??=?∑∑== 式中，n 为合成地震记录的采样序号，n =1,2,3...N ；N 为合成一道地震记录的采样点数；m =1,2,3...M ，为离散子波的采样点数；△t 为采样间隔。 这种褶积模型将地震波的实际传播过程进行了简化： 1、在合成地震记录的过程中没有考虑大地的吸收作用，所有薄层的反射波都与地震子波的形式相同，只是振幅和符号不同。 2、假设地震波垂直入射到界面上，并原路径返回。 3、假设地层横向是均匀的，在深度（纵向）方向上假设密度为常数，只是速度发生变化。 4、不考虑地震波在传播过程中的透射损失。 （二）、人工合成地震记录的方法 1、 反射系数序列 在有速度测井资料的情况下，可以用速度曲线代替波阻抗曲线，计算反射系数序列。在没有速度资料的情况下，可根据干扰波调查剖面分析的结果设计地质模型。 如设计的地质模型如图a 所示，图中H 为层厚度，V 为层速度，根据下式计算反射系数： 1 1)(--+-=N N N N N V V V V H R 式中H 为反射界面的深度，N 为反射层序号，随深度变化的反射系数序列如图b 所示。但褶积计算中需要与时间有关的反射系数，深度与时间的转换可用下列公式计算： t V H n n R H R ?=→111112),()( 1222222),()(n t V H n n R H R +?= → 蓝马Landmark R5000地震解释软件的tar包安装过程2012-03-07 23:27:44| 分类：2_工作| 标签：蓝马安装|字号大中小订阅 安装蓝马是件疼苦的事儿，下面是安装过程。 1.解压两个tar文件 tar文件在dell01的什么位置就需要把它解压到相应机子的哪个位置。 在apps目录下解压OpenWorks.tar和ora10g.tar tar –xvf OpenWorks.tar tar –xvf ora10g.tar 2.修改/etc/passwd文件 在末尾位置加两行： ora10g::105:500:Oracle Database Account:/apps/ora10g:/bin/csh owinstall::107:0:Oracle Database Account:/apps/OpenWorks:/bin/csh 3.修改/etc/group文件 在末尾位置加一行 dba:x:500:ora10g 4.在/etc/sysctl.conf中增加几行： 在末尾加下面内容 # DO NOT REMOVE: START OpenWorks Linux Configuration kernel.shmall = 2097152 kernel.shmmax = 4294967295 kernel.shmmin = 1 kernel.sem = 250 125000 200 128 fs.file-max = 65536 net.ipv4.ip_local_port_range = 1024 65000 net.core.rmem_default = 262144 net.core.rmem_max = 262144 人工合成地震记录程序设计 （一）、人工合成地震记录原理： 地震记录上看到的反射波波形是地震子波在地下各反射界面上发生反射时形成的。反射波的振幅有大有小（决定于界面反射系数的绝对值）、极性有正有负（取决于反射系数的正负）、到达时间有先有后（取决于反射界面的深度）的地震反射子波叠加的结果。 如果地震子波的波形用S (t )表示，地震剖面的反射系数为双程垂直反射时间t 的函数，用R (t )表示，那么反射波地震记录形成的物理过程在数学上就可以用S (t )的R (t )的褶积表示，即某一时刻的反射波地震记录f (t )是: )()()(t R t S t f *= 其离散形式为： ))(()()(1 t m n R t m S t n f M m ?-??=?∑= 如果大地为多层介质，在地面记录长度内可接收的反射波地震记录为： ))(()()(11t m n R t m S t n f M m N n ?-??=?∑∑== 式中，n 为合成地震记录的采样序号，n =1,2,3...N ；N 为合成一道地震记录的采样点数；m =1,2,3...M ，为离散子波的采样点数；△t 为采样间隔。 这种褶积模型将地震波的实际传播过程进行了简化： 1、在合成地震记录的过程中没有考虑大地的吸收作用，所有薄层的反射波都与地震子波的形式相同，只是振幅和符号不同。 2、假设地震波垂直入射到界面上，并原路径返回。 3、假设地层横向是均匀的，在深度（纵向）方向上假设密度为常数，只是速度发生变化。 4、不考虑地震波在传播过程中的透射损失。 （二）、人工合成地震记录的方法 1、 反射系数序列 在有速度测井资料的情况下，可以用速度曲线代替波阻抗曲线，计算反射系数序列。在没有速度资料的情况下，可根据干扰波调查剖面分析的结果设计地质模型。 如设计的地质模型如图a 所示，图中H 为层厚度，V 为层速度，根据下式计算反射系数： 1 1)(--+-=N N N N N V V V V H R 式中H 为反射界面的深度，N 为反射层序号，随深度变化的反射系数序列如图b 所示。但褶积计算中需要与时间有关的反射系数，深度与时间的转换可用下列公式计算： t V H n n R H R ?=→111112),()( 1222222),()(n t V H n n R H R +?= → 引言： 随着社会的快速发展，核电站和海洋平台迅速兴建，大型水坝、高层建筑和大跨桥梁日益增加，其中很多兴建于强震活动区。由于这些结构物的重要性，加之缺乏这类工程及相应场地的抗震经验，对其抗震性能的研究引起了社会和工程界的重视。并且相应的抗震规范都规定，在上述重要结构的设计中，应当采用地震动时程输入结构动力分析来考虑地震动时间过程影响。由于很难在天然地震中取得相应场的地峰值和反应谱，为探讨结构物在地震动反应中的耗能特性和破坏机理，必须对结构物在地震动作用下的整个过程进行模拟，用人工合成地震动方法，分析结构物及相应场地在地震动中的反应，因此对比天然地震动与人工合成地震动在相应场地反应的异同成为抗震设防的重点。 1、选取天然地震波 本论文所用的天然波取自1976年8月9日06：41唐山大地震中的一次5.7级余震记录，由于记录地点在迁安地震台，因此通常被称为“迁安波”。经过校正加速度记录信息如下：南北向记录，时间间隔0.01s，记录2320个点，持续时间23.19s，峰值为158.62gal，出现在2.37s。迁安波时程曲线如图1-1所示，其反应谱如图1-2所示。 图1-1迁安波地震记录 图1-2迁安波地震反应谱 2、人工合成地震动 在工程地震学研究中，采用多种方法来估计地震动。其中包括基于幅值和卓越周期调整的比例方法、拟合目标峰值和反应谱的数值方法、选择实际地震动记录的地震记录匹配法及半经验半理论模拟方法。 下面我们采用拟合目标峰值和反应谱的数值方法，进行人工合成地震动。随着强震动观测的发展及对地震宏观震害经验和仪器测量结果的大量分析研究发现，运用数值方法程序计算出的反应谱和加速度时程，可以通过地震动的工程特性三要素来描述即：地震动的振幅、频率和持续时间。因此我们通过控制这三要素，运用Saw软件，更改随机数200，得到图2-1和图2-2如下： 图2-1人工合成地震动时 clear; clc; %s1=input('请输入文件名: ','s'); fid=fopen('','r'); c1=fscanf(fid,'%f'); N=length(c1); for i=2:3:N k=(i-2)/3+1; deltt(k)=c1(i); vv(k)=1000000/deltt(k); rr(k)=*vv(k)^(1/4); Z(k)=vv(k)*rr(k); end n1=N/3; dp=::2303 for k=1:n1-1 R(k)=(Z(k+1)-Z(k))/(Z(k+1)+Z(k)); end figure(9); plot(dp,R); %============================================================= %对反射系数序列进行低通滤波 %============================================================== r1=fft(R); r1(1001:8716)=; figure(10); plot(abs(r1)); r2=ifft(r1); R1=real(r2); figure(11); plot(dp,R1); for i=1:n1-1 if(abs(R1(i))< R1(i)=; end end figure(12); plot(dp,R1); f=30; wl=50; t=-wl:wl; deltt=; b=(1-2*(pi*f*t*deltt).^2).*exp(-(pi*f*t*deltt).^2); figure(1); plot(b); x=conv(R1,b); xl=length(x); rl=wl; for i=rl:xl-rl-1 x1(i-rl+1)=x(i); end figure(2); plot(dp,x1); 数字信号处理实验报告 实验一、地震子波波形显示及一维地震记录合成 一、实验目的 1、认识地震子波（以雷克子波为例），对子波有直观的认识。 2、利用线性褶积公式合成一维地震记录。 二、实验内容 1、 雷克子波： ()()t f e t w m t f m πγ π2cos 2 2/2-=（零相位子波）、 ()()t f e t w m t f m πγ π2sin 2 2/2-=（最小相位子波）， 其中m f 代表子波的中心频率，γ代表子波宽度,随着γ的增大，子波能量后移，当γ=7时，最小相位子波可视为混合相位子波，这里取m f = 25 Hz ，γ= 3; 2、 根据公式编程实现零相位子波、最小相位子波的波形显示； 3、 设计反射系数)(n r (n=500)，其中0.1)100(=r ，7.0)200(-=r ，5.0)300(=r ， 4.0)400(=r ，6.0)500(=r ，其它为0; 4、 应用褶积公式∑=-= *=N m m n w m r n w n r n f 1 )()()()()(合成一维地震记录， 并图形显示； 5、 根据所学知识对实验结果进行分析。 三、实验结果： 1、零相位子波： （1）程序源代码： %编写零相位子波 t=0.002; r=3; fm=25; for n=1:51 w(n)=exp(-(2*pi*fm/r)^2*(t*n)^2)*cos(2*pi*fm*t*n); end plot(w) xlabel('n') ylabel('w') title('零相位子波') （2）图像： 2、最小相位子波： （1）程序源代码： %最小相位子波 t=0.002; r=3; fm=25; for n=1:51 w(n)=exp(-(2*pi*fm/r)^2*(t*n)^2)*sin(2*pi*fm*t*n); end plot(w) xlabel('n');ylabel('w'); title('最小相位子波') （2）图像：合成地震记录
蓝马快捷键
合成记录在地震解释中的应用
【推荐】蓝马R5000软件学习个人总结word版本 (3页)
人工合成地震记录作业
蓝马Landmark R5000地震解释软件的tar包安装过程
人工合成地震记录作业
人工合成地震动
matlab编程合成地震记录
数字信号处理(合成地震记录&fft)
相关主题
文本预览