I have modelled porosity as a discontinuty in wall.
Suppose you had a certain population of porous holes of diameter d. If they touched eachother, then the discontinuty is obviously greater, so the wall is effectively thinner at that point. From Thick Wall Pressure Vessel Theory, compute the stress on a section of wall under triaxial loading. You'll find the resulting expression to be purely a function of vessel geometry (i.e. OD and ID) and internal pressure. There is a scalar multiplier, SQRT(3).
So given the OD and internal pressure P, you can effectively compute the minimum ID given a yield strength of the material. This means that you can theoretically compute the number of porous holes of diameter "d" which when magically perfectly aligned, contributes to the change in wall. This is your answer.
Clearly in the real world, porous holes are not aligned. So what you are looking at is a statistical approach where the proximity of one porous hole would be of influence to it's nearest neighbor. You will find that in the Kepler Packing Problem, the theory lends itself very well to your particular porous hole problem.
That will then lead you to a density of porosity acceptable under Von Mises-Hencky triaxial stress, the material strength of your vessel in question and porosity density allowed in a section of wall of thickness "t". You track you allowable porous population as a limit to shrinking safety factor, i.e. the partial derivative of wall thickness to number of holes.
The result is a linear equation. Surprise!
Kenneth J Hueston, PEng
Principal
Sturni-Hueston Engineering Inc
Edmonton, Alberta Canada