Consider a two component medium . The substances filling the sets and are assumed to be electrically homogeneous and characterized by their real frequency dependent conductivity and dielectric function . The real dielectric function and conductivity of the composite are then given as
in terms of the functions and characterizing the dielectric response of the constituents.
The propagation of electromagnetic waves in the composite medium is described by the macroscopic Maxwell equations (4.4)–(4.8). In the following the magnetic permeabilities are assumed to be unity to simplify the analysis. The time variation of the fields is taken to be proportional to . Fourier transforming and inserting (4.8) into (4.4) yields where the frequency dependent displacement field
combines the free current density and the polarization current. In the quasistatic approximation one assumes that the frequency is small enough such that the inductive term on the right hand side of Faradays law (4.6) can be neglected. Introducing the complex frequency dependent dielectric function
the electric field and the displacement are found to satisfy the equations
in the quasistatic approximations. If the electric field is replaced by the potential these equations assume the same form as (5.2), and hence the methods discussed in section V.A can be employed in their analysis.
The neglect of the induced electromagnetic force is justified if the wavelength or penetration depth of the radiation is large compared to the typical linear dimension of the scatterers. If the scattering is caused by heterogeneities on the micrometer scale as in many examples of interest the approximation will be valid well into the infrared region.
The electrical conductivity of rocks fully or partially saturated with brine is an important quantity for the reconstruction of subsurface geology from borehole logs [281, 282]. The main contribution to the total conductivity of a a sample of brine filled rock comes from the electrolyte. The contribution from the rock matrix is usually negligible. 4 (This is a footnote:) 4 Nonvanishing matrix conductivity does, however, occur in veinlike ores. The electrolyte filling the pore space contributes through its intrinsic electrolytic conductivity as well as through electrochemical interactions at the interface. The total dc conductivity of the sample is written as
where is the dc conductivity of the electrolyte (usually salt water) filling the pore space and denotes the conductivity resulting from the electrochemical boundary layer at the internal surface . The surface conductivity correlates well with the specific internal surface and indirectly with other quantities related to it. The factor is called electrical formation factor. If the salinity of the pore water is high or electrochemical effects are absent the second term in (5.20) can be neglected and the formation factor becomes identical with the dimensionless resistivity of the sample normalized by the water resistivity. In the following the formation factor will be used synonymously for the dimensionless inverse dc conductivity .
The formation factor is usually correlated with the bulk porosity in a relation known as “Archie’s first equation” 
where the so called cementation index scatters widely and often obeys [281, 282]. Smaller values of are associated empirically with loosely packed media while higher values are associated with more consolidated and compacted media. Equation (5.21) implies not only an algebraic correlation between a purely geometric quantity and a transport coefficient but it also states that porous rocks do not show a conductor insulator transition at any finite porosity.
The experimental evidence for the postulated algebraic correlation (5.21) between conductivity and porosity is weak. The available range of the porosity rarely spans more than a decade. The corresponding conductivity data scatter widely for measurements on porous rocks and other media [285, 286, 254, 281, 188, 287, 288]. The most reliable tests of Archies law have been performed on artificial porous media made from sintering glass beads [285, 254, 287]. These media have a microstructure very similar to sandstone and are at the same time free from electrochemical effects. A typical experimental result for glass beads is shown in Figure 18 .
Note the small range of porosities in the figure. The existence of nontrivial power law relations in such samples is better demonstrated by correlating the conductivity with the permeability [284, 170]. In other experiments on artificial media a mixture of rubber balls and water is successively compressed while monitoring its conductivity [289, 217]. These experiments show deviations from the pure algebraic behaviour postulated by (5.21). If the cementation ”exponent” in (5.21) is assumed to depend on then it increases at low porosities in agreement with the general trend that higher values correspond to a higher degree of compaction.
A much better confirmed observation on natural and artificial porous rocks is dielectric enhancement caused by the disorder in the microstructure [290, 291, 292, 293, 294, 287, 295, 89]. Dielectric enhancement due to disorder has been studied extensively in percolation theory and experiment [296, 297, 40]. An example is shown in figure 19 for the sintered glass bead media containing thin glass plates.
In these media interfacial conductivity and other electrochemical effects can be neglected . The frequency is plotted in units of the relaxation frequency of water. Although salt water and glass are essentially dispersion free in the frequency range shown in figure 19 their mixture shows a pronounced dispersion which exceeds the values of the dielectric constants of both components. Similar results can be found in . In  the dielectric response of a large number of sandstones and carbonates is given in terms of the empirical Cole-Cole formula . Interestingly the corresponding temporal relaxation function appears within the recent theory of nonequilibrium systems  which is the same theory on which the macroscopic local porosity distributions in section III.A.5 were based.
Dielectric mixing laws express the frequency dependent dielectric function or conductivity of a two component mixture in terms of the dielectric functions of the constituents [46, 35, 40, 31]. Spectral theories express the effective dielectric function in terms of an abstract pole spectrum which is independent of the dielectric functions and of the two constitutents filling the pore and matrix space [299, 300, 301, 302, 303, 293, 304, 305, 306, 307, 308, 309]. Theoretically the effective dielectric function may be written as
The constants and are the strength and location of the poles and they reflect the influence of the microgeometry. Unfortunately these parameters do not have a direct geometrical interpretation although under the assumption of stationarity and isotropy two sum rules are known which connect integrals of the pole spectrum with the bulk porosity .
The simplest geometric theories for the effective dielectric function are mean field theories. In these approximations a small spherical cell with a randomly valued dielectric constant is embedded into a homogeneous host medium of dielectric constant . Then the electrical analogue of eq. (5.3) becomes
where the average denotes an ensemble average using the probability density of . A two component medium can be represented by the binary probability density
containing as the only geometrical input parameter. The Clausius-Mossotti approximation [310, 46] for a two component medium is obtained by setting in the limit or in the limit in (5.24). In the latter case one obtains
which will be a good approximation at low porosities. Note that the Clausius-Mossotti approximation is not symmetrical under exchanging pore and matrix. A symmetrical and also self consistent approximation is obtained from (5.24) by setting . This leads to the symmetrical effective medium approximation for a two component medium
which could have been derived also from using equation (5.25) in (5.12). The effective medium approximation is a very good approximation for microstructures consisting of a small concentration of nonoverlapping spherical grains embedded in a host. Recently, much effort has been expended to show that the EMA becomes exact for certain pathological microstructures . The so called asymmetrical or differential effective medium approximation is obtained by iterating the Clausius-Mossotti equation which gives the effective conductivity to lowest order in [285, 312, 313]. One finds the result
for spherically shaped inclusions.
The symmetric and asymmetric effective medium appoximations can be generalized to ellipsoidal inclusions because the electric field and polarization inside the ellipsoid remain uniform in an applied external field [310, 40, 312]. For aligned oblate spheroids whose quadratic form is with the effective medium theory for a two component composite results in two coupled equations
where the index denotes vertical conductivities and the index horizontal conductivities. The two equations in (5.29) are coupled through
For spheroids with identical shape but isotropically distributed orientations
Recently local porosity theory has been proposed as an alternative generalization of effective medium theories [168, 169, 170, 171, 172, 173, 174, 175]. The simplest mean field theories (5.26), (5.27) and (5.28) are based on the simplest geometric characterization theories of section III.A.1. These theories are usually interpreted geometrically in terms of grain models(see section III.B.2) with spherical grains embedded into a homogeneous host material. The generalizations (5.29), (5.34) and (5.35) are obtained by generalizing the interpretation to more general grain models. Local porosity theory on the other hand is based on generalizing the geometric characterization by using local geometry distributions (see section III.A.5) rather than simply porosity or specific surface area alone. In III.A.5 two different types of local geometry distributions were introduced: Macroscopic distributions with infinitely large measurement cells defined in (3.52), and mesoscopic distributions with measurement cells of finite volume defined in (3.33). For a mesoscopic partitioning of the sample using a simple cubic lattice with cubic unit cell the selfconsistency equation of local porosity theory for reads
where is the local percolation probability defined in section III.A.5.d, and and are the local dielectric functions of percolating or conducting (index ) and nonpercolating or blocking (index ) measuremente cells. In (5.40) it is assumed that the local dielectric response depends only on the porosity, but this may be generalized to include other geometrical characteristics.
Equation (5.36) has two interesting special cases. For a cubic measurement lattice () in the limit in which the sidelength of the cubic cells is small the one cell local porosity distribution is given by (3.31) or (3.34) if the medium is mixing. Inserting (3.31) or (3.34) into (5.36) and using , , and yields equation (5.27) for traditional effective medium theory. Note however that in the limit the local porosities become highly correlated rendering a description of the geometry in terms of the one cell function more and more inadequate. This argument does not apply in the opposite limit in which the measurement cells become very large, For stationary media the local porosities in nonoverlapping measurement cells are uncorrelated in the limit . For stationary and mixing media the local porosity distribution
which is identical to (5.27) except for the replacement of by , by , and by . Repeating the same differential replacement arguments  that lead to (5.28) for instead of gives a differential version of mesoscopic local porosity theory in the limit
As discussed in section III.A.5.e the -distribution is not the only possible macroscopic limit. Macroscopically heterogeneous media are described by the limiting macroscopic local porosity densities defined in (3.52). The -distribution is obtained in the limit or in the degenerate case arising for . The macroscopic form of local porosity theory for the effective dielectric constant is given by the integral equation
with defined in (3.52) and and . If the parameters , , and , are expressed in terms of the moments according to (3.53) the resulting effective dielectric function is found to be a function of the bulk porosity and its fluctuations. This observation indicates the possibility to study Archie’s law within the framework of local porosity theory.
Archies law (5.21) concerns the effective dc conductivity, , and it can be studied by replacing with in all the formulas of the preceding section. For notational convenience the shorthand notation will be employed in this section. Archie’s law can then be discussed by replacing with throughout and setting . From the the Clausius-Mossotti formula (5.26) one obtains the relation
for and for . Thus the symmetrical effective medium theory predicts a percolation transition at and does not agree with Archie’s law (5.21) in this respect. The same conclusion holds for the anisotropic generalization of the symmetric theory to nonspherical inclusions given in (5.29).
For the asymmetric effective medium theory in its simplest form (5.28) one finds
consistent with Archies law (5.21) with cementation exponent . This expression has found much attention because it yields [285, 315, 314, 205, 312]. There are, however, several problems with equation (5.43). Its derivation implies that the solid component is not connected . The experimentally observed behaviour is often not algebraic, and if a power law is nevertheless assumed its exponent is often very different from . Most importantly, the frequency dependent theory does not predict sufficient dielectric enhancement. The first problem can be circumvented by generalizing to a three component medium [316, 313], the second can be overcome by considering nonspherical inclusions [285, 315, 314, 205, 312]. As an example the generalization to nonspherical grains (5.34) gives
with , and the uniform spheroid model gives a similar result. The most serious problem, however, is the fact pointed out in [285, 292] that (5.28) cannot reproduce the frequency dependence of and the observed dielectric enhancement. This will be discussed further in the next section.
Local porosity theory contains geometrical information above and beyond the average porosity and cosequently it predicts more general relationships between porosity and conductivity. In its simplest form equation (5.38) leads to
which may or may not have a percolation transition depending upon whether the equation
has a solution which can be interpreted as a critical porosity. Therefore the percolation threshold can arise at any porosity including , and this allows to reconcile percolation theory with Archie’s law. Note also that the behaviour is nonuniversal5 (This is a footnote:) 5 The statement in  that local porosity theory predicts Archies law with a universal exponent is incorrect. and depends on the local percolation probability function , and the local response . Similarly the differential form (5.39) of local porosity theory yields
and is the control parameter of the percolation transition
giving the total fraction of percolating local geometries. The result (5.48)applies if for all values of , and also if at arbitrary . It holds universally as long as
where depends on the moments because depends on them through (3.53).
Compaction and consolidation processes will in general change the local porosity distributions where its dependence on or the parameters has been suppressed. Assume that it is possible to describe the consolidation process as a one parameter family of local porosity distributions depending on a parameter which characterizes the compaction process. Then the total fraction , the bulk porosity , and the integral (5.49) become functions of . If it is possible to invert the relation then the fraction becomes and equally . Therefore equations (5.48) and (5.52) become porosity-conductivity relations which depend on the consolidation process.6 (This is a footnote:) 6 The statement in  that local porosity theory predicts Archies law with a universal cementation index is incorrect. If the condition (5.51) and the asymptotic expansions and hold then these equations yield Archies law (5.21) with a nonuniversal cementation index . If the condition (5.51) does not hold and then
The measured data are the points, the solid curve represents a fit through the data. This fit was chosen to indicate that the consolidation process of sintering glass beads is expected to show a percolation transition at a small but finite threshold [173, 174]. Note however that the data of Figure 20 are consistent with the form corresponding to .
The theoretical mixing laws for the frequency dependent dielectric function discussed in section V.B.3 can be compared with experiment. Spectral theories generally give good fits to the experimental data [293, 287] but do not allow a geometrical interpretation. Geometrical theories on the other hand contain independently observable geometric characteristics, and can be falsified by experiment.
The single parameter mean field theories (5.26), (5.12) and (5.28) contain only the bulk porosity as a geometrical quantity. They are generally unable to reproduce the observed dielectric dispersion and enhancement. This is illustrated in Figure 21 which shows the experimental measurements of the real part of the frequency dependent dielectric function as solid circles .
The results were obtained for a brine saturated sample of sintered m glass spheres. The porosity of the specimen was , and the water conductivity was mS/m. A cross sectional image of the pore space has been displayed in Figure 13. The frequency in the figure is dimensionless and measured in units of the relaxation frequency of water where and are the conductivity and dielectric constant of water which are constant over the frequency range of interest.7 (This is a footnote:) 7 For water with mS/m the relaxation frequency is MHz. With the porosity known from independent measurements the simple mean field mixing laws can be tested without adjustable parameters. The prediction of the Clausius Mossotti approximation (5.26) is shown as the solid line with water as the uniform background, and as the dashed line with glass as background. The prediction of the symmetrical effective medium theory (5.27) is shown as the dotted line. The asymmetrical (differential)effective medium scheme (5.28), shown as the dash-dotted curve, appears to reproduce the high frequency behaviour correctly, but does not give the low frequency enhancement.
To compare the experimental observations with the Sen-Scala-Cohen model (5.34) or with the uniform spheroid model (5.35) the depolarization factor of the ellipsoids has to be treated as a free fit parameter  In the case of local porosity theory the local porosity distributions have been measured independently from cross sections through the pore space using image processing techniques. The resulting distributions have been displayed in Figure 15 for different sizes of the cubic measurement cells.8 (This is a footnote:) 8The sidelength of the measurement cell and the depolarization factor have been denoted by the same symbol . Their distinction should be clear from the context. The local percolation probability function on the other hand has not yet been measured directly from pore space reconstruction. Instead the result for displayed as the power law fit in Figure 20 was combined with the fact that in the limit of large measurement cells in which becomes a -distribution concentrated at (see eq. (3.32) or (3.35)). These observations and measurements motivate the Ansatz treating as a single free fit parameter. Figure 22 shows fits for the frequency dependent real dielectric function and inverse formation factor . In Figure 22 the solid circles are the experimental results.
All curves represent one parameter fits to the experimental data. The solid curve is obtained from local porosity theory (5.36) using as the best value of the fit parameter. The local porosity distributions were those of Figue 15 for measurement cells of sidelength pixels. The dashed curve corresponds to the uniform spheroid model (5.35) with a depolarization factor of as the best value of its fit parameter. The dotted curve represents the Sen-Scala-Cohen model with a best value of for the depolarization.
Similar experimental results for the dielectric dispersion have been observed in natural rock samples . Figure 8 in  compares the measurements only to the uniform spheroid model. Similar to the results of  on sintered glass beads the uniform spheroid model did not reproduce the dielectric enhancement, and required too high aspect ratios to be realistic for the observed microstructure.
Local porosity theory has also been used to estimate the broadening of the dielectric relaxation of polymers blends .