Rhombohedral Packing | Oil Authority

Rhombohedral packing (also called close-packing or face-centered cubic packing) is the tightest possible arrangement of equal spheres in three dimensions, in which each sphere is in contact with 12 neighboring spheres and the spheres occupy 74.05% of the total volume — compared to the loosest stable arrangement (cubic packing, where each sphere contacts 6 neighbors and spheres occupy 52.36% of the volume) — a theoretical model applied in petroleum engineering and petrophysics to understand the porosity of idealized granular sedimentary rock and to establish the theoretical minimum and maximum porosity bounds for unconsolidated sediment composed of uniform spherical grains.

Key Takeaways

The fundamental relationship between packing arrangement and porosity in sphere packs establishes the end-member porosity values for idealized granular sediment: cubic (orthorhombic) packing gives a theoretical maximum porosity of 47.64% (spheres in aligned layers, each sphere touching 6 neighbors); rhombohedral (close) packing gives a theoretical minimum porosity of 25.95% for a regular arrangement of equal spheres; real sediment porosity falls within this range (typically 25% to 45% for unconsolidated sands, depending on sorting, angularity, and compaction), with the actual packing geometry determined by grain size distribution, grain shape, and the depositional and burial history of the sediment.
In petroleum petrophysics, rhombohedral packing theory is used to understand the effect of grain sorting on porosity — a perfectly sorted sediment composed of identical spherical grains approaches rhombohedral packing under compaction and has lower porosity than a poorly sorted sediment where smaller grains partially fill the interstices between larger grains; this counterintuitive result (poor sorting gives lower absolute porosity than moderate sorting because bimodal size distributions allow efficient filling of pore space) is captured in the concept of the packing coordination number and is why the best reservoir sands in terms of both porosity and permeability are well-sorted, with a narrow grain size distribution that prevents fine filling of intergranular pore space.
Rhombohedral packing theory predicts a specific pore throat geometry that determines the permeability of the idealized sphere pack — the narrowest constriction between adjacent pore bodies in a close-packed sphere arrangement has a diameter that is 15.5% of the sphere diameter (for three touching spheres of equal radius), setting a theoretical minimum pore-to-throat size ratio that governs capillary entry pressure and the Washburn equation prediction of mercury injection pressure for the idealized pack; real sandstone pore throat sizes are consistently smaller than this theoretical value because angular grain contacts, cement overgrowths, and clay rims in the pore throats reduce throat sizes below the perfect sphere theoretical minimum.
The packing coordination number — the number of grain contacts per grain — is directly related to the effective stress supported by the grain framework and to the mechanical response of the sediment to burial compaction; rhombohedral packing with 12 contacts per grain provides the most stable mechanical structure and the most even stress distribution among all regular packing geometries, meaning that sediments that approach rhombohedral packing under compaction are more resistant to further porosity reduction under continued burial than loosely packed sediments with fewer grain contacts and larger void spaces that accommodate further compaction by grain rearrangement.
Cements and diagenetic modifications to the initial packing geometry fundamentally alter the applicability of theoretical packing models — quartz overgrowth cementation, carbonate cement, and clay mineral precipitation in pore throats reduce effective porosity and permeability below the values predicted by packing theory alone, while dissolution of unstable grains (feldspars, carbonates, organic matter) can create secondary porosity that increases effective porosity above the theoretical packing minimum; the theoretical packing models serve as reference end-members for understanding the maximum depth and compaction state at which the rock retains some depositional packing character before diagenesis has completely overprinted the original grain arrangement.

Fast Facts

The mathematical proof that rhombohedral (close) packing is the densest possible arrangement of equal spheres — the Kepler conjecture, stated by Johannes Kepler in 1611 — was not formally proven until 1998 when Thomas Hales published a computer-assisted proof confirming that no arrangement of equal spheres can achieve a packing density greater than 74.05%. The application of sphere packing theory to sedimentary petrology was developed by Fraser in 1935 and by Graton and Fraser in 1935, who calculated the theoretical porosity values for different packing arrangements and established the 25.95% minimum for rhombohedral packing that remains in every petrophysics textbook today. In practice, unconsolidated quartz sands range from 30% to 45% porosity, reflecting intermediate packing between the theoretical cubic and rhombohedral end-members, plus the influence of grain angularity and size sorting that deviate from the equal-sphere idealization.

What Is Rhombohedral Packing?

Imagine stacking oranges at a market — you intuitively place each row in the grooves formed by the row below, fitting the spheres together as tightly as possible. This natural stacking pattern is rhombohedral (or hexagonal close) packing, and it achieves the densest possible arrangement of equal spheres. Each sphere is in contact with 6 neighbors in its own layer, 3 neighbors in the layer below, and 3 in the layer above — a total coordination number of 12. The interstices (pore spaces) between the spheres are as small as they can be for equal spheres.

The opposite extreme is cubic packing — placing spheres in a perfectly regular grid where each sphere is directly above, below, and beside its neighbors in an aligned array. Here each sphere contacts only 6 neighbors (above, below, left, right, front, back), leaving large cubic void spaces at the centers of each unit cell. The pore space is much larger than in rhombohedral packing, giving a higher porosity (47.64% versus 25.95%).

These two idealized arrangements bracket the porosity of real granular sedimentary rocks. No arrangement of equal spheres can have porosity lower than 25.95% (without cementation reducing grain-to-grain contacts) or higher than 47.64% under gravity (a random pile of spheres). Real reservoir sands, with their angular grains, variable grain sizes, and diagenetic modifications, have porosities that fall within this range and are influenced by the same fundamental geometric principles captured by the packing theory end-members.

Rhombohedral Packing in Petrophysical Analysis

Reservoir petrophysicists use sphere packing theory as a reference framework for understanding measured core porosity values in the context of formation burial history and diagenesis. A well-sorted, clean quartz sandstone with 40% porosity at shallow depth (less than 1,000 meters) is in a packing state close to the cubic maximum; the same sandstone buried to 4,000 meters may show 15 to 20% porosity due to compaction (grain rearrangement from near-cubic toward near-rhombohedral packing), quartz overgrowth cementation (reducing inter-grain void space below the theoretical minimum for close packing), and clay mineral precipitation in pore throats. Tracking the porosity reduction pathway from initial deposition (near-maximum packing porosity) through burial (compaction-driven approach to rhombohedral packing) through diagenesis (cementation below the theoretical packing minimum) provides the quantitative framework for predicting reservoir quality at depth before drilling.

Permeability models calibrated to sphere packing geometry — particularly the Kozeny-Carman equation, which relates permeability to porosity and specific surface area through the hydraulic radius concept — use the relationship between packing coordination number, pore geometry, and pore throat size to predict permeability from porosity. The Kozeny-Carman model in its idealized form predicts that rhombohedral packing has a specific permeability (k) that scales with grain diameter squared and porosity cubed, providing a baseline permeability estimate for well-sorted reservoir sand at any given burial depth where the packing state is approximately known from the porosity measurement.

In NMR logging, the pore-to-throat size ratio of the actual reservoir rock is compared to the theoretical ratio for rhombohedral packing to assess diagenetic modification — a pore system with a pore-to-throat ratio significantly exceeding the rhombohedral theoretical value (approximately 6.5 for ideal close-packed spheres) has undergone pore throat reduction by cement or clay relative to the pore body, creating the high entry pressure and low relative permeability to hydrocarbon that is characteristic of cemented or clay-damaged sandstones even when total porosity remains relatively high.

Rhombohedral Packing Across International Jurisdictions

Canada (AER / WCSB): WCSB reservoir quality analysis for Montney, Cardium, and Mannville Formation sandstones uses sphere packing models as reference frameworks for distinguishing compaction-dominated porosity loss (approaching rhombohedral packing with depth) from cementation-dominated porosity loss (below the theoretical minimum for packing without cement) in reservoir quality prediction maps. AER requires that reserve certification reports document the petrophysical model basis for porosity estimation, and understanding the theoretical packing limits provides the geological context for explaining why certain reservoir intervals have porosity below 25% (indicating cementation has reduced porosity below the rhombohedral packing minimum) or above 40% (indicating either very shallow burial or secondary dissolution porosity).

United States (API / BSEE): US reservoir quality prediction in the Gulf Coast Miocene and Tertiary sandstones, the Permian Basin Delaware and Bone Spring sands, and the Powder River Basin Cretaceous sands uses petrophysical models grounded in sphere packing theory to establish porosity-depth trends and to distinguish primary porosity (preserved from original deposition) from secondary dissolution porosity (created by grain dissolution during diagenesis). API RP 40 (Core Analysis) and API RP 44 formation evaluation procedures implicitly use packing theory concepts in the petrophysical transformations from core measurements to reservoir properties used in SEC reserve certification.

Norway (Sodir / NORSOK): North Sea Brent Group sandstone reservoir quality prediction — a critical uncertainty in North Sea exploration economics since Brent sands lose porosity rapidly with depth due to quartz cementation — uses sphere packing compaction models combined with diagenetic cement prediction to forecast reservoir quality at depth before drilling. Equinor's reservoir quality prediction program for Jurassic sandstones uses the rhombohedral packing porosity as the compaction endpoint below which all further porosity reduction must be attributed to cementation rather than mechanical grain rearrangement, providing a mechanistic separation of compaction and cementation effects in reservoir quality models for the deep (3,000 to 5,000 meter) Brent and Cook Formation sandstones.

Middle East (Saudi Aramco): The Arab Formation carbonate reservoirs are typically analyzed using carbonate-specific pore type classification (interparticle, intraparticle, vugular, fracture) rather than sphere packing models because carbonate grains (ooids, peloids, bioclasts) have more complex geometries than the spherical model assumes; however, the underlying concept of packing efficiency and minimum theoretical porosity applies to grainstone-dominated Arab D intervals where the well-sorted, approximately spherical ooid grains provide a good first approximation to the sphere packing model. Aramco's petrophysical model for Arab D interparticle porosity in grainstone facies uses rhombohedral packing geometry as the baseline for relating grain size to pore body and pore throat dimensions used in the capillary pressure prediction model.