An analytical approach is presented to describe pressure-driven streaming current (I-str) and streaming potential (U-str) generation in geometrically complex samples, for which the classical Helmholtz-Smoluchowski (H-S) equation is known to be inaccurate. The new approach is valid under the same prerequisite conditions that are used for the development of the H-S equation, that is, the electrical double layers (EDLs) are sufficiently thin and surface conductivity and electroviscous effects are negligible. The analytical methodology is developed using linear velocity profiles to describe liquid flow inside of EDLs and using simplifying approximations to describe macroscopic flow. At first, a general expression is obtained to describe the I-str generated in different cross sections of an arbitrarily shaped sample. Thereafter, assuming that the generated Ustr varies only along the pressure-gradient direction, an expression describing the variation of generated U-str along the sample length is obtained. These expressions describing I-str and U-str generation constitute the theoretical foundation of this work, which is first applied to a set of three nonuniform cross-sectional capillaries and thereafter to a square array of cylindrical fibers (model porous media) for both parallel and transverse fiber orientation cases. Although analytical solutions cannot be obtained for real porous substrates because of their random structure, the new theory provides useful insights into the effect of important factors such as fiber orientation, sample porosity, and sample dimensions. The solutions obtained for the model porous media are used to device strategies for more accurate zeta potential determination of porous fiber plugs. The new approach could be thus useful in resolving the long-standing problem of sample geometry dependence of zeta potential measurements.