%w為角頻率，lamda，u為彈性參數(shù)，fai為孔隙度，e為裂隙密度，tao為時(shí)間尺度參數(shù)，en為流體粘滯系數(shù)，vmiu為泊松比
%r為裂隙縱橫比，xiaok為滲透率，cv為裂隙體積a為裂隙半徑，kf為流體體積模量
clear?all;
vp=2790;vs=1463;rou=2.08;fai=0.3;tao=2*10^-5;e=0.1;miu=0.3104;v=miu;xiaok=1000;
lamda=rou*vp*vp-2*rou*vs*vs;r=0.1;
u=rou*vs*vs;kf=2E9;
k=rou*（vp^2-4*v*vs/3）;
kc=?（pi*u*r*（lamda?+?u））/（kf*（lamda?+?2.0*u））;
u=rou*vs*vs;
%k=2000;
en=0.01;a=0.1;l=2;
f=5:1:104;
%以下為輸入的值；
w=2*pi*f;
cv=4*pi*a*a*a*r/3;
sigma=tao*6*xiaok*l/（en*cv*（1+kc））;
kf=（pi*u*u*r*r）/（2*（1-v）*sigma）;

%tao=4*en*a*a*a*（1-v）/（9*xiaok*l*u）;%裂隙縱橫比很小的時(shí)候

kp=4*u/（3*kf）;

y=（3*pi*（lamda+u）*（1+kp））/（4*（lamda+2*u）*（1+kc））;
yyipie=y*（lamda+2*u）/（（3*lamda+2*u）*（1+kp））;
awfenzi=（（ones（1100）+i*w*y*tao）*（lamda+2*u）/（lamda+u））*（16*e/（27*fai*（1+kp））+（lamda+u）/（3*lamda+2*miu））+i*w*tao*（1/（3*（1+kc））-yy

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件????????1753??2018-05-14?10:53??chapmann.m