地震数值模拟实验报告
本科生实验报告
实验课程数值模型模拟
学院名称地球物理学院
专业名称勘查技术与工程
学生姓名ZRY
学生学号
指导教师
实验地点624
实验成绩
二〇一五年4月二〇一五年5月
成都理工大学
《地震数值模拟》实验报告
实验二叠加地震记录的相移波动模拟实方程正演验
摘要
利用C语言编制地质模型的相移波动方程正演模拟，改变绕射点位置、速度，再做正演模拟。
关键字：地震模型；正演记录
1.1实验目的
掌握各向同性介质任意构造、水平层状速度结构地质模型的相移波动方程正演模拟基本理论、实现方法与程序编制，由正演记录初步分析地震信号的分辨率。
1.2实验内容
1、基本要求：
（1）点绕射构造和水平层状速度模型（参数如图1 所示）的正演数值模拟；
1）削波的正演；
2）无削波的震正演；
（2）计算中点和两个边界的信号位置，分析实验结果的正确性；
（3）做同样模型的褶积模型数值模拟，对比分析分析两者的异同。
（4）改变绕射点位置、速度，再做正演模拟。
2、较高要求：
（1）使用雷克子波做爆炸源，对三个不同的主频：25hz、50hz 和75hz 分别做点绕射
模型的正演模拟；
（2）设计复杂反射构造模型，再做正演模拟。
1.3实验原理
1、地震波传播的波动方程
设（x,z）为空间坐标，t 为时间，地震波传播速度为v(x,z)，则二位介质中任意位置、任意时刻的地震波场为p(z,x,t)：压缩波——纵波。则二维各向同性均匀介质中地震波传播
的遵循声波方程为
（）
2、傅里叶变换的微分性质
p(t)与其傅里叶变换的P(ω)的关系：
则有时间微分性质
ω 为频率，ω=2π/T，T 为周期。同理有空间微分性质：
k 为频率，k=2π/λ，λ为波长。
3、地震波传播的相移外推公式
令速度v 不随x 变化，只随z 变化，则利用傅里叶变换微分性质（3）和（4）式，把波动方程(1)式变换到频率-波数域，得：
或：
令：
则（5）式的解为：
包括上行波和下行波两项。正演模拟取上行波：
若和间隔为△Z ，速度v(z) 在此间隔内不随Z 变的常数,(7)式实现波场从到的延拓，即：
在深度Zj+1 开始向上延拓到Zj，若延拓深度为零，即：?Z= Z j+1-Z j=0，则
对于任意深度Z j+1 到Z j 的延拓，可得正演模拟中地震波的传播方程（延拓公式）
4、初始条件和边界条件
按照爆炸界面理论，反射界面震源在t=0 时刻同时起爆，此时刻的波场就是震源。根据不同情况，可直接使用反射系数脉冲或子波作震源。如果直接使用反射系数作震源脉冲，则初始条件可表示为：
其他
对时间t和空间x做二维傅立叶变换，则得频率-波数域的初始波场
边界条件：
其他，其他，其他
其他参数都是在范围内定义的。
5、边界处理
（1）边界反射问题
把实际无穷空间区域中求解波场的问题化为有穷区域求解时，左右两边使用零边界条件。物理上假设探区距和两个端点很远，在两个端点上收到的反射波很弱。但是，上述条件在实际中不能成立，造成零边界条件反而成为绝对阻止波通过的强反射面。在正演模拟的剖面上出现了边界假反射干涉正常界面的反射。
（2）边界强反射的处理
镶边法、削波法、吸收边界都能有效消除边界强反射。
削波法就是在波场延拓过程中，没延拓一次，在其两侧均匀衰减到零，从而消除边界强反射的影响。假设横向总长度为NX，以两边Lx 道吸波为例，有以下吸波公式：
其他
6、数字化
根据数字信号处理的采样定理，把连续的信号变为计算机能处理的数字信号，使相移法正演模拟得以实现。
频域抽样定理：一个频谱受限信号f(t)，如果时间只占据-t m—+t m 的范围，若在频域以不大于1/2t m 频率间隔?f≤1/2t m 对信号f(t)的频谱F(ω)采样，则抽样到的离散信号F1(ω)可以唯一表示原信号。
时域抽样定理：一个时间受限信号f(t)，如果频谱只占据-ωm—+ωm 的范围，则信号f(t)可以用等间隔的抽样值唯一表示出来，而时间?t 抽样间隔必须不大于1/2f m，ωm =2π f m，?t≤1/2f m。
1.4实验过程
//#include "PsFrwrdMdlParameter.h"
#include
#include
#include
#include
#include
int kkfft(float pr[], float pi[], int n, int l);
int Absorb(int Labs);
int Po2Judge(int);
int Rflct();
int Vlcty();
int WvFld0();
int PsFrwd();
int exp_ikzDz(float eikzdz[],int Ix, float Vc,int Iw,float Dw,float Dkx);
int PhaseShift();//相移
int Frqcy2Time();
#define Nx 128 //Trace Number
#define Nt 256 //Record Number
#define Nz 100 //Depth Number
#define Dx 20. //Trace Interval
#define Dt 0.004//Record Interval
#define Dz 20. //Depth Interval
#define PI 3.1415
void main()
{
int Labs=10;//Length Of Boundary Absorbing
if(Po2Judge(Nt) !=1){printf("Nt=%d is not the Power of 2\n",Nt);exit(0);} else {printf("Nt=%d is the Power of 2\n",Nt);}
if(Po2Judge(Nx) !=1){printf("Nx=%d is not the Power of 2\n",Nx);exit(0);} else {printf("Nx=%d is the Power of 2\n",Nx);}
if(Absorb(Labs) !=1){printf("Absorb is error\n"); exit(0);} else {printf("Absorb is okay\n");}
if(Rflct() !=1){printf("Reflection is error\n"); exit(0);} else {printf("Reflection is okay\n");};
if(Vlcty() !=1){printf("Vlcty is error\n"); exit(0);} else {printf("Vlcty is okay\n");}
if(WvFld0() !=1){printf("WvFld is error\n"); exit(0);} else {printf("WvFld is okay\n");}
if(PsFrwd() !=1){printf("PsFurd is error\n"); exit(0);} else {printf("PsFurd is okay\n");}
remove("Wfld0r.dat");
remove("Wfld0i.dat");
remove("PsFrwrdMdl.ncb");
remove("PsFrwrdMdl.dsp");
remove("PsFrwrdMdl");
remove("Abs");
// return( );
}
//////////////////////////////////////////////////
int Po2Judge(int N)//Judge Nt or Nx is Power of 2
{
double a=0;
for(int i=0;a { a=pow(2,i); } if(fabs(N-a)>=1)return(0); else return(1); } ////////////////////////////////////////////////// int Absorb(int Labs) //Forming Absorbing Boudary { FILE *fp_Absorb,*fp_Abs; int Ix; float Abs[Nx]; if((fp_Absorb=fopen("Absb.dat","wb"))==NULL) printf("Cannot open file ""Absorb"""); else {printf("Absb is okay\n");} if((fp_Abs=fopen("Abs.txt","w"))==NULL) printf("Cannot open file ""Abs"""); else {printf("Abs is okay\n");} for(Ix=0;Ix { Abs[Ix]=1.; } for(Ix=0;Ix { Abs[Ix]=sqrt(sin(PI*Ix/(2*(Labs-1)))); Abs[Nx-Ix-1]=Abs[Ix]; } for(Ix=0;Ix { fwrite(&Abs[Ix],sizeof(Abs[Ix]),Nx,fp_Absorb); //printf("%f\n",Abs[Ix]); fprintf(fp_Abs,"%f\n",Abs[Ix]); } fclose(fp_Absorb); fclose(fp_Abs); return(1); } ////////////////////////////////////////////////// int Rflct()//Forming Reflect Structure Model { FILE *fp_Rflct; int Ix,Iz; float Rflct[Nz]; if((fp_Rflct=fopen("Rflct.dat","wb"))==NULL) printf("Cannot open file""Reflection"""); for(Ix=0;Ix { for(Iz=0;Iz { Rflct[Iz]=0.0; if(Ix==(int)(Nx/2)&&Iz==(int)(Nz/2)) Rflct[Iz]=1.; fwrite(&Rflct[Iz],sizeof(Rflct[Iz]),1,fp_Rflct); } } fclose(fp_Rflct); return(1); } ////////////////////////////////////////////////// int Vlcty() //Forming Velociy Model { FILE *fp_Vlcty; int Ix,Iz; float Vlcty[Nz]; if((fp_Vlcty=fopen("Vlcty.dat","wb"))==NULL)printf("Cannot open file ""Vlcty"""); for(Ix=0;Ix { for(Iz=0;Iz<2*Nz/3;Iz++) { Vlcty[Iz]=3000.; // fwrite(&Vlcty[Iz],sizeof(Vlcty[Iz]),1,fp_Vlcty); } for(Iz=(2*Nz/3);Iz { Vlcty[Iz]=6000.; // fwrite(&Vlcty[Iz],sizeof(Vlcty[Iz]),1,fp_Vlcty); } fwrite(&Vlcty[0],sizeof(Vlcty[Iz]),Nz,fp_Vlcty); } fclose(fp_Vlcty); return(1); } ////////////////////////////////////////////////// int WvFld0()//Wave Field Initialization { FILE *fp_Rflct,*fp_Wfld0r,*fp_Wfld0i,*fp_Wfld0isee,*fp_Wfld0rsee,*fp_Wfld0r0see; int Ix,Iz,It; float Rflct[Nz]; float Wfld0r[Nt],Wfld0i[Nt]; if((fp_Wfld0r=fopen("Wfld0r.dat","wb"))==NULL)printf("Cannot open Wfld0r.dat"); if((fp_Wfld0i=fopen("Wfld0i.dat","wb"))==NULL)printf("Cannot open Wfld0i.dat"); if((fp_Rflct=fopen("Rflct.dat","rb"))==NULL)printf("Cannot open Rflct.dat"); //if((fp_Wfld0isee=fopen("fp_Wfld0isee.txt","wb"))==NULL)printf("Cannot open fp_Wfld0isee.dat"); //if((fp_Wfld0rsee=fopen("fp_Wfld0rsee.txt","wb"))==NULL)printf("Cannot open fp_Wfld0rsee.dat"); //if((fp_Wfld0r0see=fopen("fp_Wfld0r0see.txt","wb"))==NULL)printf("Cannot open fp_Wfld0r0see.dat"); for(Ix=0;Ix { //printf("Wavefield0_FFT:Ix=%d\n",Ix); for(Iz=0;Iz { fread(&Rflct[Iz],sizeof(Rflct[Iz]),1,fp_Rflct); for(It=0;It { Wfld0r[It]=0.; Wfld0i[It]=0.; if(It==0) Wfld0r[It]=Rflct[Iz]; } //fprintf(fp_Wfld0r0see,"%f \t",Wfld0r[Iz]); if(kkfft(Wfld0r,Wfld0i,Nt,0)!=1){printf("FFT is error");exit(0);}//原为时间域函数，变后为频率域。 for(It=0;It<(int)(Nt/2)+1;It++) { fwrite(&Wfld0r[It],sizeof(Wfld0r[It]),1,fp_Wfld0r); fwrite(&Wfld0i[It],sizeof(Wfld0i[It]),1,fp_Wfld0i); //fprintf(fp_Wfld0isee,"%f \t",Wfld0i[It]); //fprintf(fp_Wfld0rsee,"%f \t",Wfld0r[It]); } } //fprintf(fp_Wfld0isee,"\n"); //fprintf(fp_Wfld0rsee,"\n"); //fprintf(fp_Wfld0r0see,"\n"); } fclose(fp_Rflct); fclose(fp_Wfld0r); fclose(fp_Wfld0i); return(1); } ////////////////////////////////////////////////// int PsFrwd()//Wave Field PhaseShift { int PhaseShift(); int Frqcy2Time(); if( PhaseShift()!=1 ){printf("PhaseShift is error\n"); exit(0);} if( Frqcy2Time()!=1 ){printf("Frqcy2Time is error\n"); exit(0);} return(1); } ////////////////////////////////////////////////// int PhaseShift()//相移 { //1.Prepprocedure预处理 FILE *fp_Wfldr,*fp_Wfldi,*fp_Wfld0r,*fp_Wfld0i,*fp_Vlcty,*fp_Absb; float kz[2]; int Ix,Iz,Iw,Nw=Nt,Ikx,Nkx=Nx; long Mgrtn; float Vlcty[Nz]; float Wfld_r,Wfld_i; float Absb[Nx],Wfld0r[Nx],Wfld0i[Nx],Wfldr[Nx],Wfldi[Nx]; float Kxmax,Dkx,Wmax,Dw; Wmax=PI/Dt; Dw=Wmax/(float)Nt; Kxmax=PI/Dx; Dkx=Kxmax/(float)Nx; //2.Read in Velocity and Absorbing Boundary Data速度与削波数据读入 if((fp_Vlcty=fopen("Vlcty.dat","rb"))==NULL) printf("Cannot open Rflct.dat"); for(Iz=0;Iz { fread(&Vlcty[Iz],sizeof(Vlcty[Iz]),1,fp_Vlcty); } fclose(fp_Vlcty); if((fp_Absb=fopen("Absb.dat","rb"))==NULL)printf("Cannot open Absb.dat"); for(Ix=0;Ix { fread(&Absb[Ix],sizeof(Absb[Ix]),1,fp_Absb); } fclose(fp_Absb); //remove("Absb.dat"); //3.Open Initial Wave Field File and Current Wave Field File using In Wave Field Extrapolating波场文件打开 if((fp_Wfld0r=fopen("Wfld0r.dat","rb"))==NULL) printf("Cannot open Wfld0r.dat"); if((fp_Wfld0i=fopen("Wfld0i.dat","rb"))==NULL) printf("Cannot open Wfld0i.dat"); if((fp_Wfldr=fopen("Wfldr.dat","wb"))==NULL) printf("Cannot open Wfldr.dat"); if((fp_Wfldi=fopen("Wfldi.dat","wb"))==NULL) printf("Cannot open Wfldi.dat"); //4.Wave Fied Extrapolate With Every Frequency每个频率的波场延拓 for(Iw=0;Iw { printf("PhaseShift:Iw=%d\n",Iw); //4.1 Initializing Wavefield of Every Frequency: All get 0初始化当前波场 for(Ix=0;Ix { Wfldr[Ix]=0.; Wfldi[Ix]=0.; } //4.2 Wave Fied Extrapolate From Z=Zmax to Z=Zmin波场从Iz=Nz-1最深处开始，延拓到Iz=1测线深度 for(Iz=Nz-1;Iz>0;Iz--) { //4.2.1 Forming New Wavefield for Extrapolatingon Space Domain形成新波长 for(Ix=0;Ix { Mgrtn=(Ix*Nz+Iz)*(Nw/2+1)+Iw; //Take out Initial Wave Field Data With every Depth取出当前深度的初始波场 fseek(fp_Wfld0r,sizeof(Wfld0r[Ix])*Mgrtn,0); fread(&Wfld0r[Ix],sizeof(Wfld0r[Ix]),1,fp_Wfld0r); fseek(fp_Wfld0i,sizeof(Wfld0i[Ix])*Mgrtn,0); fread(&Wfld0i[Ix],sizeof(Wfld0i[Ix]),1,fp_Wfld0i); //New Wave Fied Equal to the Initial on Plus Current one新波场=初始波场+从下面延拓到此处的波场 Wfldr[Ix]=Wfldr[Ix]+Wfld0r[Ix]; Wfldi[Ix]=Wfldi[Ix]+Wfld0i[Ix]; //Boundary absorbing of New Wave Fied边界削波：新波场=新波场×削波因子 Wfldr[Ix]=Wfldr[Ix]*Absb[Ix]; Wfldi[Ix]=Wfldi[Ix]*Absb[Ix]; } //4.2.2 FFT of New Wave field From Space to WaveNumber新波场FFT到波数域 if(kkfft(Wfldr,Wfldi,Nx,0)!=1){printf("FFT is error");exit(0);} //4.2.3 Wavefield Extrapolating on Frequency-WaveNumber Domain频率-波数域波场在从Iz+1延拓到Iz for(Ikx=0;Ikx { //4.2.3.1 Computing Phaseshift data计算相移数据expikzdz（实部、虚部） if(exp_ikzDz(kz,Ikx,(float)(Vlcty[Iz]/2.),Iw,Dw,Dkx)!=1){printf("exp_ikzDz is error");exit(0);} //4.2.3.2 WaveField multiply Phaseshift data波场延拓：波场=波场×相移数据 Wfld_r=Wfldr[Ikx]*kz[0]-Wfldi[Ikx]*kz[1]; Wfld_i=Wfldi[Ikx]*kz[0]+Wfldr[Ikx]*kz[1]; Wfldr[Ikx]=Wfld_r; Wfldi[Ikx]=Wfld_i; if(Ikx!=0&&Ikx!=Nkx/2) { Wfld_r=kz[0]*Wfldr[Nkx-Ikx]-kz[1]*Wfldi[Nkx-Ikx]; Wfld_i=kz[1]*Wfldr[Nkx-Ikx]+kz[0]*Wfldi[Nkx-Ikx]; Wfldr[Nkx-Ikx]=Wfld_r; Wfldi[Nkx-Ikx]=Wfld_i; } } //4.2.4 IFFT of New Wave field From WaveNumber to Space Domain波场反FFT到空间域 if(kkfft(Wfldr,Wfldi,Nkx,1)!=1){printf("FFT is error");exit(0);} } //4.3 Storing Wave field Data on Survey Line:Z=Zmin存储延拓到了测线的波场 for(Ix=0;Ix { fwrite(&Wfldr[Ix],sizeof(Wfldr[Ix]),1,fp_Wfldr); fwrite(&Wfldi[Ix],sizeof(Wfldi[Ix]),1,fp_Wfldi); }; } //5 Close or Remove Data File关闭文件，删除中间文件 fclose(fp_Wfld0r); fclose(fp_Wfld0i); fclose(fp_Wfldr); fclose(fp_Wfldi); remove("Absb.dat"); return(1); } ////////////////////////////////////////////////// int exp_ikzDz(float eikzdz[],int Ix, float Vc,int Iw,float Dw,float Dkx) { float kz=0.; eikzdz[0]=0.; eikzdz[1]=0.; kz=sqrt(pow(Iw*Dw/Vc,2)-pow(Ix*Dkx,2)); if(kz>0.) { eikzdz[0]=(float)cos(kz*Dz); eikzdz[1]=(float)-sin(kz*Dz); } return(1); } ////////////////////////////////////////////////// int Frqcy2Time() { FILE *fp_Wfldr,*fp_Wfldi,*fp_Record; int Ix,It,Iw,Nw=Nt; float Wfldtr[Nt],Wfldti[Nt]; long AddFrmStrt; if((fp_Wfldr=fopen("Wfldr.dat","rb"))==NULL) printf("Cann't open file Wfldr.dat\n"); if((fp_Wfldi=fopen("Wfldi.dat","rb"))==NULL) printf("Cann't open file Wfldi.dat\n"); if((fp_Record=fopen("Record.dat","wb"))==NULL) printf("Cann't open file Record.dat\n"); for(Ix=0;Ix { //printf("Frqcy2Time:Ix=%d\n",Ix); for(Iw=0;Iw { AddFrmStrt=Iw*Nx+Ix; fseek(fp_Wfldr,sizeof(Wfldtr[Iw])*AddFrmStrt,0); fseek(fp_Wfldi,sizeof(Wfldti[Iw])*AddFrmStrt,0); fread(&Wfldtr[Iw],sizeof(Wfldtr[Iw]),1,fp_Wfldr); fread(&Wfldti[Iw],sizeof(Wfldti[Iw]),1,fp_Wfldi); if(Iw!=0&&Iw!=Nw/2) { Wfldtr[Nw-Iw]=Wfldtr[Iw]; Wfldti[Nw-Iw]=-Wfldti[Iw]; } } if(kkfft(Wfldtr,Wfldti,Nw,1)!=1){printf("FFT is error\n");exit(0);} for(It=0;It { fwrite(&Wfldtr[It],sizeof(Wfldtr[It]),1,fp_Record); } } fclose(fp_Wfldr); fclose(fp_Wfldi); remove("Wfldr.dat"); remove("Wfldi.dat"); fclose(fp_Record); return(1); } // Fourier Transform: FFT(i=0) or IFFT(l=1) int kkfft(float pr[], float pi[], int n, int l) { int it,m,is,i,j,nv,l0,il=0; float p,q,s,vr,vi,poddr,poddi; float fr[4096],fi[4096]; int k=0; long Ln=0; for(k=0;n-Ln>0;k++) { Ln=(long)pow(2,k); } k=k-1; for (it=0; it<=n-1; it++) { m = it; is = 0; for(i=0; i<=k-1; i++) { j = m/2; is = 2*is+(m-2*j); m = j; } fr[it] = pr[is]; fi[it] = pi[is]; } pr[0] = 1.0; pi[0] = 0.0; p = 6.283185306/(1.0*n); pr[1] = (float) cos(p); pi[1] = -(float)sin(p); if (l!=0) pi[1]=-pi[1]; for (i=2; i<=n-1; i++) { p = pr[i-1]*pr[1]; q = pi[i-1]*pi[1]; s = (pr[i-1]+pi[i-1])*(pr[1]+pi[1]); pr[i] = p-q; pi[i] = s-p-q; } for (it=0; it<=n-2; it=it+2) { vr = fr[it]; vi = fi[it]; fr[it] = vr+fr[it+1]; fi[it] = vi+fi[it+1]; fr[it+1] = vr-fr[it+1]; fi[it+1] = vi-fi[it+1]; } m = n/2; nv = 2; for (l0=k-2; l0>=0; l0--) { m = m/2; nv = 2*nv; for(it=0; it<=(m-1)*nv; it=it+nv) for (j=0; j<=(nv/2)-1; j++) { p = pr[m*j]*fr[it+j+nv/2]; q = pi[m*j]*fi[it+j+nv/2]; s = pr[m*j]+pi[m*j]; s = s*(fr[it+j+nv/2]+fi[it+j+nv/2]); poddr = p-q; poddi = s-p-q; fr[it+j+nv/2] = fr[it+j]-poddr; fi[it+j+nv/2] = fi[it+j]-poddi; fr[it+j] = fr[it+j]+poddr; fi[it+j] = fi[it+j]+poddi; } } if(l!=0) { for(i=0; i<=n-1; i++) { fr[i] = fr[i]/(1.0*n); fi[i] = fi[i]/(1.0*n); } } if(il!=0) { for(i=0; i<=n-1; i++) { pr[i] = sqrt(fr[i]*fr[i]+fi[i]*fi[i]); if(fabs(fr[i])<0.000001*fabs(fi[i])) { if ((fi[i]*fr[i])>0) pi[i] = 90.0; else pi[i] = -90.0; } else pi[i] = atan(fi[i]/fr[i])*360.0/6.283185306; } } for(i=0;i { pr[i]=fr[i]; pi[i]=fi[i]; } return ( 1 ); } 上海杉达学院 商务流程综合实训总结 单位名称：微科电子有限公司（加盖公章）姓名：陈恩娜 学院：胜祥商学院 专业：国际经济与贸易 班级： f130219 时间： 2016.11.14 工作总结历时10天的跨专业实训圆满落幕，作为国贸专业的我们参加了企业经营的模拟实训，通过这几天的实训，我也算是了解到了企业经营的基本流程和一般模式。前三天，我们基本在忙碌企业创立的事情。我们首先要做的是公司注册登记。公司注册流程共包括名称预先登记、设立登记申请书、准备申请材料、银行开户入资、验资、前置审批、报送申请材料、工商审批发照、刻制公章、开设银行帐户、办理各项登记、股东入资证明、企业机构代码、企业税务登记、企业劳动备案、社会保险登记、企业统计备案、特殊行业备案18类。我们实训时规定了公司类型为制造型企业，开始注册资金为500万元。在实训中，我们注册的公司为微科电子有限公司，地址位于上海市陆家嘴，股东为3人。申请表提交上去后，我们得到了审批，就开始了下面的企业经营规划。我们企业经营过程共分为九个相互联系又相互独立的部门。这十个部门分别为总经办、人力资源部门、物流部门、质检部门、行政部门、营销部门、生产部门、财务部门、采购部门。总经办为CEO代理，也就是我们常说的总经理，其他分别为营销总监、生产总监、物流总监、采购总监、财务总监、采购经理、人事经理、行政主管、质检经理。这九个职位分别为十一个同学完成，我作为采购部经理也参与其中。采购部的岗位职责可归纳为：依照公司生产需要及物资采购计划，全面负责公司的物料采购和供应工作；对初次进行合作的供应商进行调查了解，根据实际掌握的资料与信息做出相应的评价判定。对于符合公司要求的，方可与之开展业务往来与合作；采购工作的开展应当做到精打细算，尽力降低公司的采购成本，避免出现浪费公司资金的情况；负责对所有与公司有业务往来的供应商或供货企业进行定期的资质评价并给出明确的评价等级。针对不同的评价等级给出相应的处理意见；完成公司领导布置的其它各项工作。各岗位要各司其职，填写经营流程表，有序的完成一年的经营。 中南大学 计算机仿真与建模 实验报告 题目：理发店的服务过程仿真 姓名：XXXX 班级：计科XXXX班 学号：0909XXXX 日期：2013XXXX 理发店的服务过程仿真 1 实验案例 (2) 1.1 案例：理发店系统研究 (2) 1.1.1 问题分析 (3) 1.1.2 模型假设 (3) 1.1.3 变量说明 (3) 1.1.4 模型建立 (3) 1.1.5 系统模拟 (4) 1.1.6 计算机模拟算法设计 (5) 1.1.7 计算机模拟程序 (6) 1实验案例 1.1 案例：理发店模拟 一个理发店有两位服务员A和B顾客随机地到达该理发店,每分钟有一个顾客到达和没有顾客到达的概率均是1/2 , 其中60%的顾客理发仅用5分钟,另外40%的顾客用8分钟. 试对前10分钟的情况进行仿真。 （“排队论”，“系统模拟”，“离散系统模拟”，“事件调度法”） 1.1.1 问题分析 理发店系统包含诸多随机因素，为了对其进行评判就是要研究其运行效率， 从理发店自身利益来说，要看服务员工作负荷是否合理，是否需要增加员工等考 虑。从顾客角度讲，还要看顾客的等待时间，顾客的等待队长，如等待时间过长 或者等待的人过多，则顾客会离开。理发店系统是一个典型的排队系统，可以用 排队论有关知识来研究。 1.1.2 模型假设 1． 60％的顾客只需剪发，40％的顾客既要剪发，又要洗发； 2． 每个服务员剪发需要的时间均为5分钟，既剪发又洗发则花8分钟； 3． 顾客的到达间隔时间服从指数分布； 4． 服务中服务员不休息。 1.1.3 变量说明 u ：剪发时间（单位：分钟），u=5m ； v: 既剪发又理发花的时间（单位：分钟），v=8m ； T ： 顾客到达的间隔时间，是随机变量，服从参数为λ的指数分布，（单位： 分钟） T 0：顾客到达的平均间隔时间（单位：秒），T 0＝λ 1； 1.1.4 模型建立 由于该系统包含诸多随机因素，很难给出解析的结果，因此可以借助计算机 模拟对该系统进行模拟。 考虑一般理发店的工作模式，一般是上午9：00开始营业，晚上10：00左 右结束，且一般是连续工作的，因此一般营业时间为13小时左右。 这里以每天运行12小时为例，进行模拟。 这里假定顾客到达的平均间隔时间T 0服从均值3分钟的指数分布， 则有 3小时到达人数约为603 603=?人， 6小时到达人数约为1203 606=?人， 10小时到达人数约为2003 6010=?人， 这里模拟顾客到达数为60人的情况。 （如何选择模拟的总人数或模拟总时间） 本科生实验报告 实验课程数值模型模拟 学院名称地球物理学院 专业名称勘测技术与工程 学生 学生学号 指导教师熊高君 实验地点5417 实验成绩 2015年5月 理工大学 《地震数值模拟》实验报告 实验报告 一、实验题目： 地震记录数值模拟的褶积模型法 二、实验目的： 掌握褶积模型基本理论、实现方法与程序编制，由褶积模型初步分析地震信号的分辨率问题。 三、原理公式 1、褶积原理 地震勘探的震源往往是带宽很宽的脉冲，在地下传播、反射、绕射到测线，传播经过中高频衰减，能量被吸收。吸收过程可以看成滤波的过程，滤波可以用褶积完成。在滤波中，反射系数与震源强弱关联，吸收作用与子波关联。最简单的地震记录数值模拟，可以看成反射系数与子波的褶积。通常，反射系数是脉冲，子波取雷克子波。 （1）雷克子波： wave(t)=cos(2ft)* （2）反射系数： （3）褶积公式： 数值模拟地震记录trace(t): trace(t) =rflct(t)*wave(t)； 反射系数的参数由 z 变成了 t，怎么实现？在简单水平层介质，分垂直和非垂直入射两种实现，分别如图 1 和图 2 所示。 图1 图2 1）垂直入射： t=2h/v； 2）非垂直入射： t= 2、褶积方法 （1）离散化（数值化） 计算机数值模拟要求首先必须针对连续信号离散化处理。反射系数在空间模型中存在，不同深度反射系数不同，是深度的函数。子波是在时间记录上一延续定时间的信号，是时间的概念。在离散化时，通过深度采样完成反射系数的离散化，通过时间采样完成子波的离散化。如果记录是 Trace(t)，则记录是时间的函数，以时间采样离散化。时间采样间距以Δt 表示，深度采样间距以Δz 表示。在做多道的数值模拟时，还有横向Δx 的概念，横向采样间隔以Δx 表示。 离散化的实现：t=It×Δt；x=Ix×Δx；z=Iz×Δz； 或：It=t/Δt; Ix=x/Δx; Iz=z/Δz （2）离散序列的褶积 trace(It)= 四、实验容 1、垂直入射地震记录数值模拟的褶积模型； 2、非垂直入射地震记录数值模拟的褶积模型。 五、方法路线 中国地质大学（武汉）地空学院 地震实验报告 姓名：沈 班级：班 学号： 时间： 2015年05月 指导老师：张 一、实验目的 实验一： 1、浅层地震装备的基本组成； 2、认识GEODE96浅层地震仪的主要结构，并学会该类仪器的操作方法； 3、地震波认识。 实验二： 1、掌握浅层地震数据采集方法及注意事项 二、仪器介绍 1、仪器简介 全套美国GEOMETRICS公司生产的Geode96浅层地震仪（相当于四套独立的24道浅层地震仪）该仪器能满足折、反射地震勘探、井间勘探、面波调查等地震监测需要,应用Crystal公司的A/D转换器和高速过采样技术达到了24位地震仪的精度。频带从1.75Hz到20,000Hz,使得采样间隔可以从20毫秒到16微秒。采样到的数据叠加到32位的叠加器中，然后传回到主机的硬盘或其它介质上。内置预触发器，每道有16K的内存。用硬件相关器对震源信号进行实时相关运算。Geode包装坚固、防水、防震,有提手,重4.1公斤,用12V的外接电池可以连续工作10个小时。(如下图) 2、主要操作功能键及快捷键 注释： 1锁定与解锁；2清除界面；4检测噪声；7保存 3、操作步骤及注意事项 1、每个GEODE用数传线按规定串联，通过数传盒与笔记本电脑的USB口连 接。 2、每个GEODE接上12V电源。 3、开关接到与笔记本相连的第一个GEODE上。 4、传盒上的开关置于POWER UP处。 5、采集控制程序，并按工作需要设置好各项参数，然后进行正常数据采集工作。 6、出采集控制程序之前，应将数传盒上的开关置于POWER DOWN处。 7、卸下各连接线并清理整齐。 8、注意的是：在正常工作过程中，任何时候移动数传线与GEODE的连接头时，必须退出采集控制程序。另外Y型头上有红色标记的与GEODE的前12道相连接。而且采集控制软件运行的语言环境必须是英语(美国)。 实验课程数值模型模拟 学院名称地球物理学院 专业名称勘测技术与工程 学生姓名 学生学号 指导教师熊高君 实验地点 5417 实验成绩 2015年5月 成都理工大学 地震数值模拟》实验报告 实验报告 实验题目： 地震记录数值模拟的褶积模型法 二、实验目的： 掌握褶积模型基本理论、实现方法与程序编制，由褶积模型初步分析地震信号的分辨率问题。 三、原理公式 1、褶积原理 地震勘探的震源往往是带宽很宽的脉冲，在地下传播、反射、绕射到测线，传播经过中高频衰减，能量被吸收。吸收过程可以看成滤波的过程，滤波可以用褶积完成。在滤波中，反射系数与震源强弱关联，吸收作用与子波关联。最简单的地震记录数值模拟，可以看成反射系数与子波的褶积。通常，反射系数是脉冲，子波取雷克子波。 (1)雷克子波： wave(t)=cos(2 ft)* (2)反射系数： (3)褶积公式： 数值模拟地震记录trace(t): trace(t) =rflct(t)*wave(t)； 反射系数的参数由 z 变成了 t，怎么实现？在简单水平层介质，分垂直和非垂直 图1 图 2 1） 垂直入射： t=2h/v ； 2）非垂直入射： t= 2 、褶积方法 （1） 离散化（数值化） 计算机数值模拟要求首先必须针对连续信号离散化处理。反射系数在空间模 型中存在，不同深度反射系数不同，是深度的函数。子波是在时间记录上一延续 定时间的信号，是时间的概念。在离散化时，通过深度采样完成反射系数的离散 化，通过时间采样完成子波的离散化。如果记录是 Trace （t ），则记录是时间的 函数，以时间采样离散化。时间采样间距以Δt 表示，深度采样间距以Δz 表示。 在做多道的数值模拟时，还有横向Δx 的概念，横向采样间隔以Δx 表示。 离散化的实现：t=It×Δt；x=Ix×Δx；z=Iz×Δz； 或：I t=t/Δt; Ix=x/Δx; Iz=z/Δz （2） 离散序列的褶积 trace(It)= 宁波职业技术学院 外贸实务 II - 提高 实训报告 指导老师：江彬 班级：国贸3141 学生姓名：肖思洁 学号： 1426263133 日期： 2016-04-25 课程名称：外贸实务II-提高 1 实训目的及要求 1．1 1．2 2 实训内容及步骤（包含简要的实训步骤流程） 2．1 本人所扮演的角色 2．2 贸易资料及实训步骤 2．2．1 2．2．2 2．2．3 2．2．4 3 实训结果（包括实训项目的完成情况，代表性邮件，单据，程序或图表、结论陈述、核算表数据记录及分析等） 3．1 3．1．1 3．1．2 3．1．3 3．2 4 实训中遇到的问题及其解决方法 5 实训总结（包括心得体会、对SimTrade系统的评价、对自己实训效果的评价如实训收获不足及实训改进意见等） 6 实训评价 外贸实训报告 外贸实务实训体会总结，这次是项目过程的完成让我感觉很深刻。短短的32学时已经结束，静下心来回想这次操作模拟学习真是感受颇深。我们知道动手操作是大学教育中一个极为重要的实践性环节，通过实习，可以使我们在实践中接触与本专业相关的一些实际工作，培养和锻炼我们综合运用所学的基础理论、基本技能和专业知识，去独立分析和解决实际问题的能力，把理论和实践结合起来，提高我们的实际动手能力，为将来我们毕业后走上工作岗位打下一定的基础。通过这段时间的学习，从无知到认知，到深入了解，渐渐地我喜欢上这个专业，让我深刻的体会到学习的过程是最美的，在整个实习过程中，我每天都有很多的新的体会，新的想法。 回顾我的实习生活，感触是很深的，收获也是丰硕的。可以模拟出口商与非洲、中东等地方的外国商人做生意，他们在我公司下订单，我们再把订单下到厂里，从中赚取差额，或者作为进口商与出口商讨价还价，去除运费保险费等来赚取利益。当然对工厂的基本流程也有一定的了解。这次的实习经历我对外贸这个专业有了更加理性的认识和更深刻的体会。在这次是学习中，我学到了很多过去两年没有体会到的东西，这不仅仅只是上课模拟，也是一次对真实工作流程经历。 对实践的看法： 在操作过程中，根据本门课程的内容、特点，通过走出去、请进来等方式，精心组织方案。通过听、看、做使一些看起来繁杂的专业知识很快被我们理解和掌握。只有将理论联系实际，教学与实际相结合，才是培养我们能力的一种有效形式。 出口商+进口商+工厂，一共写了78封邮件。一共发布了8封广告和信息。 查了B2B里的多数产品信息。 银行汇率：欧元大多汇率为8.6402，美元大多为6.1463等。保险费：一切险（ALL RISKS）为0.8%，战争险（W AR RISKS）为0.08% 保险费计算方式为 (1)按CIF进口时：保险金额＝CIF货价×1.1 (2)按CFR进口时：保险金额＝CFR货价×1.1 / (1 - 1.1 ×r)，其中r为保险费率，请在"淘金网"的"保险费"页面查找，将所投险别的保险费率相加即可。 (3)按FOB进口时：保险金额＝(FOB货价+ 海运费)×1.1 / (1 - 1.1 ×r)，其中FOB 货价就是合同金额，海运费请在装船通知中查找，由出口商根据配舱通知填写，如果出口商填写错误，请其查看配舱通知。 实训的基本流程： 第一周完成了进口商，出口商，工厂，进口、出口地银行的基本资料。然后熟悉了下系统的基本轮廓，如B2B里面可以查询写什么，市场，海关等在哪个位置。根据老师的知道，试着去发广告与写邮件。 第二周确定角色，开始寻找有利信息，搜索信息，同业务伙伴建合作关系。 我先进行成本、费用、利润等的核算，若有盈利则进一步磋商合作，若亏损就跟对方进行讨价还价。过程为询盘——发盘——还盘——接受。 第三周之后进入交易准备阶段——交易磋商阶段——签订合同(T/T+FOB)与接收信用证（L/C+CIF）——履行合同阶段。 签订合同之后进行合同履行阶段。 首先作为出口商，与进口商进行磋商商定后确定的价格，之后跟工厂进行合作并进一步签订SALES CONFIRMATON。等工厂交货物发过来后，与工厂的业务就能完成。 计算机模拟手工实验 学生实验报告 学院：商学院 课程名称：计算机模拟手工实验 专业班级： 姓名： 学号： 学生实验报告 第一部分：实验概况与内容 一、实验的目的及要求 1、实验目的 本实验以模拟企业的实际会计工作为基础，按照企业会计制度和企业会计准则的要求，进行操作训练，有目的地检验和复习所学的会计理论、方法、技能和技巧通过实际的操作，使我们能够比较系统、全面地掌握工业企业会计核算的基本程序和具体方法，加强我们对会计基本理论的理解和对会计基本技能的掌握，把枯燥、抽象的书本知识转化为实际、具体的操作，使我们能够形象地掌握各种业务的处理及记账凭证的填写方法，掌握账簿的处理及登记方法，掌握成本核算方法，掌握各种报表的编制方法，掌握会计资料的整理归档方法，同时，我们可以体验在不同岗位进行不同操作，使之在实验中，培养职业道德和职业判断能力，提高职业工作能力，为我们今后从事会计实务工作打下扎实的基础。 2、实验要求 ①熟悉会计工作的基本流程，工作内容以及工作规范等基础知识； ②能够熟练的进行对实验企业所发生各经济业务的会计核算和账务处理； ③掌握实验企业建立账户、填制凭证，登记账簿、编制报表等会计实务操作流程。 二、实验内容 1、企业基本情况 津阳市永安公司是批零兼营的以零售为主的商品流通企业，主要经营五金、百货、家电等商品，分设一部四柜组，其中一部为批发部，四柜组为小百货组、五金家电组、鞋帽组、针织服装组。 开户行及账号：中国工商银行贵溪分理处 5189958。 地址：津阳市盛兴路160号。 经营规模：一般纳税人，适用增值税税率为17%。 纳税人识别号：235678902283156。 所得税税率：25% 2、内部主要财务会计制度 （1）批发商品流转业务核算的有关规定和要求： ①库存商品采用数量进价金额核算法，按商品品名开设明细账进行数量进价金额核算。 ②“商品销售收入”、“商品销售成本”账户按批发设置明细账，以便结转成本。 ③商品销售使用增值税专用发票，税率为17%。 ④商品销售成本本月末采用先进先出法，在“库存商品——批发”账户中倒算并结转成本。平时只填制出库单。 （2）零售商品流转业务核算的有关规定和要求： ①库存商品采用售价金额核算法，“库存商品”账户按零售分设小百货组、五金家电组、鞋帽组、针织服装组分户进行明细核算。 ②“商品销售收入”、“商品销售成本”账户按批发设置明细账，以便结转成本。 ③商品销售使用增值税专用发票，税率为17%。 ④商品销售实行“价税合一”、平时“商品销售收入”反映含税（增值税、下同）销售额，月末按下列公式调整为不含税销售额，以此计算冲销已销商品收入所含的增值税（进项税额）。不含税销售额=含税销售额/ （1+增值税税率） ⑤商品销售成本按含税销售额随销随转办法，注销书屋负责人的经济责任。 ⑥“商品进销差价”账户反映含税售价与不含说进价之差的数额，并按前述四柜组分别核算。 ⑦月末，按分类（柜组）差价率计算法计算并分摊已销商品实现的进销差价。（3）本公司采用的是非定额的备用金制度 3.实验过程 地震勘探实验报告记录 ————————————————————————————————作者：————————————————————————————————日期： 中国地质大学（武汉）地空学院 地震实验报告 姓名：沈 班级：班 学号： 时间： 2015年05月 指导老师：张 一、实验目的 实验一： 1、浅层地震装备的基本组成； 2、认识GEODE96浅层地震仪的主要结构，并学会该类仪器的操作方法； 3、地震波认识。 实验二： 1、掌握浅层地震数据采集方法及注意事项 二、仪器介绍 1、仪器简介 全套美国GEOMETRICS公司生产的Geode96浅层地震仪（相当于四套独立的24道浅层地震仪）该仪器能满足折、反射地震勘探、井间勘探、面波调查等地震监测需要,应用Crystal公司的A/D转换器和高速过采样技术达到了24位地震仪的精度。频带从1.75Hz到20,000Hz,使得采样间隔可以从20毫秒到16微秒。采样到的数据叠加到32位的叠加器中，然后传回到主机的硬盘或其它介质上。内置预触发器，每道有16K的内存。用硬件相关器对震源信号进行实时相关运算。Geode包装坚固、防水、防震,有提手,重4.1公斤,用12V的外接电池可以连续工作10个小时。(如下图) 2、主要操作功能键及快捷键 注释： 1锁定与解锁；2清除界面；4检测噪声；7保存 3、操作步骤及注意事项 1、每个GEODE用数传线按规定串联，通过数传盒与笔记本电脑的USB 口连接。 2、每个GEODE接上12V电源。 3、开关接到与笔记本相连的第一个GEODE上。 4、传盒上的开关置于POWER UP处。 5、采集控制程序，并按工作需要设置好各项参数，然后进行正常数据采集工作。 6、出采集控制程序之前，应将数传盒上的开关置于POWER DOWN处。 7、卸下各连接线并清理整齐。 8、注意的是：在正常工作过程中，任何时候移动数传线与GEODE的连接头时，必须退出采集控制程序。另外Y型头上有红色标记的与GEODE的前12道相连接。而且采集控制软件运行的语言环境必须是英语(美国)。 国贸专业生产实习报告 随着国际贸易的日益完善，以及中国在国际贸易的地位的不断上升，我们作为未来社会的国贸人员，为了加强社会竞争力，应培养较强的国贸工作的操作能力。于是，在结束了大三的课程后,学校给了我们一个很好的实习锻炼机会，让我们模拟国际贸易实务操作，从而从中掌握国际贸易流程。 一、实习目的 ①熟悉外贸实务的具体操作流程; ②了解、巩固与深化已经学过的理论和方法; ③增强对外贸实务的感性认识; ④提高发现问题、分析问题以及解决问题的能力。 二、实习方法： 通过进入SimTrade模拟平台，进行上机模拟操作 Simtrade外贸实习平台是一个十分成功的国际贸易模拟软件，它在很大程度上解决了学生实习难的问题。学生在网上进行国际货物买卖实务的具体操作，能很快掌握进出口的成本核算、询盘、发盘与还盘等各种基本技巧；熟悉国际贸易的物流、资金流与业务流的运作方式；切身体会到国际贸易中不同当事人面临的具体工作与他们之间的互动关系；学会外贸公司利用各种方式控制成本以达到利润最大化的思路；认识供求平衡、竞争等宏观经济现象，并且能够合理地加以利用。老师通过在网站发布新闻等行为对国际贸易环境实施宏观调控，使学生在实习中充分发挥主观能动性，真正理解并吸收课堂中所学到的知识，为将来走上工作岗位打下良好基础。 三、实习遇到的问题 1、预算错误 这是开始接触Simtrade时所最容易忽略的问题。虽然老师曾多次提醒，做贸易前一定要计算好了一切费用，选好贸易术语，最后再签定合同。但我们经常做出口商的在还没有调查进口商所在地市场的情况下就先去工厂进货了。如果工厂角色也没有做好预算，草草就签订了合同，那么可能出口商和工厂都赚不到钱。在最后交易完成后，我们经常大叫“啊，这个运费怎么比我的货物数量还多啊？”“这个保险费怎么这么贵，我要赔钱了！” 学生学号0120801080128 实验课成绩 学生实验报告书 实验课程名称材料成型数值模拟设计实验 开课学院材料学院 指导教师姓名朱春东、钱东升 学生姓名王丹丹 学生专业班级成型0801 2011-- 2012学年第一学期 实验教学管理基本规范 实验是培养学生动手能力、分析解决问题能力的重要环节；实验报告是反映实验教学水平 与质量的重要依据。为加强实验过程管理，改革实验成绩考核方法，改善实验教学效果，提高 学生质量，特制定实验教学管理基本规范。 1、本规范适用于理工科类专业实验课程，文、经、管、计算机类实验课程可根据具体情况参 照执行或暂不执行。 2、每门实验课程一般会包括许多实验项目，除非常简单的验证演示性实验项目可以不写实验 报告外，其他实验项目均应按本格式完成实验报告。 3、实验报告应由实验预习、实验过程、结果分析三大部分组成。每部分均在实验成绩中占一 定比例。各部分成绩的观测点、考核目标、所占比例可参考附表执行。各专业也可以根据具体情况，调整考核内容和评分标准。 4、学生必须在完成实验预习内容的前提下进行实验。教师要在实验过程中抽查学生预习情况， 在学生离开实验室前，检查学生实验操作和记录情况，并在实验报告第二部分教师签字栏签名，以确保实验记录的真实性。 5、教师应及时评阅学生的实验报告并给出各实验项目成绩，完整保存实验报告。在完成所有 实验项目后，教师应按学生姓名将批改好的各实验项目实验报告装订成册，构成该实验课程总报告，按班级交课程承担单位（实验中心或实验室）保管存档。 6、实验课程成绩按其类型采取百分制或优、良、中、及格和不及格五级评定。 附表：实验考核参考内容及标准 观测点考核目标成绩组成 实验预习1．预习报告 2．提问 3．对于设计型实验，着重考查设计方案的 科学性、可行性和创新性 对实验目的和基本原理 的认识程度，对实验方 案的设计能力 20% 实验过程1．是否按时参加实验 2．对实验过程的熟悉程度 3．对基本操作的规范程度 4．对突发事件的应急处理能力 5．实验原始记录的完整程度 6．同学之间的团结协作精神 着重考查学生的实验态 度、基本操作技能；严 谨的治学态度、团结协 作精神 30% 结果分析1．所分析结果是否用原始记录数据 2．计算结果是否正确 3．实验结果分析是否合理 4．对于综合实验，各项内容之间是否有分 析、比较与判断等 考查学生对实验数据处 理和现象分析的能力； 对专业知识的综合应用 能力；事实求实的精神 50% 《浅层折射波勘探》实验报告 《浅层折射波勘探》实验成绩评定表班级姓名学号 一、实验名称：浅层折射波勘探 二、实验目的 加深对地震勘探基本概念的理解，巩固已学的理论知识，了解数字地震仪的使用和仪器工作参数的选择；了解地震勘探人工震源激发，检波器的安置条件；地震折射波法野外资料的采集技术及方法，并进行资料的整理与解释；了解地震勘探野外工作施工的过程以及组织管理工作。 三、实验原理 1、折射波法基本原理 以水平界面的两层介质进行简要的说明，假设地下深度为h ，有一个水平的速度分界面R ，上、下两层的速度分别为V 1和V 2，且V 2>V 1。 如图1所示。从激发点O 至地面某一接收点D 的距离为X ，折射波旅行的路程为OK 、KE 、ED 之和，则它的旅行时t 为： 图1 水平两层介质折射波时距曲线 1 21V ED V KE V OK t ++= 式1 为了简便起见，先作如下证明：从O ，D 两点分别作界面R 的垂线，则OA ＝DG ＝h ，再自A 、G 分别作OK ，ED 的垂线，几何上不难证明∠BAK ＝∠EGF ＝i ，因 已知2 1 sin V V i = ，所以： 2 1 V V EG EF AK BK == 式2 即 21V AK V BK = 和 2 1V EG V EF = 式3 上式说明，波以速度V 1旅行BK （或EF ）路程与以速度V 2旅行AK （或EC ）路程所需的时间是相等的。将式3的关系和式1作等效置换，并经变换后可得： 2 121222122cos 2V V V V h V x V i h V x t -+=+= 式4 这就是水平两层介质的折射波时距曲线方程。它表示时距曲线是一条直线，若令x ＝0，则可得时距曲线的截距时间t 0（时距曲线延长与t 轴相交处的时间值） 2 12122102cos 2V V V V h V i h t -== 式5 式5表示出界面深度h 和截距时间t0之间的关系，当已知V 1和V 2时，可以求出界面的深度h 。 2、折射波分层解释的t 0法 折射波t 0解释法是常用的地震折射波解释方法，它是针相遇时距曲线观测系统采集发展起来的解释方法。 t 0法解释的主要原理与方法如下： t 0法又称为t 0差数时距曲线法，是解释折射波相遇时距曲线最常用的方法之一。当折射界面的曲率半径比其埋深大得很多的情况下，t 0法通常能取得很好的效果，且具有简便快速的优点。 如图2所示，设有折射波相遇的时距曲线S 1和S 2，两者的激发点分别是O 1 和O 2， 本科生实验报告 实验课程基于 Vista 系统的地震资料处理学院名称地球物理学院 专业名称勘查技术与工程（石油物探）学生姓名 学生学号 指导教师唐湘蓉 实验地点5417 实验成绩 2015年3月- 2015年5月 基于 Vista 系统的地震资料处理 一、实验目的及要求 1）认知熟悉地震资料处理软件系统--vista软件的基本功能，了解其并熟练掌握vista软件运行的基本操作； 2）了解并掌握地震数据处理的基本流程，掌握地震数据处理的流程和基本方法，选择合适的处理参数以提高地震数据处理的精度； 3）对比地震资料处理与解释的理论与实际资料处理的结果，深入理解理论，并在理论指导下提高处理解释的水平、提高资料处理的质量； 4）提高综合分析问题的能力与编写实验报告或生产报告的能力。 二、实验内容 总流程 图1 总流程图 1）加载数据 打开Vista软件后选择加入2D的SEG-Y格式的原始地震数据，本实验 所用数据为给定的SHOT-20。加载后的原始地震数据如图2： 图2 原始地震数据显示 2）道均衡 各个道由于炮检距的不同，导致的反射波的振幅的变化，因为在共反射点叠加中，要求每一个叠加道的振幅都应该相等，每一道对叠加所做的贡献是等价的，无特殊情况，一般就以记录图中间的振幅为基准，使近激发点的地震道振幅减少，增加远离激发点的地震道记录的振幅。道均衡流程模块如图3，道均衡结果如图4： 图3 道均衡流程模块 3）建立观测系统 图5 观测系统显示4）初至拾取 初至拾取结果显示如图6： 图6 初至拾取结果显示 5）初至切除 地震记录上的初至波包括直达波和浅层折射波，它们能量强且有一定延续时间，对紧接而来的浅层反射波有干涉和破坏作用。另外，动校正后会引起波形畸变，浅层尤其厉害。对这些强能量初至波和动校正畸变引起的处理办法是“切除”，即将这些波的采样值全部变为零值（充零）。初至切除流程模块如图7，初至切 除结果如图8： 图7 初至切除流程模块 国际经济与贸易专业 外贸交易模拟实验 实验报告 学号__ _ 姓名__ _____ 班级_______ 指导老师___ _______ 实验地点 一、试验目的 在Simtrade这个虚拟贸易平台中，通过扮演不同的角色，熟练掌握各种业务技巧，了解到国际贸易的物流、资金流和业务流的运作方式，增强感性认识和社会适应能力，进一步巩固、深化已学过的理论知识，提高综合运用所学知识发现问题、解决问题的能力。 二、实习时间 本次实习时间为2011-03-01至2011-03-28，共计四周 三、试验小结 本次试验运用交易方式:①L/C+CIF②L/C+CFR③T/T+CIF④T/T+CFR⑤ D/P+FOB⑥T/T+FOB 涉及到的交易产品：①洋菇罐头( 柄及碎片)CANNED STEMS AND PIECES MUSHROOMS②甜玉米罐头CANNED SWEET CORN③荔枝罐头CANNED LITCHIS④芒果罐头CANNED MANGOES⑤名牌手提包FAMOUS-BRAND HANDBAG⑥香味蜡烛SCENTED CANDLE 四、试验遇到的问题及其解决方法 1.预算与实际支出相差较大 ?没有换算成本币 在CONTRACT111中，我的进口预算表和是实际发生额都有２+汇率的差异，原因就在于我公司注册资金为欧元，而样本中为美元，而我当初为了省时省力省脑，就全部依样画弧，没有转换成欧元导致。 ?集装箱计算的问题 在CONTRACT4442中，我出口商出口预算表海运费在计算过程中集装箱数为4，而实际确实12个集装箱，结果海运费预算和实际相差甚多，之后导致FOB 价格、利润等一系列数字有出入。 ?粗心大意 学生学号实验课成绩 学生实验报告书 实验课程名称材料成型数值模拟设计实验 开课学院材料学院 指导教师姓名 学生姓名 学生专业班级成型1001班 2012-- 2013学年第二学期 实验教学管理基本规范 实验是培养学生动手能力、分析解决问题能力的重要环节；实验报告是反映实验教学水平 与质量的重要依据。为加强实验过程管理，改革实验成绩考核方法，改善实验教学效果，提高 学生质量，特制定实验教学管理基本规范。 1、本规范适用于理工科类专业实验课程，文、经、管、计算机类实验课程可根据具体情况参 照执行或暂不执行。 2、每门实验课程一般会包括许多实验项目，除非常简单的验证演示性实验项目可以不写实验 报告外，其他实验项目均应按本格式完成实验报告。 3、实验报告应由实验预习、实验过程、结果分析三大部分组成。每部分均在实验成绩中占一 定比例。各部分成绩的观测点、考核目标、所占比例可参考附表执行。各专业也可以根据具体情况，调整考核内容和评分标准。 4、学生必须在完成实验预习内容的前提下进行实验。教师要在实验过程中抽查学生预习情况， 在学生离开实验室前，检查学生实验操作和记录情况，并在实验报告第二部分教师签字栏签名，以确保实验记录的真实性。 5、教师应及时评阅学生的实验报告并给出各实验项目成绩，完整保存实验报告。在完成所有 实验项目后，教师应按学生姓名将批改好的各实验项目实验报告装订成册，构成该实验课程总报告，按班级交课程承担单位（实验中心或实验室）保管存档。 6、实验课程成绩按其类型采取百分制或优、良、中、及格和不及格五级评定。 附表：实验考核参考内容及标准 观测点考核目标成绩组成 实验预习1．预习报告 2．提问 3．对于设计型实验，着重考查设计方案的 科学性、可行性和创新性 对实验目的和基本原理 的认识程度，对实验方 案的设计能力 20% 实验过程1．是否按时参加实验 2．对实验过程的熟悉程度 3．对基本操作的规范程度 4．对突发事件的应急处理能力 5．实验原始记录的完整程度 6．同学之间的团结协作精神 着重考查学生的实验态 度、基本操作技能；严 谨的治学态度、团结协 作精神 30% 结果分析1．所分析结果是否用原始记录数据 2．计算结果是否正确 3．实验结果分析是否合理 4．对于综合实验，各项内容之间是否有分 析、比较与判断等 考查学生对实验数据处 理和现象分析的能力； 对专业知识的综合应用 能力；事实求实的精神 50% 班级：B110904 学号：B11090230 姓名：徐敬伟 一、实验目的 国际贸易实务实验，是国际贸易实务理论课程之后的一个实践性环节，目的 在于把课程中所涉及的理论知识，应用到具体的贸易实务往来中，应达到的目标 包括： 1、利用SimTrade 提供的各项资源，做好交易前的准备工作。 2、学会运用网络资源宣传企业及产品。 3、使用邮件系统进行业务磋商，掌握往来函电的书写技巧。 4、掌握不同贸易术语在海运、保险方面的差异。在询盘、发盘、还盘、接 受环节的磋商过程中，灵活使用贸易术语（CIF 、CFR 、FOB ）与结算方式（L/C 、 T/T 、D/P 、D/A ） 5 6、掌握四种主要贸易术语（L/C 、T/T 、D/P 、D/A 7 8 9 10 体会国际贸易中不 出口商、进口商和工厂。出口商分别与 进口商、工厂发生交易，其中概要可归纳于下图中： 工厂资料： 1、推销。产品制造商和出口贸易商都需要积极开发市场，寻找贸易对象，可寄送业务推广函(Sale Letter)或在计算机网络、杂志、报刊上刊登产品广告来推销自己，同时可通过参加商展等途径寻找交易对手，增进贸易机会。 2、询盘。出口商收到工厂的业务推广函或看到广告后，根据自己的需要，对有意进一步洽商的工厂予以询盘，以期达成交易。 3、发盘。工厂按买主来函要求，计算报价回函给出口商。这期间可能需要 函电多次往返接洽，最后得到关于价格条款的一致意见。 4 (Contract或 以由工厂起草。 5、生产货物。签约后，工厂即着手生产货物。 6、交货。生产完成后，工厂依合同放货给出口商。 7、支付货款。工厂放货的同时，出口商支付货款，交易完成。 8、缴税。合同完成后，工厂还需到国税局就该笔合同的收益缴付税款，增值税率与综合费用费率可在"淘金网"的"其他费用"中查到，以合同金额乘之即得税款。 出口商（L/C方式下的履约流程） 1、推销。 2 根据自己的需要，对有意进一步洽商的出口商予以询盘 易。 3 4(Contract或Agreement)。注意起草与确认合同时 目录 1. 前言 (1) 上机实践的目的及要求 (1) 主要完成的实践内容 (2) 2. 油藏特征分析 (2) 储层物性特征 (2) 流体物性特征 (2) 储层岩石物性特征 (2) 气藏数值模型建立 (2) 模型网格的划分 (2) 模型物性 (3) 模型流体性质及相渗曲线 (3) XX气藏地质储量 (3) 4. XX气藏方案优选 (3) 开发方案的优选 (3) 采速与稳产时间的关系 (4) 5. 结论认识 (4) 结论 (4) 对本实践课程的建议 (4) 1. 前言 上机实践的目的及要求 1. 掌握油藏数值模拟的上机操作流程； 2. 掌握ECLIPSE软件的数据录入、编辑和修改方法； 3. 掌握ECLIPSE软件结果输出及三维可视化方法； 4. 掌握机理模型研究方案设计的思路及方法 主要完成的实践内容 1. 油藏数值模拟数值整理； 2. 依据现有数据，应用块中心网络系统建立一个三维油藏数值模拟模型； 3. 预测单口气藏天然能量开发的最终采收率（20年）（不考虑水体能量）； 4. 预测多口气井采收率（20年）； 5. 预测不同稳产年限下，气井的合理产量（稳产5年）； 6. 水平井开发和直井开发效果对比； 2. 油藏特征分析 储层物性特征 表2-1 储层物性特征 流体物性特征 气藏数值模型建立 模型网格的划分 模型流体性质及相渗曲线 XX气藏地质储量 4. XX气藏方案优选开发方案的优选 水平井方案 水平井方案 采速与稳产时间的关系 采油速度越快，稳产时间越短。采油速度越慢，稳产时间越长。由此可见采油速度与稳产时间成反比。 5. 结论认识 结论 通过这个实验，我们了解了eclipse软件的基本操作，并且建立了一个简单的均质油藏的模型，并且成功计算了产量。这个实验然我们获益匪浅。 对本实践课程的建议 建议增加实验课的课时，其余的方面都很好。老师讲的不错，需要学习的内容都学会了。 《地震数值模拟》实验报告 一、实验题目 叠加地震记录的相移波动方程正演模拟 二、实验目的 1.掌握各向同性介质任意构造、水平层状速度结构地质模型的相移波动方程正演模拟基本理论 2.实现方法与程序编制 3.由正演记录初步分析地震信号的分辨率。 三、实验原理 1、地震波传播的波动方程 设（x,z）为空间坐标，t为时间，地震波传播速度为v(x,z)，则二位介质中任意位置、任意时刻的地震波场为p(z,x,t)：压缩波——纵波。则二维各向同性均匀介质中地震波传播的遵循声波方程为 2、傅里叶变换的微分性质 p(t)与其傅里叶变换的P(w)的关系： 3、地震波传播的相移外推公式 令速度v不随x变化，只随z变化，则利用傅里叶变换微分性质把波动方程(变换到频率-波数域，得： 4、初始条件和边界条件 按照爆炸界面理论，反射界面震源在t=0时刻同时起爆，此时刻的波场就是震源。根据不同情况，可直接使用反射系数脉冲或子波作震源。如果直接使用反射系数作震源脉冲，则初始条件可表示为： 5、边界处理 （1）边界反射问题 把实际无穷空间区域中求解波场的问题化为有穷区域求解时，左右两边使用零边界条件。物理上假设探区距Xmin与Xmax两个端点很远，在两个端点上收到的反射波很弱。但是，上述条件在实际中不能成立，造成零边界条件反而成为绝对阻止波通过的强反射面。在正演模拟的剖面上出现了边界假反射干涉正常界面的反射。 （2）边界强反射的处理 镶边法、削波法、吸收边界都能有效消除边界强反射。 削波法就是在波场延拓过程中，没延拓一次，在其两侧均匀衰减到零，从而消除边界强反射的影响。假设横向总长度为NX，以两边Lx道吸波为例，有以下吸波公式： 四、实验内容 中国地 质大学 (武汉) 地空学 院 地震实 验报告 姓名:沈 班级:班 学号: 时间: 2015年05月 指导老师: 张 一、实验目的 实验一: 1、浅层地震装备的基本组成; 2、认识GEODE96浅层地震仪的主要结构,并学会该类仪器的操作方法; 3、地震波认识。 实验二: 1、掌握浅层地震数据采集方法及注意事项 二、仪器介绍 1、仪器简介 全套美国GEOMETRICS公司生产的Geode96浅层地震仪(相当于四套独立的24道浅层地震仪)该仪器能满足折、反射地震勘探、井间勘探、面波调查等地震监测需要,应用Crystal公司的A/D转换器与高速过采样技术达到了24位地震仪的精度。频带从1、75Hz到20,000Hz,使得采样间隔可以从20毫秒到16微秒。采样到的数据叠加到32位的叠加器中,然后传回到主机的硬盘或其它介质上。内置预触发器,每道有16K的内存。用硬件相关器对震源信号进行实时相关运算。Geode包装坚固、防水、防震,有提手,重4、1公斤,用12V的外接电池可以连续工作10个小时。(如下图) 2、主要操作功能键及快捷键 注释: 1锁定与解锁;2清除界面;4检测噪声;7保存 3、操作步骤及注意事项 1、每个GEODE用数传线按规定串联,通过数传盒与笔记本电脑的USB口连接。 2、每个GEODE接上12V电源。 3、开关接到与笔记本相连的第一个GEODE上。 4、传盒上的开关置于POWER UP处。 5、采集控制程序,并按工作需要设置好各项参数,然后进行正常数据采集工作。 6、出采集控制程序之前,应将数传盒上的开关置于POWER DOWN处。 7、卸下各连接线并清理整齐。 8、注意的就是:在正常工作过程中,任何时候移动数传线与GEODE的连接头时,必须退出采集控制程序。另外Y型头上有红色标记的与GEODE的前12道相连接。而且采集控制软件运行的语言环境必须就是英语(美国)。 三、实验内容 1、浅层地震装备认识及地震波认识:第一周上午主要就是老师介绍检波器、地震仪以及实验装备,认识设备后进行采集装置的连接,全班同学轮流当做指挥员与爆破员; 2、浅层地震数据采集实验:隔一周之后的上午全体同学使用地震仪进行浅层地震数据的采集及简单的分析,并对干扰波进行识别。 simtrade实习报告4000字 实习报告频道为大家整理的simtrade实习报告4000字，供大家参考。 大四上学期期末，我们做了simtrade的实习。Simtrade是一个关于国际贸易的实验平台。通过这个软件，我们可以模拟进行国际贸易，填制单据。从而整体掌握国际贸易流程。 一、实习目的阿 四年的大学即将结束，学生已基本完成国贸专业的各基础课、专业课程的学习任务。根据教学计划，学生要进行为期2周的专业课程实验。通过实验，可以使学生熟悉外贸实务的具体操作流程，增强感性认识，并可从中进一步了解、巩固与深化已经学过的理论和方法，提高发现问题、分析问题以及解决问题的能力 二、实习方法 通过进入simtrade模拟平台，进行模拟操作。 三、实习性质 本次实验属于外贸流程模拟实验，外贸流程模拟实验是实践性教学的重要环节之一。 四、实习时间分配 本次实验时间为XX-XX学年第一学期的19-20周，共2周时间。 五、实习地点 全部安排在实验室进行 六、实习意义 学生在机房进行国际货物买卖实务的具体操作，能很快掌握进出口的成本核算、询盘、发盘与还盘等各种基本技巧;熟悉国际贸易的物流、资金流与业务流的运作方式;切身体会到国际贸易中不同当事人面临的具体工作与他们之间的互动关系;学会外贸公司利用各种方式控制成本以达到利润最大化的思路;认识供求平衡、竞争等宏观经济现象，并且能够合理地加以利用。老师通过在网站发布新闻、调整商品成本与价格、调整汇率及各项费率等方式对国际贸易环境实施宏观调控，使学生在实习中充分发挥主观能动性，真正理解并吸收课堂中所学到的知识，为将来走上工作岗位打下良好基础。 七、实习心得 (一)对Simtrade软件的认识 Simtrade是一个国际贸易模拟实习软件。这个软件囊括了几乎所有基本的国际贸易的流程及涉及的国际贸易角色和单据。它在实习角色的分配上，包括了五个实习角色。工厂、出口商、进口商、出口地银行、进口地银行。这五个角色之间都有业务要做。在这五个角色之外，软件又设置了很多虚拟角色，它们可以辅助我们完成simtrade实训总结
计算机仿真与建模实验报告
地震记录簿数值模拟的褶积模型法
-地震勘探实验报告
地震记录数值模拟的褶积模型法
simtrade外贸实务实训报告
计算机模拟手工实验报告
地震勘探实验报告记录
Simtrade实验报告
deform3D实验报告
折射波勘探实验报告全解
地震勘探资料处理
simtrade实验报告
材料成型数值模拟设计实验
simtrade实习分析报告
油藏数值模拟实验报告
叠加地震记录的相移波动方程正演模拟数值模拟实验共22页
-地震勘探实验报告
simtrade实习报告 字