l_data.m

地震、测井方面matlab代码,解释的比较详细

function wlog=l_data
% Create synthetic well log for test and demo purposes; the log curves are not
% based on any rock physics
% Written by: E. R.: October 1, 2005;
% Last updated:
%
%      wlog=l_data

param.seed=1111;
param.stdev=2;
param.first=5000;
param.last=10000;
param.step=0.5;
depth=(param.first:param.step:param.last)';
nsamp=length(depth);
param.mininc=1;
param.maxinc=10;
param.avinc=1;
param.maxindex=600;

%	Set seed for random number generator to assure that the same log 
%       is created each time the function is run
randn('state',param.seed);
rand('state',param.seed);

curves=zeros(nsamp,6);

curves(:,1)=depth;
index=index4log_sampling1(nsamp,1,param);
samples=get_samples;
temp=260-samples(index);
dt=temp.*(7500./depth).^0.1;

temp1=(temp-122)/18;
vclay=min((1-temp1),1).^0.75;
vp=1.0e6./dt;

curves(:,2)=dt;
curves(:,3)=vp;
curves(:,4)=0.23*curves(:,3).^0.25 + 0.01*randn(nsamp,1);
curves(:,5)=0.7742*vp-884/0.3 + 100*randn(nsamp,1);
curves(:,6)=vclay;

wlog=l_convert(curves,{'depth','ft','Depth';
                       'DTp','us/ft','P-sonic';
		       'Vp','ft/s','P-velocity';
		       'rho','g/cm3','Bulk density';
		       'Vs','ft/s','S-velocity';
		       'Vclay','fraction','Clay volume fraction'});


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function samples=get_samples

samples = [ ...
126.7580  131.4790  130.2920  132.6880  130.2100  130.3830
  126.6060  131.0700  130.7870  132.6880  130.8430  130.3890
  126.5780  130.9710  131.1850  132.6470  131.3070  130.6120
  126.5230  130.9590  131.3560  132.3690  131.7660  131.6500
  126.4210  130.9590  131.4890  131.9430  133.1190  132.7200
  126.1570  131.0720  131.6260  131.5700  134.1980  133.7470
  125.9010  131.2030  131.8070  131.1530  135.3200  134.3730
  125.7170  131.2330  132.1230  130.8340  135.8210  134.3730
  125.5650  131.2330  132.5230  130.5380  135.8370  134.3630
  125.5220  130.9260  133.1060  130.2070  135.8370  134.0780
  125.5220  130.4290  133.8630  129.9380  135.6860  133.2580
  125.6470  130.1760  134.6570  129.6140  135.6490  132.4390
  125.9450  130.1150  134.8390  129.4060  135.6100  131.8320
  126.0470  130.0060  134.8490  129.2500  135.5230  131.5210
  126.0470  129.9050  134.8940  129.0410  135.5210  131.2050
  125.8870  129.8600  134.9030  129.0190  135.5210  130.7620
  125.4950  129.8320  135.0460  129.0190  135.7420  130.1920
  124.7290  129.7390  135.4940  129.2880  135.7490  129.0400
  123.8910  129.6240  136.1950  129.5060  135.7490  127.9400
  122.9530  129.3830  136.9140  129.8960  135.7280  127.2380
  122.4160  129.2950  137.6440  130.2820  135.3780  126.7400
  122.3140  129.2870  137.8730  130.7060  134.9270  126.7400
  122.3140  129.2870  137.8950  131.0050  134.5500  126.7610
  122.7270  129.3640  137.8950  131.1670  134.3010  127.2440
  123.3420  129.4450  137.7670  131.2920  134.2660  128.2100
  123.7730  129.6900  137.5350  131.3040  134.1880  129.2390
  124.6810  129.8600  137.1920  131.3040  134.1610  130.2910
  125.7170  130.0130  136.7250  131.2200  134.1600  131.3010
  126.7860  130.0480  136.2790  131.2120  134.1600  131.3010
  127.8400  130.0480  136.0510  131.2120  134.1850  131.2640
  128.8590  129.9940  135.7860  131.3570  134.4230  130.5090
  129.8430  129.5980  135.4980  131.4160  134.6660  128.7100
  130.5150  129.0460  135.4320  131.4210  134.7780  127.3110
  130.6700  128.4650  135.2420  131.7010  135.0060  127.1000
  130.6700  127.9750  134.9930  132.0270  135.2640  127.0590
  130.6700  127.6460  134.7110  132.3230  135.5900  126.7570
  130.5830  127.5130  134.4180  132.6130  135.8880  126.7570
  130.5500  127.4050  134.0400  133.1320  136.2930  126.6980
  130.5500  127.3680  133.8510  133.6000  136.3460  125.6670
  130.5670  127.3650  133.8270  133.6400  136.3470  124.3410
  130.5950  127.3650  133.8270  133.6400  136.3470  123.1540
  130.6020  127.6100  134.0830  133.4070  136.3470  123.1540
  130.6020  127.9810  134.3110  133.1980  136.3470  123.1640
  130.5360  128.6520  134.3430  132.9560  136.3660  123.3450
  130.3880  129.5570  134.3430  132.8900  136.6400  123.8380
  130.2310  130.5290  134.0530  132.8860  136.9470  124.4020
  130.1620  131.1780  133.3760  132.8860  136.9720  124.9500
  130.1520  131.6500  132.7110  133.0920  136.9720  125.2910
  130.0960  131.9540  131.9920  133.4460  136.9720  125.6450
  130.0410  132.1620  131.6580  133.9350  136.6560  125.9050
  129.8520  132.2820  131.4820  134.3980  136.2130  126.5630
  129.2360  132.3110  131.4250  134.8510  135.6550  127.6620
  128.2480  132.5820  131.4150  135.1120  134.9970  129.0170
  127.1710  133.0700  131.3380  135.1310  134.3140  130.4490
  126.1590  133.6920  131.3010  135.1310  133.6120  131.8110
  125.3320  134.3030  131.1470  135.0750  133.0530  132.8950
  124.6080  134.8130  130.9210  134.7890  132.6300  134.1530
  123.9240  134.9130  130.7970  134.5710  132.3540  134.7620
  123.3900  134.9130  130.7830  134.5310  132.0790  134.7620
  123.0650  134.5470  130.7830  134.5290  131.7980  134.7490
  122.8740  133.8190  130.9500  134.5260  131.7970  134.5740
  122.2860  133.1220  131.2020  134.5250  131.7700  134.3270
  122.0670  132.6960  131.4000  134.4450  131.7290  134.0950
  122.0580  132.3230  131.5950  134.3910  131.6810  134.0400
  122.0580  132.1270  131.6270  134.3880  131.5810  134.0400
  122.4270  132.0060  131.6690  134.2530  131.1740  134.0570
  122.5580  131.9830  131.6750  134.0550  130.6770  134.2710
  122.5580  131.9830  131.6750  133.4900  130.1380  134.5960
  122.2450  132.3970  131.6100  132.9130  129.6020  135.1060
  121.6370  132.8980  131.2720  132.5510  129.0840  135.1450
  121.0730  133.0730  130.9810  132.1400  128.6030  135.1590
  120.9370  133.0910  130.7690  131.7990  128.0460  135.1670
  120.9370  133.0910  130.5120  131.3720  127.6620  135.1670
  121.8230  132.8820  130.4780  131.2880  127.3470  135.1430
  122.8680  132.5690  130.4780  131.2850  127.3470  134.7340
  123.9540  132.0310  130.5320  131.2850  127.3500  133.0280
  124.9510  131.4500  130.9420  131.2950  127.6120  133.0280
  125.7890  130.8300  131.3850  131.4470  128.0530  133.1180
  126.6410  130.3230  131.7800  131.4560  128.5300  134.1650
  127.6710  130.2320  132.0780  131.4560  129.2920  135.6210
  128.7040  130.2320  132.2180  131.3030  129.9740  136.5350
  129.5520  130.2100  132.2330  130.8200  129.9820  137.0570
  129.7520  130.2040  132.2330  129.8110  129.9820  137.0570
  129.8700  130.2040  132.1630  128.7560  129.9680  136.9310
  130.0020  130.0550  132.1530  127.7430  129.2750  135.5710
  130.1130  129.8800  132.1460  127.6900  128.6630  134.1290
  130.4200  129.6850  132.1310  127.6900  128.4590  133.7610
  130.5130  129.5360  132.0230  127.7180  128.4590  133.5790
  130.5440  129.5100  131.9870  128.5590  128.4740  133.2700
  130.5540  129.5100  131.8040  129.4410  129.1340  132.5420
  130.6390  129.5100  131.7800  130.3550  129.5830  131.7970
  131.2340  129.5100  131.5990  130.4400  129.9870  131.7010
  132.1410  129.2820  131.3900  130.4420  129.9870  131.7010
  132.9860  128.9880  131.3650  130.4420  129.9890  131.8310
  133.7560  128.7250  131.3650  129.3130  129.9890  132.9850
  133.9410  128.6790  131.4430  128.3370  130.0010  133.0630
  133.9410  128.6790  131.7690  128.2910  130.4090  133.2970
  133.8840  128.7520  132.1980  128.2910  130.4090  133.2970
  133.1810  129.2420  132.4690  128.6900  130.4080  133.1390
  132.3590  129.7320  132.6660  129.4290  130.3830  131.6610
];
samples=samples(:);