频谱激电法作为地球物理勘探的有效手段，在寻找激电异常体以及探寻深部矿产方面有重大贡献，岩石频谱参数的确定是实现该方法的重要步骤。在利用最小二乘法对复电阻率、相位参数进行反演时，各频谱参数初始值选取不同对反演效果影响较大，选取不好甚至会发散。在各频谱参数所属范围内采用随机选取初始值的方法进行反演，对复电阻率幅值曲线以及相位曲线拟合效果较好，不易陷入局部极小值的特点，能克服了较为盲目选取初始值而对结果带来的不良影响。%Spectral induced polarization (SIP)is an effective means to explore an induced polariza-tion target and deep minerals in geophysical prospecting.Induced polarization (IP)is a geophysi-cal imaging technique used to identify subsurface materials,such as ore.The method of IP is sim-ilar to electrical resistivity tomography,in that an electric current is induced into the subsurface through two electrodes and voltage is monitored through two other electrodes.IP is a measure of a delayed voltage response in earth materials.The IP effect is caused by a current-induced elec-tron-transfer reaction between electrolyte ions and metallic luster minerals and is a measurement of the electrical energy storage capacity of the earth. The IP effect can be determined by passing an induced current into the ground and measuring the change in voltage with respect to time (TDIP),or by measuring changes in phase at a given frequency with respect to a reference phase (FDIP).To produce an IP effect,fluid-filled pores must be present,since the rock is essentially an insulator.The IP effect becomes evident when these pore spaces are in contact with metallic luster minerals,graphite,clays,or other alteration products.IP effects make the apparent resistivity of the host rock change with frequency.Gener-ally,the rock resistivity decreases as the measurement frequency increases.Understanding the rock spectrum parameters,which are obtained through resistivity and phase inversion,is the most important step. There are many methods applied in rock spectrum inversion.Generally,when using the least squares algorithm,the inversion result is primarily determined by the initial value of each param-eter.In this paper,the authors will show that the initial value based on the least squares algo-rithm can be selected randomly from the scope of various parameters.In order to demonstrate the practicality of this method,the theoretical data determined by Luo Yan-zhong is compared to shale measurement data (complex resistivity and phase)measured by an SI- 1260A Impedance Analyzer.The shale was saturated before starting the experiment and the measuring conditions were room temperature and standard pressure.During the measurement,the frequency was ap-plied at a total of 61 points,varying from 0.01 Hz to 10 kHz and the measuring mode was a sym-metric quadrupole. In order to prevent the partial saturation fluid to evaporate and alter the electrical character-istics of the rock samples,both of the shale core's ends were soaked in the saturation fluid.In or-der to reduce the electrode polarization effects,the annular high-purity platinum network was used as unpolarized electrodes.In order to reduce noise disturbance,each core was measured a minimum of three times until the error of two adjacent measurements was within the permissible range.The experiment was automatically controlled by a computer during the entire process.Im-pedance and phase parameters were obtained from the experiment.Complex resistivity could be calculated through the impedance according to the geometric parameters of the shale cores.The single and double Cole-Cole models are selected to study the complex resistivity properties of shale.The inversion results show that the complex resistivity and phase curves can be described well and the spectrum parameters also can be obtained precisely.Simultaneously,relative to blindly selecting the initial value,this method has the advantages of avoiding local minima and performing at a high speed,and it overcomes the previous detrimental impact on the results.
展开▼
